1.电加热炉的遗传优化的多模型预测函数控制方法,其特征在于该方法包括以下步骤:
1建立被控对象的多模型
1.1根据工作区域,按照工作的温度范围进行i等分,i是要进行工作区域划分的个数;
1.2在每个工况区间中采集实际过程对象的实时阶跃响应数据,利用该数据建立被控对象的分数阶传递函数模型Mi,形式如下:
其中,Mi为第i个子模型,α1,i为第i个系统的微分阶次,T1,i为相应的系数,S为拉普拉斯变换算子,Km,i为模型比例增益系数,τm,i为模型的滞后时间常数;
1.3通过Oustaloup近似方法将分数阶模型数值化处理为:
其中,α为分数阶微分阶次,0<α<1,N为选定的近似阶次,Kα=Whα,Wn'=WbWu(2n-1-α)/N,Wn=WbWu(2n-1+a)/N, Wh和Wb分别为拟合频率的上限值和下限值;
1.4根据步骤1.3中Oustaloup近似方法,将步骤1.2中的分数阶模型转化为整数高阶模型,再通过采样时间TS和零阶保持器后离散化得到如下形式:
其中,fj,hj(j=1,2,…,l)均为离散近似后得到的系数,d=τ/TS为过程的滞后时间,l为离散模型的长度,y(k)为k时刻的过程模型输出,u(k-d-1)为过程模型在k-d-1时刻的输入值;
为了减少误差通过对模型进行一阶向后差分,得到如下形式:
1.5选取系统的状态变量如下:
ΔXm(k)=[Δy(k),Δy(k-1),…,Δy(k-l),Δu(k-1),…,Δu(k-l+1-d)]T (5)结合步骤1.4,得到被控对象的状态空间模型,形式如下:
其中,T为矩阵的转置符号,ΔXm(k)的维数为(2l+d-1)×1;
Bm=[0 … 0 1 0 … 0]TCm=[1 0 … 0 0 … 0]输出的误差定义为:
e(k)=y(k)-r(k) (7)其中r(k)为参考轨迹;
1.6预测函数控制的输入量是由一组相对应的基函数线性组合,表述为:
其中μi为基函数的系数,fj(i)是基函数在i时刻的值
1.7由于实际的过程存在误差,实际与模型之间的误差为:
1.8通过滚动优化,为了求取当前最优的控制量,目标函数选择如下:
2.利用遗传算法优化预测函数中基函数的个数
2.1设定种群染色体大小为N,进化的最大迭代次数为N1,随机初始化种群,通过适应度函数优化基函数个数目可得目标函数:f=J+ω·nr
其中ω为权重系数,J最优目标函数;
2.2采用十进制编码方式对染色体进行编码,第i个的染色体表示为:
其中i=1,2,…,N,N是种群的大小,nr是基函数的个数,
1≤nr≤D,D是基函数的最大个数;染色体Ci′中的元素为:ci=fmin+r(fmax-fmin),1≤i≤nr其中r是位于[0,1]之间的随机变量fmax,fmin分别为基函数的最大值和最小值;
2.3选取染色体的操作算子的具体步骤为:
2.3.1染色体交叉运算;选取交叉算子Pc,使染色体Ci′和下一个染色体Ci+1′以概率Pc进行交叉运算,产生下一代染色体;
2.3.2染色体修剪因子运算;以Pn的概率修改基函数的个数nr,进而改变染色体中的元素Ci′,从而完成染色个体变异的操作;
2.4依照步骤2.2到2.3进行循环重复优化搜索,如果达到最大的进化次数N1结束优化搜索计算,得到经过遗传算法优化后的染色体Ci′,进而得到最优的基函数个数nr';
2.5在求取最优的基函数的个数的情况下,求取目标函数的最优值从而得到控制量的参数为:
其中:
μj=-(0,…,1,0)(GTQG)-1GTQ[Yr(k)-Fu(k-1)+SΔR]
3.多模型的切换
3.1计算当前时刻每个子模型的输出yi(t),并计算此时子模型的输出与当前时刻对象的实际输出的偏差如下:ei(t)=yout(t)-yi(t),i=1,2,…N (15)其中yout(t)为对象的实际输出,ei(t)代表当前时刻第i个子模型与对象实际输出的偏差;
3.2计算切换指标函数;
其中α表示当前时刻的误差权重,β表示过去时刻的误差权重,ρ表示性能指标的记忆效应,L表示所取得过去时刻的长度,N表示模型的个数,Ji(k)表示第i个子模型的切换指标函数,根据求出模型性能指标函数中的最小的性能指标,确定当前时刻的模型的控制器:min(i)=min(Ji(k)) (17)当前控制器的输入量为当前时刻性能指标函数最小的所对应的控制器:u(t)=umin(i)(t) (18)
3.3在下一时刻到来时依照步骤2.1到3.2中的方式求解分数阶的多模型预测函数控制器的控制量,再将控制量作用于被控对象,并按照此步骤循环操作下去。