基于粒子群优化反向传播神经网络的脉动风速预测方法
【技术领域】
[0001] 本发明涉及一种采用粒子群优化算法优化神经网络学习算法对单点脉动风速进 行预测，改进传统神经网络容易陷入局部极值，导致收敛速度慢的缺陷，具体的说是一种基 于粒子群优化反向传播（Back Propagation，BP)神经网络的脉动风速预测方法。
【背景技术】
[0002] 研究风荷载时，通常把风处理为在一定时距内不随时间变化的平均风速和随时间 随机变化的脉动风速两部分，平均风速产生结构静态响应，而脉动风速产生动态响应。风作 用在高层结构时，其正负风压对结构形成风荷载，同时钝体绕流还会引起结构抖振、旋涡脱 落引起的横向振动和扭转振动。极端风荷载作用下产生的抖振和颤振会引起建筑物倒塌或 严重破坏；动态位移超限易引起墙体开裂和附属构件破坏；大幅振动会造成居住和生活的 不舒适；脉动风频繁作用也会使外墙面构件和附属物产生疲劳破坏。掌握完整的脉动风速 时程资料对于结构设计、安全具有重要意义。
[0003] 基于知识发现的风速智能模型预测方法不局限于用数学方程式描述预测自变量 与预测因子之间的关系，而是根据预测方法本身的结构机制获取预测的"知识"和"经验"。 这类方法具有自适应能力，可以较好的处理非线性、非确定性时间序列的预测问题。目前脉 动风速建模预测的方法主要有时间序列分析法、神经网络、支持向量机等方法。然而这些方 法都存在着理论或应用上的不足，如时问序列模型高阶模型参数估计难度大、低阶模型预 测精度低；支持向量机（SVM)虽然通过核函数定义的非线性变换将输入空间变换到一个高 维空间，在这个高维空间中寻找输入变量和输出变量之间的一种非线性关系，解决了 "维数 灾难"问题，但核函数的选择和参数优化决定了模型的特性；神经网络预测模型采用的是黑 箱型的学习模式，能够识别复杂性、非线性的特性。该预测模型是利用大量的历史数据来训 练神经网络，学习完成后得到输入与输出的一种映射关系。如果给定相应的输入，通过这种 映射关系就可以获得到对应的输出即预测结果。神经网络预测模型也存在着不足，比如结 构的设定较难，网络初始的权值和偏置是随机选取的，倘若这些参数的位置选择不当，则会 导致网络的收敛速度慢、陷入局部最优。
[0004] 随着智能优化算法的发展，神经网络的预测性能得到大大的提高。粒子群 (PSO-BP)是一种并行随机搜索最优化方法，具有全局搜索能力，采用粒子群优化算法来确 定神经网络初始权值和阈值，有利于神经网络网络很快收敛到全局最优解。
【发明内容】
[0005] 本发明的目的在于提供一种基于粒子群优化反向传播神经网络的脉动风速预测 方法，解决传统的传统神经网络容易陷入局部极值，导致收敛速度慢的问题。本发明利用 自回归滑动平均（Auto-Regressive and Moving Average，ARMA)方法数值有限脉动风速 作为样本数据，再建立粒子群优化算法反向传播（Particle Swarm Optimization-Back Propagati〇n，PS0-BP)神经网络预测模型，利用该模型预测所需后续时间的脉动风速，计算 实际风速与预测风速的平均绝对误差（MAE)、均方根误差（RMSE)以及相关系数（R)来评价 本方法的精度。
[0006] 根据上述发明构思，本发明采用下述技术方案：本发明粒子群优化神经网络的脉 动风速预测方法包括以下步骤：
[0007] 第一步：利用自回归滑动平均模型模拟生成垂直空间点脉动风速样本，将每一个 空间点的脉动风速样本分为训练集、测试集两部分，对其分别进行归一化处理；对数据归一 化处理后，目的是为了取消各维数据之间的数量级差别，从而避免因为输入输出数据的数 量级差别大而造成网络预测误差大的情况，取嵌入维数k = 10对进行样本数据进行相空间 重构；确定单点脉动风速样本的自回归滑动平均模型各参数，自回归滑动平均模型的自回 归阶数P = 4,滑动回归阶数q = 1 ;
[0008] 第二步：取采样时间1000s的20m、50m脉动风速样本，嵌入维数k = 10,对样本数 据进行相空间重构；给定训练样本N = ((Xi, tj IxiG Rn, Rn, i = 1，…，N}、激励函数 g(x)，设定神经网络模型输入节点m = 10、隐层节点p = 10、输出节点q = 1、网络训练次数 印ochs = 100、误差精度goal = 0· 0001及学习率Ir = 0· 1 ;将l-790s脉动风速作为训练 集，791-990S脉动风速作为测试集，归一化处理；
[0009] 第三步：设置粒子群优化算法粒子种群规模M，迭代次数K，学习因子CpC2，随机产 生粒子初始位置和初始速度，采用粒子群优化算法将神经网络的权值与偏置映射为粒子群 优化算法中的粒子，以均方根误差作为适应度函数，通过粒子的速度与位置不断更新来优 化这些参数，从而实现网络训练；
[0010] 第四步：将测试样本与利用粒子群优化算法神经网络模型和反向传播神经网络预 测的脉动风速结果对比；同时计算预测风速与实际风速的平均绝对误差、均方根误差以及 相关系数，评价本发明的有效性，同时利用未优化的神经网络对同样脉动风速进行预测，分 析比较两种方法的性能。
[0011] 优选地，上述第一步中，自回归滑动平均模型模拟m维脉动风速表示为下式：
[0013] 式中，U(t)为脉动风速;A1J,分别是mXm阶自回归模型和滑动平均模型的系数 矩阵；X (t)为mX 1阶正态分布白噪声序列；P为自回归阶数、q为滑动回归阶数，t为时间。
[0014] 优选地，所述第二步中，设定反向传播神经网络输入节点为m = 10,隐层节点为p =10,输出节点为q = 1，则反向传播神经网络完成映射f:Rm- R \隐层各节点的输入Sj为：
[0016] 反向传播神经网络激励函数为g(x)，则隐层节点的输出1^为：
[0018] 同理，输出层节点的输入L、输出y分别为： CN 105184416 A ^ ΗΠ T> 3/6 贞
[0021] 反向传播神经网络的连接权重COipVj和偏置θ ρ γ可以通过反向传播神经网络 训练求得，故y即为反向传播神经网络的预测模型。
[0022] 优选地，所述第三步中，设置神经网络输入节点为m= 10,隐层节点为p = 10,输 出节点为q = 1，激励函数g(x)为"Sigmoid" ；设定粒子群优化算法的粒子种群规模M = 30,迭代次数K =100,学习因子(^= C2= 2,将神经网络的权值与偏置映射为粒子群优化算 法中的粒子并且初始化，计算各粒子的适应度函数值，根据适应度比较不断更新粒子的速 度和位置，判断在给定的最大次数或最小误差要求的条件下，停止迭代，输出权值和偏置， 利用输出的权值和偏置对神经网络进行训练对测试集进行预测。
[0023] 优选地，所述第四步中，将粒子群优化算法神经网络模型预测结果与实际风速进 行对比，包括风速幅值、自相关函数，计算预测结果的误差指标，包括平均绝对误差、均方根 误差以及相关系数，评价预测方法的精度，同时利用未优化的神经网络对同样脉动风速进 行预测，分析比较两种方法的性能。
[0024] 本发明粒子群优化神经网络的脉动风速预测方法具有如下优点：本发明从学习规 贝1J、连接权值出发，使用粒子群优化算法对神经网络神经进行优化，将神经网络的权值和偏 置映射为粒子，并选择神经网络预测模型的均方误差作为适应度函数。数值模拟结果显示， 通过粒子群优化算法优化的神经网络预测模型的预测精度可以得到进一步提高，说明将粒 子群优化算法应用于神经网络预测模型的优化是可行的，可以作为脉动风速预测的一种有 效地方法。
【附图说明】
[0025] 图1是沿地面垂直方向20米处脉动风速模拟样本示意图；
[0026] 图2是沿地面垂直方向50米处脉动风速模拟样本示意图；
[0027] 图3为PSO-BP神经网络的预测算法的流程图；
[0028] 图4是20米PSO-BP、BP预测风速与模拟风速对比示意图；
[0029] 图5是20米PSO-BP、BP预测风速与模拟风速自相关函数对比示意图；
[0030] 图6是50米PSO-BP、BP预测风速与模拟风速对比示意图；
[0031] 图7是50米PSO-BP、BP预测风速与模拟风速自相关函数对比示意图。
【具体实施方式】
[0032] 本发明的构思如下：通过ARMA数值模拟出脉动风速作为有限样本，建立BP神经网 络预测模型，同时将权值和偏置作为PSO算法的粒子进行优化，最终建立基于粒子群优化 的BP神经网络的脉动风速预测模型。
[0033] 以下结合附图采用本发明对单点脉动风速预测作进一步详细说明，步骤如下：
[0034] 第一步，利用自回归滑动平均模型模拟生成垂直空间点脉动风速样本，将每一个 空间点的脉动风速样本分为训练集、测试集两部分，对其分别进行归一化处理；对数据归一 化处理后，目的是为了取消各维数据之间的数量级差别，从而避免因为输入输出数据的数 量级差别大而造成网络预测误差大的情况，取嵌入维数k = 10对进行样本数据进行相空间 重构；确定单点脉动风速样本的自回归滑动平均模型各参数，自回归滑动平均模型的自回 归阶数P = 4,滑动回归阶数q = 1。模拟某100米的超高层建筑，沿高度方向取每隔10米 的点作为各模拟风速点。其他相关参数见表1 :
[0035] 表1相关模拟参数表
[0037] 模拟功率谱采用Kaimal谱，只考虑高度方向的空间相关性。模拟生成20、50米脉 动风速样本分别见图1、图2。
[0038] 第一步中，自回归滑动平均模型模拟m维脉动风速表示为下式（1):
[0040] 式中，U(t)为脉动风速;A1, 分别是mXm阶自回归模型和滑动平均模型的系数 矩阵；X (t)为mX 1阶正态分布白噪声序列；P为自回归阶数、q为滑动回归阶数，t为时间。
[0041] 将得到的脉动风速按式（2)进行归一化处理：
[0043] 式中，K为归一化后脉动风速，Y1为实际脉动风速样本，y _为实际脉动风速最大 值，卩_实际脉动风速最小值。
[0044] 第二步，取采