matlab二维三维画图(2)

subplot(2,2,3) syms x3 y3

z3=(x3.^2+y3.^2)^2-2*(x3.^2-y3.^2) ezplot(z3,[-2,2,-1,1]);

title('3.(x^2+y^2)^2=2(x^2-y^2)') subplot(2,2,4) syms x4 y4

z4=(x4.^2+y4.^2)^2-4.*x4.*y4 ezplot(z4,[-1.5,1.5,-1.5,1.5]); title('4.(x^2+y^2)^2=4xy')

结果：

1.x(2/3)+y(2/3)=2(2/3)21.51y10.5000.511.52x13.(x2+y2)2=2(x2-y2)10.53y0-0.5-1-2-1012x 3

8.画出下列曲线，并分别进行标注。

? ?3t?x??1?t21) ?2??y?3t1?t2解：syms t

x=3.*t/(1+t.^2); y=3.*t.^2/(1+t.^2) ezplot(x,y,[-20,20])

结果：

2.x2+y2+2x=2sqrt(x2+y2)220-2-4-20x24.(x2+y2)2=4xy140-1-101x4yy

x = (3 t)/(t2 + 1), y = (3 t2)/(t2 + 1)2.521.5y10.5-1.5-1-0.50x0.511.5

?x?ln(1?t2)2) ?

?y?t?arctant解：syms t

x=log(1+t.^2); y=t-atan(t)

ezplot(x,y,[-2,2])

结果：

x = log(t2 + 1), y = t - atan(t)0.80.60.40.20-0.2-0.4-0.6-0.8-0.200.20.40.60.8x11.21.41.61.8y

2??x?1?t3) ? 3??y?t?t解：syms t

x=1-t.^2; y=t-t.^3

ezplot(x,y,[-1,1])

结果：

x = 1 - t2, y = t - t30.30.20.10-0.1-0.2-0.3y00.10.20.30.40.5x0.60.70.80.9

三维画图

1、在同一平面中的两个窗口分别画出心形线和马鞍面。 x2y2双曲抛物面：-?z

129

要求：1.在图形上加格栅、图例和标注

2.定制坐标

3.以不同角度观察马鞍面 解：subplot(1,2,1)

x=-5:0.1:5;y=-5:0.1:5;z=-5:0.1:5; [x,y,z]=meshgrid(x,y,z);

f=((x.^2+(9/4).*y.^2+z.^2 - 1).^3 - x.^2.*z.^3 - (9/80).*y.^2.*z.^3 ); p=patch(isosurface(x,y,z,f,0)); isonormals(x,y,z,f, p)

set(p,'FaceColor','red','EdgeColor','none'); daspect([1 1 1]); view(3)

camlight; lighting phong grid on

title('立体心') subplot(1,2,2) x=-5:0.1:5; y=-5:0.1:5;

[X,Y]=meshgrid(x,y); Z=(X.^2)./12-(Y.^2)./9; surf(X,Y,Z) shading flat

title('马鞍面')

结果：

马鞍面立体心32210-110-1-2-350-5-502010-1-25

2222222. 以不同的视角观察球面x?y?z?r和圆柱面x?y?rx所围区域。

解：r=1;

[x,y,z]=sphere(50); X=r*x;Y=r*y;Z=r*z; surf(X,Y,Z); hold on;

f=@(x,y)x.^2+y.^2-r*x;

ezsurf(f,[-1.5,2.5],[-2,1.5]); axis equal;

结果：

x2+y2-r x76543210-110-1-2y-101x2

3.画出下列图形

?x2?y2?z2?9（1）?

?y?x解：clear

t=-sqrt(6)/2:0.1:sqrt(6)/2 x=t; y=t;

z=sqrt(9-2*x.^2) plot3(x,y,z,'r') hold on

plot3(x,-y,z,'bo') axis square

结果：

3.232.82.621-2-1012-2-10

x?1yz?3xy?2z??与L2：??（2）L1： 1?412?2?1

解：clear

t=-5:0.1:5; x=t;

y=-4.*t+4; z=t-4;