资源描述：
《matlab绘制空间曲面方法简介》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、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./sqr
2、t(X.*X+Y.*Y);mesh(X,Y,Z)###################################################################################matlab画两曲面的交线悬赏分：100-解决时间：2009-10-419:37上半球面z.^2+y.^2+x.^2=1(z>=0)与平面x+y+z=1.5的交线该怎么画出来？---------------------------[X,Y,Z]=sphere(50);Z(Z<0)=-(Z(Z<0));mesh(X,Y,Z)%画上半球面holdon;ez
3、mesh('1.5-x-y',[-11])%画平面[x,y,z]=meshgrid(linspace(-1,1));contourslice(x,y,z,x+y+z-1.5,X,Y,Z,[00])%画交线axisequalview(135,20)###################################################################################两直线相交function[X,Y]=pll(X1,Y1,X2,Y2)%直线相交求交点A1=Y1(1)-Y1(2);B1=X1(2)-X1(1);C1=Y1
4、(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([-C1B1;-C2B2])/D;Y=det([A1-C1;A2,-C2])/D;调用格式：x1=[15];y1=[15];x2=[15];y2=[5,1];[x,y]=pll(x1,y1,x2,y2);plot(x1,y1,'r');holdonplot(x2,y2,'b');plot(x,y,'ko');%直线与多条直线相交xi=[12345];yi
5、=[26361];plot(xi,yi);holdonx1=[15];y1=[45];line(x1,y1);x=zeros(size(xi));y=x;fori=1:5-1x2=xi([ii+1]);y2=yi([ii+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=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;contour(Z,3);holdonc=contour(Z,3);x=[0360];y=[0400];
6、y=(y(2)-y(1))/(x(2)-x(1))*(x-x(1))+y(1);z=[00];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);holdonr0=abs(y2-sin(pi*x))
7、<=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=X.^2+Y.^2;mesh(X,Y,Z);holdonx=[-1010];y=[-103];z=[3035];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)&...(a
8、bs(Y-y(1)-(y
