注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

Matlab绘制常见三维图形--圆环面如何画  

2012-03-17 13:30:27|  分类: 绘图 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
%% 圆环面
R = 6; r = 2;
syms u v;
ezmesh((R+r*cos(u))*cos(v),(R+r*cos(u))*sin(v),r*sin(u));
axis equal;

%% 圆盘
R = 6; r = 2;
theta=linspace(0,2*pi,90);ph=linspace(r,R,30);
[t,p]=meshgrid(theta,ph);
r=t*0;
[x,y,z]=pol2cart(t,p,r);
mesh(x,y,z);
axis equal;



%% a.画椭球  
a=1;b=2;c=3;  
figure(1),title('椭球')  
ellipsoid(0,0,0,a,b,c)  
%  
shading interp  
box on;  
grid off  
colormap hsv  
%% b.单叶双曲面  
[x,y]=meshgrid(-2:0.2:2,-2:0.2:2);  
z=sqrt(c.^2*(x.^2/(a.^2)+y.^2/(b.^2)-1));  
%  
figure(2)  
subplot(121),title('单叶双曲面')  
surf(x,y,real(z))  
hold on;  
surf(x,y,-real(z))  
%  
[u,v]=meshgrid(linspace(0,2*pi,50),linspace(-1*pi/2,pi/2,50));  
x=a.*sec(v).*cos(u);  
y=a.*sec(v).*sin(u);  
z=c.*tan(v);  
subplot(122),title('单叶双曲面')  
surf(x,y,z)  
%  
shading interp  
box on;  
grid off  
colormap hsv  
%% c.双叶双曲面  
[x,y]=meshgrid(-2:0.2:2,-2:0.2:2);  
z=sqrt(c.^2*(x.^2/(a.^2)-y.^2/(b.^2)-1));  
%  
figure(3)  
subplot(121),title('双叶双曲面')  
surf(x,y,real(z))  
hold on;  
surf(x,y,-real(z))  
%  
[u,v]=meshgrid(linspace(0,2*pi,50),linspace(-1*pi/2,3*pi/2,50));  
x=a.*tan(v).*cos(u);  
y=a.*tan(v).*sin(u);  
z=c.*sec(v);  
subplot(122),title('单叶双曲面')  
mesh(x,y,z)  
%  
shading interp  
box on;  
grid off  
colormap hsv  
%% d.椭圆抛物面  
[x,y]=meshgrid(-2:0.2:2,-2:0.2:2);  
z=(x.^2/(a.^2)+y.^2/(b.^2))./2;  
%  
figure(4),title('椭圆抛物面')  
surf(x,y,z)  
%  
colormap gray  
  
%% e.马鞍面  
[x,y]=meshgrid(-2:0.2:2,-2:0.2:2);  
z=(x.^2/(a.^2) - 2.*y.^2/(b.^2));  
%  
figure(5),title('马鞍面')  
surf(x,y,z)  
%  
shading interp  
box on;  
grid off  
colormap hsv  
%% f.圆环面  
r=2;  
R=4;  
[u,v]=meshgrid(linspace(0,2*pi,50),linspace(0,2*pi,50));  
x=(R+r.*cos(v)).*cos(u);  
y=(R-r.*cos(v)).*sin(u);  
z=r.*sin(v);  
%  
figure(6)  
subplot(131),title('圆环面1')  
surf(x,y,z)  
%  
syms u v  
subplot(132),title('圆环面2')  
ezmesh('(4+2*cos(v))*cos(u)','(4-2*cos(v))*sin(u)',...  
    '2*sin(v)',[0,2*pi,0,2*pi])  
subplot(133),title('圆环面3')  
ezsurf('(4+2*cos(v))*cos(u)','(4-2*cos(v))*sin(u)',...  
    '2*sin(v)',[0,2*pi,0,2*pi])  
%  
shading interp  
box on;  
grid off  
colormap hsv  
%% g.螺旋面  
u=linspace(0,2*pi);  
v=linspace(0,6*pi);  
[U,V]=meshgrid(u,v);  
x=(4+sin(20*U)).*cos(V);  
y=(4+sin(20*U)).*sin(V);  
z=2*V;  
figure(7),title('螺旋面')  
mesh(x,y,z)  
axis equal  
%% h.螺旋卷筒  
clear all  
R=2;  
b=2;  
u=linspace(0,2);  
v=linspace(0,6*pi);  
[U,V]=meshgrid(u,v);  
x=U.*cos(V);  
y=U.*sin(V);  
z=U.*V;  
figure(8),title('螺旋卷筒')  
mesh(x,y,z)  
axis equal  
axis off  
%% i.螺旋管,俗称“弹簧”  
u=linspace(0,2*pi,100);  
v=linspace(0,6*pi,100);  
[U,V]=meshgrid(u,v);  
X=(4+sin(U)).*cos(V);  
Y=(4+sin(U)).*sin(V);  
Z=cos(U)+V;  
figure(9),title('螺旋管')  
mesh(X,Y,Z)  
axis equal  
axis off  
%% j.圆柱上的缎带  
u=linspace(0,2);  
v=linspace(0,6*pi);  
[U,V]=meshgrid(u,v);  
X=2*cos(V);  
Y=2*sin(V);  
Z=U+V;  
figure(10),title('圆柱上的缎带')  
mesh(X,Y,Z)  
axis equal  
axis off  
%% k.纽带  
u=linspace(0,2);  
v=linspace(0,6*pi);  
[U,V]=meshgrid(u,v);  
X=U.*cos(V);  
Y=U.*sin(V);  
Z=5*V;  
figure(10),title('纽带')  
mesh(X,Y,Z)  
axis equal  
axis off  
  评论这张
 
阅读(5220)| 评论(1)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016