CN100483276C - 一种基于噪声估计的自适应状态反馈预测控制方法 - Google Patents
一种基于噪声估计的自适应状态反馈预测控制方法 Download PDFInfo
- Publication number
- CN100483276C CN100483276C CNB2006100833311A CN200610083331A CN100483276C CN 100483276 C CN100483276 C CN 100483276C CN B2006100833311 A CNB2006100833311 A CN B2006100833311A CN 200610083331 A CN200610083331 A CN 200610083331A CN 100483276 C CN100483276 C CN 100483276C
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mtr
- mover
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 143
- 230000009471 action Effects 0.000 claims abstract description 30
- 230000008569 process Effects 0.000 claims description 63
- 239000011159 matrix material Substances 0.000 claims description 46
- 230000003044 adaptive effect Effects 0.000 claims description 20
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 abstract 1
- 239000000446 fuel Substances 0.000 description 13
- 238000010586 diagram Methods 0.000 description 12
- 230000001629 suppression Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 8
- 239000001257 hydrogen Substances 0.000 description 7
- 229910052739 hydrogen Inorganic materials 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000001301 oxygen Substances 0.000 description 4
- 229910052760 oxygen Inorganic materials 0.000 description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- -1 hydrogen ions Chemical class 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000012528 membrane Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 239000003570 air Substances 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000006722 reduction reaction Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于噪声估计的自适应状态反馈预测控制方法,包括以下步骤:A、利用Singular Pencil Model模型表示实际被控过程;B、采用递推的扩展卡尔曼滤波参数估计方法估计所述利用Singular PencilModel模型表示的实际被控过程的状态和模型参数;C、利用估计出来的状态和模型参数计算当前的控制作用,并将计算出来的控制作用加到实际被控过程中;D、在下一个执行周期,重复执行步骤B和C,实现对实际被控过程的自适应状态反馈预测控制。利用本发明,充分考虑了被控过程中含有的噪声信息和状态信息,大大提高了控制方法抑制噪声的能力。另外,本发明提供的控制方法计算量小,便于工业应用。
Description
技术领域
本发明涉及自动化控制技术领域,尤其涉及一种基于噪声估计的自适应状态反馈预测控制方法。
背景技术
在工业自动化控制领域,应用预测控制技术对被控过程进行控制具有广泛的应用。
现有的应用预测控制技术多是基于离线辨识的模型,或是根据物理和化学反应机理推导出来的模型,例如模型算法控制(Model AlgorithmControl)、动态矩阵控制(Dynamic Matrix Control)和广义预测控制(GeneralPredictive Control)等。另外,现有的应用预测控制技术也有将在线辨识和上述预测控制相结合的自适应预测控制方法,例如自适应的广义预测控制,模型参考自适应预测控制方法等。
在上述各种现有的应用预测控制技术中,由于噪声信息难以确定,尤其是工业上被控过程中含有不可测的噪声信息,导致现有的控制方法无法充分考虑被控过程中含有的噪声信息。即使个别现有的控制方法考虑了被控过程中含有的噪声信息,考虑的也大都是可测噪声信息,而在工业上被控过程中含有的噪声信息多是不可测噪声信息,所以现有的控制方法并没有充分考虑被控过程中含有的不可测噪声信息。
另外,由于大部分应用预测控制技术都没有引入状态反馈,所以现有的控制方法都没有充分考虑被控过程中含有的状态信息。即使个别现有的方法虽引入了状态反馈,例如状态反馈预测控制,但是受到状态是否可测的限制,如果状态是可测的,则可以直接利用测出的状态信息进行状态反馈预测控制;如果状态是不可测的,则还要单独设计状态观测器,利用设计的状态观测器对状态信息进行观测,然后再利用测出的状态信息进行状态反馈预测控制,实现过程非常复杂,实现成本也很高。
因此,基于上述两个原因,即现有的控制方法没有充分考虑被控过程中含有的噪声信息和状态信息,导致了现有的控制方法抑制噪声的能力很差,这也是目前应用预测控制技术急需解决的重要问题。
发明内容
(一)要解决的技术问题
针对上述现有技术存在的问题,本发明的主要目的在于提供一种基于噪声估计的自适应状态反馈预测控制方法,以充分考虑被控过程中含有的噪声信息和状态信息,提高控制方法抑制噪声的能力。
(二)技术方案
为达到上述目的,本发明的技术方案是这样实现的:
一种基于噪声估计的自适应状态反馈预测控制方法,该方法包括:
a、利用Singular Pencil Model模型表示实际被控过程;
假设xk是状态变量,uk是输入变量,yk是输出变量,ek为噪声变量,E*,A*,B*,C*,E0为矩阵,I为单位阵,z为前向移位算子,则建立公式
假设 则将公式 变形为公式
其中,I是单位矩阵;wk是由输入输出数据组成的矩阵;r是待估计的***模型参数向量;η是待估计的噪声模型参数向量;
b、采用递推的扩展卡尔曼滤波参数估计方法估计所述利用SingularPencil Model模型表示的实际被控过程的状态和模型参数;
所述递推的扩展卡尔曼滤波参数估计算法公式为
式中,
公式 用于修正Kk,将修正后的Kk代入公式 计算得到包含所述实际被控过程状态和模型参数的参数
c、利用估计出来的状态和模型参数计算当前的控制作用,并将计算出来的控制作用加到实际被控过程中;其中所述利用估计出来的状态和模型参数计算当前的控制作用,具体包括:
c1、利用包含实际被控过程状态和模型参数的参数 构建矩阵A、矩阵B、矩阵V和状态反馈向量x;
c2、利用构建的矩阵A、B和V构建状态反馈预测控制方法公式
c3、根据参考轨迹ys(k+p)=αy(k)+(1-α)ω、性能指标 和预测公式 计算增量形式的最优控制率;
Δu(k)=(STQS+λ)-1STQ(ys(k+p)-y(k)-ΔKX-ΔSU-ΔTY),
其中,ω是输入设定值,α为滤波系数,
c4、将Δu(k)=(STQS+λ)-1STQ(ys(k+p)-y(k)-ΔKX-ΔSU-ΔTY)代入公式u(k)=u(k-1)+Δu(k),计算得到当前的控制作用u(k);
d、在下一个执行周期,重复执行步骤b和c。
所述被控过程用二阶单入单出模型表示,rT=[a2 a1 b2 b1],ηT=[c2 c1],xT=[x2 x1],所述步骤c1包括:用rT中的ai构造矩阵 用rT中的bi构造矩阵 用ηT中的ci构造矩阵 用状态xT作状态反馈。
步骤c中所述将计算出来的控制作用加到实际被控过程中包括:将计算出来的控制作用通过数模转化,由数字量转化成模拟量,再送到执行元件中。
(三)有益效果
从上述技术方案可以看出,本发明具有以下有益效果:
1、利用本发明,通过将扩展卡尔曼滤波参数估计方法和状态反馈预测控制方法相结合,可实现同时估计被控过程中含有的状态信息和模型参数,并且估计出来的状态信息可直接用于状态反馈,估计出来的模型参数中的噪声信息可直接用于控制器设计,相当于在控制作用中引入了噪声的校正项,充分考虑了被控过程中含有的噪声信息和状态信息,使该控制方法和其它控制方法相比具有较强的抑制噪声能力,大大提高了控制方法抑制噪声的能力。
2、利用本发明,由于采用的是单值状态反馈预测控制算法,所以本发明提供的控制方法计算量小,便于工业应用。
3、利用本发明,由于采用了SPM(Singular Pencil Model)模型将状态信息和参数信息进行分离,因此可以非常方便地对状态信息和参数信息同时进行估计。
4、利用本发明,由于可以直接利用在线估计出来的状态信息做状态反馈,所以无需再设计状态观测器,使控制器的设计更简单,简化了实现过程,降低了实现成本。
5、利用本发明,由于在每个执行周期内都在线估计被控过程的状态信息和参数信息,动态更新控制器的参数,所以控制器能适应被控过程的变化。
附图说明
图1为本发明提供的自适应状态反馈预测控制总体技术方案的实现流程图;
图2为本发明提供的基于噪声估计的自适应状态反馈预测控制结构示意图;
图3为现有技术中广义预测控制(GPC)方法和本发明提供的控制方法对阶跃扰动的抑制能力仿真对比示意图;
图4为本发明实施例中燃料电池流量***控制Simulink模块示意图;
图5为本发明实施例中燃料电池流量***方法流程图;
图6为本发明实施例中燃料电池流量***的实际控制效果示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明的核心思想是:通过将扩展卡尔曼滤波参数估计方法和状态反馈预测控制方法相结合,可实现同时估计被控过程中含有的状态信息和模型参数,并且估计出来的状态信息可直接用于状态反馈,估计出来的模型参数中的噪声信息可直接用于控制器设计,相当于在控制作用中引入了噪声的校正项,充分考虑了被控过程中含有的噪声信息和状态信息,使该控制方法和其它控制方法相比具有较强的抑制噪声能力,大大提高了控制方法抑制噪声的能力。
上述扩展卡尔曼滤波参数估计方法和状态反馈预测控制方法都是目前应用预测控制技术领域广泛应用的控制方法,下面对其原理分别进行详细阐述。
卡尔曼滤波器是由美国学者卡尔曼(Rudolph Kalman)开发设计,最初是专为航天器导航而研发的,后逐渐应用到许多领域,包括各种核电站设备、人口统计建模、制造业、地层放射性探测以及模糊逻辑和神经网络学习等。卡尔曼滤波器主要用来预估那些只能被***本身间接或不精确观测的***状态。普通卡尔曼滤波器用于估计线性***的状态,对于非线性***的状态,则应使用扩展卡尔曼滤波器进行估计。事实上,在工程中几乎所有的进程都是非线性的。一些非线性***非常接近线性***,而另一些则差别甚大。因而,扩展卡尔曼滤波器更有应用价值。扩展卡尔曼滤波器只是简单地将线性卡尔曼滤波原理推广到非线性***。由于本发明方法中涉及噪声信息的估计,而噪声信息需要用输入输出信息来估计,这是一个非线性问题,因而应用扩展卡尔曼滤波进行参数估计。
状态反馈预测控制方法是将状态空间模型用于预测控制中的一种算法,由于引入了状态信息,而状态中又包含了不可测干扰的影响,相当于引入了未知干扰的前馈,大大提高了控制方法的抗干扰效果。单值的状态反馈预测控制算法还具有计算量小的优点,且和多值的算法控制效果基本相当,便于工业应用。但由于需要状态信息,往往由于状态的不可测而要另外构造状态观测器。本发明方法在参数估计出来的同时就可获得状态信息,使状态反馈预测控制方法更有工业应用价值。
本发明在将扩展卡尔曼滤波参数估计方法和状态反馈预测控制方法结合在一起的过程中,需要在状态反馈预测控制方法中考虑噪声信息,推导出带有噪声信息的状态反馈预测控制方法公式,并将在线估计出来的状态信息和参数信息用于控制方法中,形成自适应的状态反馈预测控制方法。
如图1所示,图1为本发明提供的自适应状态反馈预测控制总体技术方案的实现流程图,该方法包括以下步骤:
步骤101:利用Singular Pencil Model模型表示实际被控过程;
步骤102:采用递推的扩展卡尔曼滤波参数估计方法估计所述利用Singular Pencil Model模型表示的实际被控过程的状态和模型参数;
步骤103:利用估计出来的状态和模型参数计算当前的控制作用,并将计算出来的控制作用加到实际被控过程中;
步骤104:在下一个执行周期,重复执行步骤102和步骤103,实现对实际被控过程的自适应状态反馈预测控制。
上述步骤101中所述Singular Pencil Model模型是一种特殊的模型结构,可以用来表示线性动态***,其模型结构如公式(1)所示:
其中,是状态变量,也叫内部变量或辅助变量。是外部变量,包括***的输入、输出变量。E、F是(n+p)×n矩阵,G是(n+p)×(p+m)矩阵。D是线性算子。对于连续***,D可以是Laplace算子s,对于离散***,D可以是前向移位算子z。
以状态空间模型表示的多入多出***可以用Singular Pencil Model模型表示为公式(2):
其中,uk是输入变量,yk是输出变量,(A,B,C,D)是状态空间矩阵。滑动自回归平均模型 也可以用Singular PencilModel模型表示为公式(3):
其中,Ai和Bi是矩阵,z是前向移位算子。
由于Singular Pencil Model模型具有如下的特征,因而得到了广泛应用。
1)、统一:Singular Pencil Model模型可包含线性控制***中最常用的状态空间模型和输入输出模型。
2)、简单:可通过对公式(1)中矩阵E,F,G做一系列简单的初等行和列变换完成线性变换。
3)、通用:和状态空间描述不同的是,Singular Pencil Model模型还可以描述非正***,如PID控制器、严格正***的逆***等等。
利用Singular Pencil Model模型,应用扩展卡尔曼滤波方法可同时估计出被控过程的状态信息和参数信息。
上述步骤101中所述利用Singular Pencil Model模型表示实际被控过程包括:假设xk是状态变量,uk是输入变量,yk是输出变量,ek为噪声变量,E*,A*,B*,C*,E0为矩阵,I为单位阵,z为前向移位算子,则建立公式(4)表示实际被控过程,所述公式(4)如下:
上述在建立公式(4)表示实际被控过程之后进一步包括:假设 其中矩阵和根据r中参数的特定位置构建,矩阵根据η中参数的特定位置构建,则将用于表示实际被控过程的公式(4)变形为公式(5),所述公式(5)如下:
假设 则将公式(5)变形为公式(6),所述公式(6)如下:
由于ek不能直接测量,可用公式(7)来估计ek,所述公式(7)如下:
将上述公式(7)代入公式(6)可以估计出包含状态信息和参数信息的参数 这样,应用扩展卡尔曼滤波器可以估计sk,也就是说,在噪声不可测的情况下,状态信息和参数信息可以同时估计出来。
上述步骤102中所述采用递推的扩展卡尔曼滤波参数估计方法估计所述利用Singular Pencil Model模型表示的实际被控过程的状态信息和模型参数信息包括:
所述递推的扩展卡尔曼滤波参数估计算法公式为公式(8),所述公式(8)如下:
式中,
公式 用于修正Kk,将修正后的Kk代入公式 计算得到包含所述实际被控过程状态和模型参数的参数
另外,步骤102中所述估计实际被控过程的状态信息和模型参数信息中的模型参数信息包括***矩阵参数和噪声参数。
上述步骤103中所述利用估计出来的状态和模型参数计算当前的控制作用包括:
步骤1031:利用包含实际被控过程状态和模型参数的参数 构建A、B、V和状态反馈x;
这里,以二阶单入单出模型为例,即假设所述被控过程用二阶单入单出模型表示,在二阶单入单出模型中,rT=[a2 a1 b2 b1],ηT=[c2 c1],xT=[x2 x1],则构建A、B、V和状态反馈x包括:
用rT中的ai构造矩阵 用rT中的bi构造矩阵 用ηT中的ci构造矩阵 用状态xT作状态反馈。
步骤1032:利用构建的A、B、V构建状态反馈预测控制方法公式(9),所述公式(9)如下:
步骤1033:根据参考轨迹ys(k+p)=αy(k)+(1-α)ω、性能指标 和预测公式(9)计算增量形式的控制率公式(10),所述公式(10)如下:
Δu(k)=(STQS+λ)-1STQ(ys(k+p)-y(k)-ΔKX-ΔSU-ΔTY) (10)
其中,ω是输入设定值,α为滤波系数, Q为输出预测误差加权阵,λ为控制量的加权系数,Δu(k)为控制作用u(k)的增量;
在本步骤中,由于采用了单值状态反馈预测控制算法,该算法预测值要在一步内调整到给定值,势必造成很大的控制作用,可采用两种手段减少控制作用,一是采用公式ys(k+p)=αy(k)+(1-α)ω所述的参考轨迹,一是在优化性能指标中考虑控制作用的影响,即采用公式 所述的性能指标;其中,ω是输入设定值,α为滤波系数,α值越大,***的柔性越好,鲁棒性越强,但控制的快速性却变差,α对闭环***的动态特性和鲁棒性起着重要作用。这里,所述的鲁棒性用于表示***的抗干扰能力,鲁棒性越强,***的抗干扰能力越强。
步骤1034:将Δu(k)=(STQS+λ)-1STQ(ys(k+p)-y(k)-ΔKX-ΔSU-ΔTY)代入公式u(k)=u(k-1)+Δu(k),计算得到当前的控制作用u(k)。
上述步骤103中所述将计算出来的控制作用加到实际被控过程中包括:将计算出来的控制作用通过数模转化,由数字量转化成模拟量,如电压信号,再送到执行元件中。
在将计算出来的控制作用加到实际被控过程后,可以通过控制阀的开度大小调节实际被控过程,使被控量跟随给定值的变化。
另外,由于本发明对于每一个执行周期,都需要采用递推的扩展卡尔曼滤波参数估计方法估计所述利用Singular Pencil Model模型表示的实际被控过程的状态和模型参数,并利用估计出来的状态和模型参数计算当前的控制作用,然后将计算出来的控制作用加到实际被控过程中。所以,在下一个执行周期,需要重复执行步骤102和步骤103,这样就可以实现对实际被控过程的自适应状态反馈预测控制。
基于图1所述的本发明提供的自适应状态反馈预测控制总体技术方案的方法流程图,图2示出了本发明提供的基于噪声估计的自适应状态反馈预测控制结构示意图。
在图2中,各部分模块的功能可通过下述步骤体现:
步骤201:用Singular Pencil Model模型描述被控过程。
步骤205:给定值ω经过参考模型后得到参考轨迹ys(k+j)。
步骤207:最后根据参考轨迹ys(k+j)和预测模型的校正输出yo(k+j)的差值,计算最优控制率作用于被控过程。
步骤208:在下一个周期,重复上述步骤202至步骤207,实现对实际被控过程的自适应状态反馈预测控制。
基于图1所述的方法流程图和图2所述的结构示意图,图3示出了现有技术中GPC方法和本发明提供的控制方法对阶跃扰动的抑制能力仿真对比示意图。
在图3中,假设阶跃扰动的幅值为5,采样时间Ts为0.1秒,两种方法分别是GPC和基于噪声估计的自适应状态反馈预测控制。为衡量闭环***的控制性能,采用误差绝对值之和SAE(Sum of Absolute Value of Error)作为衡量标准,用公式表示为:
经计算,两种控制方法的SAE分别为53.5010、15.3621,由此可以看出基于噪声估计的自适应状态反馈预测控制方法抑制扰动的能力很强。
为使本发明的目的、技术方案和优点更加清楚明白,下面以质子交换膜燃料电池流量***控制为具体实施例,并参照附图,对本发明进一步详细说明。
燃料电池是一种能将化学能转换成电能的装置。质子交换膜燃料电池由于具有温度低、效率高和无污染的特点而受到广泛关注。它以氢气和氧气,或氢气和空气为原料,反应副产物是水。在氢氧燃料电池中,氢和氧在各自的电极反应。氢气在阳极催化剂的作用下电离成电子和氢离子,发生还原反应,负电子经外电路到达导电的阴极,同时余下的氢离子通过质子渗透膜被送到阴极。在阴极,氢离子与氧气发生氧化反应并从阴极吸收电子生成水。其化学反应方程式为:
阴极:1/2O2+2H++2e-→H2O
阳极:H2→2H++2e-
总反应:1/2O2+H2→H2O
在本实施例中,以中国科学院自动化研究所1KW质子交换膜燃料电池流量***为被控对象,流量***的控制目标是通过自动和快速地调节阴极和阳极进气流量使电堆的压力保持在合适的数值,同时使空气和氢气的进气流量快速地跟踪***负载的变化。
如图4所示,图4为本发明实施例中燃料电池流量***控制Simulink模块示意图。整套***包括的硬件主要有燃料电池发电装置,网络化控制器(NetCon)***,负载电动摩托车和一台上位机。这里,Simulink是Matlab软件中的一个工具箱,此外,Matlab软件中还包括实时工具箱Real-Time Workshop。
Netcon***是一个基于32位ARM核微处理器的网络化控制器。其上可运行uclinux操作***。整个***包括可视化的控制组态软件和可视化的监视组态软件,运行于上位机,上位机与控制器之间通过以太网进行通讯。控制组态软件基于Matlab、Simulink和Real-Time Workshop,能够实现和Simulink的无缝结合。
控制算法可在上位机的Matlab/Simulink环境编写,利用Simulink提供的Real-Time Workshop工具和交叉编译器可将Simulink模块图生成可以在Netcon控制器上直接运行的代码,然后通过以太网下载到高级精简指令计算机处理器(Advanced RISC Machines,ARM)上执行。核心算法通过C-Sfunction实现,具体的实现步骤如图5所示,包括以下步骤:
步骤501:利用Singular Pencil Model模型表示实际被控过程;
具体过程如下:置初值,包括将P0、协方差阵R0、预测时域P,模型阶次n、参考轨迹滤波***α、输出预测误差加权阵Q、控制量的加权系数λ和确定待估计的参数等置初值。
步骤505:计算u(k)=u(k-1)+Δu(k),加到被控过程中。
步骤506:在下一个执行周期,重复执行步骤502至步骤505,实现对实际被控过程的自适应状态反馈预测控制。
基于图4所示的燃料电池流量***控制Simulink模块示意图和图5所示的燃料电池流量***控制方法流程图,图6示出了本发明实施例中燃料电池流量***的实际控制效果示意图。
在图6中,假设预测时域取6,模型阶次n为2,参考轨迹滤波***α为0.5,输出预测误差加权阵Q为1,控制量的加权系数λ为500,采样时间取0.2秒。其控制策略是对阴极压力采用基于噪声估计的自适应状态反馈预测控制,阳极跟随阴极压力的变化采用简单的PI控制,PID参数为比例0.2,积分为0.001,阴极压力的给定值为25KPa。整个控制策略用Matlab的Simulink模块图实现。
从图6所示的实际控制效果可以看出,本发明提供的这种基于噪声估计的自适应状态反馈预测控制方法非常实用。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1、一种基于噪声估计的自适应状态反馈预测控制方法,其特征在于,该方法包括:
a、利用Singular Pencil Model模型表示实际被控过程;
假设 则将公式 变形为公式
其中,I是单位矩阵;wk是由输入输出数据组成的矩阵;r是待估计的***模型参数向量;η是待估计的噪声模型参数向量;
b、采用递推的扩展卡尔曼滤波参数估计方法估计所述利用SingularPencil Model模型表示的实际被控过程的状态和模型参数;
所述递推的扩展卡尔曼滤波参数估计算法公式为
式中,
公式 用于修正Kk,将修正后的Kk代入公式 计算得到包含所述实际被控过程状态和模型参数的参数
c、利用估计出来的状态和模型参数计算当前的控制作用,并将计算出来的控制作用加到实际被控过程中;其中所述利用估计出来的状态和模型参数计算当前的控制作用,具体包括:
c1、利用包含实际被控过程状态和模型参数的参数 构建矩阵A、矩阵B、矩阵V和状态反馈向量x;
c2、利用构建的矩阵A、B和V构建状态反馈预测控制方法公式
c3、根据参考轨迹ys(k+p)=αy(k)+(1-α(ω、性能指标 和预测公式 计算增量形式的最优控制率;
Δu(k)=(ST QS+λ)-1ST Q(ys(k+p)-y(k)-ΔKX-ΔSU-ΔTY),
其中,ω是输入设定值,α为滤波系数, Q为输出预测误差加权阵,λ为控制量的加权系数,Δu(k)为控制作用u(k)的增量;
c4、将Δu(k)=(ST QS+λ)-1ST Q(ys(k+p)-y(k)-ΔKX-ΔSU-ΔTY)代入公式u(k)=u(k-1)+Δu(k),计算得到当前的控制作用u(k);
d、在下一个执行周期,重复执行步骤b和c。
2、根据权利要求1所述的方法,其特征在于,所述被控过程用二阶单入单出模型表示,rT=[a2 a1 b2 b1],ηT=[c2 c1],xT=[x2 x1],所述步骤c1包括:
用rT中的ai构造矩阵 用rT中的bi构造矩阵 用ηT中的ci构造矩阵 用状态xT作状态反馈。
3、根据权利要求1所述的方法,其特征在于,步骤c中所述将计算出来的控制作用加到实际被控过程中包括:
将计算出来的控制作用通过数模转化,由数字量转化成模拟量,再送到执行元件中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100833311A CN100483276C (zh) | 2006-06-02 | 2006-06-02 | 一种基于噪声估计的自适应状态反馈预测控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100833311A CN100483276C (zh) | 2006-06-02 | 2006-06-02 | 一种基于噪声估计的自适应状态反馈预测控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101082811A CN101082811A (zh) | 2007-12-05 |
CN100483276C true CN100483276C (zh) | 2009-04-29 |
Family
ID=38912409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100833311A Expired - Fee Related CN100483276C (zh) | 2006-06-02 | 2006-06-02 | 一种基于噪声估计的自适应状态反馈预测控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100483276C (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2169483A1 (en) * | 2008-09-23 | 2010-03-31 | ABB Research Ltd | Method and system for controlling an industrial process |
CN101813916B (zh) * | 2009-02-19 | 2011-11-30 | 中国石油化工股份有限公司 | 非线性生产过程的自适应预测函数控制方法 |
US9400491B2 (en) * | 2013-03-15 | 2016-07-26 | Rockwell Automation Technologies, Inc. | Stabilized deteministic optimization based control system and method |
CN103247035B (zh) * | 2013-05-20 | 2017-07-11 | 重庆邮电大学 | 基于数字x光机的医学图像处理装置、方法及*** |
JP6444851B2 (ja) * | 2015-12-24 | 2018-12-26 | ファナック株式会社 | ノイズの発生原因を検出する学習機能を有する制御装置 |
CN108536016B (zh) * | 2018-04-19 | 2021-11-05 | 北京联合大学 | 一种基于模糊逆模型的网络化控制方法 |
CN108681244B (zh) * | 2018-05-11 | 2021-01-26 | 东南大学 | 基于多模型前馈的燃料电池阳极压力动态矩阵控制方法 |
CN109814387B (zh) * | 2019-01-31 | 2020-09-04 | 中国地质大学(武汉) | 一种基于动态时延的卡尔曼滤波补偿方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4320287A (en) * | 1980-01-25 | 1982-03-16 | Lockheed Electronics Co., Inc. | Target vehicle tracking apparatus |
CN1202240A (zh) * | 1995-11-15 | 1998-12-16 | 安蒂·阿米·依玛里·兰格 | 一种动态***中自适应卡尔曼滤波方法 |
CN1621981A (zh) * | 2004-12-16 | 2005-06-01 | 上海交通大学 | 基于有限元和***辨识的压电智能结构闭环***仿真方法 |
-
2006
- 2006-06-02 CN CNB2006100833311A patent/CN100483276C/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4320287A (en) * | 1980-01-25 | 1982-03-16 | Lockheed Electronics Co., Inc. | Target vehicle tracking apparatus |
CN1202240A (zh) * | 1995-11-15 | 1998-12-16 | 安蒂·阿米·依玛里·兰格 | 一种动态***中自适应卡尔曼滤波方法 |
CN1621981A (zh) * | 2004-12-16 | 2005-06-01 | 上海交通大学 | 基于有限元和***辨识的压电智能结构闭环***仿真方法 |
Non-Patent Citations (2)
Title |
---|
关于奇异束的Leverrier-Hermite算法. 王国荣.上海师范大学学报(自然科学版),第30卷第3期. 2001 * |
基于噪声方差估计的稳健融合算法. 胡士强,敬忠良.仪器仪表学报,第25卷第4增刊期. 2004 * |
Also Published As
Publication number | Publication date |
---|---|
CN101082811A (zh) | 2007-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100483276C (zh) | 一种基于噪声估计的自适应状态反馈预测控制方法 | |
Harrag et al. | Novel neural network IC-based variable step size fuel cell MPPT controller: performance, efficiency and lifetime improvement | |
Wu et al. | Predictive control of SOFC based on a GA-RBF neural network model | |
Fang et al. | Fault diagnosis of networked control systems | |
Wang et al. | Data-driven predictive control for solid oxide fuel cells | |
Damour et al. | Real-time implementation of a neural model-based self-tuning PID strategy for oxygen stoichiometry control in PEM fuel cell | |
CN107870564B (zh) | 具有解耦性能的燃料电池抗干扰控制方法 | |
CN112349936A (zh) | 基于***辨识的船用燃料电池建模及广义预测控制方法 | |
Barragán et al. | Iterative fuzzy modeling of hydrogen fuel cells by the extended Kalman filter | |
Wang et al. | Observer-based adaptive neural network control for PEMFC air-feed subsystem | |
Belmokhtar et al. | Dynamic model of an alkaline electrolyzer based an artificial neural networks | |
Han et al. | Multi-objective model predictive control with gradient eigenvector algorithm | |
CN111413938B (zh) | 一种基于折算喷氨量的scr脱硝***扰动抑制预测控制方法 | |
CN111707945A (zh) | 一种质子交换膜燃料电池内部状态观测方法 | |
CN113690874B (zh) | 一种多区域电力***的无模型自适应负荷频率控制方法 | |
Tong et al. | Hydrogen-air PEM fuel cell: integration, modeling, and control | |
Tong et al. | Real-time simplified variable domain fuzzy control of PEM fuel cell flow systems | |
CN114325398A (zh) | 一种针对质子交换膜燃料电池***的故障检测方法 | |
Li et al. | Observer-Based Control of PEMFC Air Feed Systems Via High-Order Fully Actuated System Approaches | |
Guo et al. | An improved generalized predictive control algorithm based on PID | |
CN116627042B (zh) | 执行器非对称饱和多自主体***的分布式协同跟踪方法 | |
Halledj et al. | Intelligent control based on single-input interval type-2 fuzzy iPI and super twisting nonlinear sliding mode compensator for regulating PEMFC output voltage | |
Alipouri et al. | Optimal controller design with communication delay for solid oxide fuel cell | |
Fan et al. | Improved model predictive control for a proton exchange membrane fuel cell | |
Chen et al. | Intelligent Prediction of Power Grid Project Cost Based on PSO-BP Model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090429 Termination date: 20170602 |
|
CF01 | Termination of patent right due to non-payment of annual fee |