一种废塑料裂解炉炉膛温度的分布式PID型预测函数控制方法技术

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

【技术实现步骤摘要】

本专利技术属于自动化
，涉及一种废塑料裂解炉炉膛温度的分布式PID型预测函数控制方法。
技术介绍
随着计算机网络技术的发展，控制作为信息处理手段的一种已不仅仅局限于集中式的控制结构，而是更多为分布式的控制结构所取代，给传统的控制问题带来了诸多挑战。分布式预测函数控制(DPFC)作为预测控制在分布式控制结构中的典型应用，在有效解决复杂高维大规模系统在线优化控制问题的同时，通常存在对模型要求高、稳态性能差、不能快速抑制未知干扰等不足。因此对于模型阶次，环境扰动、非线性等方面存在较大不确定性的控制系统，常规的DPFC方法往往难以达到预期的控制效果。而传统的PID控制由于其控制结构简单、操作方便、鲁棒性强等优点，至今仍被广泛的应用于实际工业过程中。如果能在实际过程中将PID控制与DPFC方法相结合，将进一步推进预测控制在分布式控制结构中的发展，在保证系统良好控制性能的同时，有效增加了控制参数设计的灵活性。
技术实现思路
本专利技术的目的是为了克服常规DPFC方法在应用于在模型阶次、环境扰动、非线性等方面存在不确定性较大的控制系统中的不足之处，提出了一种废塑料裂解炉炉膛温度的分布式PID型本文档来自技高网...

【技术保护点】
一种废塑料裂解炉炉膛温度的分布式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，采样时刻顺...

【技术特征摘要】
1.一种废塑料裂解炉炉膛温度的分布式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,…...

【专利技术属性】
技术研发人员：张日东，汪大卫，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33