objf.m

本MATLAB例程

%计算适应度函数

function [f,p]=objf(s);

inn=size(s,1);   %有inn个个体
bn=size(s,2);    %个体长度为bn

for i=1:inn
   x=n2to10(s(i,:));  %讲二进制转换为十进制
   xx=-500+x*1000/(power(2,22)-1);  %转化为[-500,500]区间的实数
   f(i)=ft(xx);  %计算函数值，即适应度
end
f=f';

%计算选择概率
ff=(f-min(f))/250000;
fsum=sum(ff.*ff);
ps=ff.*ff/fsum;

%计算累积概率
p(1)=ps(1);
for i=2:inn
   p(i)=p(i-1)+ps(i);
end
p=p';