clcclear%椎体% r=0:0.02:1;% theta=0:0.02:2*pi;% [r1,theta1]=meshgrid(r,theta);% x=r1.*cos(theta1);% y=r1.*sin(theta1);% z=r1;%球体% fai=0:0.02:pi;% theta=0:0.02:pi*2;% a=1% [fai1,theta1]=meshgrid(fai,theta);% x=a.*sin(fai1).*cos(theta1);% y=a.*sin(fai1).*sin(theta1);% z=a.*cos(fai1);%圆柱体theta=0:0.02:pi;u=0:0.02:5;[theta1,u1]=meshgrid(theta,u);x=3.*sin(2.*theta1);y=6.*sin(theta1).^2;z=u1;surf(x,y,z);axis equal

 

%柱面% x=0:0.02:1;% z=0:0.02:4;% [x1,z1]=meshgrid(x,z);% x=x1;% y=x1.^2;% z=z1;% % hold on% quiver3(1/2, 1/4, 1/2,1/10,1/10,0/10,'g') % quiver3(1/2, 1/4, 1/2,0/10,0/10,1/7,'r') % quiver3(1/2, 1/4, 1/2,1/10,-1/10,0,'y') %抛物面r=0:0.02:2;theta=0:0.02:2*pi;[r1,theta1]=meshgrid(r,theta);x=r1.*cos(theta1);y=r1.*sin(theta1);z=r1.^2;

 

 

空间椭圆曲线参数方程

x+y+z=1 切 x^2+y^2=3形成的空间椭圆

x=3 cos(theta)

y=3 sin(theta)

z=1-x-y

椭圆切面参数方程

theta=0:0.03:2*pi;r=0:0.05:3;[r1,theta1]=meshgrid(r,theta);x=r1.*cos(theta1);y=r1.*sin(theta1);z=1-x-y;surf(x,y,z)

 

 

%轮胎v=0:0.03:2*pi;u=0:0.03:2*pi;[v1,u1]=meshgrid(v,u);R=4;r=1;x=(R+r.*cos(u1)).*cos(v1);y=(R+r.*cos(u1)).*sin(v1);z=r.*sin(u1) ;surf(x,y,z)hold on quiver3(0, 0, 0,1,0,0,'g')  quiver3(0, 0, 0,0,1,0,'g')  quiver3(0, 0, 0,0,0,1,'g')

 

theta=0:0.03:2*pi;

u=0:0.03:2;

[theta1,u1]=meshgrid(theta,u);

a=4;

b=5;

c=7;

x=a.*sinh(u1).*sin(theta1);

y=b.*sinh(u1).*cos(theta1);

z=c.*cosh(u1);

mesh(x,y,z);

hold on

z=-c.*cosh(u1);

mesh(x,y,z);

grid on

axis equal

 

 

sinh->cosh  单双叶互换