CN103336437B - 基于预测函数控制优化的非自衡对象控制方法 - Google Patents
基于预测函数控制优化的非自衡对象控制方法 Download PDFInfo
- Publication number
- CN103336437B CN103336437B CN201310310130.0A CN201310310130A CN103336437B CN 103336437 B CN103336437 B CN 103336437B CN 201310310130 A CN201310310130 A CN 201310310130A CN 103336437 B CN103336437 B CN 103336437B
- Authority
- CN
- China
- Prior art keywords
- mav
- model
- controlled device
- moment
- control
- 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.)
- Active
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于预测函数控制优化的非自衡对象控制方法。本发明方法首先基于非自衡对象的阶跃响应数据建立被控对象的模型,挖掘出基本的对象特性;然后依据预测函数控制的特性去整定相应P控制器的参数;最后对被控对象实施P控制。本发明将预测函数控制的性能赋给了P控制,有效地提高了传统控制方法的性能,同时也促进了先进控制方法的应用。
Description
技术领域
本发明属于自动化技术领域,涉及一种基于预测函数控制优化的非自衡对象比例(P)控制方法。
背景技术
在实际工业对象的控制中,由于硬件、成本、实施难度等方面的限制,一些先进控制方法虽然得到了一定程度的应用,但是目前占主流的仍是PID控制,一些非自衡对象,例如储罐、锅炉汽包水位、精馏塔液位等等,通常采用P控制。预测函数控制作为先进控制方法的一种,在非自衡对象的控制中相比P控制拥有更好的控制性能,如果能将预测函数控制的性能赋给P控制,那将进一步推进先进控制方法的应用,同时也能获得更好的实际控制效果。
发明内容
本发明的目的是针对现有先进控制方法的应用不足之处,提供一种基于预测函数控制优化的非自衡对象P控制方法,以获得更好的控制性能。该方法通过结合预测函数控制和P控制,得到了一种带有预测函数控制性能的P控制方法。该方法在继承预测函数控制优良性能的同时也保证形式简单并满足实际工业过程的需要。
本发明方法首先基于非自衡对象的阶跃响应数据建立被控对象的模型,挖掘出基本的对象特性;然后依据预测函数控制的特性去整定相应P控制器的参数;最后对被控对象实施P控制。
本发明的技术方案是通过数据采集、模型建立、预测机理、优化等手段,确立了一种基于预测函数控制的P控制方法,利用该方法可有效提高控制的精度与稳定性。
本发明方法的步骤包括:
步骤(1).通过被控对象的实时阶跃响应数据建立被控对象的模型,具体方法是:
Ⅰ.给被控对象一个阶跃输入信号,记录被控对象的阶跃响应曲线。
Ⅱ.将对应的阶跃响应曲线滤波处理成一条光滑曲线,然后将数值发生变化的曲线段拟合成一条直线,计算出模型的增益参数Km:
Km=q*l
其中,q为过程模型的控制量的阶跃变化幅度,Km是建立的被控对象模型的增益系数,l是拟合的直线的斜率。
Ⅲ.记录Ⅱ中滤波处理后的光滑曲线上每个采样时刻对应的阶跃响应数据,第一个采样时刻为Ts1,相邻两个采样时刻间隔的时间为Ts1,采样时刻顺序为Ts1、2Ts1、3Ts1……;在记录的阶跃响应数据中找出数据开始上升的起始点aI,之前的数据分别记做a1,a2,…aI-1,模型的滞后时间参数τ为τ=(I-1)Ts1,最后得到的被控对象的传递函数模型为:
其中,G(s)为被控对象的传递函数,s为拉普拉斯变换算子。
步骤(2).设计被控对象的P控制器,具体方法是:
a.对得到的传递函数在采样时间Ts下加一个零阶保持器离散化,得到离散模型为
ym(k)=ym(k-1)+KmTsu(k-1-L)
ym(k)为k时刻的模型预测输出,u(k-1-L)为k-1-L时刻的控制输入,L为离散传递函数模型的时滞,L=τ/Ts。
b.计算被控对象去掉纯滞后以后在预测函数控制下的第P步预测输出,形式如下:
ymav(k)=ymav(k-1)+KmTsu(k-1)
ymav(k+P)=ymav(k)+KmPTsu(k)
其中,P为预测步长,ymav(k+P)为k时刻去掉纯滞后的被控对象在预测函数控制下的第P步预测输出,ymav(k)为k时刻去掉纯滞后的模型输出。
c.校正k时刻的实际输出得到包含未来预测信息的新的过程输出值,形式如下:
ypav(k)=yp(k)+ymav(k)-ymav(k-L)
其中,ypav(k)为校正得到的k时刻包含未来预测信息的新的过程输出值,yp(k)为k时刻的实际输出值。
d.选取预测函数控制方法的参考轨迹yr(k+P)以及目标函数J,形式如下:
yr(k+P)=βPyp(k)+(1-βP)c(k)
J=min(yr(k+P)-ymav(k+P)-e(k))2
e(k)=ypav(k)-ymav(k)+(P+L)(ypav(k)-ymav(k)-ypav(k-1)+ymav(k-1))
其中,β为参考轨迹柔化系数,c(k)为k时刻的设定值,e(k)为k时刻校正的误差值,这里采用新的误差校正方式。
e.依据步骤d中的目标函数求解P控制器中的参数,这里将控制量u(k)进行变换:
u(k)=Kpe1(k)
e1(k)=βyp(k-1)+(1-β)c(k-1)-yp(k)
其中,Kp为位置式P控制器的参数,e1(k)为k时刻参考轨迹值与实际输出值之间的误差。结合上述式子,可以求得:
f.得到P控制器的参数Kp以后构成控制量u(k)=Kpe1(k)作用于被控对象。
g.在下一时刻,依照b到f中的步骤继续求解新的Kp,依次循环。
本发明提出的一种基于预测函数控制优化的非自衡对象比例控制方法将预测函数控制的性能赋给了P控制,有效地提高了传统控制方法的性能,同时也促进了先进控制方法的应用。
具体实施方式
以锅炉汽包水位控制为例:
锅炉汽包水位是典型的带积分环节的非自衡对象,汽包水位也是表征锅炉安全运行的重要参数之一。
步骤(1).获得锅炉汽包水位的模型。
第一步:在手动模式下,在蒸汽负荷不变的情况下通过调节给水调节阀获得汽包水位的阶跃响应曲线,将对应的阶跃响应曲线滤波处理成一条光滑曲线,然后将数值发生变化的曲线段拟合成一条直线,得到直线的斜率l。
第二步:记录第一步中滤波处理后的光滑曲线上每个采样时刻对应的阶跃响应数据,第一个采样时刻为Ts1,相邻两个采样时刻间隔的时间为Ts1,采样时刻顺序为Ts1、2Ts1、3Ts1……;在记录的阶跃响应数据中找出数据开始上升的起始点aI,之前的数据分别记做a1,a2,…aI-1,模型的滞后时间参数τ为τ=(I-1)Ts1。
第三步:依据第一步得到的斜率l计算出模型的增益系数Km=q*l,这样便得到了对象的传递函数模型:
其中,G(s)为锅炉汽包水位模型的传递函数,s为拉普拉斯变换算子。
步骤(2).设计被控对象的P控制器,具体方法是:
a.对得到的传递函数在采样时间Ts下加一个零阶保持器离散化,得到锅炉汽包水位的离散模型为
ym(k)=ym(k-1)+KmTsu(k-1-L)
ym(k)为k时刻的锅炉汽包水位预测输出,u(k-1-L)为k-1-L时刻的给水阀阀门开度的输入控制量,L为离散传递函数模型的时滞,L=τ/Ts。
b.计算被控对象去掉纯滞后以后在预测函数控制下的第P步预测输出,形式如下:
ymav(k)=ymav(k-1)+KmTsu(k-1)
ymav(k+P)=ymav(k)+KmPTsu(k)
其中,P为预测步长,ymav(k+P)为k时刻去掉纯滞后的锅炉汽包水位模型在预测函数控制下的第P步预测输出,ymav(k)为k时刻去掉纯滞后的锅炉汽包水位模型输出。
c.校正k时刻的锅炉汽包水位实际输出得到包含未来预测信息的新的过程输出值,形式如下:
ypav(k)=yp(k)+ymav(k)-ymav(k-L)
其中,ypav(k)为校正得到的k时刻包含未来预测信息的的锅炉汽包水位过程输出值,yp(k)为k时刻的锅炉汽包水位实际输出值。
d.选取预测函数控制方法的参考轨迹yr(k+P)以及目标函数J,形式如下:
yr(k+P)=βPyp(k)+(1-βP)c(k)
J=min(yr(k+P)-ymav(k+P)-e(k))2
e(k)=ypav(k)-ymav(k)+(P+L)(ypav(k)-ymav(k)-ypav(k-1)+ymav(k-1))
其中,β为参考轨迹柔化系数,c(k)为k时刻锅炉汽包水位的设定值,e(k)为k时刻校正的误差值,这里采用新的误差校正方式。
e.依据步骤d中的目标函数求解P控制器中的参数,这里将给水阀阀门开度输入控制量u(k)进行变换:
u(k)=Kpe1(k)
e1(k)=βyp(k-1)+(1-β)c(k-1)-yp(k)
其中,Kp为位置式P控制器的参数,e1(k)为k时刻锅炉汽包水位参考轨迹值与实际输出值之间的误差。
结合上述式子,可以求得:
f.得到P控制器的参数Kp以后构成控制量u(k)=Kpe1(k)作用于锅炉汽包水位的给水阀阀门。
g.在下一时刻,依照b到f中的步骤继续求解锅炉汽包水位过程P控制器新的参数Kp,依次循环。
Claims (1)
1.基于预测函数控制优化的非自衡对象控制方法,其特征在于该方法的具体步骤是:
步骤(1).通过被控对象的实时阶跃响应数据建立被控对象的模型,具体是:
Ⅰ.给被控对象一个阶跃输入信号,记录被控对象的阶跃响应曲线;
Ⅱ.将对应的阶跃响应曲线滤波处理成一条光滑曲线,然后将数值发生变化的曲线段拟合成一条直线,计算出模型的增益参数Km:
Km=q*l
其中,q为过程模型的控制量的阶跃变化幅度,Km是建立的被控对象模型的增益系数,l是拟合的直线的斜率;
Ⅲ.记录Ⅱ中滤波处理后的光滑曲线上每个采样时刻对应的阶跃响应数据,第一个采样时刻为Ts1,相邻两个采样时刻间隔的时间为Ts1,采样时刻顺序为Ts1、2Ts1、3Ts1……;在记录的阶跃响应数据中找出数据开始上升的起始点aI,之前的数据分别记做a1,a2,…aI-1,模型的滞后时间参数τ为τ=(I-1)Ts1,最后得到的被控对象的传递函数模型为:
其中,G(s)为被控对象的传递函数,s为拉普拉斯变换算子;
步骤(2).设计被控对象的P控制器,具体是:
a.对得到的传递函数在采样时间Ts下加一个零阶保持器离散化,得到离散模型为:
ym(k)=ym(k-1)+KmTsu(k-1-L)
ym(k)为k时刻的模型预测输出,u(k-1-L)为k-1-L时刻的控制输入,L为离散传递函数模型的时滞,L=τ/Ts;
b.计算被控对象去掉纯滞后以后在预测函数控制下的第P步预测输出,形式如下:
ymav(k)=ymav(k-1)+KmTsu(k-1)
ymav(k+P)=ymav(k)+KmPTsu(k)
其中,P为预测步长,ymav(k+P)为k时刻去掉纯滞后的被控对象在预测函数控制下的第P步预测输出,ymav(k)为k时刻去掉纯滞后的模型输出;
c.校正k时刻的实际输出得到包含未来预测信息的新的过程输出值,形式如下:
ypav(k)=yp(k)+ymav(k)-ymav(k-L)
其中,ypav(k)为校正得到的k时刻包含未来预测信息的新的过程输出值,yp(k)为k时刻的实际输出值;
d.选取预测函数控制方法的参考轨迹yr(k+P)以及目标函数J,形式如下:
yr(k+P)=βPyp(k)+(1-βP)c(k)
J=min(yr(k+P)-ymav(k+P)-e(k))2
e(k)=ypav(k)-ymav(k)+(P+L)(ypav(k)-ymav(k)-ypav(k-1)+ymav(k-1))
其中,β为参考轨迹柔化系数,c(k)为k时刻的设定值,e(k)为k时刻校正的误差值,这里采用新的误差校正方式;
e.依据步骤d中的目标函数求解P控制器中的参数,这里将控制量u(k)进行变换:
u(k)=Kpe1(k)
e1(k)=βyp(k-1)+(1-β)c(k-1)-yp(k)
其中,Kp为位置式P控制器的参数,e1(k)为k时刻参考轨迹值与实际输出值之间的误差;
结合上述式子,可以求得:
f.得到P控制器的参数Kp以后构成控制量u(k)=Kpe1(k)作用于被控对象;
g.在下一时刻,依照b到f中的步骤继续求解新的Kp,依次循环。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310310130.0A CN103336437B (zh) | 2013-07-19 | 2013-07-19 | 基于预测函数控制优化的非自衡对象控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310310130.0A CN103336437B (zh) | 2013-07-19 | 2013-07-19 | 基于预测函数控制优化的非自衡对象控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103336437A CN103336437A (zh) | 2013-10-02 |
CN103336437B true CN103336437B (zh) | 2015-11-25 |
Family
ID=49244629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310310130.0A Active CN103336437B (zh) | 2013-07-19 | 2013-07-19 | 基于预测函数控制优化的非自衡对象控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103336437B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106338915B (zh) * | 2016-10-31 | 2019-04-26 | 浙江邦业科技股份有限公司 | 一种基于扩展状态空间预测函数控制的积分对象控制方法 |
CN107065541A (zh) * | 2017-03-22 | 2017-08-18 | 杭州电子科技大学 | 一种焦化炉炉膛压力系统模糊网络优化pid‑pfc控制方法 |
CN107728481B (zh) * | 2017-11-14 | 2021-06-04 | 江西理工大学 | 一种基于模型预测控制的闭环建模方法及装置 |
CN109856978B (zh) * | 2019-03-26 | 2022-02-15 | 广东电网有限责任公司 | 一种获取被控对象模型的方法及装置 |
CN112965365B (zh) * | 2021-02-23 | 2023-03-31 | 浙江中智达科技有限公司 | Pid控制回路的模型辨识方法、装置、系统及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329553A (zh) * | 2008-07-04 | 2008-12-24 | 东南大学 | 气化炉的预测型比例积分微分控制方法 |
CN102436178A (zh) * | 2011-11-22 | 2012-05-02 | 浙江大学 | 一种误差容忍限机制下的焦化加热炉氧含量控制方法 |
CN103116283A (zh) * | 2013-01-18 | 2013-05-22 | 杭州电子科技大学 | 一种非自衡对象的动态矩阵控制方法 |
-
2013
- 2013-07-19 CN CN201310310130.0A patent/CN103336437B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329553A (zh) * | 2008-07-04 | 2008-12-24 | 东南大学 | 气化炉的预测型比例积分微分控制方法 |
CN102436178A (zh) * | 2011-11-22 | 2012-05-02 | 浙江大学 | 一种误差容忍限机制下的焦化加热炉氧含量控制方法 |
CN103116283A (zh) * | 2013-01-18 | 2013-05-22 | 杭州电子科技大学 | 一种非自衡对象的动态矩阵控制方法 |
Non-Patent Citations (4)
Title |
---|
Predictive Functional Controller with a Similar Proportional Integral Optimal Regulator Structure: Comparison with Traditional Predictive Functional Controller and Application to Heavy Oil Coking Equipment;Zhang Ridong,Wang Shu qing;《Chinese Journal of Chemical Engineering》;20071231;第15卷(第2期);第247-253页 * |
基于阶跃响应的非自衡对象预测控制;张政江,孙优贤;《控制与决策》;20010531;第16卷(第3期);第378-379页 * |
焦化加热炉出口温度的预测函数控制;苏成利,张日东,张健明,王树青,毛国平;《化工自动化及仪表》;20070412;第34卷(第1期);第4页左栏第4.2节 * |
非自衡对象预测函数控制及其仿真研究;王国玉,韩璞,刘红军,王东风;《华北电力大学学报》;20040731;第31卷(第4期);第39-43页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103336437A (zh) | 2013-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103336437B (zh) | 基于预测函数控制优化的非自衡对象控制方法 | |
CN103389746B (zh) | 预测函数控制优化的废塑料炼油裂解炉炉膛压力控制方法 | |
CN106338915B (zh) | 一种基于扩展状态空间预测函数控制的积分对象控制方法 | |
US6688532B2 (en) | Controller, temperature controller and heat processor using same | |
CN107885259B (zh) | 一种氯化反应温度控制方法、系统及设备 | |
CN103345150B (zh) | 预测函数控制优化的废塑料炼油裂解炉炉膛温度控制方法 | |
CN100545772C (zh) | 一种燃煤锅炉系统混合控制方法 | |
CN105045233B (zh) | 火电厂热工系统中基于时间量度的pid控制器的优化设计方法 | |
CN102621883B (zh) | Pid参数整定方法及pid参数整定系统 | |
CN108489015B (zh) | 基于极点配置和帕德近似的空调系统温度控制方法 | |
CN101709869A (zh) | 燃煤锅炉过热蒸汽温度系统混合控制方法 | |
CN106647832A (zh) | 一种基于CompactLogix的模糊自整定PID的液位控制方法 | |
CN106125546A (zh) | 一种应用被控量的变化趋势的pid控制方法 | |
CN103116283A (zh) | 一种非自衡对象的动态矩阵控制方法 | |
CN105807615A (zh) | 模糊前馈反馈控制器 | |
RU2578297C1 (ru) | Способ и устройство настройки системы автоматического регулирования давления (сард) в магистральном трубопроводе для перекачивания нефтепродуктов | |
CN104528881A (zh) | 基于电动球阀调节的产水控制方法及装置、超纯水制备系统 | |
CN101709867A (zh) | 燃煤锅炉汽包水位系统混合控制方法 | |
DE102010055849A1 (de) | Verfahren zur Regelung der Heißwasserzubereitung | |
CN103760931B (zh) | 动态矩阵控制优化的油气水卧式三相分离器压力控制方法 | |
CN102520618A (zh) | 误差容忍限机制下的焦化加热炉辐射出口温度控制方法 | |
CN103760927A (zh) | 预测函数控制优化的油气水三相分离器油水液面控制方法 | |
CN102436178B (zh) | 一种误差容忍限机制下的焦化加热炉氧含量控制方法 | |
CN108132597B (zh) | 一种微分超前智能模型集pid控制器设计方法 | |
CN102880047B (zh) | 炼油工业加热炉温度过程的伴随矩阵解耦预测控制方法 |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20131002 Assignee: ZHEJIANG BONYEAR TECHNOLOGY Co.,Ltd. Assignor: HANGZHOU DIANZI University Contract record no.: X2020330000039 Denomination of invention: Predictive function control optimization-based integrating plant control method Granted publication date: 20151125 License type: Common License Record date: 20200529 |