发明名称 MIMO非最小相位CSTR反应器的非线性控制器设计方法
摘要 本发明公开了一种MIMO非最小相位CSTR反应器的非线性控制器设计方法,该控制器设计方法在保证系统外部动态满足性能要求的同时也能够保证系统内部零动态的稳定。首先根据实际的生产过程对CSTR进行系统建模,得到该控制系统的非线性数学模型;其次利用状态反馈线性化理论将得到的MIMO非线性系统进行精确反馈线性化,得到其线性化标准型;将该MIMO线性化标准型拆解为若干线性子系统以及一个SISO线性化标准型子系统;最后基于极点配置和李雅普诺夫稳定性理论,给出了具有非最小相位特性SISO线性化标准型子系统的一种非线性控制器设计方法,联合线性子系统部分的控制器,得到MIMO非最小相位CSTR最终形式的控制器。
申请公布号 CN104035336B 申请公布日期 2017.01.11
申请号 CN201410239871.9 申请日期 2014.05.30
申请人 河海大学 发明人 王万成;金晓孝;张杰;施文
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 一种MIMO非最小相位CSTR反应器的非线性控制器设计方法,其特征在于,包括以下步骤:(1)根据实际的化工生产过程对CSTR进行状态空间建模,得到系统的非线性数学模型;(2)将系统的非线性数学模型进行状态反馈线性化,得到一般形式的MIMO线性化标准型;(3)将MIMO线性化标准型进行拆分,将其分解为若干线性子系统部分和一个SISO线性化标准型子系统部分,对各子系统单独设计控制器;(4)为了确保所设计的控制器能有效控制具有非最小相位特性的SISO线性化标准型子系统,给出了该控制器的一般形式,控制器中的待定参数通过极点配置和李雅普诺夫稳定性理论来确定;所述步骤(1)中,实际的化工生产过程所选用的是利用CSTR生产环戊烯这一化工过程,建立其状态空间模型如下式所示:状态方程:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mfrac><mrow><msub><mi>dC</mi><mi>A</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>1</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>A</mi></msub><mo>-</mo><msub><mi>k</mi><mn>3</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msubsup><mi>C</mi><mi>A</mi><mn>2</mn></msubsup><mo>+</mo><mo>(</mo><msub><mi>C</mi><mrow><mi>A</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>C</mi><mi>A</mi></msub><mo>)</mo><mfrac><mi>F</mi><mi>V</mi></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>dC</mi><mi>B</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><msub><mi>k</mi><mn>1</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>A</mi></msub><mo>-</mo><msub><mi>k</mi><mn>2</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>B</mi></msub><mo>-</mo><msub><mi>C</mi><mi>B</mi></msub><mfrac><mi>F</mi><mi>V</mi></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mi>d</mi><mi>T</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mo>(</mo><mo>-</mo><msub><mi>&Delta;H</mi><mn>1</mn></msub><mo>)</mo><msub><mi>k</mi><mn>1</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>A</mi></msub><mo>+</mo><mo>(</mo><mo>-</mo><msub><mi>&Delta;H</mi><mn>2</mn></msub><mo>)</mo><msub><mi>k</mi><mn>2</mn></msub><mo>(</mo><mi>T</mi><mo>)</mo><msub><mi>C</mi><mi>B</mi></msub></mrow><mrow><msub><mi>&rho;C</mi><mi>&rho;</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>(</mo><mo>-</mo><msub><mi>&Delta;H</mi><mn>3</mn></msub><mo>)</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mi>A</mi></msub><msubsup><mi>C</mi><mi>A</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>Q</mi><mi>H</mi></msub></mrow><mrow><msub><mi>&rho;C</mi><mi>&rho;</mi></msub></mrow></mfrac><mo>+</mo><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>-</mo><mi>T</mi><mo>)</mo><mfrac><mi>F</mi><mi>V</mi></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001077070860000011.GIF" wi="1574" he="463" /></maths>系统的被控输出为:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mi>T</mi></mtd></mtr><mtr><mtd><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><msub><mi>C</mi><mi>B</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001077070860000012.GIF" wi="190" he="150" /></maths>其中,C<sub>A0</sub>为物质A'的初始浓度,T<sub>0</sub>为反应初始温度,T是反应温度,C<sub>A</sub>、C<sub>B</sub>分别是物质A'和B'的浓度,F是物质A'的输入流量,V是CSTR的常容积;控制输入u<sub>1</sub>=F/V为稀释率,u<sub>2</sub>=Q<sub>H</sub>为单位体积增加或减少的热量;C<sub>ρ</sub>是热容量,ρ是混合物的密度;ΔH<sub>i</sub>是反应热;速率系数k<sub>i</sub>(T)由下列的Arrhenius方程给出:<img file="FDA0001077070860000021.GIF" wi="598" he="119" />i=1,2,3分别对应CSTR三个反应过程,其中,各过程的反应时间常数为k<sub>i0</sub>,反应激活能量为E<sub>i</sub>,R为气体常数;控制目标是:通过调节被控输入u<sub>1</sub>=F/V和u<sub>2</sub>=Q<sub>H</sub>来使输出y<sub>1</sub>=T和y<sub>2</sub>=C<sub>B</sub>分别稳定在设定的平衡点T<sub>S</sub>、C<sub>BS</sub>上,此时物质A'的浓度相应记为C<sub>AS</sub>,(C<sub>AS</sub>,C<sub>BS</sub>,T<sub>S</sub>)为系统达到稳态时的平衡点;所述步骤(4)中,一般形式的控制器设为v=‑Gx+v<sub>NL</sub>其中,x为SISO线性化标准型子系统的状态向量,G为行增益向量,v<sub>NL</sub>是为使系统内部零动态稳定而引入的非线性补偿项。
地址 211100 江苏省南京市江宁区佛城西路8号