matlab常用的幾個(gè)適應(yīng)度評(píng)價(jià)函數(shù)

/niuyongjie/article/details/1619496

粒子群算法(6)-----幾個(gè)適應(yīng)度評(píng)價(jià)函數(shù)下面給出幾個(gè)適應(yīng)度評(píng)價(jià)函數(shù)，并給出圖形表示

頭幾天機(jī)子種了病毒，重新安裝了系統(tǒng)，不小心把程序全部格式化了，痛哭?。?！沒辦法，好多程序不見了，現(xiàn)在把這幾個(gè)典型的函數(shù)重新編寫了，把他們給出來，就算粒子群算法的一個(gè)結(jié)束吧！痛恨病毒?。。?！第一個(gè)函數(shù)：Griewank函數(shù)，圖形如下所示：適應(yīng)度函數(shù)如下：（為了求最大值，我去了所有函數(shù)值的相反數(shù)）

function

y

=

Griewank(x)

%

Griewan函數(shù)

%

輸入x,給出相應(yīng)的y值,在x

=

(

0

,

0

,…,

0

)處有全局極小點(diǎn)0.

%

編制人：

%

編制日期：

[row,col]

=

size(x);

if

row

>

1

error(

'

輸入的參數(shù)錯(cuò)誤

'

);

end

y1

=

1

/

4000

*

sum(x.

^

2

);

y2

=

1

;

for

h

=

1

:col

y2

=

y2

*

cos(x(h)

/

sqrt(h));

end

y

=

y1

-

y2

+

1

;

y

=-

y;繪制函數(shù)圖像的代碼如下：function

DrawGriewank()

%

繪制Griewank函數(shù)圖形

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)]);

end

end

surf(X,Y,z);

shading

interp第二個(gè)函數(shù)：Rastrigin函數(shù)，圖形如下所示：適應(yīng)度函數(shù)如下：（為了求最大值，我去了所有函數(shù)值的相反數(shù)）

function

y

=

Rastrigin(x)

%

Rastrigin函數(shù)

%

輸入x,給出相應(yīng)的y值,在x

=

(

0

,

0

,…,

0

)處有全局極小點(diǎn)0.

%

編制人：

%

編制日期：

[row,col]

=

size(x);

if

row

>

1

error(

'

輸入的參數(shù)錯(cuò)誤

'

);

end

y

=

sum(x.

^

2

-

10

*

cos(

2

*

pi

*

x)

+

10

);

y

=-

y;繪制函數(shù)圖像的代碼如下：function

DrawRastrigin()

%

繪制Rastrigin函數(shù)圖形

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)]);

end

end

surf(X,Y,z);

shading

interp第三個(gè)函數(shù)Schaffer函數(shù)，圖形如下所示：函數(shù)的代碼如下，因?yàn)樵摵瘮?shù)在（0,...,0）處有最大值1，因此不需要取相反數(shù)。[c-sharp]

\o"viewplain"viewplain\o"copy"copyfunction

result=Schaffer(x1)

%Schaffer

函數(shù)

%輸入x,給出相應(yīng)的y值,在x=(0,0,…,0)

處有全局極大點(diǎn)1.

%編制人：

%編制日期：

[row,col]=size(x1);

if

row>1

error('輸入的參數(shù)錯(cuò)誤');

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;

繪制函數(shù)代碼圖形的代碼如下：[c-sharp]

\o"viewplain"viewplain\o"copy"copyfunction

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

第四個(gè)函數(shù)：Ackley函數(shù)，函數(shù)圖形如下：函數(shù)的代碼如下，因?yàn)樵摵瘮?shù)在（0,...,0）處有最小值0，因此需要取相反數(shù)[c-sharp]

\o"viewplain"viewplain\o"copy"copyfunction

result=Ackley(x)

%Ackley

函數(shù)

%輸入x,給出相應(yīng)的y值,在x=(0,0,…,0)

處有全局極小點(diǎn)0,為得到最大值，返回值取相反數(shù)

%編制人：

%編制日期：

[row,col]=size(x);

if

row>1

error('輸入的參數(shù)錯(cuò)誤');

end

result=-20*exp(-0.2*sqrt((1/col)*(sum(x.^2))))-exp((1/col)*sum(cos(2*pi.*x)))+exp(1)+20;

result=-result;

繪制函數(shù)代碼圖形的代碼如下：[c-sharp]

\o"viewplain"viewplain\o"copy"copyfunction

DrawAckley()

%繪制Ackley函數(shù)圖形

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)=Ackley([X(h,l),Y(h,l)]);

end

end

surf(X,Y,z);

shading

interp

第五個(gè)函數(shù)是：Rosenbrock函數(shù)，該函數(shù)在（1,...,1）處有最小值0，為了得到最大值，取函數(shù)值的相反數(shù)。函數(shù)圖形如下所示

函數(shù)的代碼：

[c-sharp]

\o"viewplain"viewplain\o"copy"copyfunction

result=Rosenbrock(x)

%Rosenbrock

函數(shù)

%輸入x,給出相應(yīng)的y值,在x=(1,1,…,1)

處有全局極小點(diǎn)0,為得到最大值，返回值取相反數(shù)

%編制人：

%編制日期：

[row,col]=size(x);

if

row>1

error('輸入的參數(shù)錯(cuò)誤');

end

result=100*(x(1,2)-x(1,1)^2)^2+(x(1,1)-1)^2;

result=-result;

繪制函數(shù)圖形的代碼如下：[c-sharp]

\o"viewplain"viewplain\o"copy"copyfunction

DrawRosenbrock()

%繪制Rosenbrock函數(shù)圖形,