4_15.m

来自「这是一本介绍MATLAB的教程」· M 代码 · 共 18 行

x=-2:0.1:2;
y=x';
X=ones(size(y))*x;
Y=y*ones(size(x));
Z=X.^2.*exp(-(X.^2+Y.^2))
mesh(X,Y,Z)
surf(X,Y,Z)
colormap(gray)
Z1=0.05*X-0.05*Y+0.1
hold on
mesh(X,Y,Z1)
k=find(abs(Z-Z1)<0.01)
plot3(X(k),Y(k),Z1(k),'rx')
hold off