[½»Á÷] ÇóÖúbpÉñ¾­ÍøÂçÒÅ´«Ëã·¨³ö´í¡£ ÒÑÓÐ1È˲ÎÓë

% Çå¿Õ»·¾³±äÁ¿

clc

clear


%% ÍøÂç½á¹¹½¨Á¢
%½Úµã¸öÊý

inputnum=1;

hiddennum=5;

outputnum=1;



%ѵÁ·Êý¾ÝºÍÔ¤²âÊý¾Ý

input_train=[1 2 3 4 5 6 7 8 9 10];

input_test=[21 27 29 31 34 36 39 41 49 50];

output_train=[11 12 13 14 15];

output_test=[54 56 59 61 66];



%Ñ¡Á¬Ñù±¾ÊäÈëÊä³öÊý¾Ý¹éÒ»»¯

[inputn,inputps]=mapminmax(input_train);

[outputn,outputps]=mapminmax(output_train);



%¹¹½¨ÍøÂç

net=newff(inputn,outputn,hiddennum);



%% ÒÅ´«Ëã·¨²ÎÊý³õʼ»¯

maxgen=10;                         %½ø»¯´úÊý£¬¼´µü´ú´ÎÊý

sizepop=10;                        %ÖÖȺ¹æÄ£

pcross=0.3;                       %½»²æ¸ÅÂÊÑ¡Ôñ£¬0ºÍ1Ö®¼ä

pmutation=0.1;                    %±äÒì¸ÅÂÊÑ¡Ôñ£¬0ºÍ1Ö®¼ä



%½Úµã×ÜÊý

numsum=inputnum*hiddennum+hiddennum+hiddennum*outputnum+outputnum;



lenchrom=ones(1,numsum);        

bound=[-3*ones(numsum,1) 3*ones(numsum,1)];    %Êý¾Ý·¶Î§



%------------------------------------------------------ÖÖȺ³õʼ»¯--------------------------------------------------------

individuals=struct('fitness',zeros(1,sizepop), 'chrom',[]);  %½«ÖÖȺÐÅÏ¢¶¨ÒåΪһ¸ö½á¹¹Ìå

avgfitness=[];                      %ÿһ´úÖÖȺµÄƽ¾ùÊÊÓ¦¶È

bestfitness=[];                     %ÿһ´úÖÖȺµÄ×î¼ÑÊÊÓ¦¶È

bestchrom=[];                       %ÊÊÓ¦¶È×îºÃµÄȾɫÌå

%³õʼ»¯ÖÖȺ

for i=1:sizepop

    %Ëæ»ú²úÉúÒ»¸öÖÖȺ

    individuals.chrom(i,:)=Code(lenchrom,bound);    %±àÂ루binaryºÍgreyµÄ±àÂë½á¹ûΪһ¸öʵÊý£¬floatµÄ±àÂë½á¹ûΪһ¸öʵÊýÏòÁ¿£©

    x=individuals.chrom(i,:);

    %¼ÆËãÊÊÓ¦¶È

    individuals.fitness(i)=fun(x,inputnum,hiddennum,outputnum,net,inputn,outputn);   %ȾɫÌåµÄÊÊÓ¦¶È

end



%ÕÒ×îºÃµÄȾɫÌå

[bestfitness bestindex]=min(individuals.fitness);

bestchrom=individuals.chrom(bestindex,:);  %×îºÃµÄȾɫÌå

avgfitness=sum(individuals.fitness)/sizepop; %ȾɫÌåµÄƽ¾ùÊÊÓ¦¶È

% ¼Ç¼ÿһ´ú½ø»¯ÖÐ×îºÃµÄÊÊÓ¦¶ÈºÍƽ¾ùÊÊÓ¦¶È

trace=[avgfitness bestfitness];



%% µü´úÇó½â×î¼Ñ³õʼ·§ÖµºÍȨֵ

% ½ø»¯¿ªÊ¼

for i=1:maxgen

    i

    % Ñ¡Ôñ

    individuals=Select(individuals,sizepop);

    avgfitness=sum(individuals.fitness)/sizepop;

    %½»²æ

    individuals.chrom=Cross(pcross,lenchrom,individuals.chrom,sizepop,bound);

    % ±äÒì

    individuals.chrom=Mutation(pmutation,lenchrom,individuals.chrom,sizepop,i,maxgen,bound);

   

    % ¼ÆËãÊÊÓ¦¶È

    for j=1:sizepop
x=individuals.chrom(j,:);
        x=individuals.chrom(j,:); %½âÂë

        individuals.fitness(j)=fun(x,inputnum,hiddennum,outputnum,net,inputn,outputn);   

    end

   

  %ÕÒµ½×îСºÍ×î´óÊÊÓ¦¶ÈµÄȾɫÌå¼°ËüÃÇÔÚÖÖȺÖеÄλÖÃ

    [newbestfitness,newbestindex]=min(individuals.fitness);

    [worestfitness,worestindex]=max(individuals.fitness);

    % ´úÌæÉÏÒ»´Î½ø»¯ÖÐ×îºÃµÄȾɫÌå

    if bestfitness>newbestfitness

        bestfitness=newbestfitness;

        bestchrom=individuals.chrom(newbestindex,:);

    end

    individuals.chrom(worestindex,:)=bestchrom;

    individuals.fitness(worestindex)=bestfitness;

   

    avgfitness=sum(individuals.fitness)/sizepop;

   

    trace=[trace;avgfitness bestfitness]; %¼Ç¼ÿһ´ú½ø»¯ÖÐ×îºÃµÄÊÊÓ¦¶ÈºÍƽ¾ùÊÊÓ¦¶È



end

%% ÒÅ´«Ëã·¨½á¹û·ÖÎö

figure(1)

[r c]=size(trace);

plot([1:r]',trace(:,2),'b--');

title(['ÊÊÓ¦¶ÈÇúÏß  ' 'ÖÕÖ¹´úÊý£½' num2str(maxgen)]);

xlabel('½ø»¯´úÊý');ylabel('ÊÊÓ¦¶È');

legend('ƽ¾ùÊÊÓ¦¶È','×î¼ÑÊÊÓ¦¶È');

disp('ÊÊÓ¦¶È±äÁ¿');

x=bestchrom;



%% °Ñ×îÓųõʼ·§ÖµÈ¨Öµ¸³ÓèÍøÂçÔ¤²â

% %ÓÃÒÅ´«Ëã·¨ÓÅ»¯µÄBPÍøÂç½øÐÐÖµÔ¤²â

w1=x(1:inputnum*hiddennum);

B1=x(inputnum*hiddennum+1:inputnum*hiddennum+hiddennum);

w2=x(inputnum*hiddennum+hiddennum+1:inputnum*hiddennum+hiddennum+hiddennum*outputnum);

B2=x(inputnum*hiddennum+hiddennum+hiddennum*outputnum+1:inputnum*hiddennum+hiddennum+hiddennum*outputnum+outputnum);



net.iw{1,1}=reshape(w1,hiddennum,inputnum);

net.lw{2,1}=reshape(w2,outputnum,hiddennum);

net.b{1}=reshape(B1,hiddennum,1);

net.b{2}=B2;



%% BPÍøÂçѵÁ·

%ÍøÂç½ø»¯²ÎÊý

net.trainParam.epochs=100;

net.trainParam.lr=0.1;

%net.trainParam.goal=0.00001;



%ÍøÂçѵÁ·

[net,per2]=train(net,inputn,outputn);



%% BPÍøÂçÔ¤²â

%Êý¾Ý¹éÒ»»¯

inputn_test=mapminmax('apply',input_test,inputps);

an=sim(net,inputn_test);

test_simu=mapminmax('reverse',an,outputps);

error=test_simu-output_test;