发明名称 锅炉汽包水位模型辨识的方法
摘要 一种计算机应用技术领域的锅炉汽包水位模型辨识的方法。本发明在现有的化工锅炉汽包水位控制系统中添加带滞环的偏置继电辨识模块,通过对系统在继电特性作用下的周期振荡极限环进行分析,获得汽包水位模型中四个参数的辨识公式。本发明针对化工过程中常见的汽包水位问题,提出一种新型的在线精确继电辨识方法,能在不需要任何时滞和静态增益等先验信息的情况下,根据时域分析得到的准确解析表达式精确计算出四个待辨识参数,从而得到逆响应的时滞模型与实际系统相匹配。
申请公布号 CN1858489A 申请公布日期 2006.11.08
申请号 CN200610026680.X 申请日期 2006.05.18
申请人 上海交通大学 发明人 顾诞英;张卫东;王萍;蔡云泽;欧林林
分类号 F22D5/26(2006.01);F22B37/38(2006.01);G05D9/00(2006.01) 主分类号 F22D5/26(2006.01)
代理机构 上海交达专利事务所 代理人 王锡麟;王桂忠
主权项 1、一种锅炉汽包水位模型辨识的方法,其特征在于,包括如下步骤:步骤一:偏置继电器输出值μ<sub>+</sub>,μ<sub>-</sub>和滞环ε已知参数送到存储单元RAM中;步骤二:通过继电辨识测试,上位机记录辨识过程的输入输出函数u和y,读取多个稳定振荡极限环周期参数的平均值,即周期P<sub>u1</sub>,P<sub>u2</sub>,幅值A<sub>p</sub>,A<sub>d</sub>,时间间隔<img file="A2006100266800002C1.GIF" wi="58" he="55" />和<img file="A2006100266800002C2.GIF" wi="63" he="55" />,将各参数送到存储单元RAM中;步骤三:调用事先编制好的算法程序解析辨识出系统模型的四个参数,具体算法步骤如下:A)读取内存中的μ<sub>+</sub>,μ<sub>-</sub>,ε,P<sub>u1</sub>,P<sub>u2</sub>,<img file="A2006100266800002C3.GIF" wi="57" he="56" />和<img file="A2006100266800002C4.GIF" wi="63" he="57" />代入方程(1)式,计算待辨识的时间常数τ<sub>2</sub>:<maths num="001"><![CDATA[ <math><mrow><mo>-</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>&mu;</mi><mo>+</mo></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow><mrow><mo>-</mo><msub><mi>&mu;</mi><mo>-</mo></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>t</mi><mo>~</mo></mover><mi>min</mi></msub><mo>-</mo><msub><mover><mi>t</mi><mo>~</mo></mover><mi>max</mi></msub></mrow></math>]]></maths>B)将τ<sub>2</sub>代入方程组(2)、(3)和(4),读取内存中的A<sub>p</sub>和A<sub>d</sub>,计算待辨识的时间常数τ<sub>1</sub>:<maths num="002"><![CDATA[ <math><mrow><mi>&theta;</mi><mo>=</mo><msub><mover><mi>t</mi><mo>~</mo></mover><mi>min</mi></msub><mo>+</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>&mu;</mi><mo>+</mo></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="003"><![CDATA[ <math><mrow><msub><mi>K</mi><mi>p</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>A</mi><mi>p</mi></msub><mo>-</mo><msub><mi>A</mi><mi>d</mi></msub></mrow><mrow><mrow><mo>(</mo><msub><mi>&mu;</mi><mo>+</mo></msub><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>+</mo><mn>2</mn><mi>&mu;</mi><msub><mi>&tau;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>[</mo><msub><mi>&mu;</mi><mo>-</mo></msub><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><msub><mi>&mu;</mi><mo>-</mo></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><msub><mi>&mu;</mi><mo>+</mo></msub><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>&mu;</mi><mo>+</mo></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="004"><![CDATA[ <math><mrow><msub><mrow><mo>-</mo><mi>K</mi></mrow><mi>p</mi></msub><mrow><mo>(</mo><msub><mi>&mu;</mi><mo>-</mo></msub><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>+</mo><mn>2</mn><mi>&mu;&theta;</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>K</mi><mi>p</mi></msub><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>[</mo><mn>1</mn><mo>-</mo><mfrac><msup><mi>e</mi><mrow><mi>&theta;</mi><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mrow><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup></mrow></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>+</mo><msup><mi>e</mi><mrow><msub><mrow><mo>-</mo><mi>P</mi></mrow><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>-</mo><mn>2</mn><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>&tau;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow><mo>]</mo><mo>=</mo><mn>2</mn><mi>&epsiv;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>C)将τ<sub>1</sub>和τ<sub>2</sub>代入方程(2)计算出待辨识参数θ;D)将τ<sub>1</sub>和τ<sub>2</sub>代入方程(3)计算出待辨识参数K<sub>p</sub>;在以上各式中,K<sub>p</sub>——对象静态增益,τ<sub>1</sub>——对象时间常数一,τ<sub>2</sub>——对象时间常数二,θ——对象时滞,μ<sub>+</sub>——继电器输出高电平,μ<sub>-</sub>——继电器输出低电平,ε——偏置继电器的滞环,P<sub>u1</sub>——继电器输出高电平的维持时间,P<sub>u2</sub>——继电器输出低电平的维持时间,A<sub>p</sub>——对象输出振荡曲线的正峰值,A<sub>d</sub>——对象输出振荡曲线的负峰值,<img file="A2006100266800003C1.GIF" wi="61" he="54" />——从继电器跳变到对象达到正峰值的时间间隔,<img file="A2006100266800003C2.GIF" wi="58" he="55" />——从继电器跳变到对象达到负峰值的时间间隔。
地址 200240上海市闵行区东川路800号