[發明專利]基于單純形法的核電機組原動機及其調速器參數辨識的方法

有效

• 申請號: 201210252089.1
• 申請日: 2012-07-19
• 公開（公告）號: CN102799785A
• 公開（公告）日: 2012-11-28
• 發明（設計）人: 譚金;鄧少翔;金格;蔡筍;馮永新;鄧小文
• 申請（專利權）人: 廣東電網公司電力科學研究院
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 廣州知友專利商標代理有限公司 44104
• 代理人: 周克佑

1.一種基于單純形法的核電機組原動機及其調速器參數辨識的方法，包括以下步驟：

S1給定α₀,λ，μ，ε,N,K,K₁＝0，

其中:α₀為需要辨識參數的一組估計值，μ為擴張因子，μ為壓縮因子，ε為收斂標準，N等于需要辨識參數個數，K為最大迭代次數；

S2確定初始單純形，計算α_i＝α₀+h*e_i,i＝1,2,...n????（2）；

S3目標函數式（1）,根據目標函數計算各組參數對應的平均偏差，計算C_i=Q(α_i)（i=1,2,…,N）；

S4找出α_H,α_L,α_G；

S5迭代次數加1，K₁=K₁+1；

S6判斷|C_H-C_L|<εC_L?，若成立，輸出成功的α_L，C_L，結束；

S7判斷K₁>K?若成立，輸出失敗的α_L，CL，結束；或：返回步驟S1重新給定α₀,λ，μ，ε,N,K,K₁=0，重新進行參數辨識;

S8計算反射點及C_R=Q(α_R)；

S9判斷C_R<C_G?若不成立轉至步驟S15；

S10計算并判斷(1-μ)*C_H+μ*C_R<C_L?，若成立轉至步驟S13；

S11令αS=αRCS=CR;]]>

S12取代最壞點αH=αSCH=CS,]]>返回S4；

S13計算αE=(1-μ)*αH+μ*αRCE=Q(αE);]]>

S14判斷C_E<C_R?，若成立，轉至步驟S11，否則轉至步驟S12；

S15?R、H點互換；

S16計算α_s＝(1-λ)*α_H+λ*α_R及CS；

S17判斷C_S<C_G?，若成立轉至步驟S12；

S18計算α_i＝(α_L+α_i)/2（i=1,2,…,N），轉至步驟S3；

所述的目標函數定義為仿真曲線與實測曲線的平均偏差除以仿真期間實測功率與其初始值的最大偏差的絕對值，即：

e=(1T∫0T|Psim-Pmes|dt)/|ΔPmax|---(1)]]>

其中T為仿真時間長度，P_sim為仿真值，P_mes為實測值，|ΔP_max|為實測功率與其初始值的最大偏差的絕對值。

2.根據權利要求1所述的基于單純形法的核電機組原動機及其調速器參數辨識的方法，其特征是：

所述的步驟S1中的估計值取核電機組原動機及其調速器模型的典型值。