数学实验与 Matlab for k=1: 9 subplot(3, 3, k), plot(x, y), axis([0, 2*pi, -1, ID) 212多元函数的可视化与空间解析几何(三维图形) 本节通过高等数学的几个例子观察 Matlab的三维绘图功能和技巧 1.绘制二元函数 ◆观察:绘制κ=∫(x,y)=(1-x)h(x-y)的图象,作定义域的栽剪 ◆(1)观察 meshgrid指令的效果。 【a=-0.98b=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、 mesha、surf ◆(2)做團数的定义域裁剪,瘀察上述三维绘图指令的效果。 程序zxy24m clear clf a=-1;b=1;c=-15;d=15;n=20;epsl=0.01 x=linspace(a, b, n); y=linspace(c, d, n) [X,Y=meshgrid(x,y) %计算函数值z,并作 定义域裁剪 if(1-X(1j))<epsl|x(i)-Y(1,j)<epsl%i语句这样用 (1,j)=NaN; %作定义域 裁剪,定义域以外的函数值为NaN数学实验与 Matlab 10 10 for k =1:9 y=sin(k* x); subplot(3,3,k),plot(x,y),axis([0,2* pi,-1,1]) end 】 2.1.2 多元函数的可视化与空间解析几何（三维图形） 本节通过高等数学的几个例子观察 Matlab 的三维绘图功能和技巧。 1. 绘制二元函数 ◆ 观察：绘制 ( , ) (1 ) ln( ) 2 1 z = f x y = − x x − y − 的图象，作定义域的裁剪。 ◆ （1）观察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。 ◆ （2）做函数的定义域裁剪，观察上述三维绘图指令的效果。 程序 zxy2_4.m 【 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