f.m

用matlab开发的用遗传算法改进的BP算法程序

function [eval]=f(sol)
      numv=size(sol,2);
      x=sol(1:numv);
      [P,T,R,S1,S2,S]=bpinit;
      %get the value of fitness
      
for i=1:S1,
    for k=1:R,
        W1(i,k)=x(R*(i-1)+k);
    end
end

for i=1:S2,
    for k=1:S1,
        W2(i,k)=x(S1*(i-1)+k+R*S1);
    end
end

for i=1:S1,
    B1(i,1)=x((R*S1+S1*S2)+i);
end

for i=1:S2,
    B2(i,1)=x((R*S1+S1*S2+S1)+i);
end

sum=0;
SE=0;
for i=1:100
    x1=W1*P(i)+B1;
    A1=1./(1+exp(-x1));
    x2=W2*A1+B2;
    A2=1./(1+exp(-x2));
    
SE=sumpqr(T(i)-A2);
sum=sum+SE;
end
eval=10/sum;