一种焦炭炉炉膛压力的分布式PID型动态矩阵控制方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106444388A(43)申请公布日2017.02.22(21)申请号201611111853.8(22)申请日2016.12.06(71)申请人杭州电子科技大学地址310018浙江省杭州市下沙高教园区2号大街(72)发明人张日东汪大卫(74)专利代理机构杭州奥创知识产权代理有限公司33272代理人王佳健(51)Int.Cl.G05B13/04(2006.01)权利要求书4页说明书9页(54)发明名称一种焦炭炉炉膛压力的分布式PID型动态矩阵控制方法(57)摘要本发明公开了一种焦炭炉炉膛压力的分布式PID型动态矩阵控制方法。本发明首先通过采集实时阶跃响应数据建立被控对象的动态矩阵模型向量，再将大规模系统的在线优化实现问题转化为各个小规模子系统的优化求解问题，并把网络环境下的每个子系统看作一个智能体，同时各智能体之间通过网络通信完成信息交换。然后通过引入PID算子建立一种改进的性能指标，并依据纳什最优思想来设计各智能体的PID型动态矩阵控制器，再将当前时刻的即时控制律作用于每个智能体，并将时域滚动到下一时刻，重复上述优化过程，从而完成整个大规模系统的优化任务。本发明在保证系统整体控制品质的同时，有效弥补了传统DDMC方法的不足，并提高了控制参数设计的自由度。CN106444388ACN106444388A权利要求书1/4页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：Taij＝[aij(1),aij(2),…,aij(Lij)]其中T为矩阵的转置符号,aij(k)为t＝kTs时刻第j个输入对第i个输出的阶跃响应采样值,Lij为第j个输入对第i个输出的建模时域；步骤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个智能体当前k时刻的模型预测初始响应值yi,0(k)首先，在k-1时刻加入控制增量Δu1(k-1),Δu2(k-1),…,Δun(k-1)，得到第i个智能体的模型预测值yi,P(k-1)：其中，Tyi,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)],TTAii,0＝[aii(1),aii(2),…,aii(L)],Aij,0＝[aij(1),aij(2),…,aij(L)]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)表示第i2CN106444388A权利要求书2/4页个智能体k-1时刻对k,k+1,…,k+L-1时刻的初始预测值，Aii,0,Aij,0分别为第i个智能体和第j个智能体对第i个智能体的阶跃响应数据建立的矩阵，Δu1(k-1),Δu2(k-1),…,Δun(k-1)为k-1时刻各智能体的输入控制增量；然后，得到k时刻第i个智能体的模型预测误差值ei(k)：ei(k)＝yi(k)-yi,1(k|k-1)其中yi(k)表示k时刻测得的第i个智能体的实际输出值；进一步得到k时刻修正后的模型输出值yi,cor(k)：yi,cor(k)＝y