• 使用遗传算法优化BP神经网络实现非线性函数拟合

            大家好，我是带我去滑雪！

            非线性函数拟合是一种用于找到与给定数据点集合最好匹配的非线性函数的过程。非线性函数拟合通常用于以下情况：

    • 数据趋势不是线性的：当数据点之间的关系不能用线性方程来表示时，需要使用非线性函数来更好地拟合数据。
    • 复杂的数据模式：当数据包含复杂的非线性模式或曲线时，非线性函数拟合可以提供更准确的拟合结果。
    • 物理、生物、经济或其他领域的实验数据：许多自然现象和实验数据无法用简单的线性模型来解释，因此需要使用非线性模型。

            下面使用遗传算法优化的BP神经网络实现非线性函数拟合： 

    目录

    一、预备知识 

    （1）什么是遗传算法

    （2）拟合函数

    二、模型建立

    （1）遗传算法优化BP神经网络算法流程

    （2）适应度函数

    （3）选择操作

     （4）交叉

    （5） 变异

    （6） 遗传算法主函数

    （7） 将遗传算法的到的最优个体赋值进BP神经网络

     （8）结果分析

    一、预备知识 

    （1）什么是遗传算法

            遗传算法（Genetic Algorithm, GA）是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型，是一种通过模拟自然进化过程搜索最优解的方法。其主要特点是直接对结构对象进行操作，不存在求导和函数连续性的限定；具有内在的隐并行性和更好的全局寻优能力；采用概率化的寻优方法，不需要确定的规则就能自动获取和指导优化的搜索空间，自适应地调整搜索方向。遗传算法以一种群体中的所有个体为对象，并利用随机化技术指导对一个被编码的参数空间进行高效搜索。其中，选择、交叉和变异构成了遗传算法的遗传操作；参数编码、初始群体的设定、适应度函数的设计、遗传操作设计、控制参数设定五个要素组成了遗传算法的核心内容。

    • 选择操作：从旧群体中以一定的概率选择个体到新群体中，个体被选中的概率跟适应度值有关，个体适应度值越好，被选择的概率越大。 
    • 交叉操作：从个体中选择两个个体，通过染色体的交换组合，来产生新的优秀个体。交叉过程为从群体中任选两个染色体，随机选择一点或多点染色体位置进行交换。
    • 变异操作：从群体中选择一个个体，选择染色体中的一点进行变异以产生更优秀的个体。

           运行参数主要包括群体大小、遗传代数、交叉概率、变异概率 。

    （2）拟合函数

          本次需要拟合的非线性函数为：

    y=x_{1}^{2}+x_{2}^{2}

    1. clc,clear,close all
    2. x=-6:0.1:6;
    3. [X,Y]=meshgrid(x);
    4. Z=X.^2+Y.^2;
    5. mesh(X,Y,Z)

       绘制函数图像为：

    二、模型建立

    （1）遗传算法优化BP神经网络算法流程

             流程分为：BP神经网络结构确定、遗传算法优化、BP神经网络预测。由于拟合函数中有两个输入参数，1个输出参数，所以设置BP神经网络结构为2-5-1，即输入层节点为2，隐含层节点为5，输出层节点为1，总共有2x5+5x1=15个权重，5+1=6个阈值，使用遗传算法个体的编码长度为15+6=21个。从非线性函数中随机得到2000组输入输出数据，随机选择其中1900组作为训练数据，100组作为测试数据。将训练数据预测误差绝对值和作为个体适应度值，个体适应度值越小，该个体越优。

           遗传算法参数设置为：种群规模为10，进化次数为50次，交叉概率为0.4，变异概率为0.2。

    （2）适应度函数

            将训练数据预测误差作为个体适应度值。 

    1. function error = fun(x,inputnum,hiddennum,outputnum,net,inputn,outputn)
    2. w1=x(1:inputnum*hiddennum);
    3. B1=x(inputnum*hiddennum+1:inputnum*hiddennum+hiddennum);
    4. w2=x(inputnum*hiddennum+hiddennum+1:inputnum*hiddennum+hiddennum+hiddennum*outputnum);
    5. B2=x(inputnum*hiddennum+hiddennum+hiddennum*outputnum+1:inputnum*hiddennum+hiddennum+hiddennum*outputnum+outputnum);
    6. net.trainParam.epochs=20;
    7. net.trainParam.lr=0.1;
    8. net.trainParam.goal=0.00001;
    9. net.trainParam.show=100;
    10. net.trainParam.showWindow=0;
    11. net.iw{1,1}=reshape(w1,hiddennum,inputnum);
    12. net.lw{2,1}=reshape(w2,outputnum,hiddennum);
    13. net.b{1}=reshape(B1,hiddennum,1);
    14. net.b{2}=B2;
    15. net=train(net,inputn,outputn);
    16.  
    17. an=sim(net,inputn);
    18.  
    19. error=sum(abs(an-outputn));

    （3）选择操作

           选择操作采用轮盘赌法从种群中选择适应度好的个体组成新种群。

    1. function ret=select(individuals,sizepop)
    2. fitness1=10./individuals.fitness;
    3. sumfitness=sum(fitness1);
    4. sumf=fitness1./sumfitness;
    5. index=[]; 
    6. for i=1:sizepop   %转sizepop次轮盘
    7.     pick=rand;
    8.     while pick==0    
    9.         pick=rand;        
    10.     end
    11.     for j=1:sizepop    
    12.         pick=pick-sumf(i);        
    13.         if pick<0        
    14.             index=[index j];            
    15.             break;
    16.         end
    17.     end
    18. end
    19. individuals.chrom=individuals.chrom(index,:);
    20. individuals.fitness=individuals.fitness(index);
    21. ret=individuals;

     （4）交叉

    1. function ret=Cross(pcross,lenchrom,chrom,sizepop,bound)
    2.  
    3.  for i=1:sizepop 
    4.      pick=rand(1,2);
    5.      while prod(pick)==0
    6.          pick=rand(1,2);
    7.      end
    8.      index=ceil(pick.*sizepop);
    9.      pick=rand;
    10.      while pick==0
    11.          pick=rand;
    12.      end
    13.      if pick>pcross
    14.          continue;
    15.      end
    16.      flag=0;
    17.      while flag==0
    18.          pick=rand;
    19.          while pick==0
    20.              pick=rand;
    21.          end
    22.          pos=ceil(pick.*sum(lenchrom));
    23.          pick=rand;
    24.          v1=chrom(index(1),pos);
    25.          v2=chrom(index(2),pos);
    26.          chrom(index(1),pos)=pick*v2+(1-pick)*v1;
    27.          chrom(index(2),pos)=pick*v1+(1-pick)*v2; 
    28.          flag1=test(lenchrom,bound,chrom(index(1),:));  
    29.          flag2=test(lenchrom,bound,chrom(index(2),:)); 
    30.          if   flag1*flag2==0
    31.              flag=0;
    32.          else flag=1;
    33.          end    
    34.      end
    35.  end
    36. ret=chrom;

    （5） 变异

    1. function ret=Mutation(pmutation,lenchrom,chrom,sizepop,num,maxgen,bound)
    2.  
    3. for i=1:sizepop  
    4.  
    5.     pick=rand;
    6.     while pick==0
    7.         pick=rand;
    8.     end
    9.     index=ceil(pick*sizepop);
    10.     pick=rand;
    11.     if pick>pmutation
    12.         continue;
    13.     end
    14.     flag=0;
    15.     while flag==0
    16.         pick=rand;
    17.         while pick==0      
    18.             pick=rand;
    19.         end
    20.         pos=ceil(pick*sum(lenchrom)); 
    21.     
    22.         pick=rand;
    23.         fg=(rand*(1-num/maxgen))^2;
    24.         if pick>0.5
    25.             chrom(i,pos)=chrom(i,pos)+(bound(pos,2)-chrom(i,pos))*fg;
    26.         else
    27.             chrom(i,pos)=chrom(i,pos)-(chrom(i,pos)-bound(pos,1))*fg;
    28.         end  
    29.         flag=test(lenchrom,bound,chrom(i,:)); 
    30.     end
    31. end
    32. ret=chrom;

    （6） 遗传算法主函数

    1. clc
    2. clear
    3. load data input output
    4.  
    5. inputnum=2;
    6. hiddennum=5;
    7. outputnum=1;
    8.  
    9. input_train=input(1:1900,:)';
    10. input_test=input(1901:2000,:)';
    11. output_train=output(1:1900)';
    12. output_test=output(1901:2000)';
    13.  
    14. [inputn,inputps]=mapminmax(input_train);
    15. [outputn,outputps]=mapminmax(output_train);
    16.  
    17. net=newff(inputn,outputn,hiddennum);
    18.  
    19. maxgen=20;  
    20. pcross=[0.2];
    21. pmutation=[0.1];
    22.  
    23. numsum=inputnum*hiddennum+hiddennum+hiddennum*outputnum+outputnum;
    24.  
    25. lenchrom=ones(1,numsum);        
    26. bound=[-3*ones(numsum,1) 3*ones(numsum,1)]; 
    27.  
    28. individuals=struct('fitness',zeros(1,sizepop), 'chrom',[]);
    29. avgfitness=[];
    30. bestfitness=[];
    31. bestchrom=[];
    32.  
    33. for i=1:sizepop
    34.     individuals.chrom(i,:)=Code(lenchrom,bound); 
    35.     x=individuals.chrom(i,:);
    36.     individuals.fitness(i)=fun(x,inputnum,hiddennum,outputnum,net,inputn,outputn); 
    37. end
    38. FitRecord=[];
    39.  
    40. [bestfitness bestindex]=min(individuals.fitness);
    41. bestchrom=individuals.chrom(bestindex,:); 
    42. avgfitness=sum(individuals.fitness)/sizepop; 
    43. trace=[avgfitness bestfitness]; 
    44.  
    45. for i=1:maxgen
    46.     i
    47.     individuals=Select(individuals,sizepop); 
    48.     avgfitness=sum(individuals.fitness)/sizepop;
    49.     individuals.chrom=Cross(pcross,lenchrom,individuals.chrom,sizepop,bound);
    50.  
    51.     individuals.chrom=Mutation(pmutation,lenchrom,individuals.chrom,sizepop,i,maxgen,bound);
    52.     
    53.     for j=1:sizepop
    54.         x=individuals.chrom(j,:); 
    55.         individuals.fitness(j)=fun(x,inputnum,hiddennum,outputnum,net,inputn,outputn);   
    56.     end
    57.  
    58.     [newbestfitness,newbestindex]=min(individuals.fitness);
    59.     [worestfitness,worestindex]=max(individuals.fitness);
    60.     if bestfitness>newbestfitness
    61.         bestfitness=newbestfitness;
    62.         bestchrom=individuals.chrom(newbestindex,:);
    63.     end
    64.     individuals.chrom(worestindex,:)=bestchrom;
    65.     individuals.fitness(worestindex)=bestfitness;
    66.     
    67.     avgfitness=sum(individuals.fitness)/sizepop;
    68.     
    69.     trace=[trace;avgfitness bestfitness];
    70.     FitRecord=[FitRecord;individuals.fitness];
    71. end

    （7） 将遗传算法的到的最优个体赋值进BP神经网络

    1. w1=x(1:inputnum*hiddennum);
    2. B1=x(inputnum*hiddennum+1:inputnum*hiddennum+hiddennum);
    3. w2=x(inputnum*hiddennum+hiddennum+1:inputnum*hiddennum+hiddennum+hiddennum*outputnum);
    4. B2=x(inputnum*hiddennum+hiddennum+hiddennum*outputnum+1:inputnum*hiddennum+hiddennum+hiddennum*outputnum+outputnum);
    5.  
    6. net.iw{1,1}=reshape(w1,hiddennum,inputnum);
    7. net.lw{2,1}=reshape(w2,outputnum,hiddennum);
    8. net.b{1}=reshape(B1,hiddennum,1);
    9. net.b{2}=B2;
    10.  
    11. net.trainParam.epochs=100;
    12. net.trainParam.lr=0.1;
    13.  
    14.  
    15. [net,per2]=train(net,inputn,outputn);
    16.  
    17. inputn_test=mapminmax('apply',input_test,inputps);
    18. an=sim(net,inputn_test);
    19. test_simu=mapminmax('reverse',an,outputps);
    20. error=test_simu-output_test;

     （8）结果分析

          绘制遗传算法优化过程中最优个体适应度值变化图：

    1. figure(1)
    2. [r c]=size(trace);
    3. plot([1:r]',trace(:,2),'b--');
    4. title(['适应度曲线  ' '终止代数＝' num2str(maxgen)]);
    5. xlabel('进化代数');ylabel('适应度');
    6. legend('平均适应度','最佳适应度');

    输出结果：

          模型的均方误差为2.2628x10的负7次方，什么模型的效果还可以。

    更多优质内容持续发布中，请移步主页查看。

    若有问题可邮箱联系：1736732074@qq.com 

    博主的WeChat:TCB1736732074

       点赞+关注,下次不迷路！

  • 相关阅读:
    halcon脚本-边缘及骨架的识别【附源码】
    定位消耗系统资源多的查询
    vue学习之属性绑定
    Vue Router
    Vue+element-ui实现表格数据渲染+分页
    C# out参数out多个参数
    魔兽世界安装插件后进游戏闪退的一个原因。
    Apache Hive源码阅读环境搭建
    使用 Kotlin DSL 编写网络爬虫
    C++多态与虚拟：运算符重载（Operator Overloading）
  • 原文地址：https://blog.csdn.net/qq_45856698/article/details/134292075