%計算適應度函數

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

inn=size(s,1);%有inn個個體
N=50; %個體中有N個變量
bn=size(s,2)/N;    %個體中每個變量的長度為bn

for i=1:inn
   for j=1:N
     x=n2to10(s(i,(bn*(j-1)+1):(bn*j)));  %講二進制轉換為十進制
     xx=-500.0+x*1000/(power(2,bn)-1);  %轉化為[-500,500]區間的實數
     X(j)=xx; %得到N個十進制的數,帶入適應度函數求適應度
   end
  f(i)=7.03e-008-ft(X);  %計算函數值，即適應度(該適應度函數有待改進)
end
f=f';

%計算選擇概率
fsum=sum(f.*f);
ps=f.*f/fsum;

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