[發明專利]RBF粒子群優化的工業熔融指數軟測量儀表和方法

有效

• 申請號: 201310433059.5
• 申請日: 2013-09-22
• 公開（公告）號: CN103674778A
• 公開（公告）日: 2014-03-26
• 發明（設計）人: 劉興高;張明明;李見會
• 申請（專利權）人: 浙江大學
• 主分類號: G01N11/00
• 分類號: G01N11/00
• 代理公司: 杭州求是專利事務所有限公司 33200
• 代理人: 周烽

1.一種RBF粒子群優化的工業熔融指數軟測量軟測量儀表，包括用于測量易測變量的現場智能儀表、用于測量操作變量的控制站、存放數據的DCS數據庫以及熔融指數軟測量值顯示儀，所述現場智能儀表、控制站與DCS數據庫連接，其特征在于：所述軟測量儀表還包括RBF粒子群優化的工業熔融指數軟測量模型，所述DCS數據庫與所述RBF粒子群優化的工業熔融指數軟測量模型的輸入端連接，所述RBF粒子群優化的工業熔融指數軟測量模型的輸出端與熔融指數軟測量值顯示儀連接，所述RBF粒子群優化的工業熔融指數軟測量模型包括：數據預處理模塊，用于將從DCS數據庫輸入的模型訓練樣本進行預處理，對訓練樣本中心化，即減去樣本的平均值，然后對其進行標準化：

計算均值：TX‾=1NΣi=1NTXi---(1)]]>

計算方差：σx2=1N-1Σi=1N(TXi-TX‾)---(2)]]>

標準化：X=TX-TX‾σx---(3)]]>

其中，TX_i為第i個訓練樣本，N為訓練樣本數，為訓練樣本的均值，X為標準化后的訓練樣本。σ_x表示訓練樣本的標準差，σ²_x表示訓練樣本的方差。

模糊方程模塊，對從數據預處理模塊傳過來的標準化后的訓練樣本X，進行模糊化。設模糊方程系統中有c^*個模糊群，模糊群k、j的中心分別為v_k、v_j，則標準化后的第i個訓練樣本X_i對于模糊群k的隸屬度μ_ik為：

μik=(Σj=1c*(||Xi-vk||||Xi-vj||)2n-1)-1---(4)]]>

式中，n為模糊分類過程中需要的分塊矩陣指數，通常取作2，||·||為范數表達式。

使用以上隸屬度值或者它的變形以獲得新的輸入矩陣，對于模糊群k，其輸入矩陣變形為：

Φ_ik(X_i,μ_ik)=[1?func(μ_ik)?X_i]????（5）

其中func(μ_ik)為隸屬度值μ_ik的變形函數，一般取exp(μ_ik)等，Φ_ik(X_i,μ_ik)表示第i個輸入變量X_i及其模糊群k的隸屬度μ_ik所對應的新的輸入矩陣。

RBF神經網絡作為模糊方程系統的局部方程，對每個模糊群進行優化擬合。設第k個RBF神經網絡模糊方程的的輸出是：

y^ik=ΣlwlkΦlk(||Xi-Clk||)---(6)]]>

其中C_lk和w_lk是第k個RBF神經網絡模糊方程第l個結點的中心和輸出權值，Φ_lk(||X_i-C_lk||)是第k個RBF神經網絡模糊方程第l個結點的徑向基函數，由下式確定：

Φlk(||Xi-Clk||)=exp(-(||Xi-Clk||)22×σlk)---(7)]]>

其中σ_lk是對應的徑向基函數的高斯寬度，||·||為范數表達式。最后，由反模糊方法中的重心法得到最后的模糊方程系統的輸出：

y^i=Σk=1c*μiky^ikΣk=1c*μik---(8)]]>

式中，為第k個RBF神經網絡模糊方程的的輸出。

粒子群算法優化模塊，用于采用粒子群算法對模糊方程中RBF神經網絡局部方程的C_pk、σ_pk、w_pk進行優化，具體實現步驟如下：

①確定粒子群的優化參數為RBF神經網絡局部方程的C_lk、σ_lk、w_lk、粒子群個體數目popsize、最大循環尋優次數iter_max、第p個粒子的初始位置r_p、初始速度v_p、局部最優值Lbest_p以及整個粒子群的全局最優值Gbest。

②設定優化目標函數，將其轉換為適應度，對每個局部模糊方程進行評價；通過相應的誤差函數計算適應度函數，并認為誤差大的粒子適應度小，粒子p的適應度函數表示為：

f_p=1/(E_p+1)????（9）

式中，E_p是模糊方程系統的誤差函數，表示為：

Ep=1NΣi=1N(y^i-Oi)2---(10)]]>

式中，是模糊方程系統的預測輸出，O_i為模糊方程系統的目標輸出；

③按照如下公式，循環更新每個粒子的速度和位置，

v_p(iter+1)=ω×v_p(iter)+m₁a₁(Lbest_p-r_p(iter))+m₂a₂(Gbest-r_p(iter))

??????????????????????????????????（11）

r_p(iter+1)=r_p(iter)+v_p(iter+1)????（12）

式中，v_p表示更新粒子p的速度，r_p表示更新粒子p的位置，Lbest_p表示更新粒子p的個體最優值，Gbest即為整個粒子群的全局最優值，iter表示循環次數，ω是粒子群算法中的慣性權重，m₁、m₂是對應的加速系數，a₁、a₂是[0,1]之間的隨機數；

④對于粒子p，如果新的適應度大于原來的個體最優值，更新粒子的個體最優值：

Lbest_p=f_p????（13）

⑤如果粒子p的個體最優值Lbest_p大于原來的粒子群全局最優值Gbest，更新原來的粒子群全局最優值Gbest:

Gbest=Lbest_p????（14）

⑥判斷是否滿足性能要求，若是，結束尋優，得到一組優化的模糊方程的局部方程參數；否則返回步驟③，繼續迭代尋優，直至達到最大迭代次數iter_max。

所述RBF粒子群優化的工業熔融指數軟測量模型還包括：

模型更新模塊，用于模型的在線更新，定期將離線化驗數據輸入到訓練集中，更新模糊方程模型。

2.一種用如權利要求1所述的RBF粒子群優化的工業熔融指數軟測量儀表實現的軟測量方法，其特征在于：所述軟測量方法具體實現步驟如下：

1）、對丙烯聚合生產過程對象，根據工藝分析和操作分析，選擇操作變量和易測變量作為模型的輸入，操作變量和易測變量由DCS數據庫獲得；

2）、將從DCS數據庫輸入的模型訓練樣本進行預處理，對訓練樣本中心化，即減去樣本的平均值，然后對其進行標準化，使得其均值為0，方差為1。該處理采用以下算式過程來完成：

2.1）計算均值：TX‾=1NΣi=1NTXi---(1)]]>

2.2）計算方差：σx2=1N-1Σi=1N(TXi-TX‾)---(2)]]>

2.3）標準化：X=TX-TX‾σx---(3)]]>其中，TX_i為第i個訓練樣本，N為訓練樣本數，為訓練樣本的均值，X為標準化后的訓練樣本。σ_x表示訓練樣本的標準差，σ²_x表示訓練樣本的方差。

3）、對從數據預處理模塊傳過來的訓練樣本，進行模糊化。設模糊方程系統中有c^*個模糊群，模糊群k、j的中心分別為v_k、v_j，則標準化后的第i個訓練樣本X_i對于模糊群k的隸屬度μ_ik為：

μik=(Σj=1c*(||Xi-vk||||Xi-vj||)2n-1)-1---(4)]]>

式中，n為模糊分類過程中需要的分塊矩陣指數，通常取作2，||·||為范數表達式。

使用以上隸屬度值或者它的變形以獲得新的輸入矩陣，對于模糊群k，其輸入矩陣變形為：

Φ_ik(X_i,μ_ik)=[1?func(μ_ik)?X_i]????（5）

其中func(μ_ik)為隸屬度值μ_ik的變形函數，一般取exp(μ_ik)等，Φ_ik(X_i,μ_ik)表示第i個輸入變量X_i及其模糊群k的隸屬度μ_ik所對應的新的輸入矩陣。

RBF神經網絡作為模糊方程系統的局部方程，對每個模糊群進行優化擬合。設第k個RBF神經網絡模糊方程的的輸出是：

y^ik=ΣlwlkΦlk(||Xi-Clk||)---(6)]]>

其中C_lk和w_lk是第k個RBF神經網絡模糊方程第l個結點的中心和輸出權值，Φ_lk(||X_i-C_lk||)是第k個RBF神經網絡模糊方程第l個結點的徑向基函數，由下式確定：

Φlk(||Xi-Clk||)=exp(-(||Xi-Clk||)22×σlk)---(7)]]>

其中σ_lk是對應的徑向基函數的高斯寬度，||·||為范數表達式。最后，由反模糊方法中的重心法得到最后的模糊方程系統的輸出：

y^i=Σk=1c*μiky^ikΣk=1c*μik---(8)]]>

式中，為第k個RBF神經網絡模糊方程的的輸出。

4）、采用粒子群算法對模糊方程中RBF神經網絡局部方程的C_pk、σ_pk、w_pk進行優化，具體實現步驟如下：

①確定粒子群的優化參數為RBF神經網絡局部方程的C_lk、σ_lk、w_lk、粒子群個體數目popsize、最大循環尋優次數iter_max、第p個粒子的初始位置r_p、初始速度v_p、局部最優值Lbest_p以及整個粒子群的全局最優值Gbest。

②設定優化目標函數，將其轉換為適應度，對每個局部模糊方程進行評價；通過相應的誤差函數計算適應度函數，并認為誤差大的粒子適應度小，粒子p的適應度函數表示為：

f_p=1/(E_p+1)????（9）

式中，E_p是模糊方程系統的誤差函數，表示為：

Ep=1NΣi=1N(y^i-Oi)2---(10)]]>

式中，是模糊方程系統的預測輸出，O_i為模糊方程系統的目標輸出；

③按照如下公式，循環更新每個粒子的速度和位置，

v_p(iter+1)=ω×v_p(iter)+m₁a₁(Lbest_p-r_p(iter))+m₂a₂(Gbest-r_p(iter))

??????????????????????????????????（11）

r_p(iter+1)=r_p(iter)+v_p(iter+1)????（12）

式中，v_p表示更新粒子p的速度，r_p表示更新粒子p的位置，Lbest_p表示更新粒子p的個體最優值，Gbest即為整個粒子群的全局最優值，iter表示循環次數，ω是粒子群算法中的慣性權重，m₁、m₂是對應的加速系數，a₁、a₂是[0,1]之間的隨機數；

④對于粒子p，如果新的適應度大于原來的個體最優值，更新粒子的個體最優值：

Lbest_p=f_p????（13）

⑤如果粒子p的個體最優值Lbest_p大于原來的粒子群全局最優值Gbest，更新原來的粒子群全局最優值Gbest:

Gbest=Lbest_p????（14）

⑥判斷是否滿足性能要求，若是，結束尋優，得到一組優化的模糊方程的局部方程參數；否則返回步驟③，繼續迭代尋優，直至達到最大迭代次數iter_max。

所述軟測量方法還包括以下步驟：5）、定期將離線化驗數據輸入到訓練集中，更新模糊方程模型。