发明名称 基于递推最小二乘法RLS的结晶器ARX模型辨识方法
摘要 本发明涉及基于递推最小二乘法RLS的结晶器ARX模型辨识方法,即:以结晶器油缸阀开度为输入u,以结晶器位置为输出y,在采样数据基础上建立结晶器ARX模型最小平方和指标函数,利用变量u和y获得中间迭代变量Pk和Lk,针对ARX模型最小平方和指标函数,通过Pk和Lk按照最小二乘法逐步迭代计算结晶器ARX模型未知参数;该方法用于钢铁冶金行业中连铸机结晶器的在线辨识与控制器参数调整。本发明能够在线利用采样数据快速、准确逼近模型未知参数全局最优解,从而避免传统的最小二乘法需要离线辨识ARX模型的弊病;能够在模型参数发生变化时及时提供准确的辨识参数,为设计性能优良的结晶器控制系统提供了科学的数学模型。
申请公布号 CN102591208B 申请公布日期 2014.05.07
申请号 CN201210066222.4 申请日期 2012.03.14
申请人 中冶南方工程技术有限公司 发明人 张华军;蔡炜;褚学征;陈方元;尉强;周登科
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 湖北武汉永嘉专利代理有限公司 42102 代理人 王守仁
主权项 1.一种基于递推最小二乘法RLS的结晶器ARX模型辨识方法,其特征是一种利用迭代最小二乘法进行结晶器ARX模型辨识的方法,具体是:以结晶器油缸阀开度为输入u,以结晶器位置为输出y,在采样数据基础上建立结晶器ARX模型最小平方和指标函数,利用变量u和y获得中间迭代变量Pk和Lk,针对ARX模型最小平方和指标函数,通过Pk和Lk按照最小二乘法逐步迭代计算结晶器ARX模型未知参数;该方法包括以下步骤:(1)采集输入输出数据:以结晶器油缸阀开度为输入u(t),以结晶器位置为输出y(t)采集N对数据样本Z<sup>N</sup>;(2)构建结晶器白噪声干扰下的ARX模型:该ARX模型为:A(q)y(t)=B(q)u(t)+e(t),其中A(q)=1+a<sub>1</sub>q<sup>-1</sup>+a<sub>2</sub>q<sup>-2</sup>+L+a<sub>na</sub>q<sup>-na</sup>,B(q)=b<sub>1</sub>q<sup>-1</sup>+b<sub>2</sub>q<sup>-2</sup>+L+b<sub>nb</sub>q<sup>-nb</sup>,q<sup>-1</sup>为后向移动算子,q为前向移动算子,na、nb为正实数,e(t)为高斯白噪声;(3)确定ARX模型待辨识参数:令θ=[a<sub>1</sub>  a<sub>2</sub>  L  a<sub>na</sub>  b<sub>1</sub>  b<sub>2</sub>  L  b<sub>nb</sub>]为ARX模型待辨识参数;(4)计算基于参数θ的模型输出预测值:令<img file="FDA00003584069500011.GIF" wi="158" he="86" />为基于参数θ的模型输出预测值,其中预测表达式为:<img file="FDA00003584069500012.GIF" wi="986" he="93" />其中<img file="FDA00003584069500013.GIF" wi="1651" he="93" />(5)确定ARX模型辨识过程的目标函数:令带有高斯白噪声的ARX模型辨识过程的目标函数为:<img file="FDA00003584069500014.GIF" wi="1175" he="143" />(6)计算R(t),并判断R(t)是否可逆,若不可逆则执行步骤(6),否则令t=t<sub>0</sub>,执行步骤(7);计算公式为:<img file="FDA00003584069500015.GIF" wi="502" he="158" />(7)计算F(t<sub>0</sub>):计算公式为:<img file="FDA00003584069500021.GIF" wi="515" he="165" />(9)令P(t<sub>0</sub>)=R(t<sub>0</sub>)<sup>-1</sup>,<maths num="0001"><![CDATA[<math><mrow><mover><mi>&theta;</mi><mo>^</mo></mover><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>*</mo><mi>F</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>(10)令<img file="FDA00003584069500023.GIF" wi="1408" he="166" />(11)则未知参数迭代公式为<img file="FDA00003584069500024.GIF" wi="922" he="105" />经过上述步骤,实现基于递推最小二乘法RLS的结晶器ARX模型的辨识。
地址 430223 湖北省武汉市东湖新技术开发区大学园路33号