【转】MATLAB之三维作图

转自：http://hi.baidu.com/mynana/blog/item/3c79c1ce6b65cb0393457e7a.html/cmtid/5d844e08dfb79f9e0b7b8227

这几天看到Ras函数（Rastrigin's Function），这个函数对模拟退火、进化计算等算法具有很强的欺骗性，因为它有非常多的局部最小值点和局部最大值点，很容易使算法陷入局部最优，而 不能得到全局最优解。

    Ras函数为：


Matlab的Help里面有这个函数的三维图像，但是没有解释是怎样画出来的。经过请教同学，得出两种方法，现记录在此。

第一：
[x,y]=meshgrid([0:0.1:4]);
z=20+x.*x+y.*y-10*(cos(2*pi*x)+cos(2*pi*y)); 
surfc(x,y,z)
shading interp

第二：
f=@(x,y) 20+x.*x+y.*y-10*(cos(2*pi*x)+cos(2*pi*y));
ezsurf(f,[-2 2])
shading interp

三维曲线的画法

三维空间曲线要用到plot3函数，这个和plot类似。plot3函数有三个参数，x，y和z轴，比如下面的例子：

>> T = -2:0.01:2;

>> plot3(cos(2*pi*T),sin(2*pi*T),T)

如果安装了Symbolic Math Toolbox的话也可以用下面ezlpot3函数的方法：

>> ezplot3(‘cos(2*pi*T)’,’sin(2*pi*T)’,'T’,[-2 2])

三维曲面的画法

有mesh何surf两种命令来画三维曲面，它们使用的场合不同。前者是当z轴是x和y的显式函数时，后者是x,y,z中某个为其他2个的函数。

mesh函数

>> [X Y]=meshgrid(-2:.1:2, -2:.1:2);

>> Z = X.^2 – Y.^2;

>> mesh(X, Y, Z)

同理用Symbolic Math Toolbox可以直接执行

>> ezmesh(‘X.^2 – Y.^2′, [-2 2], [-2 2])

surf函数

在函数不能表示成z = f(x, y)时，需要用surf函数。比如x²+y²+z²=1.

先需要用柱面坐标或者球坐标来表示。这里用柱面坐标表示为 r²+z²=1

x = sqrt(1-z²)cosθ, x = sqrt(1-z²)sinθ;

执行matlab指令：

>> [theta, Z] = meshgrid((0:0.1:2)*pi, (-1:0.1:1));

>> X =sqrt(1 – Z.^2).*cos(theta);

>> Y =sqrt(1 – Z.^2).*sin(theta);

>> surf(X, Y, Z); axis square

同理用Symbolic Math Toolbox可以直接执行

>> ezsurf(’sqrt(1-s^2)*cos(t)’,’sqrt(1-s^2)*sin(t)’, ’s’, [-1, 1, 0, 2*pi]); axis equa

改天贴图上来