【摘要】 本发明涉及炼油工业加热炉温度过程的伴随矩阵解耦预测控制方法。传统方法控制手段与控制参数完全依赖技术人员经验，控制效果不好。本发明方法首先基于炼油工业加热炉过程模型建立伴随矩阵解耦过程模型，挖掘出基本的过程特性；然后基于该伴随矩阵解耦模型建立预测控制回路；最后通过计算预测控制器的参数，将过程对象实施预测控制。本发明方法通过数据采集、过程处理、预测机理、数据驱动、优化等手段，确立了一种炼油工业加热炉温度过程的伴随矩阵解耦预测控制方法，利用该方法可有效提高控制的精度，提高控制平稳度。本发明方法有效减少理想工艺参数与实际工艺参数的误差，保证控制装置操作在最佳状态。 【专利类型】发明授权 【申请人】杭州电子科技大学 【申请人类型】学校 【申请人地址】310018 浙江省杭州市下沙高教园区2号大街 【申请人地区】中国 【申请人城市】杭州市 【申请人区县】上城区 【申请号】CN201210358368.6 【申请日】2012-09-24 【申请年份】2012 【公开公告号】CN102880047B 【公开公告日】2014-09-17 【公开公告年份】2014 【授权公告号】CN102880047B 【授权公告日】2014-09-17 【授权公告年份】2014.0 【IPC分类号】G05B13/00 【发明人】张日东; 孔亚广; 薛安克; 王建中; 葛铭 【主权项内容】炼油工业加热炉温度过程的伴随矩阵解耦预测控制方法，其特征在于该方法的具体步骤是：(1)利用炼油工业加热炉温度过程模型建立伴随矩阵解耦模型，具体方法是：首先采集炼油工业加热炉温度过程的输入输出数据，利用该数据建立输入输出模型如下：Y(s)＝G(s)U(s)其中Y(s)、G(s)、U(s)分别为输出向量、传递函数矩阵、输入向量；Y(s)=y1(s)y2(s)· · · yN(s), ]]>U(s)=u1(s)u2(s)· · · uN(s)]]>G11(s), G12(s), …, GNN(s)表示过程的各回路传递函数，uk(s)和yk(s)分别为输入、输出变量的拉普拉斯变换，k＝1, 2, …, N，s为拉普拉斯变换算子，N为过程的输入输出变量个数，所述的输入输出数据为数据采集器中存储的数据；进一步对上述方程选取伴随矩阵解耦阵为：N(s)＝adjG(s)其中，N(s)是伴随矩阵解耦阵，adjG(s)为G(s)的伴随矩阵；将上述伴随矩阵解耦阵与过程输入输出模型合并得到：W(s)＝G(s)·N(s)＝G(s)·adjG(s)＝diagdetG(s)其中，W(s)是得到的解耦过程模型，detG(s)为G(s)的行列式，diagdetG(s)为以G(s)的行列式为元素的对角矩阵；将上述解耦过程模型通过离散化处理成过程的离散表示方式：A(z-1)yk(t)＝z-(d+1)B(z-1)uk(t)其中，A(z-1)、B(z-1)为输入变量和输出变量的多项式，d为过程的时间滞后，z-1为后移一步算子，z-(d+1)为后移d+1步算子，yk(t)、uk(t)为分别与yk(s)、uk(s)对应的离散化输出和输入变量，k＝1, 2, …, N；(2)基于该解耦过程模型设计预测控制器，具体方法是：a、定义该预测控制器多步预测输出的预测公式为：yk(t+1)=Σ i=1nA1, iyk(t+1-i)+Σ i=0mB1, iΔ uk(t-d-i)]]>其中n、m分别是输出预测和输入预测的阶次，A1, i、B1, i为输出预测和输入预测的系数，Δ为差分算子，yk(t+1)为第t+1时刻的输出预测，yk(t+1-i)、Δuk(t-d-i)为与A1, i、B1, i对应的各时刻输出变量、输入增量变量，t为当前时刻，i为第i步预测时刻；b、定义参考轨迹的组成为：yr(t+d)＝yk(t+d)yr(t+d+j)＝αyr(t+d+j-1)+(1-α)s其中，yk(t+d)是第t+d时刻的输出预测，yr(t+d)、yr(t+d+j)、yr(t+d+j-1)分别是第t+d时刻、第t+d+j时刻、第t+d+j-1时刻的参考轨迹，α是参考轨迹的柔化系数，s是过程输出设定值；c、定义输出预测和参考轨迹的向量形式为：Y＝Ym+GUYr＝(yr(t+d+1), yr(t+d+2), …, yr(t+p))T其中，Y为输出预测向量，Ym为依据a.步骤并令预测控制器多步预测输出的预测公式中输入增量变量为零时求得的数值，GU为依据步骤a并令预测控制器多步预测输出的预测公式中输入增量变量不为零时求得的数值；Yr为参考轨迹向量，yr(t+d+i), i＝1, 2, …, p为第t+d+i时刻的参考轨迹数值，p为预测步长，T为求转置算子；G、U具体表现为：G=B1, 0B2, 0B1, 00· · · · · · Bp-d, 0Bp-d-1, 0· · · B1, 0]]>U＝(Δuk(t), Δuk(t+1), …, Δuk(t+p-d-1))T其中，Bi, 0, i＝1, 2, …, p-d为相应的输入增量变量的系数，Δuk(t+i), i＝0, 1, …, p-d-1为对应的第t+i时刻输入增量变量；d.计算控制器的参数，具体是：Δuk(t)＝qT(Yr-Ym)其中qT为控制器求取系数向量；qT＝(GTG+β2I)-1GT其中β2、I为加权项系数。 【当前权利人】杭州电子科技大学 【当前专利权人地址】浙江省杭州市下沙高教园区 【专利权人类型】公立 【统一社会信用代码】12330000470009026T 【引证次数】4.0 【自引次数】1.0 【他引次数】3.0 【家族引证次数】4.0 【家族被引证次数】6

未经允许不得转载：http://www.zhongzhencnc.com/1791487711.html