1 简介

针对变压器故障的特征,结合变压器油中气体分析法以及三比值法.提出了基于海洋捕食者算法改进极限学习机的故障诊断方法.由于输入层与隐含层的权值和阈值是随机产生.传统的极限学习机可能会使隐含层节点过多,训练过程中容易产生过拟合现象.该方法运用海洋捕食者算法对极限学习机的输入层与隐含层的权值与阈值进行优化,从而提高模型的稳定性和预测精度.将诊断结果与传统的基于极限学习机故障诊断进行对比,结果表明,基于海洋捕食者算法改进极限学习机变压器故障诊断的精度更高.

基于 MPA优化 ELM 的变压器故障诊断的具体步骤为：

1)确定 ELM 的拓扑结构。即输入层神经元个数，隐含层神经元个数以及输出层神经元个数；

2)对 ELM 中输入层到隐含层的权值以及阈值进行编码，得到初始种群；

3)解码得到权值和阈值，将权值和阈值带入到ELM 的训练网络中，使用训练样本进行训练；

4)训练完成后，使用测试样本进行测试，将测试样本的期望值和预测值的误差平方和作为适应度函数；

5)对种群进行选择，交叉，变异，得到新的种群，如果满足条件，则得出了误差平方和最小的网络权值和阈值，如果不满足条件，则返回步骤 2)；

6)将优化后的权值和阈值带入到训练网络中，计算隐含层输出矩阵 H， 并求解矩阵 H 的 MoorePenrose 广义逆 H+；

7)计算输出层权值β赞=H+T；

8)将测试样本带入到模型中进行预测。

2 部分代码


          
          

           
           %_________________________________________________________________________
           
           
%  Marine Predators Algorithm source code (Developed in MATLAB R2015a)
           
           
%
           
           
function [Top_predator_fit,Top_predator_pos,Convergence_curve]=MPA(SearchAgents_no,Max_iter,lb,ub,dim,fobj)
           
           
Top_predator_pos=zeros(1,dim);
           
           
Top_predator_fit=inf; 
           
           
Convergence_curve=zeros(1,Max_iter);
           
           
stepsize=zeros(SearchAgents_no,dim);
           
           
fitness=inf(SearchAgents_no,1);
           
           
Prey=initialization(SearchAgents_no,dim,ub,lb);
           
           
Xmin=repmat(ones(1,dim).*lb,SearchAgents_no,1);
           
           
Xmax=repmat(ones(1,dim).*ub,SearchAgents_no,1);
           
           
Iter=0;
           
           
FADs=0.2;
           
           
P=0.5;
           
           
while Iter<Max_iter    
           
           
     %------------------- Detecting top predator -----------------    
           
           
 for i=1:size(Prey,1)  
           
           
    Flag4ub=Prey(i,:)>ub;
           
           
    Flag4lb=Prey(i,:)<lb;    
           
           
    Prey(i,:)=(Prey(i,:).*(~(Flag4ub+Flag4lb)))+ub.*Flag4ub+lb.*Flag4lb;                    
           
           
    fitness(i,1)=fobj(Prey(i,:));
           
           
     if fitness(i,1)<Top_predator_fit 
           
           
       Top_predator_fit=fitness(i,1); 
           
           
       Top_predator_pos=Prey(i,:);
           
           
     end          
           
           
 end
           
           
     %------------------- Marine Memory saving ------------------- 
           
           
 if Iter==0
           
           
   fit_old=fitness;    Prey_old=Prey;
           
           
 end
           
           
  Inx=(fit_old<fitness);
           
           
  Indx=repmat(Inx,1,dim);
           
           
  Prey=Indx.*Prey_old+~Indx.*Prey;
           
           
  fitness=Inx.*fit_old+~Inx.*fitness;
           
           
  fit_old=fitness;    Prey_old=Prey;
           
           
     %------------------------------------------------------------   
           
           
 Elite=repmat(Top_predator_pos,SearchAgents_no,1);  %(Eq. 10) 
           
           
 CF=(1-Iter/Max_iter)^(2*Iter/Max_iter);
           
           
 RL=0.05*levy(SearchAgents_no,dim,1.5);   %Levy random number vector
           
           
 RB=randn(SearchAgents_no,dim);          %Brownian random number vector
           
           
  for i=1:size(Prey,1)
           
           
     for j=1:size(Prey,2)        
           
           
       R=rand();
           
           
          %------------------ Phase 1 (Eq.12) ------------------- 
           
           
       if Iter<Max_iter/3 
           
           
          stepsize(i,j)=RB(i,j)*(Elite(i,j)-RB(i,j)*Prey(i,j));                    
           
           
          Prey(i,j)=Prey(i,j)+P*R*stepsize(i,j); 
           
           
          %--------------- Phase 2 (Eqs. 13 & 14)----------------
           
           
       elseif Iter>Max_iter/3 && Iter<2*Max_iter/3 
           
           
         if i>size(Prey,1)/2
           
           
            stepsize(i,j)=RB(i,j)*(RB(i,j)*Elite(i,j)-Prey(i,j));
           
           
            Prey(i,j)=Elite(i,j)+P*CF*stepsize(i,j); 
           
           
         else
           
           
            stepsize(i,j)=RL(i,j)*(Elite(i,j)-RL(i,j)*Prey(i,j));                     
           
           
            Prey(i,j)=Prey(i,j)+P*R*stepsize(i,j);  
           
           
         end  
           
           
         %----------------- Phase 3 (Eq. 15)-------------------
           
           
       else 
           
           
           stepsize(i,j)=RL(i,j)*(RL(i,j)*Elite(i,j)-Prey(i,j)); 
           
           
           Prey(i,j)=Elite(i,j)+P*CF*stepsize(i,j);  
           
           
       end  
           
           
      end                                         
           
           
  end    
           
           
     %------------------ Detecting top predator ------------------        
           
           
  for i=1:size(Prey,1)  
           
           
    Flag4ub=Prey(i,:)>ub;  
           
           
    Flag4lb=Prey(i,:)<lb;  
           
           
    Prey(i,:)=(Prey(i,:).*(~(Flag4ub+Flag4lb)))+ub.*Flag4ub+lb.*Flag4lb;
           
           
    fitness(i,1)=fobj(Prey(i,:));
           
           
      if fitness(i,1)<Top_predator_fit 
           
           
         Top_predator_fit=fitness(i,1);
           
           
         Top_predator_pos=Prey(i,:);
           
           
      end     
           
           
  end
           
           
     %---------------------- Marine Memory saving ----------------
           
           
 if Iter==0
           
           
    fit_old=fitness;    Prey_old=Prey;
           
           
 end
           
           
    Inx=(fit_old<fitness);
           
           
    Indx=repmat(Inx,1,dim);
           
           
    Prey=Indx.*Prey_old+~Indx.*Prey;
           
           
    fitness=Inx.*fit_old+~Inx.*fitness;
           
           
    fit_old=fitness;    Prey_old=Prey;
           
           
     %---------- Eddy formation and FADs? effect (Eq 16) ----------- 
           
           
  if rand()<FADs
           
           
     U=rand(SearchAgents_no,dim)<FADs;                                                                                              
           
           
     Prey=Prey+CF*((Xmin+rand(SearchAgents_no,dim).*(Xmax-Xmin)).*U);
           
           
  else
           
           
     r=rand();  Rs=size(Prey,1);
           
           
     stepsize=(FADs*(1-r)+r)*(Prey(randperm(Rs),:)-Prey(randperm(Rs),:));
           
           
     Prey=Prey+stepsize;
           
           
  end
           
           
  Iter=Iter+1;  
           
           
  Convergence_curve(Iter)=Top_predator_fit; 
           
           
end
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.
  • 7.
  • 8.
  • 9.
  • 10.
  • 11.
  • 12.
  • 13.
  • 14.
  • 15.
  • 16.
  • 17.
  • 18.
  • 19.
  • 20.
  • 21.
  • 22.
  • 23.
  • 24.
  • 25.
  • 26.
  • 27.
  • 28.
  • 29.
  • 30.
  • 31.
  • 32.
  • 33.
  • 34.
  • 35.
  • 36.
  • 37.
  • 38.
  • 39.
  • 40.
  • 41.
  • 42.
  • 43.
  • 44.
  • 45.
  • 46.
  • 47.
  • 48.
  • 49.
  • 50.
  • 51.
  • 52.
  • 53.
  • 54.
  • 55.
  • 56.
  • 57.
  • 58.
  • 59.
  • 60.
  • 61.
  • 62.
  • 63.
  • 64.
  • 65.
  • 66.
  • 67.
  • 68.
  • 69.
  • 70.
  • 71.
  • 72.
  • 73.
  • 74.
  • 75.
  • 76.
  • 77.
  • 78.
  • 79.
  • 80.
  • 81.
  • 82.
  • 83.
  • 84.
  • 85.
  • 86.
  • 87.
  • 88.
  • 89.
  • 90.
  • 91.
  • 92.
  • 93.
  • 94.
  • 95.
  • 96.

3 仿真结果

【ELM预测】基于海洋捕食者算法优化极限学习机预测Matlab代码_matlab代码

【ELM预测】基于海洋捕食者算法优化极限学习机预测Matlab代码_权值_02

​4 参考文献

[1]曹卫东, 柳坤, 倪建军,等. 海洋捕食者算法驱动的滚齿参数低碳求解方法:.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

【ELM预测】基于海洋捕食者算法优化极限学习机预测Matlab代码_故障诊断_03