[發(fā)明專利]一種中速磨制粉系統(tǒng)給煤量軟測量方法

有效

• 申請?zhí)枺?/td>: 201310662727.1
• 申請日: 2013-12-09
• 公開（公告）號: CN103699782B
• 公開（公告）日: 2017-02-01
• 發(fā)明（設計）人: 唐耀華;李文啟;郭為民;孫建華;李貴兵;高建紅;朱峰;段松濤;李炳楠;田亮;郭毅樂
• 申請（專利權）人: 國家電網公司;國網河南省電力公司電力科學研究院;大唐安陽發(fā)電有限責任公司;河南恩湃高科集團有限公司;華北電力大學
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 鄭州知己知識產權代理有限公司41132
• 代理人: 季發(fā)軍

1.一種中速磨制粉系統(tǒng)給煤量軟測量方法，其特征在于，包括以下步驟：

A：構建包含所有輔助變量的數據集，并對數據集內的數據依次進行預處理，分別計算入爐煤量與各輔助變量的相關系數；選擇相關系數中絕對值最大的n組輔助變量作為建立中速磨制粉系統(tǒng)給煤量預測模型的輔助變量；

B：根據選擇的n組輔助變量，運用多元線性回歸法建立中速磨制粉系統(tǒng)給煤量預測模型為Y=k₀+k₁X₁+k₂X₂+…+k_nX_n+ε；

其中，Y為制粉系統(tǒng)的給煤量；X₁、X₂…X_n分別為步驟A中確立的n組輔助變量；k₀、k₁、k₂…k_n是待定的回歸系數；ε為各種隨機因素對Y的影響總和。

C：利用最小二乘法參數估計法估計待定的回歸系數k₀、k₁、k₂…k_n；

D：將步驟C中求得的待定的回歸系數k₀、k₁、k₂…k_n，分別代入步驟B確定的中速磨制粉系統(tǒng)給煤量預測模型中，計算中速磨制粉系統(tǒng)給煤量。

2.根據權利要求1所述的中速磨制粉系統(tǒng)給煤量軟測量方法，其特征在于，所述的步驟A，包括以下步驟：

A1：構建含有N組數據的數據集，其中每組數據均包含所有的輔助變量，然后利用公式(1)求得每個輔助變量的數據均值；

x‾=1NΣi=1Nxi---(1);]]>

A2：利用公式（2）求得每個輔助變量的數據標準差σ；

σ2=1N-1Σi=1N(xi-x‾)2---(2);]]>

A3：利用公式（3），將數據集中的每個輔助變量的數據樣本x_i轉化為相對應的標準樣本值

x^i=xi-x‾σ---(3);]]>

A4：使用步驟A3中得到的每個輔助變量的標準樣本值計算入爐煤量與各輔助變量的相關系數，選擇相關系數中絕對值最大的n組輔助變量作為建立中速磨制粉系統(tǒng)給煤量預測模型的輔助變量；相關系數的計算公式為：

ρ=E((X-E(X))(Y-E(Y)))E(X2)-(E(X))2E(Y2)-(E(Y))2,]]>

其中，ρ為相關系數，E(·)表示求取變量的數學期望，X為輔助變量，Y為入爐煤量。

3.根據權利要求2所述的中速磨制粉系統(tǒng)給煤量軟測量方法，其特征在于，所述的步驟C，包括以下步驟：

C1：利用矩陣形式描述待定的回歸系數k₀、k₁、k₂…k_n：

Y=Y1Y2···Ym,X=1X11X12...X1n1X21X22...X2n···············1Xm1Xm2...Xmn;]]>

ϵ=ϵ1ϵ2···ϵm,K=k0k1k2···kn;]]>

C2：根據步驟C1可將步驟B中所建立的中速磨制粉系統(tǒng)給煤量預測模型Y=k₀+k₁X₁+k₂X₂+…+k_nX_n+ε簡化為Y=XK+ε，其中，Y為m維輸出向量，是制粉系統(tǒng)給煤量Y的矩陣表示形式；X為m×n維測量矩陣；K為n+1維參數向量；ε為m維噪聲向量；

C3：設表示k的最優(yōu)估值，是Y的最優(yōu)估值，則

其中，Y^=[Y^1,Y^2,...,Y^m];]]>

C4：設e_j為殘差，表示Y_j與之差，

其中，為參數向量；

C5：利用公式使殘差平方和取最小值，當滿足非奇異條件時，對J求的偏導數并令其等于0，即：

∂(J)∂(K^)=-2XTY+2XTXK^=0;]]>

化簡后可得待定的回歸系數：