一种基于粒子群优化最简回声状态网络参数的方法和系统的制作方法
【专利摘要】本发明涉及一种基于粒子群优化最简回声状态网络参数的方法和系统，其方法包括以下步骤：(1)建立最简回声状态网络模型；(2)设置最简回声状态网络模型的初始化参数；(3)建立适应度函数；(4)计算粒子目标函数；(5)更新各粒子个体最优解和结构参数；(6)判断是否达到终止条件；(7)输出粒子群最优解。本发明所引用最简回声状态网络模型包括输入层、储备池和输出层，该模型具有一个确定的储备池结构，其中的神经元以环形结构相连，增强了储备池拓扑的稳定性且简化了计算量，同时，基于粒子群优化算法运算速度快，全局搜索能力强的特点，对最简回声状态网络的多个参数进行了优化，从而得到一个粒子群最优解，提高了预测精度。
【专利说明】
-种基于粒子群优化最简回声状态网络参数的方法和系统
技术领域
[0001] 本发明设及机器学习领域，具体设及一种基于粒子群优化最简回声状态网络参数 的方法和系统。
【背景技术】
[0002] 回声状态网络ESN化Cho state化twork)是一种新型的递归神经网络，其独特的 动态储备池结构使网络具备良好的短期记忆能力，与传统的递归神经网络相比，ESN最大的 优势就是简化了网络的训练过程，解决了传统递归神经网络结构难W确定、训练算法过于 复杂的问题，同时也克服了递归网络存在的记忆渐消问题，但也存在一些问题，如：
[0003] (1)储备池内部随机拓扑，映射到未知的高维空间；
[0004] (2)ESN结构参数选取凭借经验，没有科学理论的支撑。
[0005] 运些问题阻止了ESN成为可W广泛使用的工具，为了改善运些不足，文献(Rodan A,Tino P.Minimum Complexity Echo State Network[J].iEEE Transactions on 化ural化tworks，2011，22.)提出了一种具有环形储备池拓扑结构的最简回声状态网络 MESN(Minimum ComplexityEcho Sl:ate Network),分别设置了输入权值矩阵和内部权值矩 阵的权值相同，运确保仅需调整较少的自由参数，同时该确定的储备池拓扑避免了不断重 复实验W获得性能良好的ESN结构。针对非线性系统逼近任务，该方法在保持精度不降低的 同时，极大地减少了计算量。
[0006] 相比经典回声状态网络，虽然MESN极大地简化了ESN结构，解决了第一个问题，即 储备池内部随机拓扑，映射到未知的高维空间，但是其参数选择如神经元数量、谱半径、输 入权值、内部权值和反馈权值等大多在给定的参数空间通过试凑的方法实现，或者根据经 验选择，存在很大的盲目性与不确定性，因此，如何寻找到最适合的参数解决参数选取的问 题成为了人们研究的热点方向。
【发明内容】
[0007] 本发明所要解决的技术问题是提供一种基于粒子群优化最简回声状态网络参数 的方法和系统，该方法将粒子群优化算法引入简最回声状态网络，优化其关键参数，有效地 加快了运算速度，提高了预测精度。
[000引本发明解决其技术问题所采用的技术方案是:一种基于粒子群优化最简回声状态 网络参数的方法，包括W下步骤：
[0009] 步骤1:建立一个最简回声状态网络模型，用环形拓扑结构代替原有回声状态网络 的随机拓扑结构；
[0010] 步骤2:设置最简回声状态网络模型的初始化参数，所述初始化参数包括结构参 数、粒子个体最优解、粒子群最优解、最大迭代次数及粒子群的粒子个数；
[0011] 步骤3:对最简回声状态网络模型的结构参数和最优解进行随机分配，建立最简回 声状态网络模型的适应度函数，并计算每个粒子的适应度值；
[0012] 步骤4:根据步骤3中计算的每个粒子的适应度值和步骤2中的初始化参数的粒子 个体最优解及粒子群最优解计算粒子目标函数；
[0013] 步骤5:根据粒子群算法的位置和速度更新公式，更新最简回声状态网络模型中粒 子个体最优解和结构参数，迭代次数加1;
[0014] 步骤6:判断迭代次数是否大于最大迭代次数，如果是，执行步骤7，否则，执行步骤 3至步骤5;
[0015] 步骤7:输出该粒子群最优解。
[0016] 进一步地，所述步骤1中的最简回声状态网络模型包括输入层、储备池和输出层， 所述输入层、储备池和输出层对应的输入向量、状态向量和输出向量分别为:K个输入单元U (n) =(山(n), . . .，Uk(n));N个储备池单元x(n) = (xi(n), . . .，XN(n));和L个输出单元y(n) = (yi(n),. . . ,yL(n))〇
[0017] 进一步地，所述储备池中包括多个神经元，每个神经元都有一个激励状态，且储备 池中神经元的激励状态x(n+l)，通过状态更新方程更新：
[001 引 x(n+l)=f(Winu(n+l)+Wx(n)+Wbacky(n))
[0019] 其中，u(n)为K个输入单元;x(n)为N个储备池单元;y(n)为L个输出单元;Wi嘴W分 别表示N X K输入权值矩阵和N X N储备池连接矩阵，Wbatk为一个N X L的反馈连接矩阵，分别 设置Win和W中的所有非零元素绝对值相同，等于非零值a G (-1，1 )，其权值符号通过随机方 法或者Logistic映射的方式产生，f( ?)表示储备池中神经元的激励函数。
[0020] 进一步地，所述储备池中神经元的激励函数为正弦函数或者S型函数。
[0021] 进一步地，根据所述储备池中神经元的激励状态x(n+l)，该模型的输出为：
[0022] y(n+l)=f〇ut(W〇ut(u(n+l)，x(n+l)，y(n)))
[0023] 其中，u(n)为K个输入单元;x(n)为N个储备池单元;y(n)为L个输出单元，是一 个输出函数，表示一个LX化+化L)的输出连接矩阵。
[0024] 进一步地，所述步骤3中的适应度函数包括MES的川练阶段和测试阶段的误差，可表 示为：
[0025] Fitness = JMErrortrain) +f 2 (Elrrortes t)
[0026] 训练阶段需要考虑该结构参数下ESN模型的逼近能力，用训练的误差化rortrain体 现，预测阶段考虑的是该结构参数下ESN模型的泛化能力，用训练的误差化rortest体现，其 中，fi( ? )、f2( ?)分别为Errortrain和Errortest的激励函数，Fitness为为每个粒子的适应度 值。
[0027] 进一步地，所述步骤4中计算目标函数的具体过程为：
[0028] (1)对于某一时刻的一个粒子将该粒子当前时刻的最优解/4所对应的适应 度值与该粒子个体最优解Pld的适应度值进行比较，若当前时刻的最优解括^所对应的适应 度值小于该粒子个体最优解Pld的适应度值，则用该粒子当前时刻的最优解把替换为该粒 子个体最优解Pld，否则Pld保持不变；
[0029] (2)将该粒子得到的最优解Pid所对应的适应度值与粒子群最优解Pgd所对应的适 应度值进行比较，若当前时刻的最优解Pid所对应的适应度值小于粒子群最优解Pgd的适应 度值，则将其最优解Pid的作为粒子群最优解Pgd，否则Pgd保持不变。
[0030] 进一步地，所述步骤2中结构参数包括神经元数量、输入权值矩阵、反馈权值矩阵、 储备池的权值、谱半径，其中，神经元数量初始化从10开始，变化范围在10-1000之间，其余 参数初始时刻为0~1的随机数。
[0031 ]进一步地，所述步骤5中的位置和速度更新公式如下：
[0032]
[0033]
[0034] 其中，ri，n是[0~1]区间生成的随机数，d = l，2,3,4,5,即解空间的维数为5或自 变量的个数为5，1 = 1，2-1，1是群体中粒子的个数，(31和〇2为加速因子，4表示迭代次数;1^；； 表示第i个粒子在第k次迭的第d维分量的最优解; < 表示第i个粒子在第k次迭代的第d维分 量的结构参数;=神经元数量=输入权值矩阵;=反馈权值矩阵=储备池的 权值;jfs =谱半径;Pid表示第i个粒子的个体最优解的第d维分量;Pgd表示粒子群最优解的 第d维分量。
[0035] -种基于粒子群优化最简回声状态网络参数的系统，包括建模单元、初始化单元、 建立函数单元、计算单元、更新单元、判断单元及输出单元；
[0036] 所述建模单元，用于建立一个最简回声状态网络模型，用环形拓扑结构代替原有 回声状态网络的随机拓扑结构；
[0037] 所述初始化单元，用于设置最简回声状态网络模型的初始化参数，所述初始化参 数包括结构参数、粒子个体最优解、粒子群最优解、最大迭代次数及粒子群的粒子个数；
[0038] 所述建立函数单元，用于对最简回声状态网络模型的结构参数和最优解进行随机 分配，建立最简回声状态网络模型的适应度函数，并计算每个粒子的适应度值；
[0039] 所述计算单元，用于根据计算的每个粒子的适应度值和初始化参数的粒子个体最 优解及粒子群最优解计算粒子目标函数；
[0040] 所述更新单元，用于根据粒子群算法的位置和速度更新公式，更新最简回声状态 网络模型中粒子个体最优解和结构参数.
[0041 ]所述判断单元，用于判断是否达到终止条件；
[0042] 所述输出单元，用于输出粒子群最优解。
[0043] 本发明的有益效果:本发明将粒子群优化算法引入最简回声状态网络优化其关键 参数，最简回声状态网络模型包括输入层、储备池、输出层，该模型具有一个确定的储备池 结构，其中的神经元W环形结构相连，增强了储备池拓扑的稳定性并且大大简化了计算量， 同时，基于粒子群优化算法运算速度快，全局捜索能力强的特点，对最简回声状态网络的多 个参数进行了优化，从而找寻到一个全局最优解，提高了预测精度。
【附图说明】
[0044] 图1为最简回声状态网络结构图；
[0045] 图2为本发明的方法流程图；
[0046] 图3粒子数为20时本发明逼近NARMA系统的效果图；
[0047] 图4粒子数为20本发明逼近NARMA系统的误差幅度图；
[004引图5粒子数为20时本发明逼近NARMA系统的性能评估图。
【具体实施方式】
[0049] W下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并 非用于限定本发明的范围。
[0化日]粒子群优化算法PSO(化dicle Swarm Optimization)具有原理简单，参数少，收 敛速度快，全局寻优能力强的特点，运种优化算法是根据鸟捕食的行为提出的，每个觅食的 鸟都是一个粒子，粒子在空间中W-定的速度飞行(运个速度根据它本身的飞行经验和同 伴的飞行经验来动态调整)所有粒子都有一个被目标函数决定的适应值一一运个值适用于 评价粒子的好坏程度。优化捜索正是在由运样一群随机初始化形成的粒子组成的一个种群 中，W迭代的方式进行的。
[0051] 本发明利用粒子群优化算法优化最简回声状态网络的5个关键储备池参数:神经 元数量、输入权值矩阵、反馈权值矩阵、储备池的权值及谱半径。
[0052] 如图1至图5所示，本发明提供了一种基于粒子群优化最简回声状态网络参数的方 法，包括W下步骤：
[0053] 步骤1:如图1所示，建立一个最简回声状态网络模型，用环形拓扑结构代替原有回 声状态网络的随机拓扑结构;=层最简回声状态网络模型由输入层、储备池和输出层=部 分构成，其对应的输入向量、状态向量和输出向量分别为:K个输入单元u(n) = (Ui(n)，...， uk(n)) ;N个储备池单元x(n) = (xi(n)，. . .，XN(n));和L个输出单元y(n) = (yi(n)，. . .，yL (n))，其中，该模型采用单个非线性神经元和延迟环代替传统储备池中大量随机互连的非 线性神经元，W简化其物理拓扑，运种系统更易于硬件实现，因为只包括两个元素:单个非 线性神经元和一个延迟环，储备池中的每个神经元都有一个激励状态，且储备池神经元的 激励状态x(n+l)，通过状态更新方程更新：
[0054] x(n+l)=f(Winu(n+l)+Wx(n)+Wbacky(n))
[005引其中，Wi嘴W分别表示NXK输入权值矩阵和NXN储备池连接矩阵，Wbaek为一个NXL 的反馈连接矩阵，分别设置Win和W中的所有非零元素绝对值相同，等于非零值a G (-1，1 )，其 权值符号通过随机方法或者Logistic映射的方式产生，f表示储备池神经元的激励函数，通 常为正弦函数或者S型函数。
[0056] 根据上述储备池激励状态x(n+l)，该模型的输出可W通过W下等式计算：
[0057] y(n+l)=f?t(W?t(u(n+l)，x(n+l)，y(n)))
[005引其中，u(n)为K个输入单元;x(n)为N个储备池单元;y(n)为L个输出单元，是一 个输出函数，表示表示一个LX化+化L)的输出连接矩阵。
[0059] 步骤2:设置最简回声状态网络模型的初始化参数，所述初始化参数包括（1)结构 参数:神经元数量、输入权值矩阵、反馈权值矩阵、储备池的权值、谱半径，其中，神经元数量 初始化从10开始，变化范围在10-1000之间，其余参数初始时刻为0~1的随机数；（2)最优 解:粒子个体最优解Pid和粒子群最优解Pgd;(3)粒子群的粒子个数、最大迭代次数。
[0060] 步骤3:对最简回声状态网络模型的结构参数和最优解进行随机分配，建立最简回 声状态网络模型的适应度函数，并计算每个粒子的适应度值;为确保回声状态网络的泛化 能力，需将全部训练数据分为训练集与预测集，先进行MES的川练阶段再进行预测阶段，所述 适应度函数包括MES的川练阶段和测试阶段的误差，可表示为：
[0061 ] Fitness = fi( Errortrain) +f 2 (Errortest)
[0062] 训练阶段需要考虑该结构参数下ESN模型的逼近能力，用训练的误差化rortrain体 现。预测阶段考虑的是该结构参数下ESN模型的泛化能力，用训练的误差化rortest体现，其 中，fi( ? )、f2( ?)分别为Errortrain和Errortest的激励函数，Fitness为每个粒子的适应度 值。
[0063] 步骤4:根据步骤3中计算的每个粒子的适应度值和步骤2中的初始化参数计算粒 子目标函数;计算粒子目标函数的具体过程为：
[0064] (1)对于某一时刻的一个粒子X:,,,，将该粒子当前时刻的最优解托/所对应的适应度 值与该粒子个体最优解Pid的适应度值进行比较，若当前时刻的最优解托,所对应的适应度 值小于该粒子个体最优解Pid的适应度值，则用该粒子当前时刻的最优解替换为该粒子 个体最优解Pid，否则Pid保持不变；
[0065] (2)将该粒子得到的最优解Pid所对应的适应度值与粒子群最优解Pgd所对应的适 应度值进行比较，若当前时刻的最优解Pid所对应的适应度值小于粒子群最优解Pgd的适应 度值，则将其最优解Pid的作为粒子群最优解Pgd，否则Pgd保持不变。
[0066] 步骤5:根据粒子群算法的位置和速度更新公式，更新粒子个体最优解和结构参 数，迭代次数加1;
[0067] 位晉和速度更新公式如下：
[006引
[0069]
[0070] 其中，ri，n是[0~1]区间生成的随机数，(1 = 1，2,3,4,5，即解空间的维数为5或自 变量的个数为5，1 = 1，2。礼1是群体中粒子的个数向和〇2为加速因子，4表示迭代次数;4 表示第i个粒子在第k次迭的第d维分量的最优解；表示第i个粒子在第k次迭代的第d维 分量的结构参数;=神经元数量;=输入权值矩阵；=反馈权值矩阵;xf 1 =储备池 的权值;4 =谱半径;Pid表示第i个粒子的个体最优解的第d维分量;Pgd表示粒子群最优解 的第d维分量。
[0071 ]步骤6:判断迭代次数是否大于最大迭代次数，如果是，执行步骤7，否则，执行步骤 3至步骤5;
[0072] 步骤7:输出粒子群最优解。
[0073] -种基于粒子群优化最简回声状态网络参数的系统，包括建模单元、初始化单元、 建立函数单元、计算单元、更新单元、判断单元及输出单元；
[0074] 所述建模单元，用于建立一个最简回声状态网络模型，用环形拓扑结构代替原有 回声状态网络的随机拓扑结构；
[0075] 所述初始化单元，用于设置最简回声状态网络模型的初始化参数，所述初始化参 数包括结构参数、粒子个体最优解、粒子群最优解、最大迭代次数及粒子群的粒子个数；
[0076] 所述建立函数单元，用于对最简回声状态网络模型的结构参数和最优解进行随机 分配，建立最简回声状态网络模型的适应度函数，并计算每个粒子的适应度值；
[0077] 所述计算单元，用于根据计算的每个粒子的适应度值和初始化参数的粒子个体最 优解及粒子群最优解计算粒子目标函数；
[0078] 所述更新单元，用于根据粒子群算法的位置和速度更新公式，更新最简回声状态 网络模型中粒子个体最优解和结构参数.
[0079] 所述判断单元，用于判断是否达到终止条件；
[0080] 所述输出单元，用于输出粒子群最优解。
[0081] 图3和图4分别给出了粒子数为20时本发明逼近NARMA系统的效果图和误差幅度 图，可W看出利用该方法可W有效地跟随NARMA系统的未来变化趋势。
[0082] 图5给出了粒子数为20时本发明逼近NARMA系统的性能评估图，可W看出在第16次 迭代时NRMSE最小，即MESN具有最佳的非线性逼近能力，此时优化结果为：储备池规模为 620，输入权值矩阵权值为0.55，反馈权值矩阵权值为0.21，储备池权值矩阵权值为0.48，谱 半径为0.77。
[0083] 本发明的有益效果:本发明将粒子群优化算法引入最简回声状态网络优化其关键 参数，最简回声状态网络模型包括输入层、储备池、输出层，该模型具有一个确定的储备池 结构，其中的神经元W环形结构相连，增强了储备池拓扑的稳定性并且大大简化了计算量， 同时，基于粒子群优化算法运算速度快，全局捜索能力强的特点，对最简回声状态网络的多 个参数进行了优化，从而找寻到一个全局最优解，提高了预测精度。
[0084] W上所述仅为本发明的较佳实施例，并不用W限制本发明，凡在本发明的精神和 原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种基于粒子群优化最简回声状态网络参数的方法，其特征在于，包括W下步骤： 步骤1:建立一个最简回声状态网络模型，用环形拓扑结构代替原有回声状态网络的随 机拓扑结构； 步骤2:设置最简回声状态网络模型的初始化参数，所述初始化参数包括结构参数、粒 子个体最优解、粒子群最优解、最大迭代次数及粒子群的粒子个数； 步骤3:对最简回声状态网络模型的结构参数和最优解进行随机分配，建立最简回声状 态网络模型的适应度函数，并计算每个粒子的适应度值； 步骤4:根据步骤3中计算的每个粒子的适应度值和步骤2中的初始化参数的粒子个体 最优解及粒子群最优解计算粒子目标函数； 步骤5:根据粒子群算法的位置和速度更新公式，更新最简回声状态网络模型中粒子个 体最优解和结构参数，迭代次数加1; 步骤6:判断迭代次数是否大于最大迭代次数，如果是，执行步骤7，否则，执行步骤3至 步骤5; 步骤7:输出粒子群最优解。2. 如权利要求1所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，所述步骤1中的最简回声状态网络模型包括输入层、储备池和输出层，所述输入层、储备 池和输出层对应的输入向量、状态向量和输出向量分别为：K个输入单元u(n) = (ui (η), . . . ,uk(n)) ;N个储备池单元x(n) = (xi(n), . . . ,ΧΝ(η));和L个输出单元y(n) = (yi (η)，...，yL(n))。3. 如权利要求2所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，所述储备池中包括多个神经元，每个神经元都有一个激励状态，且储备池中神经元的激 励状态x(n+l)为： x(n+l)=f(Winu(n+l)+Wx(n)+Wbacky(n)) 其中，U (η)为K个输入单元;X (η)为N个储备池单元;y (η)为L个输出单元;win和W分别表 示N X K输入权值矩阵和N X N储备池连接矩阵，Wbatk为一个N X L的反馈连接矩阵，分别设置 Win和W中的所有非零元素绝对值相同，等于非零值ae (-1，1)，其权值符号通过随机方法或 者Logistic映射的方式产生，f( ·)表示储备池中神经元的激励函数。4. 如权利要求3所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，所述储备池中神经元的激励函数为正弦函数或者S型函数。5. 如权利要求3所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，根据所述储备池中神经元的激励状态x(n+l)，该模型的输出为： y(n+l)=严(rut(u(n+l)，x(n+l)，y(n))) 其中，u(n)为K个输入单元;x(n)为N个储备池单元;y(n)为L个输出单元，rut是一个输 出函数，表示一个LX化+化L)的输出连接矩阵。6. 如权利要求1所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，所述步骤3中的适应度函数包括最简回声状态网络模型训练阶段和测试阶段的误差，表 示为： Pi tne S S - f 1 (Errortra in)+f2(Errortest) 训练阶段需要考虑该结构参数下ESN模型的逼近能力，用训练的误差Errortrain体现;预 测阶段考虑的是该结构参数下ESN模型的泛化能力，用训练的误差化rortest体现，其中，fi (· )、f2( ·)分别为Elrrortrain和Elrrortest的激励函数，Fitness为每个粒子的适应度值。7. 如权利要求1所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，所述步骤4中计算粒子目标函数的具体过程为： (1) 对于某一时刻的一个粒子，将该粒子当前时刻的最优解所对应的适应度值与 该粒子个体最优解Pid的适应度值进行比较，若当前时刻的最优解所对应的适应度值小 于该粒子个体最优解Pid的适应度值，则用该粒子当前时刻的最优解托，替换为该粒子个体 最优解Pid，否则Pid保持不变； (2) 将该粒子得到的最优解pid所对应的适应度值与粒子群最优解pgd所对应的适应度值 进行比较，若当前时刻的最优解Pid所对应的适应度值小于粒子群最优解Pgd的适应度值，贝U 将其最优解Pid的作为粒子群最优解Pgd,否则Pgd保持不变。8. 如权利要求1所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，所述步骤2中结构参数包括神经元数量、输入权值矩阵、反馈权值矩阵、储备池的权值、 谱半径，其中，神经元数量初始化从10开始，变化范围在10-1000之间，其余参数初始时刻为 0~1的随机数。9. 如权利要求8所述的一种基于粒子群优化最简回声状态网络参数的方法，其特征在 于，所述步骤5中的位置和速度更新公式如下：其中，η，η是[0~U区间生成的随机数，(1=1，2,3,4,5，即解空间的维数为5或自变量的个 数为5，1 = 1，2-礼1是群体中粒子的个数，〇和〇2为加速因子古表示迭代次数^表示第1个粒 子在第k次迭的第d维分量的最优解;X；；表示第i个粒子在第k次迭代的第d维分量的结构参数； 却=神经元救量：4崎资入权值化阵：鸣=反馈权植化阵；4=储备池的权值：4=谱半径： Pid表示第i个粒子的个体最优解的第d维分量;Pgd表示粒子群最优解的第d维分量。10. -种基于粒子群优化最简回声状态网络参数的系统，其特征在于，包括建模单元、 初始化单元、建立函数单元、计算单元、更新单元、判断单元及输出单元； 所述建模单元，用于建立一个最简回声状态网络模型，用环形拓扑结构代替原有回声 状态网络的随机拓扑结构； 所述初始化单元，用于设置最简回声状态网络模型的初始化参数，所述初始化参数包 括结构参数、粒子个体最优解、粒子群最优解、最大迭代次数及粒子群的粒子个数； 所述建立函数单元，用于对最简回声状态网络模型的结构参数和最优解进行随机分 配，建立最简回声状态网络模型的适应度函数，并计算每个粒子的适应度值； 所述计算单元，用于根据计算的每个粒子的适应度值和初始化参数的粒子个体最优解 及粒子群最优解计算粒子目标函数； 所述更新单元，用于根据粒子群算法的位置和速度更新公式，更新最简回声状态网络 模型中粒子个体最优解和结构参数； 所述判断单元，用于判断是否达到终止条件； 所述输出单元，用于输出粒子群最优解。
【文档编号】H04L12/24GK105978732SQ201610478359
【公开日】2016年9月28日
【申请日】2016年6月27日
【发明人】孙晓川, 陈扬, 张明辉, 李莹琦
【申请人】华北理工大学