matlab第9-10讲byXR浅析.ppt

例5.14 已知6x30，15y36，求不定方程2x+5y=126的整数解。 程序如下： x=7:29; y=16:35; [x,y]=meshgrid(x,y); %在[5,29]×[14,35]区域生成网格坐标 z=2*x+5*y; k=find(z==126); %找出解的位置 x(k),y(k) %输出对应位置的x,y即方程的解 2. 绘制三维曲面的函数 surf函数和mesh函数的调用格式为： mesh(x,y,z,c)绘制三维网格图 surf(x,y,z,c)绘制曲面图，各线条之间的补面用颜色填充。 例5.15 用三维曲面图表现 函数z=sin(y)cos(x)。 程序1： x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos(x); mesh(x,y,z);xlabel(x-axis),ylabel(y-axis),zlabel(z-axis);title(mesh); 程序2： x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos(x); surf(x,y,z);xlabel(x-axis),ylabel(y-axis),zlabel(z-axis);title(surf); 程序3： x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos(x); plot3(x,y,z);xlabel(x-axis),ylabel(y-axis),zlabel(z-axis);title(plot3-1);grid; 例5.16 绘制两个直径相等的圆管的相交图形。 程序如下: %两个等直径圆管的交线 m=30; z=1.2*(0:m)/m; r=ones(size(z)); theta=(0:m)/m*2*pi; x1=r*cos(theta);y1=r*sin(theta); %生成第一个圆管的坐标矩阵 z1=z*ones(1,m+1); x=(-m:2:m)/m; x2=x*ones(1,m+1);y2=r*cos(theta); %生成第二个圆管的坐标矩阵 z2=r*sin(theta); surf(x1,y1,z1); %绘制竖立的圆管 axis equal,axis off hold on surf(x2,y2,z2); %绘制平放的圆管 axis equal,axis off title(两个等直径圆管的交线); hold off 例5.17 分析由函数z=x2-2y2构成的曲面形状及与平面z=a的交线。 程序如下： [x,y]=meshgrid(-10:0.2:10); z1=(x.^2-2*y.^2)+eps; %第1个曲面 a=input(a=?); z2=a*ones(size(x)); %第2个曲面 subplot(1,2,1);mesh(x,y,z1);hold on;mesh(x,y,z2); %分别画出两个曲面 v=[-10,10,-10,10,-100,100];axis(v);grid; %第1子图的坐标设置 hold off; r0=abs(z1-z2)=1; %求两曲面z坐标差小于1的点 xx=r0.*x; yy=r0.*y; zz=r0.*z2; %求这些点上的x,y,z坐标，即交线坐标 subplot(1,2,2); plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),*); %在第2子图画出交线 axis(v);grid; %第2子图的坐标设置 例5.18 在xy平面内选择区域[-8,8]×[-8,8]，绘制函数的4种三维曲面图。 程序如下： [x,y]=meshgrid(-8:0.5:8); z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps); subplot(2,2,1); meshc(x,y,z);%等高线 title(meshc(x,y,z)) subplot(2,2,2); meshz(x,y,z);%带底座 title(meshz(x,y,z)) subplot(2,2,3); surfc(x,y,z);%等