【预测模型-ELM预测】基于蝙蝠算法优化极限学习机预测matlab代码

1 简介

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

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

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

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

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

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

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

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

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

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

2 部分代码

clc,clear,close all

warning off

% BA算法参数

maxiter = 200;  % 迭代次数

sizepop = 10;  % 种群数量

% 频率范围

popmin1 = -1;  popmax1 = 1; % x1  频率

popmin2 = -1;  popmax2 = 1; % x2  频率

Qmin = 0.1;    % 最小频率

Qmax = 0.5;    % 最大频率

A = 0.1;          % 音量 (不变或者减小)

impluse = 0.85;   % 脉冲率 (不变或增加)

Vmin = -1; % 最小速度

Vmax = 1;  % 最大速度

%% 初始化种群

            V(j,1)=Vmax;

        end

        if V(j,1)<Vmin

            V(j,1)=Vmin;

        end

        % V--x2

        if V(j,2)>Vmax

            V(j,2)=Vmax;

        end

        if V(j,2)<Vmin

            V(j,2)=Vmin;

        end

        

        pop(j,:) = pop(j,:) + 0.5*V(j,:);

        

        % 脉冲率

        if rand>impluse

            pop(j,:) = zbest + A * randn(1,2);

%             x1 = popmin1 + (popmax1-popmin1)*rand;

%             x2 = popmin2 + (popmax2-popmin2)*rand;

%             pop(j,:) = [x1,x2];

        end

        

        % x1  越界限制

        if pop(j,1)>popmax1

            pop(j,1)=popmax1;

        end

        if pop(j,1)<popmin1

            pop(j,1)=popmin1;

        end

        % x2  越界限制

        if pop(j,2)>popmax2

            pop(j,2)=popmax2;

        end

        if pop(j,2)<popmin2

            pop(j,2)=popmin2;

        end

        

        % 适应度更新

        fitness(j) = fun(pop(j,:));

        

        % 比较  个体间比较

        if fitness(j)<fitnessgbest(j)

            fitnessgbest(j) = fitness(j);

            gbest(j,:) = pop(j,:);

        end

        if fitness(j)<bestfitness

            bestfitness = fitness(j);

            zbest =  pop(j,:);

        end

        

    end

    fitness_iter(i) = bestfitness;

end

disp('最优解')

disp(zbest)

fprintf('\n')

figure('color',[1,1,1])

plot(fitness_iter,'ro-','linewidth',2)

figure('color',[1,1,1])

loglog(fitness_iter,'ro-','linewidth',2)

axis tight

3 仿真结果

4 参考文献

[1]郑嘉利, & 王哲. (2019). 一种基于蝗虫算法和极限学习机的RFID室内定位方法.  CN109598320A.

[2]陈绍炜等. "基于蝙蝠算法优化ELM的模拟电路故障诊断研究." 电子测量技术 2(2015):138-141.

​​​​​​​