基于麻雀搜索算法的极限学习机(ELM)分类算法-附代码

智能算法研学社（Jack旭）

已于 2024-01-14 14:59:42 修改

阅读量6.8k

点赞数 5

分类专栏：智能优化算法应用机器学习 # ELM极限学习机文章标签：神经网络机器学习算法深度学习

于 2020-12-14 16:47:27 首次发布

本文链接：https://blog.csdn.net/u011835903/article/details/111177850

版权

智能优化算法应用同时被 3 个专栏收录

4161 篇文章 380 订阅

订阅专栏

机器学习

1950 篇文章 40 订阅

订阅专栏

ELM极限学习机

276 篇文章 2 订阅

订阅专栏

基于麻雀搜索算法的极限学习机(ELM)分类算法

文章目录

基于麻雀搜索算法的极限学习机(ELM)分类算法

摘要：本文利用麻雀搜索算法对极限学习机进行优化，并用于分类问题

1.极限学习机原理概述

典型的单隐含层前馈神经网络结构如图1 所示，由输入层、隐含层和输出层组成，输入层与隐含层、隐含层与输出层神经元间全连接。其中，输入层有 n 个神经元，对应 n 个输入变量，隐含层有 l个神经元；输出层有 m 个神经元，对应 m 个输出变量。为不失一般性，设输入层与隐含层间的连接权值 w 为:
$=\left[\begin{matrix}w_{11}&w_{12}&...&w_{1,n}\\ w_{21}&w_{22}&...&w_{2n}\\ ...\\ w_{l1}&w_{l2}&...&w_{ln} \end{matrix}\right]\tag{1}$
其中， $w_n$ 表示输入层第 $i$ 个神经元与隐含层第 $j$ 个神经元间的连接权值。

在这里插入图片描述

设隐含层与输出层间的连接权值，为 $\beta$ :
$\beta =\left[\begin{matrix} \beta_{11}&\beta_{12}&...&\beta_{1m}\\ \beta_{21}&\beta_{22}&...&\beta_{2m}\\ ...\\ \beta_{l1}&\beta_{l2}&...&\beta_{lm} \end{matrix}\right] \tag{2}$
其中，自 $\beta_{jk}$ 表示隐含层第 j 个神经元与输出层第 k个神经元间的连接权值。

设隐含层神经元的阈值值 b 为：
$=\left[\begin{matrix}b_1\\ b_2\\ ...\\ b_l \end{matrix}\right]\tag{3}$
设具有 Q 个样本的训练集输入矩阵 X 和输出矩阵 Y 分别为
$=\left[\begin{matrix}x_{11}&x_{12}&...&x_{1Q}\\ x_{21}&x_{22}&...&x_{2Q}\\ ...\\ x_{n1}&x_{n2}&...&x_{nQ} \end{matrix}\right]\tag{4}$

$=\left[\begin{matrix}y_{11}&y_{12}&...&y_{1Q}\\ y_{21}&y_{22}&...&y_{2Q}\\ ...\\ y_{m1}&y_{m2}&...&y_{mQ} \end{matrix}\right]\tag{5}$

设隐含层神经元的激活函数为 g(x)，则由图1 可得，网络的输出 T 为:
$[t_1,..,t_Q]_{m*Q},t_j = [t_{1j},...,t_{mj}]^T =\left[\begin{matrix}\sum_{i=1}^t\beta_{i1}g(w_ix_j + b_i)\\ \sum_{i=1}^t\beta_{i2}g(w_ix_j + b_i)\\ ...\\ \sum_{i=1}^t\beta_{im}g(w_ix_j + b_i) \end{matrix}\right]_{m*1},(j=1,2,...,Q)\tag{6}$
式(6)可表示为：
$H\beta = T’ \tag{7}$
其中， T’为矩阵 T 的转置； H 称为神经网络的隐含层输出矩阵，具体形式如下：
$H(w_1,...,w_i,b_1,...,b_l,x_1,...,x_Q) =\left[\begin{matrix} g(w_1*x_1 + b_1)&g(w_2*x_1 + b_2)&...&g(w_l*x_1 + b_l)\\ g(w_1*x_2 + b_1)&g(w_2*x_2 + b_2)&...&g(w_l*x_2 + b_l)\\ ...\\ g(w_1*x_Q + b_1)&g(w_2*x_Q + b_2)&...&g(w_l*x_Q + b_l) \end{matrix}\right]_{Q*l}$

2.ELM学习算法

由前文分析可知，ELM在训练之前可以随机产生 w 和 b ，只需确定隐含层神经元个数及隐含层和神经元的激活函数（无限可微），即可计算出 $\beta$ 。具体地， ELM 的学习算法主要有以下几个步骤：

（1）确定隐含层神经元个数，随机设定输入层与隐含层间的连接权值 w 和隐含层神经元的偏置 b ;

（2）选择一个无限可微的函数作为隐含层神经元的激活函数，进而计算隐含层输出矩阵 H ;

（3）计算输出层权值： $\beta = H^+T'$

值得一提的是，相关研究结果表明，在 ELM 中不仅许多非线性激活函数都可以使用（如 S 型函数、正弦函数和复合函数等)，还可以使用不可微函数，甚至可以使用不连续的函数作为激活函数。

3.分类问题

本文对乳腺肿瘤数据进行分类。采用随机法产生训练集和测试集，其中训练集包含 500 个样本，测试集包含 69 个样本。

4.基于麻雀搜索算法优化的ELM

麻雀搜索算法的具体原理参考博客：https://blog.csdn.net/u011835903/article/details/108830958。

由前文可知，ELM的初始权值和阈值都是随机产生。每次产生的初始权值和阈值具有满目性。本文利用麻雀搜索算法对初始权值和阈值进行优化。适应度函数设计为训练集的错误率与测试集的错误率的和，以期望使训练得到的网络在测试集和训练集上均有较好的结果:
$f i t n ess = a r g min (T r ain E rror R a t e + T es tE rror R a t e) 。$

5.测试结果

麻雀算法相关参数如下：

%训练数据相关尺寸
R = size(Pn_train,1);
S = size(Tn_train,1);
N = 20;%隐含层个数
%% 定义麻雀优化参数
pop=20; %种群数量
Max_iteration=50; %  设定最大迭代次数
dim = N*R + N*S;%维度，即权值与阈值的个数
lb = [-1.*ones(1,N*R),zeros(1,N*S)];%下边界
ub = [ones(1,N*R),ones(1,N*S)];%上边界

将经过麻雀优化后的SSA-ELM与基础ELM进行对比。

预测结果如下图

麻雀收敛曲线如下：

在这里插入图片描述

数据结果如下：

麻雀优化ELM结果展示：----------------
训练集正确率Accuracy = 99%(495/500)
测试集正确率Accuracy = 95.6522%(66/69)
病例总数：569 良性：357 恶性：212
训练集病例总数：500 良性：318 恶性：182
测试集病例总数：69 良性：39 恶性：30
良性乳腺肿瘤确诊：39 误诊：0 确诊率p1=100%
恶性乳腺肿瘤确诊：27 误诊：3 确诊率p2=90%

传统ELM结果展示：----------------
训练集正确率Accuracy = 92.8%(464/500)
测试集正确率Accuracy = 84.058%(58/69)
病例总数：569 良性：357 恶性：212
训练集病例总数：500 良性：318 恶性：182
测试集病例总数：69 良性：39 恶性：30
良性乳腺肿瘤确诊：36 误诊：3 确诊率p1=92.3077%
恶性乳腺肿瘤确诊：22 误诊：8 确诊率p2=73.3333%

从上述数据可以看出，SSA-ELM训练得到的网络，无论是在测试集和训练集上的正确率均高于基础ELM训练得到的网络。SSA-ELM具有较好的性能。