Matlab计算与仿真 6.绘制二元函数 绘制二=f(x,y)=(1-x)hn(x-y)的图象,作定义域的裁剪。 观察 meshgrid指令的效果 a=-0.98:b=0.98:c=-1;d=1:n=10 x=linspace(a, b,n): y=linspace(c, d, n) [X, y]=meshgrid(x, y) plot(X, y,+) 三维绘图指令mesh、 mesh、surf 做函数的定义域裁剪,观察上述三维绘图指令的效果。 clear clf a=-:b1:c=-15d=15m20eps1=001 xElinspace(a, b, n): y=linspace(c, d, n): [X, y]meshgrid(x, y): fori=1:n%计算函数值z,并作定义域裁剪 or J= f(1-×(j)eps1×(Jj)-Y()eps1%f语句这样用 z(j)=№aN%作定义域裁剪,定义域以外的函数值为NaN else z(小)=1000*sq+(1-X(j)^-1*og(x(小)Y(j) end Zz=-20*ones(1, n): plot3(×X,zz) grid off: hold on %画定义域的边界线 mesh(Xy,z)%绘图,读者可用 mesh,surf, mesh在此替换之 xlabel(x') ylabel(y), zlabel(z). box on%把三维图形封闭在箱体里 运行上述程序,有关向量存储在工作空间中,在此基础上,观察上述等值线绘制指令 的运行效果 [cs, h]contour(X,y, z, 15) clabel(s, h, labelspacing, 244)Matlab计算与仿真 6. 绘制二元函数 绘制 ( , ) (1 ) ln( ) 2 1 z = f x y = − x x − y − 的图象，作定义域的裁剪。 - 观察meshgrid指令的效果。 a=-0.98;b=0.98;c=-1;d=1;n=10; x=linspace(a,b, n); y=linspace(c,d, n); [X,Y]=meshgrid(x, y); plot(X,Y,'+') 三维绘图指令mesh、meshc、surf。 - 做函数的定义域裁剪，观察上述三维绘图指令的效果。 clear,clf, a=-1;b=1;c=-15;d=15;n=20;eps1=0.01; x=linspace(a,b,n);y=linspace(c,d,n); [X,Y]=meshgrid(x,y); for i=1:n %计算函数值 z ，并作定义域裁剪 for j=1:n if (1-X(i,j))<eps1|X(i,j)-Y(i,j)<eps1 %if 语句这样用 z(i,j)=NaN; %作定义域裁剪，定义域以外的函数值为 NaN else z(i,j)=1000*sqrt(1-X(i,j))^-1.*log(X(i,j)-Y(i,j)); end end end zz=-20*ones(1,n); plot3(x,x,zz); grid off; hold on %画定义域的边界线 mesh(X,Y,z) %绘图，读者可用 meshz, surf, meshc 在此替换之 xlabel('x'),ylabel('y'),zlabel('z'), box on %把三维图形封闭在箱体里 运行上述程序，有关向量存储在工作空间中，在此基础上，观察上述等值线绘制指令 的运行效果。 [cs,h]=contour(X,Y,z,15); clabel(cs,h,'labelspacing',244)