转自: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函数。比如x2+y2+z2=1.
先需要用柱面坐标或者球坐标来表示。这里用柱面坐标表示为 r2+z2=1
x = sqrt(1-z2)cosθ, x = sqrt(1-z2)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
改天贴图上来