gabp.m

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

clear all
tic,
[P,T,R,S1,S2,S]=bpinit;
bounds=ones(S,1)*[-10 10];
num=60;
pop=initializega(num,bounds,'fitness');

gen=500;
[x endPop bPop trace]=ga(bounds,'fitness',[],pop,[1e-6 1 1],'maxGenTerm',gen,...'normGeomSelect',[0.009],['arithXover'],[2],'nonUnifMutation',[2 gen 3]);
for i=1:S
    x(i)=endPop(1,i);
end;
sum=0;

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

for i=1:100
    x1=W1*P(i)+B1;
    A1=1./(1+exp(-x1));
    x2=W2*A1+B2;
    A2=1./(1+exp(-x2));
    
    YY(i)=A2;
end

i=1:1:100
plot(i,YY(i),'r+',i,T(i),'b-');
toc