%%[Matalb code]
x=1:0.1:2;
y=1:0.1:2;
[X,Y]=meshgrid(x,y);
z=f(X,Y)
mesh(z)
假设z=exp(-x.^2-y.^2),则绘图代码如下:
f=@(x,y)exp(-x.^2-y.^2)
x=1:0.1:2;
y=1:0.1:2;
[X,Y]=meshgrid(x,y);
z=f(X,Y);
mesh(z)
axis('tight')
图片绘制为:
追问定义域条件中有一个x<y,应该如何表示呢?
追答不好意思。刚没看清楚。只需要在绘图前添加z(X>=Y)=nan即可。
如果想让图片平滑,可以点的取值密集一些。
f=@(x,y)exp(-x.^2-y.^2)
x=1:0.01:2;
y=1:0.01:2;
[X,Y]=meshgrid(x,y);
z=f(X,Y);
z(X>=Y)=nan;
mesh(X,Y,z)