基于Morlet小波核的LSSVM脉动风速预测方法
【技术领域】
[0001] 本发明涉及一种采用Morlet小波函数构造小波核的最小二乘支持向量机的单点 脉动风速预测方法，具体的说是一种基于Morlet小波核的LSSVM脉动风速预测方法。
【背景技术】
[0002] 支持向量机出色的学习性能，尤其针对于小样本问题，一直是机器学习以及数据 挖据算法的研究热点。核函数在支持向量机中是至关重要的，它的引入极大地提高了学习 机器的非线性处理能力，保持了学习机器在高维空间中的内在线性，使得学习的过程容易 得到控制。显然支持向量机的性能在很大程度上取决于核函数的好坏，因此近年来关于支 持向量机的研究大部分都集中在支持向量机核函数的研究。目前国内外关于核函数的研究 主要可概括为以下几个方面：（1)组合现有的核函数构造出新的核函数；(2)改进现有的核 函数设计新的核函数；（3)核函数参数优化问题；(4)针对参数优化的模型选择问题。
[0003] 小波分析是近些年来发展起来的一种数学分析方法，它的特点是小波变换可以同 时在时域和频域都局有较好的分析效果，而且具备稀疏变换和多分辨分析的功能，而稀疏 变化的核函数有助于提高模型的精度和迭代的收敛速度；同时如果对平滑函数缺乏先验知 识，多尺度插值方法是最好的。这些小波理论独有的特点使其成为了解决非线性问题的研 究热点，而且已经出现了一些使用小波分析进行非线性建模的方法。常用高斯（RBF)核函 数对边界处信号逼近和多尺度信号逼近性能不是很好，而基于小波理论构造的核函数能够 弥补传统核函数在逼近性能方面的不足，有效地提升支持向量机的推广能力。由于小波核 具有局部化、多层次、多分辨的优点，同时可以进行多尺度分析，因此可以改进RBF核函数 在回归估计方面的不足。但是多尺度分析中的尺度参数，当所选的尺度参数较大时，会引起 算法的复杂化，影响算法的执行效率；反之如果尺度参数较小，必然会影响多尺度分析的效 果，继而影响算法的性能。
[0004] Morlet子波核函数不仅具有非线性映射的特征而且也继承了小波分析的多尺度 分析和核稀疏变化特性，能够对信号的逐级精细描述。因而采用子波核函数的小波支持向 量机能够以较高的精度逼近任意函数，这是传统的核函数所不具备的。
【发明内容】
[0005] 本发明的目的在于提供一种基于Morlet小波核的LSSVM脉动风速预测方法，其利 用ARM (Auto-Regressive and Moving Average，自回归滑动平均）模拟脉动风速样本，基 于Morlet母小波函数构造Morlet核函数，使该核函数具备多尺度分析和系数变换的特性， 建立Morlet核的最小二乘支持向量机（LSSVM)的模型，改进传统核函数对边界处信号逼近 和多尺度信号逼近性能不是很好的缺点，利用该模型对单点风速进行预测。计算实际风速 与预测风速的平均绝对误差（ME)、均方根误差（RMSE)以及相关系数（R)评价本方法的有 效性。
[0006] 本发明采用下述技术方案：一种基于Morlet小波核的LSSVM脉动风速预测方法， 其特征在于，其包括以下步骤：
[0007] 第一步：利用ARMA模型模拟生成垂直空间点脉动风速样本，将每个空间点的脉动 风速样本分为训练集、测试集两部分，采用Matlab对样本归一化处理；
[0008] 第二步：根据一维Morlet母小波函数，构造满足Mercer定理的Morlet小波核 函数，使该核函数具备多尺度分析和系数变换的特性，建立基于Morlet核函数的LSSVM模 型；
[0009] 第三步：引入PSO优化方法，对Morlet核函数的参数优化：小波核函数参数1、惩 罚参数c进行寻优，确定最优模型参数，利用PSO优化后的Morlet核函数将脉动风速训练 样本变换成为核函数矩阵，映射到高维特征空间，得到脉动风速训练样本的非线性模型，利 用此模型对脉动风速测试样本进行预测；
[0010] 第四步：将测试样本和预测的脉动风速结果对比，计算预测风速与实际风速的平 均绝对误差、均方根误差以及相关系数。
[0011] 优选地，所述第一步中，ARMA模型模拟m维脉动风速表示为下式：
[0013] 式中，U(t)为脉动风速；别是mXm阶AR和MA模型的系数矩阵；X(t)为 mX 1阶正态分布白噪声序列；P为自回归阶数、q为滑动回归阶数。
[0014] 优选地，所述第二步中，给定一维母小波函数一Mor I et小波，并证明该小波满足允 许支持向量机核函数条件，如下式：
[0016] 因为该小波为复数形式，所以取实数Morlet小波函数为如下式：
[0018] 定义Morlet小波核函数为如下式：
[0020] 优选地，所述第三步中，设置粒子群规模m = 30,随机产生核参数的初始位置，确 定待优化参数的范围，并设置最大迭代速度；通过迭代优化不断更新粒子的速度和位置，最 终根据终止迭代次数或适应度条件确定最优参数，建立Morlet小波核的LSSVM模型。
[0021] 本发明LSSVM脉动风速预测方法具有如下优点：小波函数具有稀疏变化和多尺度 性质，而稀疏变化的核函数有助于提高模型的精度和迭代的收敛速度；同时如果对平滑函 数缺乏先验知识，多尺度插值方法是最好的。小波核函数能够对信号的逐级精细描述，因而 采用子波核函数的小波支持向量机能够以较高的精度逼近任意函数，这是传统的核函数所 不具备的。根据运行结果表明，基于多核的LSSVM方法预测得到的脉动风速与实际脉动风 速吻合很好，可以作为脉动风速预测的一种有效方法。
【附图说明】
[0022] 图1是30米高度处脉动风速模拟样本示意图；
[0023] 图2是基于优化Morlet核的LSSVM脉动风速预测方法设计框架图示意图；
[0024] 图3是30米优化Morlet小波核LSSVM预测风速与实际风速对比示意图；
[0025] 图4是30米优化Morlet小波核LSSVM预测风速与实际风速自相关函数对比示意 图；
[0026] 图5是30米优化Morlet小波核LSSVM预测风速与实际风速功率谱密度函数对比 示意图。
【具体实施方式】
[0027] 本发明的构思如下：考虑到小波具有稀疏变化和多尺度分析的特征，而稀疏变化 的核函数有助于提高模型的精度和迭代的收敛速度；同时如果对平滑函数缺乏先验知识， 多尺度插值方法是最好的。因此在小波核函数的基础上又提出了多尺度小波核函数，进一 步提高了核函数的性能，同时针对多尺度核函数有尺度选取的问题，本发明又提出了 PSO 来优化多尺度核函数中尺度的选择。根据Mercer定理构造出新的Morlet小波核函数。改 进常用RBF核函数对边界处信号逼近和多尺度信号逼近性能不是很好的缺点。而基于小波 理论构造的核函数能够弥补传统核函数在逼近性能方面的不足，有效地提升支持向量机的 推广能力。由于小波核具有局部化、多层次、多分辨的优点，同时可以进行多尺度分析，因此 可以改进RBF核函数在回归估计方面的不足。利用Morlet核函数LSSVM模型对脉动风速 进行预测，验证结果的有效性。
[0028] 以下结合附图采用本发明对单点脉动风速预测作进一步详细说明，本发明基于 Morlet小波核的LSSVM脉动风速预测方法包括如下步骤：
[0029] 第一步，利用ARMA(自回归滑动）模型模拟生成空间点脉动风速样本，将每个空间 点的脉动风速样本分为训练集、测试集两部分，采用Matlab对样本归一化处理；
[0030] 所述第一步中，ARMA模型模拟m维脉动风速表示为下式（1):
[0032] 式⑴中，U(t)为脉动风速；A1JAv别是mXm阶AR和MA模型的系数矩阵；X(t) 为mX 1阶正态分布白噪声序列；P为自回归阶数，q为滑动回归阶数。
[0033] ARMA模型的自回归阶数p = 4,滑动回归阶数q = 1。模拟某200米的超高层建 筑，沿高度方向取每隔10米的点作为各模拟风速点，其他相关参数见表1 :
[0034] 表1相关模拟参数表
[0036] 模拟功率谱采用Kaimal谱，只考虑高度方向的空间相关性。模拟生成30米脉动 风速样本分别见图1。
[0037] 原始样本数据进行相应的预处理有利于加快模型的样本训练速度和收敛速度，提 高预测精度。本发明采用归一化处理方法，如下式（2):
[0039] 式中，为归一化后脉动风速，Y1为实际脉动风速样本，y_为实际脉动风速最大 值，卩_实际脉动风速最小值。
[0040] 第二步，根据一维Morlet母小波函数，构造满足Mercer定理的Morlet小波核函 数，使该核函数具备多尺度分析和系数变换的特性，建立基于Morlet核函数的LSSVM模型。
[0041] 所述第二步中，针对SVM中的不等式约束，Suykens根据正则化理论改变标准SVM 的约束条件和风险函数：将不敏感损失函数被误差的二次平方项代替作为损失函数，不等 式约束条件转变成等式约束条件。因此LSSVM将求解二次规划问题转化成求