发明名称 一种基于LMIs的汽包水位控制系统设计方法
摘要 一种基于LMIs的汽包水位控制系统设计方法,包括以下步骤:a.将给水流量和蒸汽流量作为汽包水位的影响因素,确定汽包水位平衡函数;b.将汽包水位平衡的传递函数转化成状态空间表达式的形式;c.设置汽包水位设定点,引入一个新的状态向量,由线性矩阵不等式求得矩阵M和P;d.利用矩阵M和P构造用于控制给水流量的状态反馈控制器。本发明针对锅炉汽包水位的动态特性,提出了一种基于LMIs的汽包水位控制系统的设计方法,该方法不仅可以保证系统的稳定性,而且能够确保控制动作满足饱和约束条件,从而提高了汽包水位的控制精度,保证了锅炉系统的安全运行。
申请公布号 CN106125781A 申请公布日期 2016.11.16
申请号 CN201610742294.4 申请日期 2016.08.28
申请人 华北电力大学(保定) 发明人 魏乐;秦腾腾;侯君虹
分类号 G05D9/12(2006.01)I 主分类号 G05D9/12(2006.01)I
代理机构 石家庄冀科专利商标事务所有限公司 13108 代理人 李羡民;高锡明
主权项 一种基于LMIs的汽包水位控制系统设计方法,其特征是,所述方法包括以下步骤:a.将给水流量和蒸汽流量作为汽包水位的影响因素,确定汽包水位平衡函数:<maths num="0001"><math><![CDATA[<mrow><mi>h</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>g</mi><mn>1</mn></msub><mi>s</mi></mfrac><mo>&lsqb;</mo><msub><mi>q</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>q</mi><mi>D</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mfrac><msub><mi>g</mi><mn>2</mn></msub><mrow><msub><mi>&tau;</mi><mn>0</mn></msub><mi>s</mi><mo>+</mo><mn>1</mn></mrow></mfrac><msub><mi>q</mi><mi>D</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>g</mi><mn>3</mn></msub><mfrac><mrow><mi>&beta;</mi><mi>&tau;</mi><mi>s</mi><mo>+</mo><mn>1</mn></mrow><mrow><msup><mi>&tau;</mi><mn>2</mn></msup><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><mi>&zeta;</mi><mi>&tau;</mi><mi>s</mi><mo>+</mo><mn>1</mn></mrow></mfrac><msub><mi>q</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001095626330000011.GIF" wi="1286" he="127" /></maths>其中:h(s)代表锅炉汽包水位(%);q<sub>W</sub>(s)是给水流量的变化(%);q<sub>D</sub>(s)是蒸汽流量的变化(%);g<sub>1</sub>、g<sub>2</sub>、g<sub>3</sub>、β、τ和τ<sub>0</sub>是在一定负载下的常数;s是传递函数中的复变量;b.将汽包水位平衡的传递函数转化成状态空间表达式的形式:<maths num="0002"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>q</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>q</mi><mi>D</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001095626330000017.GIF" wi="798" he="63" /></maths>h(t)=Cx(t)X(t)=[X<sub>1</sub>(t),X<sub>2</sub>(t)]<sup>T</sup>代表了系统的状态变量,A(t)、B<sub>1</sub>(t)、B<sub>2</sub>(t)和C为参数矩阵:<img file="FDA0001095626330000012.GIF" wi="1044" he="215" />C=[1 ‑1];c.设置h<sub>r</sub>(t)为汽包水位设定点,引入一个新的状态向量<img file="FDA0001095626330000013.GIF" wi="312" he="143" />其中,<img file="FDA0001095626330000014.GIF" wi="532" he="94" />令<img file="FDA0001095626330000015.GIF" wi="1271" he="143" /><img file="FDA0001095626330000016.GIF" wi="406" he="79" />由线性矩阵不等式:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mrow><mo>&lsqb;</mo><mover><mi>A</mi><mo>&OverBar;</mo></mover><mi>P</mi><mo>&rsqb;</mo></mrow><mi>s</mi></msub><mo>+</mo><msub><mrow><mo>&lsqb;</mo><msub><mover><mi>B</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mi>M</mi><mo>&rsqb;</mo></mrow><mi>s</mi></msub></mrow></mtd><mtd><msub><mover><mi>B</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub></mtd><mtd><mrow><mi>P</mi><msup><mover><mi>C</mi><mo>&OverBar;</mo></mover><mi>T</mi></msup></mrow></mtd><mtd><msup><mi>M</mi><mi>T</mi></msup></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mrow><mo>-</mo><msup><mi>&gamma;</mi><mn>2</mn></msup><mi>I</mi></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mrow><mo>-</mo><mi>I</mi></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mrow><mo>-</mo><mi>R</mi></mrow></mtd></mtr></mtable></mfenced><mo>&lt;</mo><mn>0</mn></mrow>]]></math><img file="FDA0001095626330000021.GIF" wi="876" he="358" /></maths>求得矩阵M和P,其中,γ是正标量,Q和R是x<sub>e</sub>(t)和q<sub>W</sub>(t)的加权因子;d.利用矩阵M和P构造用于控制给水流量q<sub>W</sub>(s)的状态反馈控制器:q<sub>W</sub>(t)=KX(t),其中矩阵K=MP<sup>‑1</sup>。
地址 071003 河北省保定市北市区永华北大街619号