发明名称 多变量时滞系统解耦控制器解析设计方法
摘要 一种工业过程控制技术领域的多变量时滞系统解耦控制器解析设计方法,步骤如下:1)首先对化工多变量时滞过程的传递函数矩阵辨识模型;2)设计用于保证控制器能实现的两个调节因子;3)设计解耦的控制器C为如上式形式,如果对象传递函数矩阵分解后G<SUB>O</SUB>(s)中仍然含有时滞项,将导致设计的控制器是无穷高阶,这时使用控制器降阶技术对该控制器进行降阶。本发明采用输入输出设计方法,设计方法容易理解,使用方便,不需要状态信息。本发明设计的控制器可以实现系统输出解耦,并且只需要调节有限的几个参数,就可以方便有效的定量调节系统性能和鲁棒性。
申请公布号 CN1794117A 申请公布日期 2006.06.28
申请号 CN200510112229.5 申请日期 2005.12.29
申请人 上海交通大学 发明人 张卫东;欧林林;刘媛媛;张彬;王萍
分类号 G05B13/00(2006.01);G05B13/04(2006.01);G05B15/02(2006.01) 主分类号 G05B13/00(2006.01)
代理机构 上海交达专利事务所 代理人 王锡麟;王桂忠
主权项 1、一种多变量时滞系统解耦控制器解析设计方法,其特征在于,具体步骤如下:1)首先对化工多变量时滞过程的传递函数矩阵辨识模型<img file="A2005101122290002C1.GIF" wi="601" he="226" />其中<maths num="001"><![CDATA[ <math><mrow><msub><mi>g</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>g</mi><mi>oij</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&theta;</mi><mi>ij</mi></msub><mi>s</mi></mrow></msup><mo>,</mo></mrow></math>]]></maths>它是指从被控过程的第i个输入到第j个输出的传递函数,g<sub>0ij</sub>(s)是其稳定正则的有理传递函数部分,θ<sub>ij</sub>是其相应的过程传输时滞,i,j=1,…p,进行分解,分解形式如下:G(s)=G<sub>D</sub>(s)G<sub>O</sub>(s),其中,<maths num="002"><![CDATA[ <math><mrow><msub><mi>G</mi><mi>D</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mi>diag</mi><mo>{</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mi>s</mi></mrow></msup><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&theta;</mi><mi>p</mi></msub><mi>s</mi></mrow></msup><mo>}</mo><mo>,</mo><msub><mi>G</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>G</mi><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>G</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中θ<sub>i</sub>取为G<sup>-1</sup>(s)的第i列中最大预估值;2)设计用于保证控制器能实现的两个调节因子:<maths num="003"><![CDATA[ <math><mrow><msub><mi>G</mi><mi>N</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mi>diag</mi><mo>{</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>r</mi><mi>z</mi></msub></munderover><msup><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><mi>s</mi><mo>+</mo><msub><mi>z</mi><mi>j</mi></msub></mrow><mrow><mi>s</mi><mo>+</mo><msubsup><mi>z</mi><mi>j</mi><mo>*</mo></msubsup></mrow></mfrac><mo>)</mo></mrow><msub><mi>k</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub></msup><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>r</mi><mi>z</mi></msub></munderover><msup><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><mi>s</mi><mo>+</mo><msub><mi>z</mi><mi>j</mi></msub></mrow><mrow><mi>s</mi><mo>+</mo><msubsup><mi>z</mi><mi>j</mi><mo>*</mo></msubsup></mrow></mfrac><mo>)</mo></mrow><msub><mi>k</mi><mi>pj</mi></msub></msup><mo>}</mo></mrow></math>]]></maths><maths num="004"><![CDATA[ <math><mrow><mi>J</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mi>diag</mi><mo>{</mo><mfrac><mn>1</mn><msup><mrow><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mi>s</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><msub><mi>p</mi><mn>1</mn></msub></msup></mfrac><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mfrac><mn>1</mn><msup><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>n</mi></msub><mi>s</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>p</mi><mi>n</mi></msub></msup></mfrac><mo>}</mo></mrow></math>]]></maths>其中z<sub>j</sub>是G<sup>-1</sup>(s)第i列中不稳定的极点,k<sub>ij</sub>是z<sub>j</sub>的最大个数,λ<sub>i</sub>为控制器可调参数,p<sub>i</sub>取G(s)中第i列元素的最大相对阶次;3)设计解耦的控制器C为如下形式:<maths num="005"><![CDATA[ <math><mrow><mi>C</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>G</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>G</mi><mi>D</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>G</mi><mi>N</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>J</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>I</mi><mo>-</mo><msub><mi>G</mi><mi>D</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>G</mi><mi>N</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>J</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>在上述过程中,如果对象传递函数矩阵分解后G<sub>O</sub>(s)中仍然含有时滞项,将导致设计的控制器是无穷高阶,这时使用控制器降阶技术对该控制器进行降阶。
地址 200240上海市闵行区东川路800号