exm030302.m

来自「matlab电子书籍」· M 代码 · 共 27 行

clear
x=-5:0.1:5;
y=(-2.5:0.1:2.5)';
N=length(x);
M=length(y);
for ii=1:M
    for jj=1:N
        X0(ii,jj)=x(jj);
        Y0(ii,jj)=y(ii);
        Z0(ii,jj)=sin(abs(x(jj)*y(ii)));
    end
end  

%（2）
[X,Y]=meshgrid(x,y);
Z=sin(abs(X.*Y));  

%（3）
norm(Z-Z0)  

%（4）
surf(X,Y,Z)
xlabel('x')
ylabel('y')
shading interp
view([190,70])