|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
这是本人收集高手们的作品,络加修改
- x=0:0.5:10;
- z=x.^2;
- y=zeros(size(x));
- theta=pi/20;
- xx=x;yy=y;zz=z;
- for i=1:40
- M=[cos(i*theta) sin(i*theta) 0;-sin(i*theta) cos(i*theta) 0;0 0 1];
- temp=M*[x;y;z];
-
- xx(i+1,:)=temp(1,:);yy(i+1,:)=temp(2,:);zz(i+1,:)=temp(3,:);
- axis([-12 12 -12 12 -10 100])
- mesh(xx,yy,zz);
- axis off
- pause(0.5)
- if i==40
- break
- end
- clf
- end
复制代码
- clear;
- z=0:100;
- for n=1:100
- t= 2*pi*(n-1)/100:1/2*pi/100:2*pi*n/100;
- X=sqrt(z)'*cos(t);
- Y=sqrt(z)'*sin(t);
- Z=z'*ones(1,5);
- surf(X,Y,Z);
- colormap(spring)
- shading interp
- axis square
- box on
- axis([-11 11 -11 11 -1 101])
- hold on
- axis off
- pause(0)
- end
复制代码
- x=0:0.5:10;
- z=x.^2;
- y=zeros(size(x));
- theta=pi/20;
- xx=x;yy=y;zz=z;
- for i=1:40
- M=[cos(i*theta) sin(i*theta) 0;-sin(i*theta) cos(i*theta) 0;0 0 1];
- temp=M*[x;y;z];
-
- xx(i+1,:)=temp(1,:);yy(i+1,:)=temp(2,:);zz(i+1,:)=temp(3,:);
- axis([-12 12 -12 12 -10 100])
- mesh(xx,yy,zz);
- axis off
- pause(0.5)
- if i==40
- break
- end
- clf
- end
复制代码 |
|