CN107908106A - 二次再热机组再热汽温的自降阶多回路集中预估控制*** - Google Patents
二次再热机组再热汽温的自降阶多回路集中预估控制*** Download PDFInfo
- Publication number
- CN107908106A CN107908106A CN201711057680.0A CN201711057680A CN107908106A CN 107908106 A CN107908106 A CN 107908106A CN 201711057680 A CN201711057680 A CN 201711057680A CN 107908106 A CN107908106 A CN 107908106A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- msup
- 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
-
- 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/048—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 using a predictor
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种二次再热机组再热汽温的自降阶多回路集中预估控制***,该控制***由烟气再循环风机转速控制回路和烟气挡板开度控制回路构成,将一次再热汽温和二次再热汽温与设定值的偏差送入预估控制器,计算得出烟气再循环风机转速和烟气挡板开度的优化值,保证机组一次再热汽温和二次再热汽温维持在合理安全的范围内,本发明采用预估控制能够方便地处理大惯性多变量***的优化问题,计算过程清晰、简单,工程应用时,编程实施方便,对未来输出偏差的预测能及时调节相应的控制量,一次再热汽温和二次再热汽温能稳定在合理的范围之内,保证了机组的稳定性和安全性,控制效果较传统PID控制好。
Description
技术领域
本发明涉及一种二次再热机组再热汽温的自降阶多回路集中预估控制***,属于热 能动力工程和自动控制领域。
背景技术
随着对火电机组经济性要求进一步的提高,机组不断向更高参数发展,并采用二次 再热技术,进一步提高机组效率。与此同时,机组的复杂性进一步提高,受控对象的惯性增大,对再热汽温的控制要求更高。二次再热机组再热汽温的控制一般通过调节烟气 再循环风机转速和烟气挡板开度,减温水调节一般用于事故喷水情况。超临界机组*** 的复杂性要求控制精度更高,控制动作更快。***的惯性大,对象模型阶次高,若不进 行降阶处理,对采样周期的要求较为苛刻。若采用传统PID控制,达不到理想的效果, 这就需要我们探索其他控制方案。
发明内容
本发明二次再热超超临界机组的再热汽温控制过程具有大惯性、大滞后的特点提出 了一种二次再热机组再热汽温的自降阶多回路集中预估控制***,预估控制能通过预测 未来的输出计算当前时刻的最优控制量,***的鲁棒性较好,同时控制效果良好,其原理清晰,易于在线计算,控制效果良好。
为解决上述技术问题,本发明提供了一种二次再热机组再热汽温的自降阶多回路集 中预估控制***,以烟气再循环风机转速和烟气挡板开度作为输入,以一次再热汽温和二次再热汽温作为输出,包括烟气再循环风机转速控制回路、烟气挡板开度控制回路, 该预估控制***所采用的控制器为预估控制器,分别将一次再热汽温的设定值和二次再 热汽温的设定值与所述预估控制器预测的输出的偏差送入所述预估控制器,该预估控制 器输出烟气再循环风机转速和烟气挡板开度的优化控制序列,所述优化控制序列中取当 前时刻的控制作用作用于经降阶的实际对象模型,得到一次再热汽温、二次再热汽温的 实际输出,下一时刻继续相同计算,实现滚动优化。
所述预估控制器包括:
预测模块,用于预测未来各个采样时刻的输出;
优化性能指标计算模块,用于根据设定的性能指标计算控制范围内最优的控制增量;
控制实施模块,用于将计算所得最优控制序列应用于***。
进一步的,预估控制器所依赖的对象模型由实验数据拟合得出,通过在多个负荷点 做阶跃响应试验,建立各负荷点上的线性传递函数模型,中间负荷的模型通过已建立的相邻负荷点上的线性传递函数模型通过插值的方法计算得出,模型如下:
对对象模型进行Pade近似法降阶,对
其幂级数展开式为
其中,
代入
pij0=Cij0
pij1=Cij1+Cij0qij1
pijl=Cijl+Cijl-1qij0+…Cij0qijl
Cijl+1+Cijlqij1+…+Cijl-k+1qijk=0
Cijl+2+Cijl+1qij1+…+Cijl-k+2qijk=0
Cijl+k+Cijl+k-1qij1+…+Cijlqijk=0 (5)
解得pijs(s=1,2,…l),qijt(t=1,2…k),降阶模型如下:
根据降阶后的模型推导集中式预估控制器所依据的可控自回归积分滑动平均模型, 即CARIMA模型,模型如下:
其中,
模型可转化成:
通过预估控制器的预测模型为:
Y=F1ΔU+F2ΔU(k-j)+GY(k) (10)
其中,
Y=[y1(k+1)…y1(k+N)y2(k+1)…y2(k+N)]T
ΔU=[Δu1(k) … Δu1(k+Nu-1) Δu2(k) … Δu2(k+Nu-1)]T
ΔU(k-j)=[Δu1(k-1) … Δu1(k-nb1) Δu2(k-1) … Δu2(k-nb2)]T
Y(k)=[y1(k) … y1(k-na1) y2(k) … y2(k-na2)]T
nb1=max(nb11,nb21),nb2=max(nb12,nb22),na1=na11,na2=na22
N为预测时域,Nu为控制时域,Yp=F2ΔU(k-j)+GY(k)为基于过去输入输出的输 出预测响应;
F1、F2、G的求解通过求解如下Diophantine方程:
i,j=1,2,l为预测步数
其中,
进一步的,预估控制器的优化性能指标为:
J=[F1ΔU+F2ΔU(k-j)+GY(k)-Yr]T[F1ΔU+F2ΔU(k-j)+GY(k)-Yr]+ΔUTΓΔU
(11)
其中,Yr=[y1r(k+1) … y1r(k+N) y2r(k+1) … y2r(k+N)]T为设定值,Γ为 控制权矩阵;
时,所述预估控制器的最优控制增量为:
ΔU(k)=(F1 TF1+Γ)-1F1 T[Yr-F2ΔU(k-j)-GY(k)] (12)。
进一步的,控制实施模块取计算得到的最优控制增量序列中当前时刻k的控制增量 作用于***:
uj(k)=uj(k-1)+Δuj(k),j=1,2 (13)
再以k+1时刻为基点进行下一时刻的最优控制增量序列计算,实现滚动优化。
有益效果:本发明与现有技术相比,本发明提出的一种二次再热机组再热汽温的自降阶多回路集中预估控制***,具有以下优点:
1、预估控制器的算法计算过程清晰、简单,工程应用时,编程实施非常方便;
2、能够适用于被控过程是线性模型或非线性模型情况,对象具有大惯性、大滞后等特点时也能保持良好的控制效果;
3、能保证一次再热汽温和二次再热汽温维持在稳定安全的范围内;
4、控制器计算得到的烟气再循环风机转速和烟气挡板开度在最优范围内,且变化幅度合理,不会使得一次再热汽温和二次再热汽温有较大的波动,提高***经济性的同 时保证了安全性。
附图说明
图1为本发明实施例控制***示意图。
具体实施方式
下面结合实施例进一步阐述该发明方法。
如图1所示为一种二次再热机组再热汽温的自降阶多回路集中预估控制***原理图,再热汽温控制***包括烟气再循环风机转速控制回路和烟气挡板开度控制回路,该 再热汽温控制***的输入量为烟气再循环风机转速和烟气挡板开度,其输出量为一次再 热汽温和二次再热汽温,输入、输出量之间存在耦合,对象惯性较大,采用可降阶的预 估控制算法,图1中:输入数据即烟气再循环风机转速、烟气挡板开度;输出数据为一 次再热汽温、二次再热汽温。根据历史输入输出数据拟合出对象模型并进行降阶处理。 y1r,y2r分别为一次再热汽温、二次再热汽温的设定值,为根据过去及当前时刻的 输入和输出以及未来的输入计算得到的未来输出预测值,u1,u2分别为根据优化性能指 标计算得到的烟气再循环风机、烟气挡板开度的最优控制序列,y1,y2为一次再热汽温、 二次再热汽温的模型实际输出。
再热汽温控制***对象模型由实验数据拟合得出,通过在多个负荷点做阶跃响应试 验,建立各负荷点上的线性传递函数模型,中间负荷的模型通过已建立的相邻负荷点上的线性传递函数模型通过插值的方法计算得出,模型如下:
对象模型采用Pade近似法进行降阶:对
将其幂级数展开,
其中,
代入
pij0=Cij0
pij1=Cij1+Cij0qij1
pijl=Cijl+Cijl-1qij0+…Cij0qijl
Cijl+1+Cijlqij1+…+Cijl-k+1qijk=0
Cijl+2+Cijl+1qij1+…+Cijl-k+2qijk=0
Cijl+k+Cijl+k-1qij1+…+Cijlqijk=0 (5)
解得pijs(s=1,2,…l),qijt(t=1,2…k),降阶模型如下:
根据降阶后的模型推导集中式预估控制器所依据的可控自回归积分滑动平均模型, 即CARIMA模型,模型如下:
其中,
模型可转化成:
通过预估控制器的预测模块,预测未来时刻的输出,预测模型为:
Y=F1ΔU+F2ΔU(k-j)+GY(k) (10)
其中,
Y=[y1(k+1) … y1(k+N) y2(k+1) … y2(k+N)]T
ΔU=[Δu1(k) … Δu1(k+Nu-1) Δu2(k) … Δu2(k+Nu-1)]T
ΔU(k-j)=[Δu1(k-1) … Δu1(k-nb1) Δu2(k-1) … Δu2(k-nb2)]T
Y(k)=[y1(k) … y1(k-na1) y2(k) … y2(k-na2)]T
nb1=max(nb11,nb21),nb2=max(nb12,nb22),na1=na11,na2=na22
N为预测时域,Nu为控制时域,Yp=F2ΔU(k-j)+GY(k)为基于过去输入输出的输出预 测响应。
F1、F2、G的求解通过求解如下Diophantine方程:
i,j=1,2,l为预测步数
其中,
确定所述控制器的优化性能指标计算模块,并计算控制范围内最优的控制增量,优 化性能指标为:
J=[F1ΔU+F2ΔU(k-j)+GY(k)-Yr]T[F1ΔU+F2ΔU(k-j)+GY(k)-Yr]+ΔUTΓΔU
(11)
其中,Yr=[y1r(k+1) … y1r(k+N) y2r(k+1) … y2r(k+N)]T为设定值,Γ为控制权矩阵。根据控制量所允许的波动范围确定控制加权系数。
时,所得控制增量即为最优控制增量:
ΔU(k)=(F1 TF1+Γ)-1F1 T[Yr-F2ΔU(k-j)-GY(k)] (12)
通过预估控制器的控制实施模块,将计算所得最优控制律应用于***,计算得到的 最优控制增量序列中当前时刻k的控制增量为:
uj(k)=uj(k-1)+Δuj(k),j=1,2 (13)
再以k+1时刻为基点进行下一时刻的最优控制增量序列计算,实现滚动优化。
下面以某电厂660MW超超临界二次再热机组采用本发明的优化控制***为例,详细说明本发明内容。
采样周期取为10s,预测时域N取750,控制时域Nu取10,控制权矩阵取其中 在650MW负荷附近,一次再热汽温、二次再热汽温分 别为601.2℃、601.5℃,烟气再循环风机转速为72%,烟气挡板开度为50%,将一次再 热汽温、二次再热汽温设定值分别阶跃增加4℃,结果显示,一次再热汽温、二次再热 汽温能快速跟随设定值,其中一次再热汽温最大偏差为0.5℃,二次再热汽温最大偏差 为0.7℃,偏差很小。
以上实例表明:本实施例的二次再热机组再热汽温的自降阶多回路集中预估控制系 统,能有效改善二次再热机组再热汽温控制***的控制性能,一次再热汽温、二次再热汽温快速响应设定值变化,波动小,并维持在安全范围内,机组经济性和安全性均得到 保障。
Claims (4)
1.一种二次再热机组再热汽温的自降阶多回路集中预估控制***,其特征在于:以烟气再循环风机转速和烟气挡板开度作为输入,以一次再热汽温和二次再热汽温作为输出,包括烟气再循环风机转速控制回路、烟气挡板开度控制回路,该预估控制***所采用的控制器为预估控制器,分别将一次再热汽温的设定值和二次再热汽温的设定值与所述预估控制器预测的输出的偏差送入所述预估控制器,该预估控制器输出烟气再循环风机转速和烟气挡板开度的优化控制序列,所述优化控制序列中取当前时刻的控制作用作用于经降阶的实际对象模型,得到一次再热汽温、二次再热汽温的实际输出,下一时刻继续相同计算,实现滚动优化;
所述预估控制器包括:
预测模块,用于预测未来各个采样时刻的输出;
优化性能指标计算模块,用于根据设定的性能指标计算控制范围内最优的控制增量;
控制实施模块,用于将计算所得最优控制序列应用于***。
2.根据权利要求1所述的一种二次再热机组再热汽温的自降阶多回路集中预估控制***,其特征在于:所述预估控制器所依赖的对象模型由实验数据拟合得出,通过在多个负荷点做阶跃响应试验,建立各负荷点上的线性传递函数模型,中间负荷的模型通过已建立的相邻负荷点上的线性传递函数模型通过插值的方法计算得出,模型如下:
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
对对象模型进行Pade近似法降阶,对
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>s</mi>
<mo>+</mo>
<mn>...</mn>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>m</mi>
</mrow>
</msub>
<msup>
<mi>s</mi>
<mi>m</mi>
</msup>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>s</mi>
<mo>+</mo>
<mn>...</mn>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>n</mi>
</mrow>
</msub>
<msup>
<mi>s</mi>
<mi>n</mi>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>&le;</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其幂级数展开式为
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>p</mi>
</mrow>
</msub>
<msup>
<mi>S</mi>
<mi>p</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>p</mi>
<mo>!</mo>
</mrow>
</mfrac>
<msup>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
代入
pij0=Cij0
pij1=Cij1+Cij0qij1
pijl=Cijl+Cijl-1qij0+…Cij0qijl
Cijl+1+Cijlqij1+…+Cijl-k+1qijk=0
Cijl+2+Cijl+1qij1+…+Cijl-k+2qijk=0
Cijl+k+Cijl+k-1qij1+…+Cijlqijk=0 (5)
解得pijs(s=1,2,…l),qijt(t=1,2…k),降阶模型如下:
<mrow>
<msubsup>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>s</mi>
<mo>+</mo>
<mo>...</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>l</mi>
</mrow>
</msub>
<msup>
<mi>s</mi>
<mi>l</mi>
</msup>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>s</mi>
<mo>+</mo>
<mo>...</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<msup>
<mi>s</mi>
<mi>k</mi>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>&le;</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
根据降阶后的模型推导集中式预估控制器所依据的可控自回归积分滑动平均模型,即CARIMA模型,模型如下:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&xi;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&xi;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>/</mo>
<mi>&Delta;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>...</mn>
<msub>
<mi>a</mi>
<mrow>
<msub>
<mi>ijn</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>...</mn>
<msub>
<mi>b</mi>
<mrow>
<msub>
<mi>ijn</mi>
<mrow>
<mi>b</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>b</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
模型可转化成:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>&Delta;</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>&Delta;</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&xi;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&xi;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
通过预估控制器的预测模型为:
Y=F1ΔU+F2ΔU(k-j)+GY(k) (10)
其中,
Y=[y1(k+1) … y1(k+N) y2(k+1) … y2(k+N)]T
ΔU=[Δu1(k) … Δu1(k+Nu-1) Δu2(k) … Δu2(k+Nu-1)]T
ΔU(k-j)=[Δu1(k-1) … Δu1(k-nb1) Δu2(k-1) … Δu2(k-nb2)]T
Y(k)=[y1(k) … y1(k-na1) y2(k) … y2(k-na2)]T
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mn>111</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>F</mi>
<mn>121</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mn>211</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>F</mi>
<mn>221</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mn>112</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>F</mi>
<mn>122</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mn>212</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>F</mi>
<mn>222</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>G</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>G</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>G</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
nb1=max(nb11,nb21),nb2=max(nb12,nb22),na1=na11,na2=na22
N为预测时域,Nu为控制时域,Yp=F2ΔU(k-j)+GY(k)为基于过去输入输出的输出预测响应;
F1、F2、G的求解通过求解如下Diophantine方程:
i,j=1,2,l为预测步数
其中,
<mrow>
<msub>
<mi>E</mi>
<mrow>
<msub>
<mi>ii</mi>
<mi>l</mi>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>e</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>...</mn>
<msub>
<mi>e</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<msub>
<mi>ii</mi>
<mi>l</mi>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>...</mn>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>g</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>g</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<msub>
<mi>ij</mi>
<mi>l</mi>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>...</mn>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>f</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</msub>
<msup>
<mi>q</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>f</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>u</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>u</mi>
</msub>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<msub>
<mi>ij</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>e</mi>
<mrow>
<msub>
<mi>ii</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>N</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的一种二次再热机组再热汽温的自降阶多回路集中预估控制***,其特征在于:所述预估控制器的优化性能指标为:
J=[F1ΔU+F2ΔU(k-j)+GY(k)-Yr]T[F1ΔU+F2ΔU(k-j)+GY(k)-Yr]+ΔUTΓΔU (11)
其中,Yr=[y1r(k+1) … y1r(k+N) y2r(k+1) … y2r(k+N)]T为设定值,Γ为控制权矩阵;
时,所述预估控制器的最优控制增量为:
ΔU(k)=(F1 TF1+Γ)-1F1 T[Yr-F2ΔU(k-j)-GY(k)] (12)。
4.根据权利要求3所述的一种二次再热机组再热汽温的自降阶多回路集中预估控制***,其特征在于:所述控制实施模块取计算得到的最优控制增量序列中当前时刻k的控制增量作用于***:
uj(k)=uj(k-1)+Δuj(k),j=1,2 (13)
再以k+1时刻为基点进行下一时刻的最优控制增量序列计算,实现滚动优化。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711057680.0A CN107908106B (zh) | 2017-11-01 | 2017-11-01 | 二次再热机组再热汽温的自降阶多回路集中预估控制*** |
PCT/CN2018/088141 WO2019085446A1 (zh) | 2017-11-01 | 2018-05-24 | 二次再热机组再热汽温的自降阶多回路集中预估控制*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711057680.0A CN107908106B (zh) | 2017-11-01 | 2017-11-01 | 二次再热机组再热汽温的自降阶多回路集中预估控制*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107908106A true CN107908106A (zh) | 2018-04-13 |
CN107908106B CN107908106B (zh) | 2019-11-12 |
Family
ID=61842417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711057680.0A Active CN107908106B (zh) | 2017-11-01 | 2017-11-01 | 二次再热机组再热汽温的自降阶多回路集中预估控制*** |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN107908106B (zh) |
WO (1) | WO2019085446A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108762086A (zh) * | 2018-06-19 | 2018-11-06 | 中国大唐集团科学技术研究院有限公司华东分公司 | 基于模型预测控制的二次再热蒸汽温度控制装置及控制*** |
CN109654475A (zh) * | 2018-11-26 | 2019-04-19 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种二次再热机组再热汽温控制方法 |
WO2019085446A1 (zh) * | 2017-11-01 | 2019-05-09 | 东南大学 | 二次再热机组再热汽温的自降阶多回路集中预估控制*** |
CN111123704A (zh) * | 2019-12-09 | 2020-05-08 | 国家能源集团泰州发电有限公司 | 一种二次再热机组再热汽温的节能型预测控制方法 |
CN111561693A (zh) * | 2019-11-29 | 2020-08-21 | 上海圆曦电力科技有限公司 | 一种电站锅炉状态空间汽温控制的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102387311A (zh) * | 2010-09-02 | 2012-03-21 | 深圳Tcl新技术有限公司 | 一种合成视频图像的装置以及合成视频图像的方法 |
CN102566433A (zh) * | 2012-02-08 | 2012-07-11 | 上海锅炉厂有限公司 | 二次再热直流锅炉-汽轮机dbc协调控制*** |
US20120259437A1 (en) * | 2011-04-08 | 2012-10-11 | General Cybernation Group Inc. | Model-free adaptive control of advanced power plants |
CN104075309A (zh) * | 2014-06-30 | 2014-10-01 | 章礼道 | 燃用准东煤的二次再热电站锅炉 |
CN106019939A (zh) * | 2016-06-20 | 2016-10-12 | 东南大学 | 基于非线性智能预测控制技术的超临界机组再热汽温优化控制*** |
CN106169108A (zh) * | 2016-07-14 | 2016-11-30 | 河海大学 | 含电池储能***的主动配电网短期有功优化方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10301916A (ja) * | 1997-04-23 | 1998-11-13 | Ishikawajima Harima Heavy Ind Co Ltd | 再熱器ガスダンパ開度制御装置 |
JP2000257801A (ja) * | 1999-03-10 | 2000-09-22 | Babcock Hitachi Kk | プロセスの制御方法と流動層ボイラ及びその燃料供給方法 |
CN106642069B (zh) * | 2016-10-08 | 2018-09-25 | 西安西热控制技术有限公司 | 超超临界二次再热发电机组再热汽温复合控制策略 |
CN106871097A (zh) * | 2017-01-10 | 2017-06-20 | 上海电力学院 | 电站两次再热机组的锅炉烟气余热回收*** |
CN107908106B (zh) * | 2017-11-01 | 2019-11-12 | 东南大学 | 二次再热机组再热汽温的自降阶多回路集中预估控制*** |
-
2017
- 2017-11-01 CN CN201711057680.0A patent/CN107908106B/zh active Active
-
2018
- 2018-05-24 WO PCT/CN2018/088141 patent/WO2019085446A1/zh active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102387311A (zh) * | 2010-09-02 | 2012-03-21 | 深圳Tcl新技术有限公司 | 一种合成视频图像的装置以及合成视频图像的方法 |
US20120259437A1 (en) * | 2011-04-08 | 2012-10-11 | General Cybernation Group Inc. | Model-free adaptive control of advanced power plants |
CN102566433A (zh) * | 2012-02-08 | 2012-07-11 | 上海锅炉厂有限公司 | 二次再热直流锅炉-汽轮机dbc协调控制*** |
CN104075309A (zh) * | 2014-06-30 | 2014-10-01 | 章礼道 | 燃用准东煤的二次再热电站锅炉 |
CN106019939A (zh) * | 2016-06-20 | 2016-10-12 | 东南大学 | 基于非线性智能预测控制技术的超临界机组再热汽温优化控制*** |
CN106169108A (zh) * | 2016-07-14 | 2016-11-30 | 河海大学 | 含电池储能***的主动配电网短期有功优化方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019085446A1 (zh) * | 2017-11-01 | 2019-05-09 | 东南大学 | 二次再热机组再热汽温的自降阶多回路集中预估控制*** |
CN108762086A (zh) * | 2018-06-19 | 2018-11-06 | 中国大唐集团科学技术研究院有限公司华东分公司 | 基于模型预测控制的二次再热蒸汽温度控制装置及控制*** |
CN108762086B (zh) * | 2018-06-19 | 2021-03-23 | 中国大唐集团科学技术研究院有限公司华东分公司 | 基于模型预测控制的二次再热蒸汽温度控制装置及控制*** |
CN109654475A (zh) * | 2018-11-26 | 2019-04-19 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种二次再热机组再热汽温控制方法 |
CN109654475B (zh) * | 2018-11-26 | 2020-09-04 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种二次再热机组再热汽温控制方法 |
CN111561693A (zh) * | 2019-11-29 | 2020-08-21 | 上海圆曦电力科技有限公司 | 一种电站锅炉状态空间汽温控制的方法 |
CN111561693B (zh) * | 2019-11-29 | 2021-12-17 | 上海圆曦电力科技有限公司 | 一种电站锅炉状态空间汽温控制的方法 |
CN111123704A (zh) * | 2019-12-09 | 2020-05-08 | 国家能源集团泰州发电有限公司 | 一种二次再热机组再热汽温的节能型预测控制方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2019085446A1 (zh) | 2019-05-09 |
CN107908106B (zh) | 2019-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107908106A (zh) | 二次再热机组再热汽温的自降阶多回路集中预估控制*** | |
CN104122797B (zh) | 一种新型火电机组负荷多变量预测控制方法 | |
CN107045289A (zh) | 一种电加热炉温度的非线性神经网络优化pid控制方法 | |
CN107515598A (zh) | 基于多参数动态矩阵控制的火电机组分布式协调控制*** | |
Zhou et al. | Adaptive NN control for nonlinear systems with uncertainty based on dynamic surface control | |
CN101763035A (zh) | Rbf神经网络整定pid与模糊免疫控制方法 | |
CN106019939A (zh) | 基于非线性智能预测控制技术的超临界机组再热汽温优化控制*** | |
CN108167802B (zh) | 低负荷下锅炉负荷的多模型智能寻优预测控制方法 | |
CN110501909A (zh) | 基于扰动观测器的增强鲁棒特性的模糊预测控制方法 | |
Wang et al. | An improved coordinated control strategy for boiler-turbine units supplemented by cold source flow adjustment | |
Chen et al. | Robust adaptive control of maximum power point tracking for wind power system | |
CN113091088B (zh) | 一种基于双阶段神经网络的锅炉燃烧广义预测控制方法 | |
Zhou et al. | Study on PID parameters tuning based on particle swarm optimization | |
CN114326387A (zh) | 一种火电机组脱硝控制装置及控制方法 | |
Hou et al. | Parameters optimization of ADRC based on adaptive CPSO algorithm and its application in main-steam temperature control system | |
Guolian et al. | Multiple-model predictive control based on fuzzy adaptive weights and its application to main-steam temperature in power plant | |
CN107918282B (zh) | 用于互联电网自动发电控制的可变域mpc方法 | |
Wang et al. | Neuro-PID control of heat exchanger in an Organic Rankine Cycle system for waste heat recovery | |
El Mjabber et al. | An adaptive control for a variable speed wind turbine using RBF neural network | |
Bin et al. | Pitch angle control based on renforcement learning | |
Ma et al. | ANN-based real-time parameter optimization via GA for superheater model in power plant simulator | |
CN110442028B (zh) | 基于模糊预测的反分叉控制方法 | |
Chen et al. | Experimental analysis of control loops with different delay times in the supply air system of a radiator test rig | |
Wang et al. | Study on Fuzzy immune PID Controller for main steam pressure control system in marine steam power plant | |
Dolezel et al. | Piecewise-linear neural model for helicopter elevation control |
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 |