粒子群算法-----几个适应度评价函数.docVIP

粒子群算法-----几个适应度评价函数 ? 第一个函数：Griewank函数，图形如下所示： 适应度函数如下：（为了求最大值，我去了所有函数值的相反数） ? function?y = Griewank(x)% Griewan函数% 输入x,给出相应的y值,在x = ( 0 , 0 ,…, 0 )处有全局极小点0.% 编制人：% 编制日期：[row,col] = size(x);if ?row 1 ????error( 输入的参数错误 );endy1 = 1 / 4000 * sum(x. ^ 2 );y2 = 1 ;for ?h = 1 :col????y2 = y2 * cos(x(h) / sqrt(h));endy = y1 - y2 + 1 ;y =- y; 绘制函数图像的代码如下： function?DrawGriewank()% 绘制Griewank函数图形x = [ - 8 : 0.1 : 8 ];y = x;[X,Y] = meshgrid(x,y);[row,col] = size(X);for ?l = 1 :col???? for ?h = 1 :row????????z(h,l) = Griewank([X(h,l),Y(h,l)]);????endendsurf(X,Y,z);shading?interp 第二个函数：Rastrigin函数，图形如下所示： 适应度函数如下：（为了求最大值，我去了所有函数值的相反数） ? function?y = Rastrigin(x)% Rastrigin函数% 输入x,给出相应的y值,在x = ( 0 , 0 ,…, 0 )处有全局极小点0.% 编制人：% 编制日期：[row,col] = size(x);if ?row 1 ????error( 输入的参数错误 );endy = sum(x. ^ 2 - 10 * cos( 2 * pi * x) + 10 );y =- y; 绘制函数图像的代码如下： function?DrawRastrigin()% 绘制Rastrigin函数图形x = [ - 5 : 0.05 : 5 ];y = x;[X,Y] = meshgrid(x,y);[row,col] = size(X);for ?l = 1 :col???? for ?h = 1 :row????????z(h,l) = Rastrigin([X(h,l),Y(h,l)]);????endendsurf(X,Y,z);shading?interp 第三个函数Schaffer函数，图形如下所示： 函数的代码如下，因为该函数在（0,...,0）处有最大值1，因此不需要取相反数。 [c-sharp] \o view plain view plain \o copy copy \o print print \o ? ? function?result=Schaffer(x1)?? %Schaffer?函数?? %输入x,给出相应的y值,在x=(0,0,…,0)?处有全局极大点1.?? %编制人：?? %编制日期：?? [row,col]=size(x1);?? if?row1?? ????error(输入的参数错误);?? end?? x=x1(1,1);?? y=x1(1,2);?? temp=x^2+y^2;?? result=0.5-(sin(sqrt(temp))^2-0.5)/(1+0.001*temp)^2;?? ?? 绘制函数代码图形的代码如下： [c-sharp] \o view plain view plain \o copy copy \o print print \o ? ? function?DrawSchaffer()?? x=[-5:0.05:5];?? y=x;?? [X,Y]=meshgrid(x,y);?? [row,col]=size(X);?? for?l=1:col?? for?h=1:row?? z(h,l)=Schaffer([X(h,l),Y(h,l)]);?? end?? end?? surf(X,Y,z);?? shading?interp?? ?? 第四个函数：Ackley函数，函数图形如下： 函数的代码如下，因为该函数在（0,...,0）处有最小值0，因此需要取相反数 [c-sharp] \o view plain view plain \o copy copy \o print print \o ? ? function?result=A