tch45.m

来自「MATLAB实现简单图象处理的程序,在matlab2006a上通过」· M 代码 · 共 41 行

t = -pi:pi/200:8*pi;
clf,set(gcf,'position',[348 31 450 270])
clf,plot3(cos(t),sin(t),t)
clf,comet3(cos(t),sin(t),t)
comet(t,tan(sin(t))-sin(tan(t)))

t_final=100; x0=[0; 0; 1e-10];
[t,x] = ode45('lorenzeq',[0, t_final], x0);
plot3(x(:,1), x(:,2), x(:,3)); 
comet3(x(:,1), x(:,2), x(:,3)); 

subplot(221),sphere(3), title('n=3'),axis equal
subplot(222),sphere(6), title('n=6'),axis equal
subplot(223),sphere,    title('n=20'),axis equal
subplot(224),sphere(50),title('n=50'),axis equal

t=0:.2:2*pi; R=cos(t);
subplot(221),cylinder(R,3), title('n=3'),axis equal
subplot(222),cylinder(R,6), title('n=6'),axis equal
subplot(223),cylinder(R),   title('n=20'),axis equal
subplot(224),cylinder(R,50),title('n=50'),axis equal

t=0:.2:2*pi; R=[0 5];
subplot(221),cylinder(R,3), title('n=3')
subplot(222),cylinder(R,6), title('n=6')
subplot(223),cylinder(R),   title('n=20')
subplot(224),cylinder(R,50),title('n=50')

[x,y]=meshgrid(-3:0.1:3,-2:0.1:2);
z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y);
subplot(122)
mesh(x,y,z)
axis([-3 3 -2 2 -0.7 1.5])
surf(x,y,z)
colorbar
colormap('cool')

x=0:0.2:pi;  R=exp(-x.^2/2).*sin(x);
[x,y,z]=cylinder(R,20);
surf(x,y,z*pi),axis([-0.6,0.6,-0.6,0.6,0,pi])