1.一种基于变分异方差高斯过程回归的短期风速预测方法，其特征在于，包括如下步骤：

步骤一：收集风电场历史实测短期风速数据，根据历史实测风速数据建立风速时间序列，通过完全集合经验模态分解将其分解为K个本征模态函数分量IMFk(k＝1,2,3,...,K)和一个残差变量Rk；

步骤二：对于步骤一中每个本征模态函数分量IMFk(k＝1,2,3,...,K)和残差Rk，计算其偏自相关函数值，选择95％置信水平下显著的滞后时间序列作为输入变量，目标向量为输出变量，并将全部样本数据的前70％作为训练样本，剩余30％作为检验样本；

步骤三：对于步骤一中每个本征模态函数分量IMFk(k＝1,2,3,...,K)和残差Rk，建立变分异方差高斯过程回归模型，输出最优变分参数Λ和模型最优超参数所述变分异方差高斯过程回归模型是引入边缘变分方法估计异方差高斯过程回归的预测分布，所述异方差高斯过程记为G3，其由两个方差相同的高斯过程组成：第一个高斯过程模型称为G1，对无噪声目标输出进行建模；第二个高斯过程模型称为G2，对数噪声进行建模，即，z(xi)＝lg(r(xi))，该z过程由协方差函数kz控制，其超参数为θz；所述变分异方差高斯过程回归的步骤如下：

3.1)对于给定的训练集D＝{(xi,yi)(i＝1,2,...,n)}，应用第一个高斯过程G1预测输出向量y，构建GPR模型，其中xi表示d维输入向量，yi表示对应的输出向量，n表示训练样本的数量；

3.2)根据G1估计训练数据的噪声水平，即zi＝lg{var[yi,G(xi,X,y)]}，形成新的数据集D'＝{(xi,zi)(i＝1,2,...,n)}；

3.3)运用第二个高斯过程G2对数据集D'建模；

3.4)运用异方差高斯过程G3对数据集D建模，其中，用第二个高斯过程G2预测噪声水平ri，采用边缘变分方法估计异方差高斯过程回归的预测分布，称为变分异方差高斯过程；

3.5)若G3不收敛，令G1等于G3并返回3.2)，否则输出变分参数Λ和模型超参数步骤四：将步骤三所得的最优变分参数Λ和最优模型超参数 以及步骤二中的检验样本代入变分异方差高斯过程回归预测模型，得到检验阶段的预测值，并对所有子模态的预测结果进行求和重构，得到原始风速时间序列的预测值；

步骤五：计算预测值和实测值的均方根误差RMSE、平均绝对误差MAE、相关系数R和平均绝对百分比误差MAPE，评价所提出的风速预测模型的性能。

2.根据权利要求1中所述的基于变分异方差高斯过程回归的短期风速预测方法，其特征在于，在异方差高斯过程回归模型中，所述3.1)中输入向量X*的预测分布p(y*|X*,X,y)服从高斯分布，即：

p(y*|X*,X,y)～N(m*,∑*)其中，m*和∑*分别表示y*的期望和协方差，N表示正态分布，样本集的大小，p(y*|X*,X,y)为输入向量X*的预测分布。

3.根据权利要求2所述的基于变分异方差高斯过程回归的短期风速预测方法，其特征在于，所述异方差高斯过程回归模型相应的变分，即变分异方差高斯过程回归模型为：F(q(f),q(z))＝logp(y)‑KL(q(f)q(z)||p(f,z|y))式中，KL(.||.)表示kullback‑leibler散度，q(f)和q(z)为任意的变概率密度，p(y)表示异方差高斯过程回归的边际对数似然，p(f,z|y)表示潜在变量f(x)和误差噪声z(x)的联合概率分布。

4.根据权利要求3所述的基于变分异方差高斯过程回归的短期风速预测方法，其特征在于，将q(z)的分布限制为多元正态分布，即q(z)＝N(z|u,∑)，则异方差高斯过程回归模型的边缘变分MV边界可以写为：

式中，Kz表示过程z的核矩阵，其超参数为θz，0为N维全零矩阵，由于∫N(y|f,R)N(f|0,Kf)df＝N(y|0,Kf+R),其中R为主对角元素满足[R]ii＝exp([μ]i‑[∑]ii/2)的对角矩阵，则上式可进一步简化为：

式中，tr(·)表示矩阵的迹。

5.根据权利要求1所述的基于变分异方差高斯过程回归的短期风速预测方法，其特征在于，所述步骤五中四项统计指标计算公式如下：s o

式中，Yi和Yi分别表示第i个样本的模拟值和实测值； 和 分别表示样本数据集的平均模拟值和平均实测值；N表示样本集的大小。