File:Total variation.gif
页面内容不支持其他语言。
Total_variation.gif (179 × 236像素,文件大小:57 KB,MIME类型:image/gif、循环、15帧、3.0秒)
摘要
描述Total variation.gif |
English: Illustration of total variation |
日期 | (UTC) |
来源 | 自己的作品 |
作者 | Oleg Alexandrov |
许可协议
Public domainPublic domainfalsefalse |
![]() |
我,本作品著作权人,释出本作品至公有领域。这适用于全世界。 在一些国家这可能不合法;如果是这样的话,那么: 我无条件地授予任何人以任何目的使用本作品的权利,除非这些条件是法律规定所必需的。 |
Source code (MATLAB)
% Illustration of total variation
function main()
% KSmrq's colors
red = [0.867 0.06 0.14];
blue = [0, 129, 205]/256;
green = [0, 200, 70]/256;
yellow = [254, 194, 0]/256;
white = 0.99*[1, 1, 1];
gray = 0.5*[1, 1, 1];
black = [0, 0, 0];
% Set up the grid and other parameters
N = 300;
A = 0.2; B = 2.8+A;
C=-1; D = 3;
X = linspace(A, B, N);
s = 1.4 + A;
Y = (X-s).^3 - (X - s);
Y = (Y-min(Y))/(max(Y)-min(Y))*(D-C)+C;
lw = 3.4; % linewidth
fs = 28; % font size
smallrad = 0.07;
numP = 15;
for p=1:numP
t = (p-1)/(numP-1);
x = A+t*(B-A);
y = interp1(X, Y, x);
figure(1); clf;
set(gca, 'fontsize', fs);
set(gca, 'linewidth', 0.4*lw)
hold on;
shifty1 = 1.2*smallrad;
shifty2 = 0.5;
plot_axes (-0.4, B+0.1, C-shifty1, D + shifty2, lw/1.5, fs, gray);
axis([A-0.5, B, C-shifty1, D + shifty2]);
plot([0, x], [y, y], 'linewidth', lw/2, 'linestyle', '--', 'color', gray);
plot(X, Y, 'color', blue, 'linewidth', lw);
ball(x, y, smallrad, green);
ball(0, y, smallrad, red);
axis equal; axis off;
% save to disk
file = sprintf('Frame%d.eps', 1000+p);
disp(file);
saveas(gcf, file, 'psc2')
pause(0.1);
end
% Converted to gif with the command
% convert -antialias -loop 10000 -delay 20 -scale 65% -compress LZW Frame10* Total_variation.gif
function plot_axes (A, B, C, D, lw, fs, color)
arrow_size = 0.3;
sharpness = pi/7;
arrow_type = 1;
arrow([A, 0], [B, 0], lw, arrow_size, sharpness, arrow_type, color);
arrow([0, C], [0, D], lw, arrow_size, sharpness, arrow_type, color);
small = 0.4;
text(B-0.5*small, -0.7*small, 'x', 'fontsize', fs);
text(-small, D - 0.3*small, 'y', 'fontsize', fs);
function arrow(start, stop, thickness, arrow_size, sharpness, arrow_type, color)
% Function arguments:
% start, stop: start and end coordinates of arrow, vectors of size 2
% thickness: thickness of arrow stick
% arrow_size: the size of the two sides of the angle in this picture ->
% sharpness: angle between the arrow stick and arrow side, in radians
% arrow_type: 1 for filled arrow, otherwise the arrow will be just two segments
% color: arrow color, a vector of length three with values in [0, 1]
% convert to complex numbers
i=sqrt(-1);
start=start(1)+i*start(2); stop=stop(1)+i*stop(2);
rotate_angle=exp(i*sharpness);
% points making up the arrow tip (besides the "stop" point)
point1 = stop - (arrow_size*rotate_angle)*(stop-start)/abs(stop-start);
point2 = stop - (arrow_size/rotate_angle)*(stop-start)/abs(stop-start);
if arrow_type==1 % filled arrow
% plot the stick, but not till the end, looks bad
t=0.5*arrow_size*cos(sharpness)/abs(stop-start); stop1=t*start+(1-t)*stop;
plot(real([start, stop1]), imag([start, stop1]), 'LineWidth', thickness, 'Color', color);
% fill the arrow
H=fill(real([stop, point1, point2]), imag([stop, point1, point2]), color);
set(H, 'EdgeColor', 'none')
else % two-segment arrow
plot(real([start, stop]), imag([start, stop]), 'LineWidth', thickness, 'Color', color);
plot(real([stop, point1]), imag([stop, point1]), 'LineWidth', thickness, 'Color', color);
plot(real([stop, point2]), imag([stop, point2]), 'LineWidth', thickness, 'Color', color);
end
function ball(x, y, radius, color) % draw a ball of given uniform color
Theta=0:0.1:2*pi;
X=radius*cos(Theta)+x;
Y=radius*sin(Theta)+y;
H=fill(X, Y, color);
set(H, 'EdgeColor', color);
说明
添加一行文字以描述该文件所表现的内容
此文件中描述的项目
描绘内容
创作作者 简体中文(已转写)
某些值没有维基数据项目
版权状态 简体中文(已转写)
受版权保护,并由版权所有者释出于公有领域 简体中文(已转写)
版权所有 简体中文(已转写)
著作权持有者释出至公有领域 简体中文(已转写)
25 9 2008
文件历史
点击某个日期/时间查看对应时刻的文件。
日期/时间 | 缩略图 | 大小 | 用户 | 备注 | |
---|---|---|---|---|---|
当前 | 2008年9月25日 (四) 04:46 | ![]() | 179 × 236(57 KB) | Oleg Alexandrov | {{Information |Description={{en|1=Illustration of total variation}} |Source=Own work by uploader |Author=Oleg Alexandrov |Date=~~~~~ |Permission= |other_versions= }} <!--{{ImageUpload|full}}--> |
文件用途
以下页面使用本文件:
全域文件用途
以下其他wiki使用此文件:
- fr.wikipedia.org上的用途
- nl.wikipedia.org上的用途