[0003] 本发明的目的是为了有效的处理系统过程建模不够准确和控制特性偏移极端的问题,而提出了一种水泥熟料冷却过程的分数阶模型广义二自由度控制方法。
[0004] 本发明首先建立一个分数阶广义过程模型,定义新型的广义输入得到离散系统的广义过程高阶模型;然后,引入闭环等效的广义开环控制算法,结合二自由度多性能控制思想设计广义二自由度控制器。最后以水泥生产过程中的水泥熟料冷却环节为例,说明本发明能够有效的处理系统过程建模不够准确和控制特性偏移极端的问题。
[0005] 本发明方法的步骤包括:
[0006] 步骤1、建立系统的分数阶广义过程模型,具体步骤如下:
[0007] 1.1假设系统被控对象的分数阶传递函数模型,表示如下:
[0008]
[0009] 其中,K为分数阶模型增益常数,T是分数阶模型惯性时间常数,s为拉普拉斯变换算子,θ为模型滞后时间,0<ω≤1为分数阶微分阶次。
[0010] 1.2分数阶微分项的高阶近似,采用方法如下:
[0011]
[0012] 其中,N为选定的近似阶次,Wh和Wb分别为分数阶拟合频率的上限值和下限值。
[0013] 1.3根据步骤1.2,并通过Riemann-Liouville定理可以将步骤1.1中的分数阶模型转换为如下的高阶模型:
[0014]
[0015] 其中,i1、i2、...、in为高阶模型分母的阶次,且为整数;j1、j2、...、jn为高阶模型分子的阶次,且为整数;n1、n2、...、nn和m1、m2、...、mn分别为方程多项式的系数。
[0016] 1.4忽略系统扰动,根据步骤1.3可以得到系统离散过程线性模型输入输出关系,表示如下:
[0017] y(z)=Gp(z)u(z)
[0018] 其中,z是离散变换因子,y(z)、u(z)分别是过程输出、输入的z变换;为过程被控对象的脉冲传递函数。
[0019] 1.5根据内模控制原理,建立系统的控制过程模型,如下;
[0020]
[0021] 其中, 是系统内部模型脉冲传递函数,一般设为Gu+(z)是全通滤波器脉冲传递函数,包含时滞和有半平面零点,Gu-(z)是具有最小相位系统脉冲传递函数;GI(z)为内模控制脉冲传递函数模型,且GI(z)=Gu-(z)-1F(z),λ为滤波系数,F(z)为滤波器脉冲传递函数,Gu-(z)-1为Gu-(z)的倒数,z-1为差分算子。
[0022] 1.6根据系统的模型对偶控制,将步骤1.4中的过程线性模型改写为如下形式:
[0023] A(z)y(z)=B+(z)B-(z)u(z)
[0024] 其中,
[0025] a0=1;
[0026] β0≠0; b0≠0。
[0027] A(z)为Gp(z)的最小公分母多项式,B+(z)为Gp(z)传递函数带z-1项和时间延迟的零点多项式,B-(z)为所有在开单位圆内零点的多项式;z-i为z-1的i次方;ai、βi、bi分别为对应多项式的过程系数。
[0028] 1.7根据步骤1.6,定义一个新型的广义输入,形式如下:
[0029] V(z)=Q(z)B-(z)u(z)+P(z)y(z)
[0030] 其中, q0=1;V(z)是定义的新型广义输入z变换形式,pi,qi是对应过程系数。
[0031] 1.8根据步骤1.7,可以将步骤1.6中的过程模型改写为,如下形式:
[0032] T(z)y(z)=B+(z)V(z)
[0033] 其中,T(z)是过程极点多项式,且T(z)=A(z)Q(z)+B+(z)P(z)。将上式等效为闭环过程,可得到: 为广义被控对象脉冲传递函数,e(z)为等效误差,GI(z)=T(z)M(z),M(z)=F(z)/B+(z);与步骤1.6类比可得到T(z)的多项式形式类比于A(z)的形式,采用收缩原开环极点配置方法可得到: 将其带入
式T(z)=A(z)Q(z)+B+(z)P(z)中,通过丢番图方程求解可得Q(z)和P(z),0≤α≤1为模型广义调节参数。
[0034] 步骤2、设计分数阶广义二自由度过程控制器,具体步骤如下:
[0035] 2.1将闭环反馈响应等效为广义开环形式,如下:
[0036]
[0037] 其中,ysp(z)为系统设定值的z变换形式,Ku=B+(1)/T(1),为系统的稳态增益;D(z)为系统的反馈控制器。
[0038] 2.2根据步骤2.1可以求得系统控制器的形式,如下所示:
[0039]
[0040] 2.3根据步骤2.2,可以得到分数阶系统广义控制量展开形式:
[0041]
[0042] 为了加快系统响应速度,将上式中的偏差项系数用系数ψ替换,且ψ>1/KU,即可以得到新型广义控制律为:
[0043]
[0044] 2.4根据步骤1.7、步骤1.8和步骤2.3可以得到系统过程控制律,如下:
[0045] Q(z)B-(z)[KuT(z)-B+(z)]u(z)=(ψ-Ku)T(z)e(z)-P(z)[KuT(z)-B+(z)]y(z)[0046] 2.5根据步骤2.1到步骤2.4,依次做循环求解得到系统的分数阶模型广义响应过程控制量u(z),再将其作用于被控对象。
[0047] 本发明通过数据采集、模型建立、预测机理、优化等手段,提出了一种水泥熟料冷却过程的分数阶模型广义响应控制方法,首先该方法通过分数阶处理广义模型多项式提高了系统的预测精度,增强了系统准确性;其次采用二自由度分别控制系统的响应速度和系统超调稳定性,可有效提高系统的控制性能。