遗传算法优化加热炉温度的分数阶预测函数控制方法
【专利摘要】本发明公开了一种遗传算法优化加热炉温度的分数阶预测函数控制方法。本发明首先采用Oustaloup近似方法将分数阶系统近似为整数阶系统，基于Oustaloup近似模型建立预测输出模型，然后将整数阶预测函数控制方法扩展到分数阶预测函数控制方法中，将分数阶微分算子引入目标函数，并采用遗传算法来优化微分算子，从而使经优化后获得更加合理的控制效果。本发明方法可有效提高系统的控制性能。
【专利说明】
遗传算法优化加热妒溫度的分数阶预测函数控制方法
技术领域
[0001] 本发明属于自动化技术领域，设及一种遗传算法优化加热炉溫度的分数阶预测函 数控制方法。
【背景技术】
[0002] 化工过程是我国流程工业过程的重要组成部分，其要求是供给合格的工业产品， W满足我国工业的需要。对于许多复杂的对象使用整数阶微分方程无法精确描述，用分数 阶微分方程能更精确地描述对象特征和评估产品性能。预测函数(PFC)作为先进控制方法 的一种，具有计算量小、鲁棒性强、控制性能好等特点，在实际过程控制中获得了大量成功 的应用。对于分数阶系统，传统PID控制方法和整数阶预测函数控制方法对运类对象的控制 效果并不是很好，运就需要我们研究具备良好控制性能的控制器来控制运类用分数阶模型 描述的实际被控对象。虽然分数阶预测函数方法增加自由度，但对于调整参数也增加了难 度，针对运类问题提出了采用遗传算法来优化运些参数，从而使经优化后获得更加合理的 控制效果，保证了跟踪精度。
【发明内容】
[0003] 本发明的目的是针对分数阶系统描述的加热炉溫度过程，提供一种遗传算法优化 加热炉溫度的分数阶预测函数控制方法，W保证分数阶系统跟踪性能并保障良好的控制效 果。
[0004] 本发明方法首先采用Oustaloup近似方法将分数阶系统近似为整数阶系统，基于 Oustaloup近似模型建立预测输出模型，然后将整数阶预测函数控制方法扩展到分数阶预 测函数控制方法中，将分数阶微分算子引入目标函数，并采用遗传算法来优化微分算子，从 而使经优化后获得更加合理的控制效果。
[0005] 本发明方法的步骤包括：
[0006] 步骤1、建立工业加热炉对象的分数阶线性模型，具体方法是：
[0007] 1.1采集加热炉实际过程的实时输入输出数据，利用该数据建立被控对象在时刻t 的分数阶微分方程模型，形式如下：
[000引
[0009] 其中，曰1，曰2为微分阶次，c〇，ci，C2为相应的系数，y(t)，u(t)分别为过程的输出和输 入。
[0010] 1.2根据分数阶微积分定义，对步骤1.1中的模型进行拉氏变换，得到被控对象的 传递函掛化击血下，
[0011]
[001^ 其中，S为复变量。
[001引 1.3由Ous化Ioup近似方法得到微分算子S。的近似表达形式如下：
[0014]
[001引其中，a为分数阶微分阶次，0 < a < I，N为选定的近似阶次：
，Wb和Wh分别为选定的拟合频率的下限和上 限。
[0016] 1.4根据步骤1.3中的方法，将步骤1.2中的分数阶系统近似为整数阶高阶系统模 型，对得到的高阶模型在采样时间Ts下加零阶保持器离散化，得到如下形式的模型：
[0017]
[001引其中，門，田（j = l，2，.''Ls)均为离散近似后得到的系数，Ls为离散模型的长度。
[0019] 步骤2、设计加热炉对象的分数阶预测函数控制器，具体方法如下：
[0020] 2.1计算被控对象在预测函数控制下的i巧预测输出，形式如下：
[0021]
[0022] 其中，P为预测时域，y化+i)为k+i时刻过程的预测模型输出，i = l，2,…P。
[0023] 2.2对步骤2.1中的式子进行整理变换，得到如下模型：
[0024]
[0025]
[0026]
[0027]
[002引
[0029]
[0030] 1 其中，T为转置符号;结合上述式子，得到被控对象的预测输出模型为：
[0032]
[0033]
[0034] 2.3修正当前时刻被控对象的预测输出模型，得到修正后的对象模型，形式如下：
[0035]
[0036] E = [e化+1) ,e(k+2),,e(k+P) ]T
[0037] e(k+i) =yp化）-y(k),i = 1,2,...P
[0038] 其中，yp化)是k时刻被控对象的实际输出值，e化+i)为k+i时刻被控对象的实际输 出值与模型预测输出的差值；
[0039] 2.4选取预测函数控制方法的参考轨迹yr化+i)和目标函数Jfpfc,其形式如下：
[0040] yr(k+i)=入Iyp化）+ (1-入i)c(k),i = 0,1,...P
[0041]
[0042] 其中，yr化+i)为k+i时刻的参考轨迹，A为参考轨迹的柔化系数，C化)为k时刻的设 定值/巧/(〇表示函数f(t)在[tl，t2让的丫次积分，D为微分符号；
[00创依化GV//m.w/J - Lw?斯n，分数阶微积分定义，对上述目标函数在采样时间Ts进行 离散化，得到：
[0044]
[0045] 其中，
[0046] Yr= [yr化+1)，yr化+2),…yr(k+P)]T
[0047] A (Ts, 丫）=Tsdiag(邮-1,邮-2,...,mi,m〇)
[004引
[0049]
[0050] 2.5依据步骤2.4中的目标函数求解过程输入的最优值，即最优控制律，形式如下：
[0化1 ]
[0052] 2.6在k+L时刻，L= 1，2，3，…，依照步骤2.巧Ij步骤2.5中的步骤依次循环求解分数 阶预测函数控制器的控制量U化+1)，再将其作用于被控对象。
[0053] 步骤3、基于遗传算法优化最优控制律中的A，具体步骤是：
[0054] 3.1首先使用四位二进制进行二进制编码，得到第一代染色体。
[0055] 3.2选取遗传算法的适应度函数，并计算个体的适应度值，形式如下：
[0化6]
[0057]其中，F是个体的适应度函数。当适应度函数值大于适应度预设值fz时，遗传算法 终止。
[005引3.3利用轮转法来确定选择算子，形式如下：
[0化9]
[0060] 其中，P(Ci)是个体Cl的选择概率，f (Cl)是个体Cl的适应度。
[0061] 3.4利用步骤3.3中的选择算子将染色体适应度较高的个体选择出来W交叉概率 Pc进行交叉操作，产生下一代个体。
[0062] 3.5选取合适的变异算子，形式如下：
[0063]
[0064] 其中，a日表示变异概率Pm的初始值，b日是变异概率的程度，g式进化的代数，g日是变 异概率改变很大的进化代数，a是变异速率。
[0(?日]!V fi#个化狱 曰十丰釉掛細趙N时.得到剪接算子，形式如下；
[0066]
[0067] 其中，P取正数。
[0068] 3.7依照步骤3.2计算每个个体的适应度值，并判断是否满足终止条件，如果满足， 则为参数的最优解，进行下一步操作，如果不满足，则执行步骤3.2到步骤3.6,直到找到满 足终止条件为止。
[0069] !V S对垫化化讲巧傭招.化古加下.
[0070]
[0071]
[0072] 其中，Q是长度为L的二进制解码产生的整数，Xj,min和Xj,max分别为输入变量的最小 值和最大值，Wmax是高斯函数的基宽的最大值。
[0073] 通过步骤3优化后带入步骤2中调整微分算子参数，在下一时刻依照步骤1到步骤3 的步骤继续对实际过程进行控制，依次循环。
[0074] 本发明的有益效果:本发明的技术方案是通过数据采集、模型建立、预测机理、优 化等手段，确立了一种工业加热炉系统的分数阶预测函数控制优化方法，该方法可有效提 高系统的控制性能。
【具体实施方式】
[0075] 下面W工业加热炉溫度控制过程为例，对本发明作进一步说明。
[0076] 在工业加热炉溫度控制过程中，为评估模型的性能，由工业生产的实际应用系统 得到溫度的实时数据。
[0077] 遗传算法优化工业加热炉系统的分数阶预测函数控制方法的具体步骤包括：
[0078] 步骤1、建立实际过程中被控对象的分数阶线性模型，具体方法是：
[0079] 1.1采集实际过程对象的实时输入输出数据，利用该数据建立被控对象在时刻t的 分数阶微A方巧横巧，形式化下：
[0080]
[0081] 其中，曰1，曰2为微分阶次，co，ci，C2为相应的系数，y(t)，u(t)分别为过程的输出和输 入。
[0082] 1.2根据分数阶微积分定义，对步骤1.1中的模型进行拉氏变换，得到被控对象的 传递函才*n 了，
[0083]
[0084] 其中，S为复变量。
[0085] 1.3由Ous化Ioup近似方法得到微分算子S*"的近似表达形式如下：
[0086]
[0087] 其中，a为分数阶微分阶次，0 < a < 1，N为选定的近似阶次，
Wb和Wh分别为选定的拟合频率的下限和上 限。
[0088] 1.4根据步骤1.3中的方法，将步骤1.2中的分数阶系统近似为整数阶高阶系统模 型，对得面I的鳥阶摇巧#要睹时巧T。下巾n塞阶化桂柴唐惭化，得到如下形式的模型；
[0089]
[0090] 其中，門，田（j = l，2，'''Ls)均为离散近似后得到的系数，Ls为离散模型的长度。
[0091 ]步骤2、设计被控对象的分数阶预测函数控制器，具体方法如下：
[0092] 2.1计算被控对象在预测函数控制下的i步预测输出，形式如下：
[0093]
[0094] 其中，P为预测时域，y化+i)为k+i时刻过程的预测模型输出，i = l，2,…P。
[00M] 2.2对步骤2.1中的式子进行整理变换，得到如下模型：
[0096] AY = BYpast+CuA)+DUpast
[0097] 其中，
[009引 Y= [y化+1) ,;y(k+2),,;y(k+P) ]T
[0099] Ypast=Wk)，y(k-l)，...，y(；k-Ls+l)]T
[0100] Upast=[u化-I)，u(k-2),u(；k-Ls+l)]T
[0101]
[0102]
[0103] 其中，T为转置符号;结合上述式子，得到被控对象的预测输出模型为：
[0104]
[0105]
[0106] :模型，得到修正后的对象模型，形式如下：
[0107]
[010 引 E=[eA+l)，e(k+2)，...，e(k+P)]T
[0109] e(k+i)=yp(^k)-y(k),i = l,2,.-.P
[0110] 其中，yp化)是k时刻被控对象的实际输出值，e化+i)为k+ i时刻被控对象的实际输 出值与模型预测输出的差值；
[0111] 2.4选取预测函数控制方法的参考轨迹yr化+i)和目标函数Jfpfc，其形式如下：
[01 12] yr(k+i )=入Iyp化）+ ( 1-入I)c(k) , i = 0,1, ...P
[0113]
[0114] 其中，yr化+i)为k+i时刻的参考轨迹，A为参考轨迹的柔化系数，C化)为k时刻的设 定值八；：\化)表示函数f(t)在[tl，t2让的丫次积分，D为微分符号；
[011引依据伤麻m'aW-五绍《淡视分数阶微积分定义，对上述目标函数在采样时间Ts进行 离散化，得到：
[0116]
[0117] 其中，
[0118] Yr= [yr化+1) ,yr(k+2), ...yr(k+P) ]T
[0119] A (Ts, 丫）=Tsdiag(邮-I,邮-2,...,mi,m〇)
[0120] Wg 二 田扣_口
[01別]阳;"=l，Vg>0 时：
对 q<0,阳 ^i=O;
[0122] 2.5依据巧骤2.4中的目标函数求解过括输入的最优値，即最优控制律，形式如下：
[0123]
[0124] 2.6在k+L时刻，L= 1，2，3，…，依照步骤2.巧Ij步骤2.5中的步骤依次循环求解分数 阶预测函数控制器的控制量U化+1)，再将其作用于被控对象。
[0125] 步骤3、利用遗传算法优化最优控制律中的A，具体步骤是：
[01%] 3.1首先使用四位二进制进行二进制编码，得到第一代染色体。
[0127] 3.2洗取溃佑尊法的适应度函数，并计算个体的适应度值，形式如下：
[012 引
[0129] 其中，F是个体的适应度函数。当适应度函数值大于适应度预设值fz时，遗传算法 终止。
[0130] 3.3刑巧轮掉法夹确定选择算子，形式如下：
[0131]
[0132] 其中，P(Ci)是个体Cl的选择概率，f (Cl)是个体Cl的适应度。
[0133] 3.4利用步骤3.3中的选择算子将染色体适应度较高的个体选择出来W交叉概率 Pc进行交叉操作，产生下一代个体。
[0134] 3.5选取合适的巧异算子，形式如下：
[0135]
[0136] 其中，a日表示变异概率Pm的初始值，b日是变异概率的程度，g式进化的代数，g日是变 异概率改变很大的进化代数，a是变异速率。
[0137] 3.6在个化敬目大干种趙规横N时，得到剪接算子，形式如下：
[013 引
[0139] 其中，P取正数。
[0140] 3.7依照步骤3.2计算每个个体的适应度值，并判断是否满足终止条件，如果满足， 则为参数的最优解，进行下一步操作，如果不满足，则执行步骤3.2到步骤3.6,直到找到满 足终止条件为止。
[0141] 3.8对染色体进行解码，形式如下：
[0144] 其中，Q是长度为L的二进制解码产生的整数，Xj,min和Xj,max分别为输入变量的最小
[0142]
[0143] 值和最大值，Wmax是高斯函数的基宽的最大值。
[0145]通过步骤3优化后带入步骤2中调整微分算子参数，在下一时刻依照步骤1到步骤3 的步骤继续对实际过程进行控制，依次循环。
【主权项】
1.遗传算法优化加热炉温度的分数阶预测函数控制方法，其特征在于该方法包括以下 步骤： 步骤1、建立工业加热炉对象的分数阶线性模型，具体是： 步骤1.1采集加热炉实际过程的实时输入输出数据，利用该数据建立被控对象在时刻t 的分数阶微分方程模型，形式如下：其中，s为复变量；步骤1.3由Oustaloup近似方法得到微分算子8°的近似表达形式如下： 其中，，α2为微分阶次，CQ，C1，C2为相应的系数，y (t)，u(t)分别为过程的输出和输入； 步骤1.2根据分数阶微积分定义，对步骤1.1中的模型进行拉氏变换，得到被控对象的 传递函数形式如下：其中，α为分数阶微分阶次，0〈α〈 1，N为选定的近似阶次，〃< w::= wb和wh分别为选定的拟合频率的下限和上限； 步骤1.4根据步骤1.3中的方法，将步骤1.2中的分数阶系统近似为整数阶高阶系统模 型，对得到的高阶模型在采样时间1下加零阶保持器离散化，得到如下形式的模型：其中，Fj，^(j = l，2,…Ls)均为离散近似后得到的系数，LS为离散模型的长度； 步骤2、设计加热炉对象的分数阶预测函数控制器，具体如下： 步骤2.1计算被控对象在预测函数控制下的i步预测输出，形式如下：其中，P为预测时域，y(k+i)为k+i时刻过程的预测模型输出，i = l，2,…P; 步骤2.2对步骤2.1中的式子进行整理变换，得到如下模型： AY = BYpast+Cu(k)+DUpast 其中， Y=[y(k+l),y(k+2),---,y(k+P)]T Ypast=[y(k)，y(k-l)，H.，y(k-Ls+l)]T UPast= [u(k-l) ,u(k-2), ··· ,u(k-Ls+l)]T 1_. Δ其中，τ为转置符号;结合上述式子，得到被控对象的预测输出模型为： 其中步骤2.3修正当前时刻被控对象的预测输出模型，得到修正后的对象模型，形式如下：E=[e(k+1) ,e(k+2),··· ,e(k+P)]τ e(k+i)=yP(k)-y(k) ,i = l,2, ·-·Ρ 其中，yP (k)是k时刻被控对象的实际输出值，e (k+i)为k+i时刻被控对象的实际输出值 与模型预测输出的差值； 步骤2.4选取预测函数控制方法的参考轨迹yr (k+i)和目标函数Jfpfc，其形式如下： yr(k+i)=A1yP(k) + (l-A1)c(k) ,i = 〇,l,-··Ρ其中，(k+i)为k+i时刻的参考轨迹，λ为参考轨迹的柔化系数，c (k)为k时刻的设定值， γ$/(〇表示函数f(t)在[tits]上的γ次积分，D为微分符号； 依据GW/mm/^-Letnikov分数阶微积分定义，对上述目标函数在采样时间Ts进行离散 化，得到： 其中，Yr=[yr(k+1) ,yr(k+2) ,---yr(k+P)]T Λ ( Τ? . Υ ) = T^d i ae (mp-ι . ττιρ-9 .. mi . mn ) 步骤2.5依据步骤2.4中的目标函数求解过程输入的最优值，即最优控制律，形式如下：步骤2.6在k+L时刻，L= 1，2，3，…，依照步骤2.1到步骤2.5中的步骤依次循环求解分数 阶预测函数控制器的控制量u(k+l)，再将其作用于被控对象； 步骤3、基于遗传算法优化最优控制律中的λ，具体是： 步骤3.1首先使用四位二进制进行二进制编码，得到第一代染色体； 步骤3.2选取遗传算法的适应度函数，并计算个体的适应度值，形式如下：1 1 其中，F是个体的适应度函数;当适应度函数值大于适应度预设值fz时，遗传算法终止； 步骤3.3利用轮转法来确定选择算子，形式如下：其中，P (Ci)是个体Ci的选择概率，f (Ci)是个体Ci的适应度； 步骤3.4利用步骤3.3中的选择算子将染色体适应度较高的个体选择出来以交叉概率p。 进行交叉操作，产生下一代个体； 步骤3.5选取合适的变异算子，形式如下：其中，ao表示变异概率pm的初始值，bo是变异概率的程度，g式进化的代数，go是变异概率 改变很大的进化代数，a是变异速率； 步骤3.6在个体数目大于种群规模N时，得到剪接算子，形式如下：其中，P取正数； 步骤3.7依照步骤3.2计算每个个体的适应度值，并判断是否满足终止条件，如果满足， 则为参数的最优解，进行下一步操作，如果不满足，则执行步骤3.2到步骤3.6,直到找到满 足终止条件为止； 步骤3.8对染色体进行解码，形式如下：其中，Q是长度为L的二进制解码产生的整数，Xj,min和Xj, max分别为输入变量的最小值和 最大值，Wmax是高斯函数的基宽的最大值； 通过步骤3优化后带入步骤2中调整微分算子参数，在下一时刻依照步骤1到步骤3的步 骤继续对实际过程进行控制，依次循环。
【文档编号】G05B13/04GK105955350SQ201610534739
【公开日】2016年9月21日
【申请日】2016年7月5日
【发明人】张日东, 张俊锋, 王玉中
【申请人】杭州电子科技大学