参数方程
>> t=0:0.01:2*pi;
x=t.^3.*sin(3*t).*exp(-t);y=t.^3.*cos(3*t).*exp(-t);z=t.^2;plot3(x,y,z),grid
>> [x,y]=meshgrid(-3:0.1:3,-2:0.1:2);
z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y);
>> mesh(x,y,z)
用点运算的形式来构造Z矩阵
surf(x,y,z)
> [x,y]=meshgrid(-2:.1:2);z=1./(sqrt((1-x).^2+y.^2))+.1./(sqrt((1+x).^2+y.^2));surf(x,y,z),shading flat
空洞为奇点
>> [x1,x2]=meshgrid(-1.5:.1:1.5,-2:.1:2);p=0.5457*exp(-0.75*x2.^2-3.75*x1.^2-1.5*x1).*(x1+x2>1)+0.7575*exp(-x2.^2-6*x1.^2).*((x1+x2>-1)&(x1+x2<=1))+0.5457*exp(-0.75*x2.^2-3.75*x1.^2+1.5*x1).*(x1+x2<=-1);surf(x1,x2,p),xlim([-1.5 1.5]);
这段代码使用meshgrid函数创建了一个二维的网格,然后根据一定的数学公式计算了一个二维的函数p。最后使用surf函数绘制了这个函数的三维图像。 使用surf函数绘制了函数p的三维图像,并通过xlim([-1.5 1.5])设置了x轴的显示范围为-1.5到1.5。
>> syms u v;x=cos(u)+v*cos(u)*cos(u/2);y=sin(u)+v*sin(u)*cos(u/2);z=v*sin(u/2);
ezsurf(x,y,z,[0,2*pi,-0.5,0.5])
[x,y,z]=sphere(50);surf(x,y,z),hold on
x1=0.3*x+0.9;y1=0.3*y-0.8;z1=0.3*z+0.6;surf(x1,y1,z1)