基于极限学习机的脉动风速快速预测方法
【技术领域】
[0001] 本发明涉及一种采用单隐层前馈神经网络学习算法对单点脉动风速进行快速预 测，改进传统数据驱动技术预测需要参数寻优和模型选择导致算法时间长的缺陷，具体的 说是一种基于极限学习机（ExtremeLearningMachine，ELM)的脉动风速快速预测方法。
【背景技术】
[0002] 研究风荷载时，通常把风处理为在一定时距内不随时间变化的平均风速和随时间 随机变化的脉动风速两部分，平均风速产生结构静态响应，而脉动风速产生动态响应。风作 用在高层结构时，其正负风压对结构形成风荷载，同时钝体绕流还会引起结构抖振、旋涡脱 落引起的横向振动和扭转振动。极端风荷载作用下产生的抖振和颤振会引起建筑物倒塌或 严重破坏；动态位移超限易引起墙体开裂和附属构件破坏；大幅振动会造成居住和生活的 不舒适；脉动风频繁作用也会使外墙面构件和附属物产生疲劳破坏。掌握完整的脉动风速 时程资料对于结构设计、安全具有重要意义。
[0003] 基于数据驱动的样本学习训练为脉动风速速测提供可行的方法。目前脉动风速建 模预测的方法主要有时间序列分析法、人工神经网络、支持向量机等方法。然而这些方法 都存在着理论或应用上的不足，如时问序列模型高阶模型参数估计难度大、低阶模型预测 精度低；支持向量机虽然通过核函数定义的非线性变换将输入空间变换到一个高维空间， 在这个高维空间中寻找输入变量和输出变量之间的一种非线性关系，解决了"维数灾难"问 题，但核函数的选择和参数优化决定了模型的特性；人工神经网络模型应用较为成熟，作为 一种数据驱动算法，具有逼近任意非线性函数的能力，可以映射出序列问复杂的非线性关 系，从而在风速和风功率预测中得到广泛应用，然而传统的人工神经网络方法存在一些问 题，如算法运行时问长，容易陷入局部极小等。
[0004] 2004年由南洋理工大学黄广斌副教授提出极限学习机是一种简单易用、有效的单 隐层前馈神经网络学习算法。该算法在随机选择输入层权值和隐层神经元的前提下，仅通 过一步计算即可求得网络输出权值，同传统神经网络相比，极限学习机极大地提高了网络 的泛化能力和学习速度，具有较强的非线性拟合能力。极限学习机只需要设置网络的隐层 节点个数，在算法执行过程中不需要调整网络的输入权值以及隐元的偏置，并且产生唯一 的最优解，因此具有学习速度快且泛化性能好的优点。
【发明内容】
[0005] 本发明所要解决的技术问题是提供一种基于极限学习机的脉动风速快速预测方 法，解决传统的传统数据驱动技术预测模耗时严重等问题，通过数值模拟出脉动风速样本 与新型数据驱动技术ELM结合，利用数值模拟为数据驱动模拟提供样本数据，再通过数据 驱动技术预测所需后续时间的脉动风速，为抗风设计提供所需的完整风速时程曲线的预测 方法，节约了大量的时间成本，并计算实际风速与预测风速的平均绝对误差（MAE)、均方根 误差（RMSE)以及相关系数（R)评价本方法的精度。
[0006] 根据上述发明构思，本发明采用下述技术方案：一种基于极限学习机的脉动风速 快速预测方法，其特征在于，其包括以下步骤：
[0007] 第一步，利用ARMA模型模拟生成垂直空间点脉动风速样本，将每一个空间点的脉 动风速样本分为训练集、测试集两部分，对其分别进行归一化处理，取嵌入维数k= 10对进 行样本数据进行相空间重构；
[0008] 第二步，给定训练样本N= {(χ〇tjRn,t#Rn,i= 1，…，N};激励函数g(x) 和隐含层节点数M，由于激活函数可以无穷可微，只需要在训练之前设置适当的隐含层节点 数，为输入权和隐含层偏差进行随机赋值，建立单隐层前馈神经网络学习数学模型，计算隐 藏层输出矩阵和隐藏节点与输出神经元的连接权向量；
[0009] 第三步，利用第二步中计算的隐藏层输出矩阵和隐藏节点与输出神经元的连接权 向量建立回归数学模型，利用该模型对脉动风速测试集样本进行预测；
[0010] 第四步：将测试样本与分别利用极限学习机和PS0-MK-LSSVM预测的脉动风速结 果对比，同时计算预测风速与实际风速的平均绝对误差、均方根误差以及相关系数，评价本 发明的有效性，同时利用粒子群PS0优化组合核的最小二乘支持向量机对同样脉动风速进 行预测，分析比较两种方法的性能。
[0011] 优选地，所述第一步中ARMA模型模拟m维脉动风速表示为下式：
[0013] 式中，U(t)为脉动风速；AyB#别是mXm阶AR和MA模型的系数矩阵；X(t)为 mX1阶正态分布白噪声序列；P为自回归阶数、q为滑动回归阶数，t为时间。
[0014] 优选地，所述第二步中单隐层前馈神经网络学习数学模型，对于N个不同样本 (Xi,ti)，隐含层节点数目是M，激励函数为g(x)，其中Xi=[Xxi2,…，xin]TeRn，ti= [L，tl2，…，tinJTeR'其单隐层前馈神经网络学习数学模型SLFN的数学模型为：
[0016] 式中，a;=[aai2，…，ain]l连接第i个隐含层结点的输入权值；b;是第i个隐 含层结点的偏差；β[au，al2，…，ainJT是连接第i个隐含层结点的输出权值；激励函数 g(x)可以是"Sigmoid"、"Sine"等；tj为第j个节点的输出值。
[0017] 优选地，所述第三步中，设置隐层结点数Μ= 20、激励函数g(x)为"Sigmoid"，计 算的隐层输出矩阵和隐层结点与输出神经元的连接权向量，建立回归数学模型Ηβ=T，对 测试集进行预测，式中Τ为输出节点的输出矩阵。
[0018] 本发明极限学习机的脉动风速快速预测方法具有如下优点：相比于传统的神经网 络算法，ELM在训练的过程中不需要调整输入权值和偏置，只需根据相应算法来调整输出权 值β，便可获得一个全局最优解；相对于PS0-MK-LSSVM，其参数选择较为容易，训练速度显 著提升，且不会陷入局部最优，在时间消耗很少的前提下能达到很好的准确率。
【附图说明】
[0019] 图1是沿地面垂直方向20米处脉动风速模拟样本示意图。
[0020] 图2是沿地面垂直方向50米处脉动风速模拟样本示意图。
[0021] 图3是20米ELM与PS0-LSSVM预测风速与模拟风速对比示意图。
[0022] 图4是20米ELM与PS0-LSSVM预测风速与模拟风速自相关函数对比示意图。
[0023] 图5是50米ELM与PS0-LSSVM预测风速与模拟风速对比示意图。
[0024] 图6是50米ELM与PS0-LSSVM预测风速与模拟风速自相关函数对比示意图。
【具体实施方式】
[0025] 下面结合附图给出本发明较佳实施例，以详细说明本发明的技术方案。
[0026] 本发明极限学习机的脉动风速快速预测方法包括以下步骤：
[0027]第一步，利用ARMA(Auto-RegressiveandMovingAverageModel，自回归滑动 平均）模型模拟生成垂直空间点脉动风速样本，将每一个空间点的脉动风速样本分为训练 集、测试集两部分，对其分别进行归一化处理，取嵌入维数k= 10对进行样本数据进行相空 间重构；确定单点脉动风速样本的ARMA模型各参数，ARMA模型的自回归阶数p= 4,滑动回 归阶数q= 1。模拟某100米的超高层建筑，沿高度方向取每隔10米的点作为各模拟风速 点。其他相关参数见表1 :
[0028] 表1相关模拟参数表
[0029]
[0030] 模拟功率谱采用Kaimal谱，只考虑高度方向的空间相关性。模拟生成20、50米脉 动风速样本分别见图1、图2。
[0031] 第一步中，ARMA模型模拟m维脉动风速表示为下式（1):
[0033] 式中，U(t)为脉动风速；AuB#别是mXm阶AR和MA模型的系数矩阵；X(t)为 mX1阶正态分布白噪声序列；P为自回归阶数、q为滑动回归阶数，t为时间。
[0034] 将得到的脉动风速按式（2)进行归一化处理：
[0036] 式中，为归一化后脉动风速，yi为实际脉动风速样本，y_为实际脉动风速最大 值，卩_实际脉动风速最小值。
[0037] 第二步，给定训练样本N= {(Xi,tjRn,t#Rn,i= 1，…，N};激励函数g(x) 和隐含层节点数M，由于激活函数可以无穷可微，只需要在训练之前设置适当的隐含层节点 数，为输入权和隐含层偏差进行随机赋值，建立单隐层前馈神经网络学习数学模型，计算隐 藏层输出矩阵和隐藏节点与输出神经元的连接权向量；根据极限学习机ELM原理，激励函 数为"Sigmoid"和隐含层结点数为Μ= 20,给定训练样本N= {(Xi,tjRn,t#Rn,i =1，…，N}，得到ELM的算法为：
[0038] (1)随机选取权值ai，偏置h(i= 1，…，Μ);
[0040] (3)计算输出权重β:β=Η2Τ。
[0041] 第二步中，单隐层前馈神经网络学习数学模型建立：
[0042] 对于N个不同样本（Xl，，隐含层节点数目是M，激励函数为g(x)，其中Xl = [Xll，xl2，…，xin]TeRn，