m=50;
n=60;
x=zeros(m,1);%生成長度爲m的一維數組
y=zeros(n,1);%生成長度爲n的一維數組
for i=1:m
x(i)=-4.0+i*8.0/m;%對x進行初始化
end
for j=1:n
y(j)=-4.0+8.0*j/n;%對y進行初始化
end
%利用x,y得到z=f(x,y),這裏z爲一個二維數組
z=zeros(m,n);%z爲二維數組,行列數可以不一樣
for i=1:m
for j=1:n
tempx=x(i);
tempy=y(j);
z(i,j)=tempx^3+tempy^2;%利用z=f(x,y)來計算z的值
end
end
surf(z,'EdgeColor','None');%繪製z的3D圖
shading interp;
生成的效果如下: