基于单纯形法的核电机组原动机及其调速器参数辨识的方法技术

技术编号:8022588 阅读:202 留言:0更新日期:2012-11-29 04:48
本发明专利技术提供了一种基于单纯形法的核电机组原动机及其调速器参数辨识方法,其特点在于将单纯形法用于核电机组原动机及其调速器参数辨识,可减少计算工作量,较快收敛,并可避开对目标函数求导等问题,得到目标函数的解析形式。所述的基于单纯形法的核电机组原动机及其调速器参数辨识方法对控制参数寻优收敛比较快,计算工作量小,能有效辨识实际核电机组原动机及其调速器的参数,具有重要的工程应用价值。

【技术实现步骤摘要】

本专利技术涉及一种核电机组原动机及其调速器参数辨识的方法,尤其是涉及一种。
技术介绍
分析电力系统的基本工具是数字仿真软件,主流的电力系统分析软件都支持多种核电机组的原动机及其调速器模型,但在参数方面,均是未经实际验证的“典型值”,在研究机电暂态稳定时,可能不会带来很大的误差,但是,对于特高压电网组建成的大电网,原动机及其调速器的影响将逐步显现出来,因此有必要通过实测数据辨识出实际参数值代替目前的“典型值”。最小二乘法是一种经典有效的参数辨识方法,但是当实测数据量很大、待辨识参 数个数较多时,计算量将会非常大。最速下降法和共轭梯度法是单变量参数辨识的有效方法,但是这两种方法都要求计算目标函数的梯度,而实际上在许多寻优问题中常常得不到目标函数的解析形式,另外最速下降法在最优点附近往往收敛速度很慢,为获得较高的寻优精度,往往要在最优点附近搜索很多次。
技术实现思路
本专利技术所要解决的技术问题,就是提供一种,本专利技术不需要得到目标函数的导数,即便当实测数据量较大、待辨识参数个数较多时,计算量也较小,收敛速度较快,搜索次数较少,辨识出的参数能保证原动机及其调速系统模型能准确用于大电网。解决上述问题,本专利技术采用如下的技术方案一种,包括以下步骤SI 给定 a 0,X,ii,e,N,K,K1=O,其中为需要辨识参数的一组估计值,U为扩张因子,U为压缩因子,e为收敛标准,N等于需要辨识参数个数,K为最大迭代次数;S2确定初始单纯形,计算a i = ai = I, 2,…n (2);S3目标函数式(1),根据目标函数计算各组参数对应的平均偏差,计算Ci=Q(CIi)(i=l,2,...,N);S4 找出 aH,aL,aG;S5 迭代次数加 I,K1=Kfl ;S6判断I Ch-ClI〈 e Cl,若成立,输出成功的4,Q,结束;S7判断IPK 若成立,输出失败的a y CL,结束;或返回步骤SI重新给定a ^入,ii , e,N,K,K1=O,重新进行参数辨识;S8 计算反射点权利要求1.一种,包括以下步骤 51给定 α。,λ,μ,ε ,NjKjK1 = O, 其中为需要辨识参数的一组估计值,μ为扩张因子,μ为压缩因子,ε为收敛标准,N等于需要辨识参数个数,K为最大迭代次数; 52确定初始单纯形,计算 Cti=Ct 0+h*ei, i = I, 2,... η (2); S3目标函数式(I),根据目标函数计算各组参数对应的平均偏差,计算Ci=QUi)(i=l,2, -,N);S4 找出 α H, a L, a G ; S5迭代次数加1,K1=KJl ; S6判断I Ch-ClI <ε Cl ,若成立,输出成功的a L,CL,结束; S7判断K1X 若成立,输出失败的4,CL,结束;或返回步骤SI重新给定%,λ,μ,ε,N, K, K1=O,重新进行参数辨识; S8 计算反射点2.根据权利要求I所述的,其特征是 所述的步骤Si中的估计值取核电机组原动机及其调速器模型的典型值。全文摘要本专利技术提供了一种基于单纯形法的核电机组原动机及其调速器参数辨识方法,其特点在于将单纯形法用于核电机组原动机及其调速器参数辨识,可减少计算工作量,较快收敛,并可避开对目标函数求导等问题,得到目标函数的解析形式。所述的基于单纯形法的核电机组原动机及其调速器参数辨识方法对控制参数寻优收敛比较快,计算工作量小,能有效辨识实际核电机组原动机及其调速器的参数,具有重要的工程应用价值。文档编号G06F19/00GK102799785SQ201210252089公开日2012年11月28日 申请日期2012年7月19日 优先权日2012年7月19日专利技术者谭金, 邓少翔, 金格, 蔡笋, 冯永新, 邓小文 申请人:广东电网公司电力科学研究院本文档来自技高网...

【技术保护点】
一种基于单纯形法的核电机组原动机及其调速器参数辨识的方法,包括以下步骤:S1给定α0,λ,μ,ε,N,K,K1=0,其中:α0为需要辨识参数的一组估计值,μ为扩张因子,μ为压缩因子,ε为收敛标准,N等于需要辨识参数个数,K为最大迭代次数;S2确定初始单纯形,计算αi=α0+h*ei,i=1,2,...n????(2);S3目标函数式(1),根据目标函数计算各组参数对应的平均偏差,计算Ci=Q(αi)(i=1,2,…,N);S4找出αH,αL,αG;S5迭代次数加1,K1=K1+1;S6判断|CH?CL|K?若成立,输出失败的αL,CL,结束;或:返回步骤S1重新给定α0,λ,μ,ε,N,K,K1=0,重新进行参数辨识;S8计算反射点及CR=Q(αR);S9判断CRαS=αRCS=CR;S12取代最坏点αH=αSCH=CS,返回S4;S13计算αE=(1-μ)*αH+μ*αRCE=Q(αE);S14判断CEe=(1T∫0T|Psim-Pmes|dt)/|&Delta;Pmax|---(1)其中T为仿真时间长度,Psim为仿真值,Pmes为实测值,|ΔPmax|为实测功率与其初始值的最大偏差的绝对值。FDA00001907276700011.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员:谭金邓少翔金格蔡笋冯永新邓小文
申请(专利权)人:广东电网公司电力科学研究院
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1