现代控制理论大作业

1、现代控制理论(主汽温对象模型)班级:学号:姓名:目录一. 背景及模型建立1. 火电厂主汽温研究背景及意义2. 主汽温对象的特性3. 主汽温对象的数学模型二.分析1. 状态空间表达2. 化为约当标准型状态空间表达式并进行分析3. 系统状态空间表达式的求解4. 系统的能控性和能观性5. 系统的输入输出传递函数6. 分析系统的开环稳定性7. 闭环系统的极点配置8. 全维状态观测器的设计9. 带状态观测器的状态反馈控制系统的状态变量图10. 带状态观测器的闭环状态反馈控制系统的分析三.结束语1. 主要内容2. 问题及分析3. 评价一.背景及模型建立1.火电厂主汽温研究背景及意义火电厂锅炉主汽温控制决定

2、着机组生产的经济性和安全性。由于锅炉的蒸汽容量非常大、过热汽管道很长，主汽温调节对象往往具有大惯性和大延迟，导致锅炉主汽温控制存在很多方面的问题，影响机组的整个工作效率。主汽温系统是表征锅炉特性的重要指标之一，主汽温的稳定对于机组的安全运行至关重要。其重要性主要表现在以下几个方面：(1) 汽温过高会加速锅炉受热面以及蒸汽管道金属的蠕变，缩短其使用寿命。例如，12CrMoV 钢在585环境下可保证其应用强度的时间约为10万小时，而在 595时，其保证应用强度的时间可能仅仅是 3 万小时。而且一旦受热面严重超温，管道材料的强度将会急剧下降，最终可能会导致爆管。再者，汽温过高也会严重影响汽轮机的汽缸

3、、汽门、前几级喷嘴和叶片、高压缸前轴承等部件的机械强度，从而导致设备损坏或者使用年限缩短。(2) 汽温过低，会使得机组循环热效率降低，增大煤耗。根据理论估计可知：过热汽温每降低10，会使得煤耗平均增加0.2%。同时，汽温降低还会造成汽轮机尾部的蒸汽湿度增大，其后果是，不仅汽轮机内部热效率降低，而且会加速汽轮机末几级叶片的侵蚀。此外，汽温过低会增大汽轮机所受的轴向推力，不利于汽轮机的安全运行。(3) 汽温变化过大会使得管材及有关部件产生疲劳，此外还将引起汽轮机汽缸的转子与汽缸的胀差变化，甚至产生剧烈振动，危及机组安全运行。据以上所述，工艺上对汽温控制系统的质量要求非常严格，一般控制误差范围在&#

4、177;5。主汽温太高会缩短管道的使用寿命，太低又会降低机组效率。所以必须实现汽温系统的良好控制。而汽温被控对象往往具有大惯性、大延时、非线性，时变一系列的特性，造成对象的复杂性，增加了控制的难度。现代控制系统中有很多关于主汽温的控制方案，本文我们着重研究带状态观测器的状态反馈控制对主汽温的控制2.主汽温对象的特性2.1主汽温对象的静态特性主汽温被控对象的静态特性是指汽温随锅炉负荷变化的静态关系。过热器的传热形式、结构和布置将直接影响过热器的静态特性。现代大容量锅炉多采用对流过热器、辐射过热器和屏式过热器。对流过热器布置在4501000烟气温度的烟道中，受烟气的横向和纵向冲刷，烟气以对流方式将

5、热量传给管道。而辐射过热器则是直接吸收火焰和高温烟气的辐射能。屏式过热器布置在炉膛内上部或出口处，属于辐射或半辐射过热器。对于辐射过热器，当锅炉负荷增加时，必须增加燃料量和风量。随着炉膛温度的提高，辐射传热量也将增加。但炉膛温度提高得并不多，而蒸汽量与燃料量是成正比的，这样就造成辐射传热量的增加赶不上蒸发量的增加。再者来说，当负荷增加时，强化燃烧后炉膛出口烟气温度将升高，这表明每公斤燃料产生的烟气带出炉膛的热量增多，这也说明了炉膛辐射吸热量的相对减少。所以，辐射受热面的汽温是随着锅炉负荷的增加而降低的。对于对流式受热面，当锅炉负荷增加和燃料消耗量增大时，流过对流受热面的烟气流速也增加，从而使对

6、流的放热系数增大。另外，因为炉膛出口烟温升高，即进入对流受热面的吸热量增加值，超过了流过对流受热面的蒸汽流量的增加值。所以，对流受热面的汽温是随锅炉负荷的增加而升高的。2.2 主汽温对象的动态特性2.2.1蒸汽流量扰动下主汽温的动态特性在蒸汽流量D产生阶跃扰动时，主汽温变化的响应曲线如图2-2-1所示：图2-2-1 蒸汽流量扰动下主汽温响应曲线蒸汽流量D扰动下，主汽温调节对象动态特性的特点是: 有延迟、有惯性、有自衡能力。当锅炉负荷增加时，通过对流式过热器的烟气温度和流速都增加，因此对流式过热器出口汽温升高；但对于辐射式过热器，炉膛内烟温升高增加的辐射传热量小于蒸汽流量增加所需的吸热量，因此辐

7、射式过热器出口汽温下降。2.2.2烟气热量扰动下汽温的动态特性在烟气热量产生阶跃扰动下，汽温变化的响应曲线如图2-2-2所示：图2-2-2 烟气热量扰动下汽温响应曲线烟气热量扰动下，汽温调节对象动态特性的特点是:有迟延、有惯性、有自平能力。由于烟气热量变化时，沿过热器长度使烟气和过热蒸汽之间的传热量同时变化，因此，汽温反应较快，其时间常数Tc和迟延时间均较小。2.2.3减温水量扰动下汽温的动态特性减温水量扰动下，汽温调节对象动态特性的特点是: 有迟延、有惯性、有自衡能力。由于现代大型锅炉过热器管路很长，因此减温水量变化时汽温反应较慢，其时间常数Tc和迟延时间均较大。在减温水量产生阶跃扰动下，汽

8、温变化的响应曲线如图2-2-3所示：图2-2-3 减温水量扰动下汽温响应曲线3.主汽温对象的数学模型3.1主汽温对象模型的建立分析和设计自动控制系统的一个首要任务是建立系统的数学模型，因为不论要了解的是简单系统的特性还是复杂系统的特性，都必须掌握系统中各变量之间的相互动态关系。尽管，电厂汽温控制系统的动态特性复杂，具有多变量、非线性和分布复杂的特点，难以建立精确的数学模型，但是运用多种知识，建立能相对反映汽温控制系统动态性能的数学模型，还是对研究相应的自动控制方法大有益处。锅炉过热器是由辐射过热器、对流过热器和减温器等组成，其任务是将从汽包出来的饱和蒸汽加热到一定的数值，然后送往汽轮机去做功。

9、通常称减温器前的过热器为前级过热器，减温器后的过热器为后级过热器。图2-3-1 过热蒸汽喷水减温系统示意图过热器布置在高温烟道中，大型锅炉的过热器往往分为若干段，在各段之间设置喷水减温器，即采用汽温的分段控制，温度调节用减温水由锅炉的给水系统提供。其示意图如图2-3-1。图中，s为过热器出口蒸汽温度，它是控制系统的被调量；Ws1、Ws2是减温器的喷水量，它是控制系统的调节量。3.2主汽温对象的模型目前电厂中广泛使用的汽温调节系统主要是基于常规控制器构成的串级控制系统。因为常规PID控制器的参数是根据被控对象的数学模型来整定的。而汽温调节对象的动态特性随机组负荷的变化而有较大变化，因此，当机组负

10、荷变化时，该控制系统的工作点建立多固定模型，设计相应的控制器，在运行过程中根据系统工况变化选择相应的控制器输出以获得满意的控制效果，通过控制器在不同工况下的切换使系统达到全局稳定性的要求。某机组导前区模型：，惰性区模型：，总体模型：。对应的一种状态空间模型为二.分析1.状态空间表达式2.化为约当标准型状态空间表达式并进行分析程序:>>clear;>>A=-1/30 1/30 0;0 -1/30 1/30;0 0 -1/12;>>B=0 0 1.1*2.8/12'>>C=1 0 0;>>D=0;>>P,J=jorda

11、n(A);>>A1=inv(P)*A*P>>B1=inv(P)*B>>C1=C*P结果:变换矩阵特征值 因此约当标准型状态空间表达式为 系统矩阵A的特征根有重根,因此,经线性变换,将A化为约当标准型.3.系统状态空间表达式的求解线性定常非齐次状态方程的解可表示为矩阵A有重特征根时,当输入脉冲函数时,则系统状态方程的解 因此4.系统的能控性和能观性(1) 能控性代数判据:线性定常系统状态完全能控的充要条件为:矩阵是满秩的,或表示成 模态判据1:若系统矩阵A为对角型,则系统能控的充要条件是输入矩阵B没有任何一行的元素全部为零 模态判据2:若系统矩阵为约当型,则系

12、统能控的充要条件是:输入矩阵B中对应于互异的特征值的各行,没有一行的元素全为零;输入矩阵B中与每个约当块最后一行相对应的各行,没有一行的元素全为零.>>Qc=ctrb(A,B)>>rank(Qc) 因此系统能控(2) 能观性代数判据:系统定常系统状态能观测的充要条件为:的秩为n 模态判据1:若系统A为对角型,则系统能观测的充要条件是输出矩阵C没有任何一列的元素全部为零 模态判据2:若系统A为约当型,则系统能观测的充要条件是:输出矩阵C中对应于互异的特征值的各列,没有一列的元素全为零;输出矩阵C中与每个约当块的第一列相对应的各列,没有一列的元素全为零>>Qo=

13、obsv(A,C)>>rank(Qo)因此系统能观5.系统的输入输出传递函数>>G=ss(A,B,C,D);>>Gc=tf(G)Transfer function: 0.0002852-s3 + 0.15 s2 + 0.006667 s + 9.259e-0056.分别用特征根方法和Lyapunov第二法分析系统的开环稳定性(1) 特征根方法:特征根分别为-0.0333,-0.0333,-0.0833,均在左半平面,因此系统开环稳定(2) Lyapunov第二法:1) 由平衡点方程 解出唯一平衡点为坐标原点2) 设李氏函数为 取,则矩阵由下式确定: 则 解出

14、 用塞尔维斯判据检验P的正定性: 可知,正定,系统在原点处的平衡状态是渐近稳定的.7.采用状态反馈将系统的闭环极点配置到合适的值,使得闭环系统阶跃响应的上升时间比开环系统阶跃响应的上升时间缩短3倍左右>>step(G)开环阶跃响应如图一,上升时间约为200s图一:>>M=-0.15 -0.15 -0.4;>>K=acker(A,B,M)>>eig(A-B*K)>>G1=ss(A-B*K,B,C,D);>>G2=tf(G1)>>step(G2)期望反馈闭环极点为-0.15,-0.15,-0.4状态反馈矩阵为验证配

15、置结果ans=-0.1500,-0.1500,-0.4000a = x1 x2 x3 x1 -0.03333 0.03333 0 x2 0 -0.03333 0.03333 x3 -4.492 -2.975 -0.6333b = u1 x1 0 x2 0 x3 0.2567c = x1 x2 x3 y1 1 0 0d = u1 y1 0Continuous-time model.经过状态反馈后,系统的传递函数阵为Transfer function: 0.0002852-s3 + 0.7 s2 + 0.1425 s + 0.009并且经过状态反馈后,系统的阶跃响应为图二,上升时间约为50s图二:

16、8.为系统设计全维状态观测器，将观测器的极点配置到合适的值(1)配置状态观测器的极点>>N=-3 -3 -8;>>L=acker(A',C',N)'>>eig(A-L*C)配置观测器期望极点为-3,-3,-8全维状态观测器的反馈阵为验证配置结果,ans=-8.0000,-3.0000,-3.0000(2)生成系统状态估计器>>est=estim(G,L)a = x1 x2 x3 x1 -13.88 0.03333 0 x2 -1661 -0.03333 0.03333 x3 -6.061e+004 0 -0.08333b

17、= u1 x1 13.85 x2 1661 x3 6.061e+004c = x1 x2 x3 y1 1 0 0 y2 1 0 0 y3 0 1 0 y4 0 0 1d = u1 y1 0 y2 0 y3 0 y4 0Input groups: Name Channels Measurement 1 Output groups: Name Channels OutputEstimate 1 StateEstimate 2,3,4 Continuous-time model.(3)状态观测器的阶跃响应>>step(4)输出和观测状态的传递函数>>tf(est)Transf

18、er function from input to output. 13.85 s2 + 56.99 s + 72 #1: - s3 + 14 s2 + 57 s + 72 13.85 s2 + 56.99 s + 72 #2: - s3 + 14 s2 + 57 s + 72 1661 s2 + 2214 s + 71.96 #3: - s3 + 14 s2 + 57 s + 72 6.061e004 s2 + 4041 s + 67.35 #4: - s3 + 14 s2 + 57 s + 72Input groups: Name Channels Measurement 1 Output

19、 groups: Name Channels OutputEstimate 1 StateEstimate 2,3,4 9.画出带状态观测器的状态反馈控制系统的状态变量图>>T=1:300;plot(T/50,y(1:300,:);ylabel('y');legend('系统','观测');>>T=1:300;plot(T/50,x1(1:300,:);ylabel('x1');legend('系统','观测');>>T=1:300;plot(T/50,x2(1:

20、300,:);ylabel('x2');legend('系统','观测');>>T=1:300;plot(T/50,x3(1:300,:);ylabel('x3');legend('系统','观测');(1) 跟踪正弦波正弦x1正弦x2正弦x3正弦y(2) 跟踪阶跃信号阶跃x1阶跃x2阶跃x3阶跃y10.求带状态观测器的闭环状态反馈控制系统的状态空间表达式，并求闭环传递函数、分析闭环系统的稳定性(1)带状态观测器的闭环状态反馈控制系统的状态空间表达式设能控能观测的受控系统为 状态反馈控制

21、率为 状态观测器方程为 由以上三式可得整个闭环系统的状态空间表达式为 将它写成分块矩阵的形状 即(2)系统的闭环传递函数带观测器状态反馈闭环系统的传递函数阵等于直接状态反馈闭环系统的传递函数阵,它玉是否采用观测器反馈无关,即 (3)闭环系统的稳定性由观测器构成状态反馈的闭环系统,其特征多项式等于状态反馈部分的特征多项式和观测器部分的特征多项式的乘积,而且相互独立.这个闭环系统的特征值为-0.15,-0.15,-0.4,-3,-3,-8,因此闭环系统稳定.三.结束语1.主要内容本次大作业是以火电厂主汽温为研究背景,通过分析主汽温的静态以及动态特性(蒸汽流量扰动下主汽温的动态特性, 烟气热量扰动下

22、汽温的动态特性,减温水量扰动下汽温的动态特性)进行主汽温对象模型的建立,最后得到状态空间模型为:.以上的模型是老师建立好的,但是通过阅读这些资料并且对比这些模型的建立过程,我了解了建模的基本步骤,以及如何对模型进行简化.我做的工作主要包括将系统状态空间表达式化为标准型(约当标准型),线性定常齐次状态方程的求解,控制系统的能控性及能观性的判别,控制系统的稳定性的判断,闭环系统的极点配置,状态观测器的设计.2.问题及分析在熟读课本的基础上,借助MATLAB进行辅助计算及对部分环节进行仿真,对控制系统有较为全面的分析.同时也遇到了很多问题.开始的时候由于盲目的追求进度,忽略了这个系统矩阵的特征根有重根,就想把它化为对角线标准型.后来通过查阅资料用语句求出了矩阵特征根及特征向量.在每一章的学习中,针对同一类问题都有若干种不同的解法,比如关于线性定常非齐次状态方程转移矩阵的计算就有直接法,拉普拉斯变换法,化矩