[發明專利]一種農產品價格預測方法及裝置

在審

• 申請號: 201611256179.2
• 申請日: 2016-12-30
• 公開（公告）號: CN106651464A
• 公開（公告）日: 2017-05-10
• 發明（設計）人: 魏帥;王德運;郭海湘
• 申請（專利權）人: 中國地質大學(武漢)
• 主分類號: G06Q30/02
• 分類號: G06Q30/02;G06Q50/02
• 代理公司: 北京華沛德權律師事務所11302
• 代理人: 房德權

1.一種農產品價格預測方法，其特征在于，包括：

獲得第一類型農產品的M個價格預測樣本，M至少為百數量級；

通過VMD變分模態分解將所述M個價格預測樣本分解為k個具有不同振動頻率特征的價格VM變分模態分量，其中，k為大于1的整數；

輸入所述k個價格VM分量至經過螢火蟲算法優化的SVM預測模型中進行預測，以預測到與所述k個價格VM分量對應的k個VM分量預測值；

將所述k個VM分量預測值線性疊加，以得到對所述第一類型農產品的價格預測值。

2.如權利要求1所述的農產品價格預測方法，其特征在于，所述通過VMD變分模態分解將所述M個價格預測樣本分解為k個具有不同振動頻率特征的價格VM分量，包括：

步驟A1：初始化和n，其中，為所對應的傅里葉轉換，為{λ¹}所對應的傅里葉轉換，為初始化后中心頻率，u_k表示模態函數，λ為拉格朗日乘法算子，n代表迭代次數；

步驟A2：根據如下公式更新模態函數u_k：

${\widehat{\mathrm{\μ}}}_k^{n+1}\left(\mathrm{\ω}\right)=\frac{\hat{f}\left(\mathrm{\ω}\right)-\underset{i\≠k}{\mathrm{\Σ}}\widehat{\mathrm{\μ}}\left(\mathrm{\ω}\right)+\frac{\widehat{\mathrm{\λ}}\left(\mathrm{\ω}\right)}{2}}{1+2\mathrm{\α}{(\mathrm{\ω}-{\mathrm{\ω}}_k)}^2},$

以及根據如下公式更新中心頻率ω_k：

${\mathrm{\ω}}_k^{n+1}=\frac{{\∫}_0^{\mathrm{\∞}}\mathrm{\ω}{\left|{\hat{y}}_k^{n+1}\left(\mathrm{\ω}\right)\right|}^{2}d\mathrm{\ω}}{{\∫}_0^{\mathrm{\∞}}{\left|{\hat{y}}_k^{n+1}\left(\mathrm{\ω}\right)\right|}^{2}d\mathrm{\ω}},$

其中，為當前剩余量的維納濾波，為當前模態函數功率譜的重心，對應為f(t)、u_i(t)、λ(t)、的傅里葉轉換，f(t)為所述價格預測樣本，u_i(t)為第i個模態函數，λ(t)表示拉格朗日乘法算子，為u_k經過n+1次迭代后的值，1≤i≤k；

步驟A3：根據如下公式更新λ；

${\widehat{\mathrm{\λ}}}^{n+1}\left(\mathrm{\ω}\right)\←{\widehat{\mathrm{\λ}}}^n\left(\mathrm{\ω}\right)+\mathrm{\τ}\[\hat{f}\left(\mathrm{\ω}\right)-\underset{1}{\overset{k}{\mathrm{\Σ}}}{\widehat{\mathrm{\μ}}}_k^{n+1}\left(\mathrm{\ω}\right)\],$

其中，為λ(ω)經過n次迭代后的傅里葉轉換，τ為更新系數；

步驟A4：若則停止迭代，并將u_k的輸出作為k個價格VM分量，否則返回步驟A2。