CN107065538A - 汽包锅炉‑汽轮机组的模糊跟踪控制方法 - Google Patents
汽包锅炉‑汽轮机组的模糊跟踪控制方法 Download PDFInfo
- Publication number
- CN107065538A CN107065538A CN201710116314.1A CN201710116314A CN107065538A CN 107065538 A CN107065538 A CN 107065538A CN 201710116314 A CN201710116314 A CN 201710116314A CN 107065538 A CN107065538 A CN 107065538A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- mtd
- mtr
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Abstract
本发明公开了一种汽包锅炉‑汽轮机组的模糊跟踪控制方法。针对于汽包锅炉‑汽轮机系统的强耦合、非线性、严格的输入限制等问题,为了使整个系统能够快速跟踪设定值,提高控制系统的性能,并减少燃料的损耗,本发明基于前馈加反馈的控制思想,设计了一种模糊跟踪控制方法。前馈控制输入信号仅由跟踪目标决定,形式简单,易于实现,反馈控制输入环节加入一个处理器,得到了带时滞的模糊状态反馈控制器,提高了控制系统的实时性。该控制方法能够有效地提高机组的动态响应品质,快速地跟踪设定值,跟踪误差小,控制精度高,提高了经济效益。
Description
技术领域
本发明属于热能动力工程与自动控制技术领域,具体涉及一种汽包锅炉-汽轮机组的新型模糊跟踪控制方法。
背景技术
锅炉-汽轮机系统是现代火力发电厂的重要组成部分,它的性能好坏直接影响了整个电力系统的性能。随着国民经济的迅猛发展,高参数、大容量的单元机组在电网中占有主导地位。大机组要根据电网的频率偏差以及中央调度来参与电网的调峰和调频,这就要求单元机组有良好的负荷适应能力,也就是它的主要运行参数在负荷变化过程中相对稳定。此外,由于国家‘节能减排’政策的实施,迫使电厂将降低成本、减少污染排放等因素考虑进来,基于上述考虑,研究设计快速负荷响应的可靠性强的控制方法具有非常重要的意义。
汽包锅炉-汽轮机系统是一个复杂的控制系统,具有多变量耦合、非线性、严格的输入限制、需调节参数多,计算量大等特点。目前,对该系统的研究受到了广泛的关注,出现了一系列的控制方法。传统的PID控制方法是建立在局部工作点线性化的模型基础上,仅对系统运行在平衡点附近的小范围有效,在负荷变化范围较大的情况下,有导致工作点偏离或者稳定性变差的问题。基于非线性理论的控制方法主要有反步控制方法、反馈线性化方法、自适应控制方法等,能从根本上解决系统大范围负荷变化的问题,但算法理论复杂,方案不易实现。滑模变结构控制给汽包锅炉-汽轮机系统提供了另一种有效的控制策略,它的优点是响应迅速,易于实现,缺点是控制输入信号变化大,在控制过程中会发生抖振现象。
Takagi-Sugeno(T-S)模糊控制,作为一种智能控制方法,为解决非线性系统的控制问题提供了新的途径。T-S模糊模型已被证明能够以任意精度逼近定义在紧集上的非线性模型。该模型具有结构简单、逼近能力强等特点。模糊控制技术具有控制器设计简便,适应于许多非线性系统、鲁棒性强等特点。将模糊控制技术应用在汽包锅炉-汽轮机组的跟踪控制上,在研究中已经验证该方法用于状态反馈控制中可以使跟踪误差渐进趋于零,可以很好的解决机组模型的多变量耦合、非线性等问题。因此,该方法很适合用于汽包锅炉-汽轮机系统的跟踪控制中。
发明内容
针对现有技术中汽包锅炉-汽轮机控制系统中存在非线性等的不足,本发明提出了一种汽包锅炉-汽轮机系统的前馈加带有时滞的模糊状态反馈跟踪控制方法,该控制方法提出的前馈控制仅由跟踪目标决定,使用模糊逻辑系统来逼近系统中的非线性项,在输出端加入一个处理器来实时地监控系统的各种状态,并计算出带有时滞的模糊反馈控制器。
本发明为实现上述目的采用的技术方案是:一种汽包锅炉-汽轮机组的模糊跟踪控制方法,包括以下步骤:
步骤1建立汽包锅炉-汽轮机160MW燃油机组的数学模型:
式中,P,E,ρf为系统状态变量,分别表示蒸汽压力(kg/cm2),输出功率(MW)和锅炉流体密度(kg/cm3);三个控制输入u1(t)、u2(t)、u3(t)分别代表燃料阀执行开度、节流阀执行开度和给水阀执行开度,单位为pu;P,E和L代表系统的输出,L表示锅炉水位;
αcs(t)和qe(t)分别为质量系数和蒸发率,且满足
控制输入为阀门,其上下限和变化率须满足如下条件:
锅炉汽轮机的动态模型转化为:
其中,
步骤2设计前馈控制输入信号:
定义跟踪误差变量:
其中,y1ref(t),y2ref(t),y3ref(t)为给定的跟踪信号,且y3ref(t)≡0;则x1ref(t)=y1ref(t),x2ref(t)=y2ref(t),x3ref(t)为待求量;
设:
将(6)代入(4)可得:
其中,
为反馈控制输入信号;
为前馈控制输入信号,且满足比较(7)和(4),求出前馈控制输入信号为:
步骤3设计反馈控制输入信号:
1)假设在锅炉流体密度达到稳态值时,系统输入信号u(t),输出信号y(t)都达到稳态值,也就是
将(8),(9)代入(7)的最后一个表达式,并整理,可得关于的一个一元二次方程:
λ1x3ref(t)2+λ2x3ref(t)+λ3=0, (10)
其中,
λ1=0.1307;
求解方程(10),并舍弃较小解,可求得
同样的,在处理器中采样y1(t),y2(t),y3(t)我们可求出x3(t)
其中,
通过处理器计算,可输出任意时刻的状态
2)设计带有时滞的模糊状态反馈控制器,
控制器设为:
使得
其中,0≤τ(t)≤τ;
模型(14)是一个非线性模型,将其模糊化,可得含有两个线性子系统的T-S模糊模型:
IFTHEN
IFTHEN
其中,
隶属度函数为:
模糊状态反馈控制器可表示为:
其中,Ki,i=1,2,为待求量,综上模糊闭环系统可表示为:
基于上述模型,定义如下Lyapunov-Krasovskii泛函:
其中P,Q,R为正定矩阵,对(18)求导,可得:
根据Wirtinger-based积分不等式和倒立凸不等式,将上述公式中的积分项扩大为
其中
定义:
并假设ei=[0n×(i-1)n In×n 0n×(i+1)n],i=1,2,...,6;
由于:
根据(19),(20),(21),得
其中,
Ξ<0可满足此时闭环系统是稳定的;
令G1=εG2,X=G2 -1,(23)式左右两端同乘以diag{X,X,X,X,X,X}T及其转置,得:
其中:
若要满足条件(24),只需求解如下条件:
其中,
若条件(25)成立,则状态反馈控制增益矩阵可由Ki=YiX-1求出;
3)根据不同的给定的采样周期T,得到不同的时滞上限τ,求出控制器增益矩阵Ki,i=1,2,从而得到状态反馈控制信号:
步骤4求出最终的控制输入信号u(t):
根据步骤2和步骤3所述的前馈控制信号和状态反馈控制信号根据公式最终求出的控制信号。
进一步,所述的方法还包括基于李雅普诺夫稳定理论对控制器设计的数值仿真验证;
应用数值仿真来验证由于采样周期的不同导致的不同的时滞输入的控制方法的有效性,通过数值分析可见,采样周期越小,则跟踪误差越小,控制效果越好。
进一步,所述的仿真验证,三个状态蒸汽压力、输出功率和锅炉流体密度的初始值选择为[P L ρf]=[75.6 15.27 299.6]。
本发明的优点是:
针对汽包锅炉-汽轮机系统的非线性、变量的强耦合性、输入信号的强制限制等问题,为了使该系统能够快速响应负荷变化,达到新的平衡点,设计了一种汽包锅炉-汽轮机系统的前馈加反馈控制方法。本发明提供的前馈输入信号仅由跟踪目标决定,形式简单,易于实现。考虑到系统的部分状态,例如锅炉流体密度,不能直接测量的特点,设计出基于处理器的时滞模糊反馈控制器。两类控制输入信号结合在一起,有效地解决了汽包锅炉-汽轮机系统的非线性控制问题,控制精度高,响应速度快,输入信号变化平缓,有效的节约了能源。仿真结果表明这种控制方法的有效性,实现了理想的控制效果,具有实际意义。
附图说明
图1本发明涉及的汽包锅炉-汽轮机系统结构示意图;
图2本发明涉及的汽包锅炉-汽轮机系统的控制策略;
图3a、图3b代表蒸汽压力的控制效果和跟踪误差曲线;
图3c、图3d代表输出功率的控制效果和跟踪误差曲线;
图3e代表锅炉水位的控制效果和跟踪误差曲线;
图3f、图g和图3h分别代表燃料阀执行开度、节流阀执行开度和给水阀执行开度;
图4a、图4b和图4c分别代表当τ=5时系统的蒸汽压力的跟踪误差曲线、输出功率的跟踪误差曲线和锅炉水位的跟踪误差曲线;
图5a、图5b和图5c分别代表当τ=50时系统的蒸汽压力的跟踪误差曲线、输出功率的跟踪误差曲线和锅炉水位的跟踪误差曲线。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明:
汽包锅炉-汽轮机组的先进的协调控制方法,包括如下步骤:一种汽包锅炉-汽轮机组的模糊跟踪控制方法,包括以下步骤:
步骤1建立数学模型
建立汽包锅炉-汽轮机燃油机组的数学模型:
式中,P,E,ρf为系统状态变量,分别表示蒸汽压力(Kg/cm2),输出功率(MW)和锅炉流体密度(Kg/cm3);三个控制输入,u1(t),u2(t),u3(t)分别代表燃料阀执行开度、节流阀执行开度和给水阀执行开度,单位为pu;P,E和L代表系统的输出,L表示锅炉水位;
动态模型是基于160MW燃油机组建立的,是具有三输入、三输出的三阶非线性模型。
αcs(t)和qe(t)分别为质量系数和蒸发率,且满足
控制输入为阀门,其上下限和变化率须满足如下条件:
为了简化分析,设令:
x(t)=(x1(t) x2(t) x3(t))T=(P E ρf)T,
y(t)=(y1(t) y2(t) y3(t))T=(P E L)T,
锅炉汽轮机的动态模型转化为:
其中,
步骤2设计前馈控制输入信号:
本发明要设计一种前馈加反馈的控制方法,以实现三个输出变量的跟踪控制。
定义跟踪误差变量:
其中y1ref(t),y2ref(t),y3ref(t)为给定的跟踪信号,且y3ref(t)≡0;则x1ref(t)=y1ref(t),x2ref(t)=y2ref(t),x3ref(t)为待求量;
设
将(6)代入(4)可得:
其中,
为反馈控制输入信号;
为前馈控制输入信号,且满比较(7)和(4),求出前馈控制输入信号为:
步骤3设计反馈控制输入信号:
由于公式(7)中有一个未知量x3ref(t),要求解状态反馈控制器,我们需要知道三个状态变量。和可通过输出测得。但是未知的。因为= x3(t)-x3ref(t),先求x3ref(t),也就是系统达到稳态时,锅炉流体密度的跟踪值。
1)假设在锅炉流体密度达到稳态值时,系统输入信号u(t),输出信号y(t)都达到稳态值,也就是:
将(8),(9)代入(7)的最后一个表达式,并整理,可得关于x3ref(t)的一个一元二次方程:
λ1x3ref(t)2+λ2x3ref(t)+λ3=0 (10)
其中,
λ1=0.1307
求解方程(10),并舍弃较小解,可求得
同样的,在处理器中采样y1(t),y2(t),y3(t),我们可求出x3(t)
其中,
通过处理器计算,可输出任意时刻的状态
2)设计带有时滞的模糊状态反馈控制器
在系统输出端加入处理器,以获得不同时刻的系统的状态。
由于需计算求出,在实际应用中,在系统输出端加入一个处理器,用来实时地计算由于我们预设的处理器采集输出信号的时间间隔T,以及处理器处理数据需要耗费一定的时间,反馈输入信号必然会有一定的时滞。其大小主要是由采集信号时间间隔T所决定。因此,对系统(7)的状态反馈控制问题变成了基于输入时滞的状态反馈控制问题。控制器设为:
使得
其中,0≤τ(t)≤τ;
模型(14)是一个非线性模型,将其模糊化,可得含有两个线性子系统的T-S 模糊模型:
IFTHEN
IFTHEN
其中,
隶属度函数为:
模糊状态反馈控制器可表示为:
其中,Ki,i=1,2,为待求量,综上模糊闭环系统可表示为:
3)设计控制器:
基于上述模型,定义如下Lyapunov-Krasovskii泛函:
其中,P,Q,R为正定矩阵,对(18)求导,可得:
根据Wirtinger-based积分不等式和倒立凸不等式,将上述公式中的积分项扩大为
其中,
定义:并假设ei=[0n×(i-1)n In×n 0n×(i+1)n],i=1,2,...,6;
由于:
根据(19),(20),(21),得
其中,
Ξ<0可满足此时闭环系统是稳定的;
令G1=εG2,X=G2 -1,(23)式左右两端同乘以diag{X,X,X,X,X,X}T及其转置,
其中:
若要满足条件(24),只需求解如下条件:
其中,
若条件(25)成立,则状态反馈控制增益矩阵可由Ki=YiX-1求出;
4)根据不同的给定的采样周期T,得到不同的时滞上限τ,求出控制器增益矩阵Ki,i=1,2,从而得到状态反馈控制信号:
应用Matlab的LMI Toolbox,令ε=10,τ=1,求解不等式(25),得到一组状态反馈控制增益矩阵:
令ε=10,τ=5,得到一组状态反馈控制增益矩阵:
令ε=10,τ=50,得到一组状态反馈控制增益矩阵:
步骤4求出最终的控制输入信号:
根据步骤2和步骤3所述的前馈控制信号和状态反馈控制信号根据公式最终求出的控制信号。
仿真是根据最大容许的时滞τ分别等于1,5,50三种情况下进行的,令表达式(25)中的预设参数ε=10。对汽包锅炉-汽轮机系统控制方法的仿真结果如图1-5所示。
图1是汽包锅炉-汽轮机系统结构示意图;图2是汽包锅炉-汽轮机系统的控制策略图;图3a-图3h分别代表当τ=1时系统的三个输出的控制情况、跟踪误差以及三个控制输入信号。从图3b、图3d和图3e可见系统的输出可以很好的跟踪给定的参考信号。图4a-图4c分别代表当τ=5时系统的蒸汽压力的跟踪误差曲线、输出功率的跟踪误差曲线和锅炉水位的跟踪误差曲线;图5a-图4c分别代表当τ=50时系统的蒸汽压力的跟踪误差曲线、输出功率的跟踪误差曲线和锅炉水位的跟踪误差曲线。
仿真结果表明,汽包锅炉-汽轮机的前馈加带有时滞的模糊反馈控制方法可以保证跟踪误差快速趋近于零,证明本发明控制方法的有效性和鲁棒性。随着反馈输入信号时滞上界的增大,跟踪误差,特别是蒸汽压力的跟踪误差在一段时间内会明显增大,可见,对系统输出的采样时间间隔越小,控制效果越好。
当然,以上说明仅仅为本发明的部分实施例,本发明并不限于列举上述实施例,应当说明的是,任何熟悉本领域的技术人员在本说明书的教导下,所做出的所有等同替代、明显变形形式,均落在本说明书的实质范围之内,理应受到本发明的保护。
Claims (3)
1.一种汽包锅炉-汽轮机组的模糊跟踪控制方法,其特征在于,包括以下步骤:
步骤1建立汽包锅炉-汽轮机燃油机组的数学模型:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>P</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mn>0.0018</mn>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>P</mi>
<mfrac>
<mn>9</mn>
<mn>8</mn>
</mfrac>
</msup>
<mo>+</mo>
<mn>0.9</mn>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>0.15</mn>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>E</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.073</mn>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>0.016</mn>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>P</mi>
<mfrac>
<mn>9</mn>
<mn>8</mn>
</mfrac>
</msup>
<mo>-</mo>
<mn>0.1</mn>
<mi>E</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>d&rho;</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>85</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mn>141</mn>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1.1</mn>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>0.19</mn>
</mrow>
<mo>)</mo>
</mrow>
<mi>P</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>L</mi>
<mo>=</mo>
<mn>0.05</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.1307</mn>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mn>100</mn>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>9</mn>
</mfrac>
<msub>
<mi>q</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>67.975</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,P,E,ρf为系统状态变量,分别表示蒸汽压力(kg/cm2),输出功率(MW)和锅炉流体密度(kg/cm3);三个控制输入,u1(t),u2(t),u3(t)分别代表燃料阀执行开度、节流阀执行开度和给水阀执行开度,单位为pu;P,E和L代表系统的输出,L表示锅炉水位;
αcs(t)和qe(t)分别为质量系数和蒸发率,且满足
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mn>0.001538</mn>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<mn>0.8</mn>
<mi>P</mi>
<mo>-</mo>
<mn>25.6</mn>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1.0394</mn>
<mo>-</mo>
<mn>0.0012304</mn>
<mi>P</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.854</mn>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>0.147</mn>
</mrow>
<mo>)</mo>
</mrow>
<mi>P</mi>
<mo>+</mo>
<mn>45.59</mn>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2.514</mn>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2.096</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
控制输入为阀门,其上下限和变化率须满足如下条件:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.007</mn>
<mo>&le;</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>0.007</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mo>&le;</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>0.02</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.05</mn>
<mo>&le;</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>0.05</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
锅炉汽轮机的动态模型转化为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mi>y</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0.05</mn>
<mrow>
<mo>(</mo>
<mn>0.1307</mn>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mn>100</mn>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>q</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>9</mn>
</mfrac>
<mo>-</mo>
<mn>67.975</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
步骤2设计前馈控制输入信号:
定义跟踪误差变量:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>y</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>3</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,y1ref(t),y2ref(t),y3ref(t)为给定的跟踪信号,且y3ref(t)≡0;则x1ref(t)=y1ref(t),x2ref(t)=y2ref(t),x3ref(t)为待求量;
设
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
将(6)代入(4)可得:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0.05</mn>
<mrow>
<mo>(</mo>
<mn>0.1307</mn>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mn>100</mn>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>q</mi>
<mo>^</mo>
</mover>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>9</mn>
</mfrac>
<mo>-</mo>
<mn>67.975</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>q</mi>
<mo>^</mo>
</mover>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.854</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>0.147</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>45.59</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2.514</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2.096</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
为反馈控制输入信号;
为前馈控制输入信号,且满足
比较(7)和(4),求出前馈控制输入信号为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.0018</mn>
<mo>&times;</mo>
<mn>0.016</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
<mrow>
<mn>0.073</mn>
<mo>&times;</mo>
<mn>0.9</mn>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.0018</mn>
<mo>&times;</mo>
<mn>0.1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>0.073</mn>
<mo>&times;</mo>
<mn>0.9</mn>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.15</mn>
<mo>&times;</mo>
<mn>1.1</mn>
<mo>&times;</mo>
<mn>0.016</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>141</mn>
<mo>&times;</mo>
<mn>0.073</mn>
<mo>&times;</mo>
<mn>0.9</mn>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.15</mn>
<mo>&times;</mo>
<mn>1.1</mn>
<mo>&times;</mo>
<mn>0.1</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</mrow>
</msup>
</mrow>
<mrow>
<mn>141</mn>
<mo>&times;</mo>
<mn>0.073</mn>
<mo>&times;</mo>
<mn>0.9</mn>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.15</mn>
<mo>&times;</mo>
<mn>0.19</mn>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>141</mn>
<mo>&times;</mo>
<mn>0.9</mn>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>0.016</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>0.073</mn>
</mfrac>
<mo>+</mo>
<mn>0.1</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</mrow>
</msup>
<mo>/</mo>
<mn>0.073</mn>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1.1</mn>
<mo>&times;</mo>
<mn>0.016</mn>
<mo>)</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>141</mn>
<mo>&times;</mo>
<mn>0.073</mn>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1.1</mn>
<mo>&times;</mo>
<mn>0.1</mn>
<mo>)</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</mrow>
</msup>
</mrow>
<mrow>
<mn>141</mn>
<mo>&times;</mo>
<mn>0.073</mn>
</mrow>
</mfrac>
<mo>-</mo>
<mn>0.19</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>141</mn>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3设计反馈控制输入信号:
1)假设在锅炉流体密度达到稳态值时,系统输入信号u(t),输出信号y(t)都达到稳态值,也就是
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
将(8),(9)代入(7)的最后一个表达式,并整理,可得关于x3ref(t)的一个一元二次方程:
λ1x3ref(t)2+λ2x3ref(t)+λ3=0, (10)
其中,
λ1=0.1307,
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>100</mn>
<mo>&times;</mo>
<mn>0.001538</mn>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.8</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>25.6</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1.0394</mn>
<mo>-</mo>
<mn>0.0012304</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>0.854</mn>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mn>0.147</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>45.59</mn>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mrow>
<mn>9</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2.514</mn>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mn>2.096</mn>
</mrow>
<mn>9</mn>
</mfrac>
<mo>-</mo>
<mn>67.975</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>100</mn>
<mrow>
<mo>(</mo>
<mn>0.8</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<mn>25.6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1.0394</mn>
<mo>-</mo>
<mn>0.0012304</mn>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
求解方程(10),并舍弃较小解,可求得
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>4</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
同样的,在处理器中采样y1(t),y2(t),y3(t),我们可求出x3(t)
<mrow>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msup>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>-</mo>
<mn>4</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<msup>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&lambda;</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>100</mn>
<mo>&times;</mo>
<mn>0.001538</mn>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>0.8</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>25.6</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1.0394</mn>
<mo>-</mo>
<mn>0.0012304</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>0.854</mn>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>0.147</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mn>45.59</mn>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mn>9</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2.514</mn>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2.096</mn>
</mrow>
<mn>9</mn>
</mfrac>
<mo>-</mo>
<mn>67.975</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>0.05</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
<mrow>
<msup>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mn>100</mn>
<mrow>
<mo>(</mo>
<mn>0.8</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<mn>25.6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1.0394</mn>
<mo>-</mo>
<mn>0.0012304</mn>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
3
通过处理器计算,可输出任意时刻的状态
2)设计带有时滞的模糊状态反馈控制器
控制器设为:
<mrow>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>K</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
使得
其中,0≤τ(t)≤τ;
模型(14)是一个非线性模型,将其模糊化,可得含有两个线性子系统的T-S模糊模型:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>F</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>50</mn>
<mo>,</mo>
<mi>T</mi>
<mi>H</mi>
<mi>E</mi>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>F</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>150</mn>
<mo>,</mo>
<mi>T</mi>
<mi>H</mi>
<mi>E</mi>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.016</mn>
<mo>&times;</mo>
<msup>
<mn>50</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mn>0.19</mn>
<mn>85</mn>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0.9</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.0018</mn>
<mo>&times;</mo>
<msup>
<mn>50</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.15</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mn>0.073</mn>
<mo>&times;</mo>
<msup>
<mn>50</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1.1</mn>
<mn>85</mn>
</mfrac>
</mrow>
</mtd>
<mtd>
<mfrac>
<mn>141</mn>
<mn>85</mn>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.016</mn>
<mo>&times;</mo>
<msup>
<mn>150</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mn>0.19</mn>
<mn>85</mn>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0.9</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.0018</mn>
<mo>&times;</mo>
<msup>
<mn>150</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.15</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mn>0.073</mn>
<mo>&times;</mo>
<msup>
<mn>150</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1.1</mn>
<mn>85</mn>
</mfrac>
</mrow>
</mtd>
<mtd>
<mfrac>
<mn>141</mn>
<mn>85</mn>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
隶属度函数为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mn>150</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
<mrow>
<msup>
<mn>150</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
<mo>-</mo>
<msup>
<mn>50</mn>
<mfrac>
<mn>1</mn>
<mn>8</mn>
</mfrac>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
模糊状态反馈控制器可表示为:
<mrow>
<mover>
<mi>u</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Ki,i=1,2,为待求量,综上模糊闭环系统可表示为:
<mrow>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<msub>
<mi>K</mi>
<mi>j</mi>
</msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
基于上述模型,定义如下Lyapunov-Krasovskii泛函:
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>P</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
<mn>0</mn>
</munderover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>Q</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
<mo>+</mo>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
<mn>0</mn>
</munderover>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>&theta;</mi>
</mrow>
<mi>t</mi>
</munderover>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>R</mi>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
<mi>d</mi>
<mi>&theta;</mi>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,P,Q,R为正定矩阵,对(18)求导,可得:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mn>2</mn>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>P</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>Q</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>Q</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msup>
<mi>&tau;</mi>
<mn>2</mn>
</msup>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>R</mi>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>&tau;</mi>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
<mn>0</mn>
</munderover>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>R</mi>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mi>d</mi>
<mi>s</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
根据Wirtinger-based积分不等式和倒立凸不等式,将上述公式中的积分项扩大为
其中,
定义:
<mrow>
<mi>&xi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mfrac>
<mn>1</mn>
<mrow>
<mi>&tau;</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</munderover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
<mfrac>
<mn>1</mn>
<mrow>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>t</mi>
</munderover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>,</mo>
</mrow>
并假设ei=[0n×(i-1)n In×n 0n×(i+1)n],i=1,2,...,6;
由于:
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mi>B</mi>
<mi>K</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
根据(19),(20),(21),得
<mrow>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>&xi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>&Xi;</mi>
<mi>&xi;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
Ξ<0可满足此时闭环系统是稳定的;
令G1=εG2,X=G2 -1,(23)式左右两端同乘以diag{X,X,X,X,X,X}T及其转置,得:
其中,
<mrow>
<mi>Y</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>Y</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mi>X</mi>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<mi>P</mi>
<mi>X</mi>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msup>
<mi>X</mi>
<mi>T</mi>
</msup>
<mi>Q</mi>
<mi>X</mi>
<mo>;</mo>
</mrow>
若要满足条件(24),只需求解如下条件:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&Xi;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo><</mo>
<mn>0</mn>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&Xi;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>&Xi;</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>j</mi>
<mi>i</mi>
</mrow>
</msub>
<mo><</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mn>2</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
若条件(25)成立,则状态反馈控制增益矩阵可由Ki=YiX-1求出;
3)根据不同的给定的采样周期T,得到不同的时滞上限τ,求出控制器增益矩阵Ki,i=1,2,从而得到状态反馈控制信号:
步骤4求出最终的控制输入信号u(t):
根据步骤2和步骤3所述的前馈控制信号和状态反馈控制信号根据公式最终求出的控制信号。
2.根据权利要求1所述的汽包锅炉-汽轮机组的控制方法,其特征在于:所述的方法还包括基于李雅普诺夫稳定理论对控制器设计的数值仿真验证;
应用数值仿真来验证由于采样周期的不同导致的不同的时滞输入的控制方法的有效性,通过数值分析可见,采样周期越小,则跟踪误差越小,控制效果越好。
3.根据权利要求1所述的汽包锅炉-汽轮机组的控制方法,其特征在于:所述的仿真验证,三个状态蒸汽压力、输出功率和锅炉流体密度的初始值选择为[P L ρf]=[75.6 15.27299.6]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710116314.1A CN107065538B (zh) | 2017-03-01 | 2017-03-01 | 汽包锅炉-汽轮机组的模糊跟踪控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710116314.1A CN107065538B (zh) | 2017-03-01 | 2017-03-01 | 汽包锅炉-汽轮机组的模糊跟踪控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107065538A true CN107065538A (zh) | 2017-08-18 |
CN107065538B CN107065538B (zh) | 2019-11-08 |
Family
ID=59622878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710116314.1A Expired - Fee Related CN107065538B (zh) | 2017-03-01 | 2017-03-01 | 汽包锅炉-汽轮机组的模糊跟踪控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107065538B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113110357A (zh) * | 2021-05-10 | 2021-07-13 | 山东大学 | 一种火力发电单元机组协调控制方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101725999A (zh) * | 2009-12-11 | 2010-06-09 | 中国电力科学研究院 | 一种生物质振动炉排锅炉的负荷控制方法 |
CN101922708A (zh) * | 2010-06-01 | 2010-12-22 | 广东电网公司电力科学研究院 | 一种基于强化燃烧的大型循环流化床机组协调控制装置 |
US20150017591A1 (en) * | 2013-07-02 | 2015-01-15 | General Electric Company | Systems and methods for advanced closed loop control and improvement of combustion system operation |
-
2017
- 2017-03-01 CN CN201710116314.1A patent/CN107065538B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101725999A (zh) * | 2009-12-11 | 2010-06-09 | 中国电力科学研究院 | 一种生物质振动炉排锅炉的负荷控制方法 |
CN101922708A (zh) * | 2010-06-01 | 2010-12-22 | 广东电网公司电力科学研究院 | 一种基于强化燃烧的大型循环流化床机组协调控制装置 |
US20150017591A1 (en) * | 2013-07-02 | 2015-01-15 | General Electric Company | Systems and methods for advanced closed loop control and improvement of combustion system operation |
Non-Patent Citations (4)
Title |
---|
CHEN JIAN,等: "Stabilisation of T-S Fuzzy Systems via Static Output Feedback:An Iterative Method", 《PROCEEDINGS OF THE 35TH CHINESE CONTROL CONFERENCE》 * |
CHONG LIN,等: "Fuzzy Normalization and Stabilization for a Class of Nonlinear Rectangular Descriptor Systems", 《NEUROCOMPUTING》 * |
刘臻臻,等: "基于自由矩阵型不等式的分布时滞系统稳定性分析", 《青岛大学学报(工程技术版)》 * |
张爱芳, 林崇, 孙立荣: "基于不同模糊权的时滞T-S模糊系统的能稳控制器设计", 《青岛大学学报(工程技术版)》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113110357A (zh) * | 2021-05-10 | 2021-07-13 | 山东大学 | 一种火力发电单元机组协调控制方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107065538B (zh) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Extended state observer based fuzzy model predictive control for ultra-supercritical boiler-turbine unit | |
CN104122797B (zh) | 一种新型火电机组负荷多变量预测控制方法 | |
Kong et al. | Nonlinear multivariable hierarchical model predictive control for boiler-turbine system | |
Fang et al. | Backstepping-based nonlinear adaptive control for coal-fired utility boiler–turbine units | |
Daren et al. | Nonlinear coordinated control of drum boiler power unit based on feedback linearization | |
Prasad et al. | Plant-wide predictive control for a thermal power plant based on a physical plant model | |
CN104714526B (zh) | 基于凝结水节流调节预估的负荷控制系统及方法 | |
CN107831652A (zh) | 一种基于冷端系统储能的机组负荷智能优化控制方法 | |
Wang et al. | An improved coordinated control strategy for boiler-turbine units supplemented by cold source flow adjustment | |
Wang et al. | Flexible electric power control for coal-fired units by incorporating feedwater bypass | |
CN107065538A (zh) | 汽包锅炉‑汽轮机组的模糊跟踪控制方法 | |
Carrasco et al. | Novel modelling for a steam boiler under fast load dynamics with implications to control | |
Sheng et al. | Sliding mode controller with AGA for drum water level of ship boiler | |
Ghabraei et al. | Multivariable robust regulation of an industrial boiler-turbine with model uncertainties | |
Liu et al. | Adaptive Fuzzy PID Control Applied in Boiler System of Power Plant | |
Wang et al. | Design and realization of automatic control system for boiler based on model free adaptive control | |
Wu et al. | Fuzzy modeling and predictive control of power plant steam temperature system | |
HABBI et al. | Design of a fuzzy model-based controller for a drum boiler-turbine system | |
Moradi et al. | Multivariable control of an industrial boiler-turbine unit with nonlinear model: A comparison between gain scheduling & feedback linearization approaches | |
Ma et al. | Neural network inverse control for the coordinated system of a 600MW supercritical boiler unit | |
CN106125781A (zh) | 一种基于LMIs的汽包水位控制系统设计方法 | |
Li et al. | MIMO CCS model identification based on GA for closed-loop field data | |
Hou et al. | Nonlinear predictive control of a boiler-turbine system based on TS fuzzy model | |
Ma et al. | Research on Control Strategy of Main Steam Temperature System of Thermal Power Plant Boilers | |
Granberg | Control of a process with large time constants and significant time delay |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191108 |
|
CF01 | Termination of patent right due to non-payment of annual fee |