[发明专利]基于单纯形法的核电机组原动机及其调速器参数辨识的方法

摘要

本发明提供了一种基于单纯形法的核电机组原动机及其调速器参数辨识方法，其特点在于将单纯形法用于核电机组原动机及其调速器参数辨识，可减少计算工作量，较快收敛，并可避开对目标函数求导等问题，得到目标函数的解析形式。所述的基于单纯形法的核电机组原动机及其调速器参数辨识方法对控制参数寻优收敛比较快，计算工作量小，能有效辨识实际核电机组原动机及其调速器的参数，具有重要的工程应用价值。

主权项

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 = α R C S = C R ; ]]>S12取代最坏点 α H = α S C H = C S , ]]>返回S4；S13计算 α E = ( 1 - μ ) * α H + μ * α R C E = Q ( α E ) ; ]]>S14判断C_E<C_R，若成立，转至步骤S11，否则转至步骤S12；S15R、H点互换；S16计算α_s＝(1-λ)*α_H+λ*α_R及CS；S17判断C_S<C_G，若成立转至步骤S12；S18计算α_i＝(α_L+α_i)/2（i=1,2,…,N），转至步骤S3；所述的目标函数定义为仿真曲线与实测曲线的平均偏差除以仿真期间实测功率与其初始值的最大偏差的绝对值，即： e = ( 1 T ∫ 0 T | P sim - P mes | dt ) / | ΔP max | - - - ( 1 ) ]]>其中T为仿真时间长度，P_sim为仿真值，P_mes为实测值，|ΔP_max|为实测功率与其初始值的最大偏差的绝对值。