[发明专利]一种废塑料裂解炉炉膛温度的分布式PID型预测函数控制方法

摘要

本发明公开了一种废塑料裂解炉炉膛温度的分布式PID型预测函数控制方法。本发明首先通过采集控制系统的实时阶跃响应数据建立过程对象的阶跃响应模型向量，再将大规模系统的在线优化问题转化为各个小规模子系统的优化求解问题，并把网络环境下的每个子系统看作一个智能体，同时各智能体之间通过网络通信实现信息共享。然后通过引入PID算子对DPFC方法的性能指标进行改进，并设计各智能体的PID型预测函数控制器。再将当前时刻所得即时控制律作用于每个智能体，并将时域滚动至下一时刻，最后通过依次迭代循环完成整个大规模系统的优化任务。本发明有效弥补了传统 DPFC方法的不足，并提高了控制参数设计的自由度。

主权项

一种废塑料裂解炉炉膛温度的分布式PID型预测函数控制方法，其特征在于该方法包括以下步骤：步骤1.通过废塑料裂解炉炉膛温度的实时阶跃响应数据建立被控对象的阶跃响应模型向量，具体是：1.1根据分布式控制的思想，将一个N输入N输出的大规模系统分散为N个智能体子系统；1.2在稳态工况下，以第j个智能体控制量为输入对第i个智能体输出量进行阶跃响应实验，分别记录第j(1≤j≤N)个输入对第i(1≤i≤N)个输出的阶跃响应曲线；1.3将步骤1.2得到的阶跃响应曲线进行滤波处理，拟合成一条光滑曲线，记录光滑曲线上每个采样时刻对应的阶跃响应数据，第一个采样时刻为Ts，相邻两个采样时刻的间隔时间为Ts，采样时刻顺序为Ts、2Ts、3Ts……；被控对象的阶跃响应将在某一个时刻tL＝LijTs后趋于平稳，当aij(t)(t>Lij)与aij(Lij)的误差和测量误差有相同的数量级时，即可认为aij(Lij)近似等于阶跃响应的稳态值；建立第j个输入对第i个输出之间的阶跃响应模型向量aij：aij＝[aij(1),aij(2),…,aij(Lij)]T其中aij(k)为t＝kTs时刻第j个输入对第i个输出的阶跃响应采样值,Lij为第j个输入对第i个输出的建模时域，T为矩阵的转置符号；步骤2.设计第i个智能体的PID型预测函数控制器，具体如下：2.1利用步骤1获得的阶跃响应模型向量aij建立被控对象的阶跃响应矩阵，其形式如下：其中Aij为第j个输入对第i个输出的P×M阶阶跃响应矩阵，P为预测控制的优化时域长度，M为预测控制的控制时域长度，且Lij＝L(1≤i≤3,1≤j≤3),M<P<L，L为系统的统一建模时域，N＝3为输入输出个数；2.2根据过程对象的期望值及跟踪特性选择基函数，并建立第i个智能体的控制量：ui(k+j)=Σn=1Eμi,n(k)×fi,kn(j),j=0,1,...,P-1]]>其中ui(k+j)表示第i个智能体在k+j时刻的控制量，E表示为基函数的个数，μi,n表示基函数的线性加权系数，fi,kn(n＝1,2,…,E)表示第i个智能体的基函数，fi,kn(j)表示第i个智能体的基函数fi,kn在t＝jTs时刻的值，Ts表示采样周期；2.3获取第i个智能体在当前k时刻的模型预测初始响应值yi,0(k)首先，在k‑1时刻加入各智能体的控制增量△u1(k‑1),△u2(k‑1),…,△un(k‑1)，得到第i个智能体的模型预测值yi,P(k‑1)：yi,P(k-1)=yi,0(k-1)+Aii,0Δui(k-1)+Σj=1,j≠inAij,0Δuj(k-1)]]>△ui(k+j)＝ui(k+j)‑ui(k+j‑1)结合步骤2.2进一步推导可得yi,P(k-1)=yi,0(k-1)-Aii,0ui,0+Gii,0(k-1)μi(k-1)+Σj=1,j≠inGij,0(k-1)μj(k-1)]]>其中，yi,P(k‑1)＝[yi,1(k|k‑1),yi,1(k+1|k‑1),…,yi,1(k+L‑1|k‑1)]Tyi,0(k‑1)＝[yi,0(k|k‑1),yi,0(k+1|k‑1),…,yi,0(k+L‑1|k‑1)]T,Aii,0＝[aii(1),aii(2),…,aii(L)]T,Aij,0＝[aij(1),aij(2),…,aij(L)]TFi,0E＝[fi,k1(‑1),fi,k2(‑1),…,fi,kE(‑1)],μi(k)＝[μi,1(‑1),μi,2(‑1),…,μi,E(‑1)]TGij,0＝Aij,0Fj,0E,ui,0＝ui(k‑2)yi,1(k|k‑1),yi,1(k+1|k‑1),…,yi,1(k+L‑1|k‑1)分别表示第i个智能体在k‑1时刻对k,k+1,…,k+L‑1时刻的模型预测值，yi,0(k|k‑1),yi,0(k+1|k‑1),…,yi,0(k+L‑1|k‑1)分别表示第i个智能体k‑1时刻对k,k+1,…,k+L‑1时刻的初始预测值，Aii,0,Aij,0分别为第i个智能体和第j个智能体对第i个智能体的阶跃响应数据建立的矩阵，ui,0为第i个智能体k‑2时刻的控制输入；然后得到第i个智能体在k时刻的模型预测误差值ei(k)：ei(k)＝yi(k)‑yi,1(k|k‑1)其中yi(k)表示在k时刻测得的第i个智能体实际输出值；进一步得到k时刻修正后的模型输出值yi,cor(k)：yi,cor(k)＝yi,0(k‑1)+h*ei(k)其中，yi,cor(k)＝[yi,cor(k|k),yi,cor(k+1|k),…,yi,cor(k+L‑1|k)]T,h＝[1,α,…,α]Tyi,cor(k|k),yi,cor(k+1|k),…,yi,cor(k+L‑1|k)分别表示第i个智能体k时刻对k,k+1,…,k+L‑1时刻预测模型的修正值，h为误差补偿的权向量，α为误差校正系数；最后得到第i个智能体在k时刻的模型预测的初始响应值yi,0(k)：yi,0(k)＝Syi,cor(k)其中，S为L×L阶的状态转移矩阵，2.4获取第i个智能体在M个连续的控制增量作用下的预测输出值yi,PM，具体是：yi,PM(k)=yi,P0(k)-A0i,P0uP0+Gii(k)μi(k)+Σj=1,j≠inGij(k)μj(k)]]>其中，yi,PM(k)＝[yi,M(k+1|k),yi,M(k+2|k),…,yi,M(k+P|k)]Tyi,P0(k)＝[yi,0(k+1|k),yi,0(k+2|k),…,yi,0(k+P|k)]TA0i,P0＝[A0i1,P0,A0i2,P0,…,A0iN,P0]T,A0ij,P0＝[aij(1),aij(2),…,aij(M),…,aij(P)]TuP0＝[u1,P0,u2,P0,…,uN,P0]T,μi(k)＝[μi,1(k),μi,2(k),…,μi,E(k)]TGij(k)＝AijFj,Eyi,P0(k)是yi,0(k)的前P项，yi,0(k+1|k),yi,0(k+2|k),…,yi,0(k+P|k)分别表示第i个智能体k时刻对k+1,k+2,…,k+P时刻的初始预测输出值；2.5选取第i个智能体的性能指标Ji(k)，形式如下：minJi(k)=(wi(k)-yi,PM(k))TKIi(wi(k)-yi,PM(k))+(Δwi(k)-Δyi,PM(k))TKpi(Δwi(k)-Δyi,PM(k))+(Δ2wi(k)-Δ2yi,PM(k))TKdi(Δ2wi(k)-Δ2yi,PM(k))+Δui,M(k)TRiΔui,M(k)]]>其中，wi(k)＝[ωi(k+1),ωi(k+2),…,ωi(k+P)]Tωi(k+ε)＝λεyi(k)+(1‑λε)c(k)(ε＝1,2,…,P)△wi(k)＝[△ωi(k+1),△ωi(k+2),…,ωi(k+P)]T△yi,PM(k)＝[△yi,M(k+1|k),△yi,M(k+2|k),…,△yi,M(k+P|k)]T△2wi(k)＝[△2ωi(k+1),△2ωi(k+2),…,△2ωi(k+P)]T△2yi,PM(k)＝[△2yi,M(k+1|k),△2yi,M(k+2|k),…,△2yi,M(k+P|k)]T△ωi(k+ε)＝ωi(k+ε)‑ωi(k+ε‑1)△yi,M(k+ε|k)＝yi,M(k+ε|k)‑yi,M(k+ε‑1|k)△2ωi(k+ε)＝△ωi(k+ε)‑△ωi(k+ε‑1)△2yi,M(k+ε|k)＝△yi,M(k+ε|k)‑△yi,M(k+ε‑1|k)分别为第i个智能体的比例系数矩阵、积分系数矩阵、微分系数矩阵，为第i个智能体的控制加权系数矩阵，ωi(k+ε)为第i个智能体在k+ε时刻的参考轨迹，yi(k)为k时刻第i个智能体的实际输出，c(k)为k时刻第i个智能体的期望输出，λ为参考轨迹的柔化因子；2.6对步骤2.5中的性能指标进行转换，形式如下：minJi(k)=E0i(k)TKIiE0i(k)+ΔE0i(k)TKpiΔE0i(k)+Δ2E0i(k)TKdiΔ2E0i(k)+Δui,M(k)TRiΔui,M(k)]]>E0i(k)=wi(k)-yi,PM(k)=[e0i(k+1),e0i(k+2),...,e0i(k+P)]T]]>ΔE0i(k)=Δwi(k)-Δyi,PM(k)=[Δe0i(k+1),Δe0i(k+2),...,Δe0i(k+P)]T]]>Δ2E0i=Δ2wi(k)-Δ2yi,PM(k)=[Δ2e0i(k+1),Δ2e0i(k+2),...,Δ2e0i(k+P)]T]]>进一步得到Δe0i(k+ϵ)=Δωi(k+ϵ)-Δyi,M(k+ϵ|k)=ωi(k+ϵ)-yi,M(k+ϵ|k)-(ωi(k+ϵ-1)-yi,M(k+ϵ-1|k))=e0i(k+ϵ)-e0i(k+ϵ-1)]]>同理可得Δ2e0i(k+ϵ)=Δe0i(k+ϵ)-Δe0i(k+ϵ-1)]]>其中引入矩阵进而有ΔE0i(k)=S1E0i(k)Δ2E0i(k)=S1ΔE0i(k)=S12E0i(k)]]>进一步可将性能指标变换为minJi(k)=E0i(k)TKIiE0i(k)+E0i(k)TS1TKpiS1E0i(k)+E0i(k)T(S12)TKdi(S12)E0i(k)+Δui,M(k)TRiΔui,M(k)=E0i(k)TQiE0i(k)+Δui,M(k)TRiΔui,M(k)]]>其中，2.7依据纳什最优的概念，对性能指标求极值，得到形式如下的纳什最优解：μi*(k)=Dii[wi(k)-yi,P0(k)+A0i,P0uP0-Σj=1,j≠inGijμj*(k)]]]>其中，2.8重复步骤2.2至步骤2.7，得到第i个智能体在k时刻的新一轮迭代最优解为：μil+1(k)=Dii[wi(k)-yi,P0(k)+A0i,P0uP0-Σj=1,j≠inGijμjl(k)]]]>进一步得到k时刻整个系统的纳什最优解：μl+1(k)＝D1[w(k)‑YP0(k)+A0i,P0uP0]+D0μl(k)其中：μl+1(k)=[μ1l+1(k),μ2l+1(k),...,μnl+1(k)]T]]>μl(k)=[μ1l(k),μ2l(k),...,μnl(k)]T]]>ω(k)＝[ω1(k),ω2(k),…,ωn(k)]T，YP0(k)＝[y1,P0(k),y2,P0(k),…,yn,P0(k)]T2.9由第i个智能体k时刻的纳什最优解得到最优控制量ui(k)，并将其作用于第i个智能体；2.10在下一时刻，重复步骤2.2到2.9继续求解第i个智能体的纳什最优解进而得到整个大规模系统的最优解μ*(k+1)，并依次循环。