发明名称 一种短期风速预测方法
摘要 本发明是一种风电场风速时间序列预测方法,其特点是,包括的步骤有:利用风速采集仪器每隔一小时记录一次同一地区的风速数据,整理采集的原始风速数据,形成风速时间序列用于分析预测;运用局域均值分解算法对风速时间序列进行多尺度分解,分解得到多个独立分量;通过对各个独立分量的延迟时间和嵌入维数的计算,运用相空间重构理论对各个独立分量进行相空间重构;利用多核最小二乘支持向量机回归模型对相空间重构后的各个独立分量进行建模预测;将预测结果进行叠加得到最终的风速时间序列的预测结果。具有方法科学、风速时间序列的预测准确、可靠,应用价值高等优点。
申请公布号 CN102539822B 申请公布日期 2014.10.29
申请号 CN201110459050.2 申请日期 2011.12.31
申请人 东北电力大学 发明人 孙斌;姚海涛;李洪伟;周云龙
分类号 G01P5/00(2006.01)I;G06F19/00(2011.01)I 主分类号 G01P5/00(2006.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 一种短期风速预测方法,其特征是:它包括以下步骤:(1)利用风速采集仪器每隔一小时记录一次同一地区的风速数据,整理采集的原始风速数据,形成风速时间序列x用于分析预测;(2)运用快速独立分量分析算法对指定时间范围内的风速时间序列x进行多尺度分解,运用自相关法和假邻近法分别求取风速时间序列x的延迟时间和嵌入维数,运用相空间重构理论对指定时间范围内的风速时间序列x的分离信号分别进行相空间重构,得到K个重构后的分离信号y<sub>k</sub>,k=1,2,…,K;设x=(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>n</sub>)<sup>T</sup>为n个点的风速时间序列,每个风速值都是m个独立源信号的线性组合s=(s<sub>1</sub>,s<sub>2</sub>,…,s<sub>m</sub>)<sup>T</sup>,即<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mi>Hs</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>h</mi><mi>j</mi></msub><msub><mi>s</mi><mi>j</mi></msub><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000513694770000011.GIF" wi="1232" he="133" /></maths>式中:H=[h<sub>1</sub>,h<sub>2</sub>,…,h<sub>n</sub>]为m×n阶满秩源信号混合矩阵,其中n≥m,式(1)即为快速独立分量分析算法的信号混合模型,快速独立分量分析算法就是在源信号s和混合矩阵H都未知的情况下,构建一个分离矩阵W,从风速时间序列x中分离出相互统计独立的分量y,使各分量逼近源信号s,这样快速独立分量分析算法问题的求解就可表示为y=Wx=(y<sub>1</sub>,y<sub>2</sub>,…,y<sub>n</sub>)<sup>T</sup>≈(s<sub>1</sub>,s<sub>2</sub>,…,s<sub>n</sub>)<sup>T</sup>  (2)采用快速独立分量分析算法算法寻找最优分离矩阵W,具体步骤如下(a)对风速时间序列x做均值和白化处理,得到新的信号<img file="FDA0000513694770000012.GIF" wi="92" he="62" /><img file="FDA0000513694770000013.GIF" wi="52" he="62" />满足<img file="FDA0000513694770000014.GIF" wi="269" he="76" />I为单位矩阵,(b)任意选择具有单位方差的初始分离矩阵w<sub>i</sub>(0),(c)计算<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>E</mi><mo>[</mo><mover><mi>x</mi><mo>~</mo></mover><mi>G</mi><mrow><mo>(</mo><msubsup><mi>w</mi><mi>i</mi><mi>T</mi></msubsup><mover><mi>x</mi><mo>~</mo></mover><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mi>E</mi><mo>[</mo><mi>G</mi><mrow><mo>(</mo><msubsup><mi>w</mi><mi>i</mi><mi>T</mi></msubsup><mover><mi>x</mi><mo>~</mo></mover><mo>)</mo></mrow><mo>]</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000513694770000015.GIF" wi="902" he="76" /></maths>(d)归一化处理分离矩阵w<sub>i</sub>(n+1)=w<sub>i</sub>(n+1)/||w<sub>i</sub>(n+1)||<sub>2</sub>,(e)判定w<sub>i</sub>(n+1)是否收敛,若收敛则分离出一个独立分量<img file="FDA0000513694770000016.GIF" wi="266" he="76" />否则返回步骤(c),(f)判定信号<img file="FDA0000513694770000017.GIF" wi="45" he="61" />中的多个独立分量是否已经全部分量完毕,若没有则返回(b),否则分离过程结束;(3)利用最小二乘支持向量机回归模型对相空间重构后的各个分离信号进行建模预测;<img file="FDA0000513694770000018.GIF" wi="1101" he="74" />式中,<img file="FDA0000513694770000019.GIF" wi="78" he="64" />为非线性映射函数,将输入样本数据y<sub>k</sub>映射到高维特征空间,ω为权向量,b为偏置量,T表示矩阵的转置;(4)将各分离信号的预测值<img file="FDA0000513694770000021.GIF" wi="398" he="74" />进行叠加得到最终的风速时间序列x的预测结果<img file="FDA0000513694770000022.GIF" wi="67" he="66" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>y</mi><mo>~</mo></mover><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mover><mi>y</mi><mo>~</mo></mover><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000513694770000023.GIF" wi="1293" he="134" /></maths>
地址 132012 吉林省吉林市船营区长春路169号