CN113359467A - 一种工业过程中基于分数阶模型的改进内模控制方法 - Google Patents
一种工业过程中基于分数阶模型的改进内模控制方法 Download PDFInfo
- Publication number
- CN113359467A CN113359467A CN202110741822.5A CN202110741822A CN113359467A CN 113359467 A CN113359467 A CN 113359467A CN 202110741822 A CN202110741822 A CN 202110741822A CN 113359467 A CN113359467 A CN 113359467A
- Authority
- CN
- China
- Prior art keywords
- model
- generalized
- fractional order
- order
- internal model
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
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)
Abstract
本发明公开了一种工业过程中基于分数阶模型的改进内模控制方法,包括如下步骤:步骤1:建立***的分数阶广义过程模型;步骤2:设计内模前馈控制器。本发明首先建立***的分数阶过程模型,利用Oustaloup近似方法将该模型近似为一个高阶的整数阶模型,然后定义新型的广义输入得到离散***的广义过程高阶模型。接着引入闭环等效的内模控制算法,结合广义过程模型设计内模控制器。最后以加热炉炉温控制为例,说明该算法有效的处理了***的纯滞后和模型失配问题,相对于传统的内模控制来说也提高了模型的精准度、明显的加快了扰动回复速度。
Description
技术领域
本发明属于自动化技术领域,涉及一种工业过程中的数字内模控制技术。
背景技术
在现阶段的工业自动化控制领域中,PID控制因为结构简单,调节方便仍是大多数工业控制的首选算法。但是对于一些模型失配明显、大时滞***来说调节效果一般。后提出的史密斯预估控制和内模控制较好的解决了这些问题。而目前的内模控制器大多采用S域传递函数形式的数学模型,这就需要一个由连续到离散模型的转换。转换过程不可避免地会出现一定程度上的偏差。因此设计一个可以应对大时滞、失配和扰动的数字化控制方法是很有必要的。
发明内容
本发明提出了一种改进的二自由度数字内模控制方案,目的是为了有效地处理常规控制方法难以应对的***过程纯滞后、扰动等不利因素。
该方法首先建立***的分数阶过程模型,利用Oustaloup近似方法将该模型近似为一个高阶的整数阶模型,然后定义新型的广义输入得到离散***的广义过程高阶模型。接着引入闭环等效的内模控制算法,结合广义过程模型设计内模控制器。最后以加热炉炉温控制为例,说明该算法有效的处理了***的纯滞后和模型失配问题,相对于传统的内模控制来说也提高了模型的精准度、明显的加快了扰动回复速度。
本发明方法的步骤包括:
步骤1:建立***的分数阶广义过程模型。
1.1假设被控对象的分数阶模型表示如下:
其中,Gp(s)是被控对象,s是拉普拉斯变换算子,μ为微分项的微分阶次,且有0<μ≤1。K为分数阶模型增益常数,T是分数阶模型惯性时间常数,θ是模型的滞后时间。
1.2为了将其转变为一组有效的离散多项式,用Oustaloup近似方法对其进行近似:
1.3根据步骤1.2,并通过Riemann-Liouville(RL)定理将步骤1.1中的分数阶模型转换为如下的高阶连续模型:
其中,αm>αm-1>…>α1>α0>0,βn>βn-1>…>β1>β0>0分别表示分数阶传递函数分子、分母的微分阶次,且βn>αm。
1.4根据步骤1.3,并考虑到***可能遭遇扰动将线性***过程模型的离散化表达式写成如下形式:
其中,y(z)、u(z)、di(z)分别是过程输出、输入和第i个扰动的z变换;Gp(z)、Gi(z)是过程被控对象和第i个扰动的脉冲传递函数。
1.5将步骤1.4中过程离散模型Gp(z)进行分解:
其中
Gp+(z)=B+(z-1)包含Gp(z)的非最小相位零点和纯滞后;B-(z-1)的所有零点在z平面开环单位圆内,A(z-1)是过程模型Gp(z)极点表达式,ai、βi、bi分别是A(z-1)、B+(z-1)、B-(z-1)的过程系数。
1.6根据步骤1.5,步骤1.4可写成:
1.7定义新型的广义输入:
V(z)=Q(z-1)B-(z-1)u(z)+P(z-1)y(z)
1.8由步骤1.7,步骤1.6可转化为:
该式表示了广义开环***,它的极点是多项式G(z-1)的零点。对其进行极点配置使G(z-1)具有A(z-1)多项式的形式:
1.9使用Diophantine方程对G(z-1)进行求解:
为了获取最小阶解,令nq=k+n+,np=na-1,可以求出P(z-1)、Q(z-1)的最小阶解。1.10由步骤1.8还可以得到广义输出和输入之间的等效模型:
步骤2:设计内模前馈控制器。
2.1对步骤1.10中的等效广义模型进行内模分解:
2.2设计滤波器f(z):
其中ε是可调参数,k为延迟时间系数,T为初始***的时间常数。
2.3由步骤2.3和2.4可得内模控制器Gimc(z):
即
2.4由步骤2.3可求出内模前馈控制器Gc(z):
2.5又因为有
e(z)=r(z)-y(z)
其中e(z)为***误差;r(z)为设定值。结合步骤1.7和步骤2.4可得控制量表达式:
具体实施方式
以工业过程中常见的加热炉温度控制为例对该算法进行进一步说明。
其中被控对象是加热炉炉内的温度,通过继电器的通断比例来控制加热的占空比,最终达到控制升温速度的效果,因此将继电器的占空比作为控制量。
步骤1:建立加热炉***的分数阶过程模型,收缩该模型的开环极点最终将其等效为一个广义过程模型。
1.1首先建立加热炉***的分数阶模型,表示如下:
其中,s是拉普拉斯变换算子,μ为微分项的微分阶次,且0<μ≤1。K为分数阶模型增益常数,T是分数阶模型惯性时间常数,θ是加热炉加热的滞后时间。
1.2为了将加热炉模型转变为一组有效的离散多项式,用Oustaloup近似方法对其进行近似:
1.3根据步骤1.2,并通过Riemann-Liouville(RL)定理将步骤1.1中的分数阶模型转换为如下的高阶连续模型:
其中,αm>αm-1>…>α1>α0>0,βn>βn-1>…>β1>β0>0分别表示分数阶传递函数分子、分母的微分阶次,且βn>αm。
1.4根据步骤1.3,并考虑到***可能遭遇扰动将加热炉过程模型的离散化形式写成如下形式:
其中,y(z)、u(z)、di(z)分别是过程输出、输入和第i个扰动的z变换;Gp(z)、Gi(z)是过程被控对象和第i个扰动的脉冲传递函数。
1.5将步骤1.4中加热炉过程模型Gp(z)进行分解:
其中
Gp+(z)=B+(z-1)包含Gp(z)的非最小相位零点和加热炉纯滞后; 其中B-(z-1)的所有零点在z平面开环单位圆内,A(z-1)是过程模型Gp(z)极点表达式。ai、βi、bi分别是A(z-1)、B+(z-1)、B-(z-1)的过程系数。
1.6结合步骤1.5步骤1.4中的加热炉离散***可转化为:
1.7定义一个新型的电加热炉温度控制***的广义输入:
V(z)=Q(z-1)B-(z-1)u(z)+P(z-1)y(z)
1.8由步骤1.7,步骤1.6可写成:
该式表示了广义加热炉开环***,它的极点是多项式G(z-1)的零点。对其进行极点配置使G(z-1)具有A(z-1)多项式的形式:
1.9使用Diophantine方程对G(z-1)进行求解:
为了获取最小阶解,令nq=k+n+,np=na-1,可以求出P(z-1)、Q(z-1)的最小阶解。1.10由步骤1.8还可以得到广义输出和输入之间的等效模型:
步骤2:设计加热炉控制***的前馈控制器。
2.1对步骤1.10中的等效广义加热炉模型进行分解:
2.2引入滤波器f(z):
其中ε是可调参数,k为延迟时间系数,T为初始***的时间常数。
2.3由步骤2.1和2.2可得内模控制器Gimc(z):
即
2.4由步骤2.3可求得电加热炉炉内温度控制***的继电器通断时间占空比反馈控制器的形式,如下所示Gc(z):
2.5又因为有
e(z)=r(z)-y(z)
其中e(z)为***误差;r(z)为温度设定值;y(z)是传感器测得的炉温。结合步骤1.4和步骤2.4可得加热炉温度控制律:
2.6根据步骤2.1到步骤2.6,依次做循环求解基于加热炉炉温控制的占空比u(z),再将其作用于加热炉***进行温度控制。
Claims (3)
1.一种工业过程中基于分数阶模型的改进内模控制方法,其特征在于包括如下步骤:
步骤1:建立***的分数阶广义过程模型;
步骤2:设计内模前馈控制器。
2.如权利要求1所述的工业过程中基于分数阶模型的改进内模控制方法,其特征在于:所述步骤1包括如下步骤:
1.1假设被控对象的分数阶模型表示如下:
其中,Gp(s)是被控对象,s是拉普拉斯变换算子,μ为微分项的微分阶次,且有0<μ≤1;K为分数阶模型增益常数,T是分数阶模型惯性时间常数,θ是模型的滞后时间;
1.2用Oustaloup近似方法对其进行近似:
1.3根据步骤1.2,并通过Riemann-Liouville(RL)定理将步骤1.1中的分数阶模型转换为如下的高阶连续模型:
其中,αm>αm-1>…>α1>α0>0,βn>βn-1>…>β1>β0>0分别表示分数阶传递函数分子、分母的微分阶次,且βn>αm;
1.4根据步骤1.3,并考虑到***可能遭遇扰动将线性***过程模型的离散化表达式写成如下形式:
其中,y(z)、u(z)、di(z)分别是过程输出、输入和第i个扰动的z变换;Gp(z)、Gi(z)是过程被控对象和第i个扰动的脉冲传递函数;
1.5将步骤1.4中过程离散模型Gp(z)进行分解:
其中
Gp+(z)=B+(z-1)包含Gp(z)的非最小相位零点和纯滞后;B-(z-1)的所有零点在z平面开环单位圆内,A(z-1)是过程模型Gp(z)极点表达式,ai、βi、bi分别是A(z-1)、B+(z-1)、B-(z-1)的过程系数;
1.6根据步骤1.5,步骤1.4可写成:
1.7定义新型的广义输入:
V(z)=Q(z-1)B-(z-1)u(z)+P(z-1)y(z)
1.8由步骤1.7,步骤1.6可转化为:
该式表示了广义开环***,它的极点是多项式G(z-1)的零点;其进行极点配置使G(z-1)具有A(z-1)多项式的形式:
1.9使用Diophantine方程对G(z-1)进行求解:
为了获取最小阶解,令nq=k+n+,np=na-1,可以求出P(z-1)、Q(z-1)的最小阶解;
1.10由步骤1.8还可以得到广义输出和输入之间的等效模型:
3.如权利要求1所述的工业过程中基于分数阶模型的改进内模控制方法,其特征在于:所述步骤2包括如下步骤:
2.1对步骤1.10中的等效广义模型进行内模分解:
2.2设计滤波器f(z):
其中ε是可调参数,k为延迟时间系数,T为初始***的时间常数;
2.3由步骤2.3和2.4可得内模控制器Gimc(z):
即
2.4由步骤2.3可求出内模前馈控制器Gc(z):
2.5又因为有
e(z)=r(z)-y(z)
其中e(z)为***误差;r(z)为设定值;结合步骤1.7和步骤2.4可得控制量表达式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110741822.5A CN113359467B (zh) | 2021-06-30 | 2021-06-30 | 一种工业过程中基于分数阶模型的改进内模控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110741822.5A CN113359467B (zh) | 2021-06-30 | 2021-06-30 | 一种工业过程中基于分数阶模型的改进内模控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113359467A true CN113359467A (zh) | 2021-09-07 |
CN113359467B CN113359467B (zh) | 2023-05-02 |
Family
ID=77537644
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110741822.5A Active CN113359467B (zh) | 2021-06-30 | 2021-06-30 | 一种工业过程中基于分数阶模型的改进内模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113359467B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974798A (zh) * | 2016-07-05 | 2016-09-28 | 杭州电子科技大学 | 电加热炉的多模型分数阶加权预测函数控制方法 |
CN108897223A (zh) * | 2018-08-02 | 2018-11-27 | 杭州电子科技大学 | 一种工业加热炉分数阶预测控制方法 |
CN109001975A (zh) * | 2018-08-02 | 2018-12-14 | 杭州电子科技大学 | 一种工业加热炉多模型分数阶控制方法 |
CN109116738A (zh) * | 2018-09-27 | 2019-01-01 | 杭州电子科技大学 | 一种工业加热炉的二自由度内模控制分析方法 |
CN109557810A (zh) * | 2018-11-29 | 2019-04-02 | 杭州电子科技大学 | 一种基于新型二自由度内模pid的加热炉温度控制方法 |
CN110109360A (zh) * | 2019-05-22 | 2019-08-09 | 杭州电子科技大学 | 一种工业过程广义大琳响应控制方法 |
CN110879576A (zh) * | 2019-10-15 | 2020-03-13 | 杭州电子科技大学 | 水泥熟料冷却过程的分数阶模型广义二自由度控制方法 |
-
2021
- 2021-06-30 CN CN202110741822.5A patent/CN113359467B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974798A (zh) * | 2016-07-05 | 2016-09-28 | 杭州电子科技大学 | 电加热炉的多模型分数阶加权预测函数控制方法 |
CN108897223A (zh) * | 2018-08-02 | 2018-11-27 | 杭州电子科技大学 | 一种工业加热炉分数阶预测控制方法 |
CN109001975A (zh) * | 2018-08-02 | 2018-12-14 | 杭州电子科技大学 | 一种工业加热炉多模型分数阶控制方法 |
CN109116738A (zh) * | 2018-09-27 | 2019-01-01 | 杭州电子科技大学 | 一种工业加热炉的二自由度内模控制分析方法 |
CN109557810A (zh) * | 2018-11-29 | 2019-04-02 | 杭州电子科技大学 | 一种基于新型二自由度内模pid的加热炉温度控制方法 |
CN110109360A (zh) * | 2019-05-22 | 2019-08-09 | 杭州电子科技大学 | 一种工业过程广义大琳响应控制方法 |
CN110879576A (zh) * | 2019-10-15 | 2020-03-13 | 杭州电子科技大学 | 水泥熟料冷却过程的分数阶模型广义二自由度控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113359467B (zh) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Borase et al. | A review of PID control, tuning methods and applications | |
KR100407203B1 (ko) | 다변수예측제어기에서독립적인피드포워드제어의결합방법 | |
CN109557810B (zh) | 一种基于新型二自由度内模pid的加热炉温度控制方法 | |
CN105892296B (zh) | 一种工业加热炉***的分数阶动态矩阵控制方法 | |
CN112462614B (zh) | 改进的ladrc线性自抗扰控制***及参数整定方法 | |
Padhan et al. | Modified Smith predictor and controller for time delay processes | |
CN105159065B (zh) | 一种非线性自抗扰控制***稳定性判断方法 | |
Ahamad et al. | A comparative study of PID controller tuning techniques for time delay processes | |
CN110109360B (zh) | 一种工业过程广义大琳响应控制方法 | |
Seghiri et al. | Fractional order adaptive MRAC controller design for high-accuracy position control of an industrial robot arm | |
CN109001975B (zh) | 一种工业加热炉多模型分数阶控制方法 | |
CN113359467A (zh) | 一种工业过程中基于分数阶模型的改进内模控制方法 | |
CN113377008A (zh) | 一种pid控制***及参数整定方法 | |
CN113467340A (zh) | 改进的一阶线性自抗扰控制***及其参数整定方法 | |
CN108897223B (zh) | 一种工业加热炉分数阶预测控制方法 | |
Ren et al. | Design and implementation of the PI-type active disturbance rejection generalized predictive control | |
Sladka et al. | Position control with 2dof pi controller | |
Sethi et al. | Design of fractional order MRAPIDC for inverted pendulum system | |
CN114019786A (zh) | 一种pi在线切换到pid的控制***及参数整定方法 | |
Chowdary et al. | Robust controller design for first order plus time delay systems using kharitonov theorem | |
Neçaibia et al. | Optimal auto-tuning of fractional order PI λ D μ controller for a DC motor speed using Extremum seeking | |
Reyes et al. | A PID using a non-singleton fuzzy logic system type 1 to control a second-order system | |
Motorga et al. | Adaptive Control using AI for the Induction Motors which Drive the Hearth of an Industrial Furnace | |
Oprzędkiewicz et al. | The fractional order PID control of the forced air heating system | |
Sivaramakrishnan et al. | Design of Hybrid control for Isothermal Continuous stirred tank Reactor |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Zhou Baoguo Inventor after: Wu Feng Inventor after: Zhang Ridong Inventor before: Zhou Baoguo Inventor before: Wu Feng Inventor before: Zhang Ridong |
|
GR01 | Patent grant | ||
GR01 | Patent grant |