参考http://huzhyi21.blog.163.com/blog/static/1007396201061052214302/
x=linspace(-2, 2, 25); % 在x轴上取25点
y=linspace(-2, 2, 25); % 在y轴上取25点
[xx,yy]=meshgrid(x, y); % xx和yy都是21x21的矩阵
zz=xx.*exp(-xx.^2-yy.^2); % 计算函数值，zz也是21x21的矩阵
surf(xx, yy, zz); % 画出立体曲面图
colorbar;
Matlab中绘制等高线：contour
>> contour(xx,yy,zz)
>> colorbar
网状图与等高线同时绘制：surfc
>> surfc(xx,yy,zz)
>> colorbar
以上代码在octave-3.6.4-gcc4.6.2中可运行，运行结果：