发明名称 基于C‑C与ELM的短期风速预测方法
摘要 本发明提供了一种基于C‑C与ELM的短期风速预测方法,该方法考虑到原始的单变量时间序列风速数据间的最大动态演化信息,利用C‑C法进行相空间重构以确定风速时间序列间的相关关系,得到风速变化特征的最佳嵌入维数和延迟时间,接着采用学速度快、收敛精度高的极限学机方法进行风速的快速预测,通过与神经网络和支持向量机方法预测结果的对比可知,本发明提出的预测方法不仅提高了2小时的短期风速预测精度,且所用方法的学及预测时间短,适用于风电场短期风速的现场快速预测。
申请公布号 CN106126906A 申请公布日期 2016.11.16
申请号 CN201610456924.1 申请日期 2016.06.22
申请人 重庆科技学院 发明人 苏盈盈;刘兴华;刘君;李太福;黄利军
分类号 G06F19/00(2011.01)I;G06N3/08(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 重庆为信知识产权代理事务所(普通合伙) 50216 代理人 龙玉洪
主权项 一种基于C‑C与ELM的短期风速预测方法,其特征在于,包括如下步骤:S1:采集同一地区的风速数据,形成风速时间序列;S2:利用小波降噪方法对采集到的风速数据进行滤波降噪;S3:利用C‑C法对滤波降噪后的风速数据进行相空间重构:即采用C‑C法确定风速时间序列间的相关关系,得到风速变化特征的最佳嵌入维数m和延迟时间τ:设x={x<sub>i</sub>|i=1,2,…,N}为风速时间序列,以延迟时间τ,嵌入维数m重构相空间X={X<sub>j</sub>},其中X<sub>j</sub>={x<sub>j</sub>,x<sub>j+τ</sub>,…,x<sub>j+(m‑1)τ</sub>}(j=1,2,…,M)为重构后相空间中的第j个点,M为相空间重构后总的点数;S4:运用极限学习机、BP神经网络及支持向量机对重构后的风速数据建模,并进行短期风速预测:其中,运用极限学习机对重构后的风速数据建模的具体方法为:设C‑C法得到的M个新样本为(X<sub>j</sub>,T<sub>j</sub>),其中,X<sub>j</sub>=[x<sub>j1</sub>,x<sub>j2</sub>,...,x<sub>jn</sub>]<sup>T</sup>∈R<sup>n</sup>为输入变量,T<sub>j</sub>=[t<sub>j1</sub>,t<sub>j2</sub>,...,t<sub>jm</sub>]<sup>T</sup>∈R<sup>m</sup>为输出变量,j=1,2,…,M,对于一个带有L个隐层节点的标准前馈神经网,极限学习机的数学模型为:<img file="FDA0001025322570000011.GIF" wi="1261" he="150" />式中:w<sub>i</sub>=[w<sub>i1</sub>,w<sub>i2</sub>,…,w<sub>in</sub>]<sup>T</sup>为第i个输入节点和隐层节点的连接权重,β<sub>i</sub>=[β<sub>i1</sub>,β<sub>i2</sub>,…,β<sub>im</sub>]<sup>T</sup>为第i个隐层节点和输出节点的连接权重,b<sub>i</sub>为第i隐层节点的阈值,w<sub>i</sub>·X<sub>j</sub>为w<sub>i</sub>与X<sub>j</sub>的内积,g(X<sub>j</sub>)为激活函数,y<sub>j</sub>为第j个输入样本的输出值,对于一个给定的训练集{(X<sub>j</sub>,T<sub>j</sub>)|X<sub>j</sub>∈R<sup>n</sup>,T<sub>j</sub>∈R<sup>m</sup>,j=1,…,M}:S41:随机生成输入权重w<sub>i</sub>和阈值b<sub>i</sub>,i=1,…,L;S42:计算隐层输出矩阵H,其中H=g(w<sub>i</sub>·X<sub>j</sub>+b<sub>i</sub>),g(X<sub>j</sub>)为激活函数;S43:计算输出权重<img file="FDA0001025322570000021.GIF" wi="739" he="167" />式中,T为实际的风速历史数据,H<sup>T</sup>为输出矩阵H的转置,<img file="FDA0001025322570000022.GIF" wi="75" he="62" />为输出矩阵H的广义逆矩阵,<img file="FDA0001025322570000023.GIF" wi="54" he="134" />为小正数;S44:判断残差||E||是否大于等于预设精度ε,且隐层节点个数L小于隐层节点个数的最大值L<sub>max</sub>,如果是,则进入步骤S45,否则,进入步骤S5:S45:隐层节点个数L=L+1;S46:对于新隐层节点个数L,随机生成输入权重<img file="FDA0001025322570000024.GIF" wi="91" he="63" />阈值b<sub>L</sub>;S47:计算新隐层节点的输出权重:<img file="FDA0001025322570000025.GIF" wi="374" he="158" />S48:计算新隐层节点L的残差:E=E‑β<sub>L</sub>·H<sub>L</sub>,并进入步骤S44;S5:对比步骤S4中三种预测模型的性能指标,找到最佳短期风速预测模型。
地址 401331 重庆市沙坪坝区虎溪大学城重庆科技学院
您可能感兴趣的专利