构造旋转体

aspen

1. % 构造旋转体
   
2. % 圆柱
   
3. subplot(221);
   
4. t=ones(1,40)*5;
   
5. [x,y,z]=cylinder(t);
   
6. c=ones(size(z));
   
7. surf(x,y,z,c)
   
8. shading faceted
   
9. title('圆柱')
   
10. axis off
    
11. 
    
12. % 圆锥
    
13. subplot(222);
    
14. t=0:.1:2;
    
15. [x,y,z]=cylinder(fliplr(t));
    
16. c=ones(size(z));
    
17. surf(x,y,z,c)
    
18. shading faceted
    
19. title('圆椎')
    
20. axis off
    
21. 
    
22. % 椭球
    
23. subplot(223)
    
24. t=linspace(0,pi,100);
    
25. [x,y,z]=cylinder(sin(t)/4);
    
26. c=ones(size(z));
    
27. mesh(x,y,z,c)
    
28. shading faceted
    
29. title('椭球')
    
30. set(gca,'position',[0.13 0 0.327023 0.44])
    
31. axis off
    
32. 
    
33. % 一般旋转体
    
34. subplot(224)
    
35. t=linspace(0,pi,100);
    
36. [x,y,z]=cylinder(2+cos(t).^2);
    
37. c=ones(size(z));
    
38. mesh(x,y,z,c)
    
39. shading faceted
    
40. title('一般旋转体')
    
41. set(gca,'position',[0.56 0 0.327023 0.44])
    
42. axis off

复制代码

aspen

1. t=[1:10000];
   
2. tr=randn(size(t));
   
3. r=cos(t).*cos(t*6)+cos(t*9).*cos(tr*9);
   
4. x=r.*cos(t+pi/2);
   
5. y=r.*sin(t+pi/2);
   
6. H=plot(x,y);
   
7. k=1;
   
8. xlabel('Please press "space" key and stop this program!',...
   
9.   'fontsize',12,'color','r');
   
10. set(gcf,'doublebuffer','on');
    
11. d=0;
    
12. while k
    
13.    s=get(gcf,'currentkey');
    
14.    if strcmp(s,'space');
    
15.        clc;k=0;
    
16.    end
    
17.    t=t+cos(d+1)*10;
    
18.    r=0.2*sin(t)+sin(3*t).*sin(5*t)+1.9*sin(7*t).*sin(9*t);
    
19.    x=r.*cos(t);
    
20.    y=r.*sin(t);
    
21.    set(H,'xdata',x,'ydata',y);
    
22.    pause(0.1);
    
23. end
    
24. figure(gcf);

复制代码

mabiyang001

yanhaoguo

顶一下，最好再能做个动画的

happy

yanhaoguo 发表于 2010-10-26 16:17

                               
登录/注册后可看大图

顶一下，最好再能做个动画的

参考论坛动画的相关程序修改就能够实现相关的旋转动画效果

yanhaoguo

回复 happy 的帖子

好，有时间去学习下

happy

yanhaoguo 发表于 2010-10-27 10:05

                               
登录/注册后可看大图

回复 happy 的帖子

好，有时间去学习下

改出来后分享一下，

zhouyang664

献上旋绕动画一个：曲线xy=1绕z轴旋转一周
figure;
view(3);
axis([-3,3,-3,3,0,3]);
axis off
hold all;
plot3([0,0],[0,0],[-3,3],'k')
plot3([-3,3],[0,0],[0,0],'k')
plot3([0,0],[-3,3],[0,0],'k')
ax=[3,2.7,2.7;0,.2,-.2;0,0,0];
fill3(ax(1,:),ax(2,:),ax(3,:),'k')
ax=[0,.2,-.2;0,0,0;3,2.7,2.7];
fill3(ax(1,:),ax(2,:),ax(3,:),'k')
ax=[0,.2,-.2;3,2.7,2.7;0,0,0];
fill3(ax(1,:),ax(2,:),ax(3,:),'k')%绘制坐标系
t=1/2:.001:2;
for theta=0:pi/100:2*pi
z=t;x=sin(theta)*1./t;y=cos(theta)*1./t;
line(x,y,z,'color',rand(1,3));
pause(.005)
end

通过效果图可以看出，line绘制的是曲线图，而mesh绘制的才是曲面图，
不知有没有高手可以奉献一个曲面的动画？或者是mesh就不能产生这种类型的动画？