matlab绘制空间曲面方法简介

matlab 曲面绘制%plot3(x,y,z)用来绘制3维曲线图，而不能绘制曲面图！就是把所有的 %(x,y,z)点连接在一起。而画曲面图必须用surf和mesh函数，而这两个 %函数都需要知道对应x，y向量交叉点内所有点处的z值，所以得用 %x=-200:0.3:200;y=-200:0.1:200;%X,Y=meshgrid(x,y);%Z=X.*Y/sqrt(X.*X+Y.*Y)%获得z值，再用mesh(X, Y, Z);绘曲面图，下面为代码 x=-200:0.3:200;y=-200:0.1:200;X,Y=meshgrid(x,y);Z=X.*Y./sqrt(X.*X+Y.*Y);mesh(X,Y,Z)#matlab画两曲面的交线悬赏分：100 - 解决时间：2009-10-4 19:37上半球面z2+y.2+x2=l(z=0)与平面x+y+z=1.5的交线该怎么画出来？ X,Y,Z=sphere(50);Z(Z0)=-(Z(Z0);mesh(X,Y,Z)%画上半球面hold on;ezmesh(1.5-x-y,-1 1)%画平面x,y,z=meshgrid(linspace(-1,1);contourslice(x,y,z,x+y+z-1.5,X,Y,Z,0 0)%画交线axis equalview(135,20)# #两直线相交function X,Y=pll(X1,Y1,X2,Y2) % 直线相交求交点A1=Y1(1)-Y1(2);B1=X1(2)-X1(1);C1=Y1(2)*X1(1)-Y1(1)*X1(2);A2=Y2(1)-Y2(2);B2=X2(2)-X2(1); C2=Y2(2)*X2(1)-Y2(1)*X2(2);D=det(A1,B1;A2,B2);X=det(-C1 B1;-C2 B2)/D;Y=det(A1 -C1;A2,-C2)/D;调用格式：x1=1 5;y1=1 5;x2=1 5;y2=5,1; x,y=pll(x1,y1,x2,y2);plot(x1,y1,r);hold onplot(x2,y2,b); plot(x,y,ko);%直线与多条直线相交xi=1 2 3 4 5;yi=2 6 3 6 1; plot(xi,yi);hold onx1=1 5;y1=4 5;line(x1,y1); x=zeros(size(xi);y=x;for i=1:5-1x2=xi(i i+1);y2=yi(i i+1); x,y=pll(x1,y1,x2,y2);plot(x,y,ro)end%直线与曲线相交x=-8:0.1:8;y=x;X,Y=meshgrid(x,y);R=sqr t(X.2+Y2)+eps;Z二sin(R)./R;contour(Z,3);hold onc=contour(Z,3);x=0 360;y=0 400;y=(y(2)-y(1)/(x(2)-x(1)*(x-x(1)+y(1);z=0 0; line(x,y,z);c=c;X=c(:,1);Y=c(:,2); r0=abs(Y-(y(2)-y(1)/(x(2)-x(1)*(X-x(1)+y(1)=.93;zz=0;yy=r0.*Y;xx=r0.*X; plot(xx(r0=0),yy(r0=0),r)%曲线与曲线相交 x=0:pi/400:2*pi;，x=x;y1=sin(pi*x);y2=cos(pi*x);plot(x,y1,x,y2);hold on r0=abs(y2-sin(pi*x)=0.02;yy=r0.*y1;xx=r0.*x;plot(xx(r0=0),yy(r0=0),r.)直线与曲面相交x=-8:0.3:8;y=x;X,Y=meshgrid(x,y);Z=X2+Y.2;mesh(X,Y,Z);hold onx=-10 10;y二-10 3;z=30 35;line(x,y,z);r0=(abs(Y-y(1)-(y(2)-y(1)/(x(2)-x(1)*(X-x(1)=0.45)&. (abs(Z-z(1)-(z(2)-z(1)/(x(2)-x(1)*(X-x(1)0.45)&.(abs(Y-y(1)-(y(2)-y(1)/(z(2)-z(1)*(Z-z(1)=0.45); zz=r0.*Z;yy=r0.*Y;xx=r0.*X;plot3(xx(r0=0),yy(r0=0),zz(r0=0),r*)平面与曲面相交x=-8:0.1:8;y=x;X,Y=meshgrid(x,y);Z1=2*ones(size(X);Z2=X.2-Y2;mesh(X,Y,Z1);hold onmesh(X,Y,Z2);r0=(abs(Z1-Z2)=.65); zz=r0.*Z1;yy=r0.*Y;xx=r0.*X;plot3(xx(r0=0),yy(r0=0),zz(r0=0),k*)clcdisp(观察曲面后，按任意键画交线)；pauseclf plot3(xx(r0=0),yy(r0=0),zz(r0=0),k*);%曲面与多个截平面相交y=-10:0.5:10;z=y;Z,Y=meshgrid(z,y);X=Z;X1=0*ones(size(Z);X2=3*ones(size(Z);X3=-3*ones(size(Z);Z4=(X2-Y.2)/10;mesh(X1,Y,Z);hold onmesh(X2,Y,Z)mesh(X3,Y,Z);mesh(X,Y,Z4);r1=(abs(X1-X)0.05);r2=(abs(X2-X)0.05);r3=(abs(X3-X)0.05);zz1=r1.*Z4;yy1=r1.*Y;xx1=r1.*X;zz2=r2.*Z4;yy2=r1.*Y;xx2=r1.*X;zz3=r3.*Z4;yy3=r1.*Y;xx3=r1.*X;plot3(xx1(r1=0),yy1(r1=0),zz1(r1=0),k*); plot3(xx2(r2=0),yy2(r2=0),zz2(r2=0),k*); plot3(xx3(r3=0),yy3(r3=0),zz3(r3=0),k*); colormap(hsv)clc;disp(观察曲面后，按任意键画交线)；hold off平面与曲面相交y=-8:0.4:8;z=y;Z,Y=meshgrid(z,y);X=Z;X1=zeros(size(Z);Z2=zeros(size(Z);Z3=(X2-Y.2)/10;mesh(X1,Y,Z);hold onmesh(X,Y,Z2);mesh(X,Y,Z3);r1=(abs(X1-X)0.05);r2=(abs(Z3-Z2)0.05);r3=(abs(X1-X)0.05)&(abs(Z-Z2)=0.05);zz1=r1.*Z3;yy1=r1.*Y;xx1=r1.*X;zz2=r2.*Z3;yy2=r2.*Y;xx2=r2.*X;zz3=r3.*Z;yy1=r3.*Y;xx1=r3.*X1;plot3(xx1(r1=0),yy1(r1=0),zz1(r1=0),k*); plot3(xx2(r2=0),yy2(r2=0),zz2(r2=0),k*); plot3(xx3(r3=0),yy3(r3=0),zz3(r3=0),k*);colormap(hsv);引自：