多变量过程的蒸馏塔模型预测控制优化PID控制方法技术

本发明专利技术公开了一种多变量过程的蒸馏塔模型预测控制优化PID控制方法，本发明专利技术结合扩展非最小状态空间模型预测控制和传统PID控制算法。首先建立多变量过程的蒸馏塔状态空间模型，挖掘出基本的过程特性；然后结合蒸馏塔内多变量过程的状态过程和输出误差，建立扩展非最小状态空间模型。在模型的基础上，基于优化思想，利用模型预测控制方法来整定PID控制器的参数，最后对蒸馏塔内多变量过程实现PID控制，有效提高了多变量系统的控制性能，明显减小系统超调量，改善了控制性能，又保证了控制结构简单；弥补了传统PID控制的不足，又保证了良好的控制性能。

【技术实现步骤摘要】

本专利技术属于自动化
，涉及一种多变量过程的蒸馏塔模型预测控制优化 PID控制方法。
技术介绍
在实际工业过程中，由于PID控制结构简单，操作方便，已经广泛地应用于实际过 程控制系统中。但是，随着工业的发展，人们对产品质量、能源消耗和环境保护的要求越来 越高，对工业过程的控制精度要求也日益严格，传统的PID控制方法难以进一步提升控制 水平，控制过程变得更加复杂。针对蒸馏塔这一多变量对象，传统的PID控制的参数整定 过程比较复杂，可能往往达不到所要求的控制性能指标和控制精度。模型预测控制（MPC) 作为先进控制技术的一种，具有对模型精度要求低、在线计算方便、控制效果好等特点，在 实际复杂工业过程中已经得到非常成功的应用。基于扩展非最小状态空间模型预测控制 (EWSSMPC)在多变量过程控制中比PID控制拥有更好的控制性能，如果能够在实际过程中 将E匪SSMPC和PID控制技术相结合，系统的控制性能将得到进一步的改善，同时又能保证 控制结构的形式比较简单。
技术实现思路
本专利技术的目的是针对传统PID控制在多变量过程控制的不足之处，提出了一种多 变量过程的蒸馏塔模型预测控制优化PID控制方法，该方法结合扩展非最小状态空间模型 预测控制和传统PID控制算法，弥补了传统PID控制的不足，又保证了良好的控制性能。 本专利技术方法首先建立多变量过程的蒸馏塔状态空间模型，挖掘出基本的过程特 性；然后结合蒸馏塔内多变量过程的状态过程和输出误差，建立扩展非最小状态空间模型。 在模型的基础上，基于优化思想，利用模型预测控制方法来整定PID控制器的参数，最后对 蒸馏塔内多变量过程实现PID控制，有效提高了多变量系统的控制性能。 本专利技术的技术方案是通过数据采集、预测机理、滚动优化等手段，确立了在蒸馏塔 内多变量过程的状态空间模型预测控制优化PID控制方法，利用该方法可明显减小系统超 调量，改善了控制性能，又保证了控制结构简单。 本专利技术方法的步骤如下： 步骤（1).建立被控对象的扩展非最小状态空间模型，具体步骤如下： 1-1.通过采集被控对象的实时数据，用最小二乘法建立多输入多输出系统的模型其中， CN 105159095 A m 2/7 页 J Y (k)为P维输出，U (k)为q维输入，氏…Hn山，Iv" Ln是系统需要辨识的系数。 通过最小二乘法辨识系统模型Ir，结果如下： 其中， 1-2.将步骤1-1中辨识得到的模型进一步处理为如下形式： 其中，y (k)和u (k)分别是k时刻输出和输入，Λ为后移算子。 1-3.选取非最小状态空间变量Λ X (k)，形式如下： 其中，Λ x(k)的维数是 m = pX (n-l)+qXn。 将步骤1-2中的模型经过转换后，可得状态空间模型： Ax(k+l)、Ay(k+l)分别是第k+Ι时刻的状态和输出，Au(k)是第k时刻的输入 变量增量值，A、B、C分别对应的是状态矩阵、输入矩阵和输出矩阵。一为p维的单位矩阵， Iq为q维的单位矩阵。 1-4定义预期的输出r(k)，则跟踪误差表示为： e (k) = y (k) ~r (k) 结合步骤1-3中的状态空间模型和定义的跟踪误差得到： 其中Λ u (k)，Λ r(k+l)分别是经过后移算子后的输入和预期的输出。 为了获得扩展非最小状态空间模型，构造一个新的状态变量如下： 进一步扩展状态空间模型，形式如下： 0是一个mX q维的零矩阵，Iq是一个q维单位矩阵。 步骤（2).设计多变量过程的PID控制器，具体步骤是： 2-1.为了预测未来时刻的状态，定义如下形式的矩阵： r (k+i) = a xy (k) + (1- a x) c (k) 其中，P是预测时域，α是参考轨迹柔化系数矩阵，c(k)是在k时刻的设定值。 将未来时刻的状态写成如下形式： 5 2-1.选择性能指标函数J，形式如下： 其中，Q和γ分别是多变量过程的状态和输入的权重矩阵。 2-3.结合扩展非最小状态空间预测控制（EWSSMPC)优化思想，设计PID控制器， 将控制量写成增量形式，用如下式子表述：CN 105159095 A IX m "ti 4/7 页 其中，kp(k)，kjk)，kd(k)分别代表的是在k时刻比例系数矩阵、积分系数矩阵和 微分系数矩阵， ei (k)是在k时刻设定值矩阵和实际输出矩阵的差值。 进一步将控制量简化为如下形式： 求解上述式子，得到PID控制器的参数向量： 则得到的多变量控制系统的PID控制参数如下： 2-4.得到PID控制器的参数Mkhk1 (k)，kd(k)以后构成控制量u (k)作用于被控 对象，u (k) = u (k-1) +kp (k) (e! (k) -e! (k-1)) +Ici (k) e! (k) +kd (k) (e! (k) -2e! (k-1) +e! (k-2))。 2-5.在k+1时刻，重复步骤2-1至2-4,继续求解新的控制量u (k+1)，并依次循环 求解 kp (k+1)，Ici (k+1)，kd (k+1)，1 = 1，2, 3，…。 本专利技术有益效果如下： 本专利技术方法首先基于实时多输入多输出数据建立对象的状态空间模型，然后结合 多输入多输出对象的状态过程和输出误差建立扩展的非最小状态空间模型。在此模型的基 础上，依据E匪SSMPC的方法来优化PID控制器的参数，最后对多变量过程对象实现PID控 制，该方法明显地减小了系统超调量，获得了良好的控制效果。【具体实施方式】 以蒸馏塔这一多变量过程为例： 这里以蒸馏塔内多个变量控制为例加以描述，蒸馏塔内是一个复杂的情形，同时 可能受到多种不同的控制量的影响。 步骤（1).建立被控对象的扩展非最小状态空间模型，具体步骤如下： 1-1.通过采集被控对象的实时数据，用最小二乘法建立多输入多输出系统的模型 其中， Y(k)为p维输出，U(k)为q维输入，氏…Hn山，LyLn是系统需要辨识的系数。 通过最小二乘法辨识系统模型^，结果如下： 其中， 1-2.将步骤1-1中辨识得到的模型进一步处理为如下形式： Δ y (k+1) +H1A y (k) +Η2Δ y (k-1) +*··+ΗηΔ y (k-n+1) = L1Au (k) +L2A u (k-1) +··· +LnA u (k-n+1) 其中，y (k)和u (k)分别是k时刻输出和输入，Λ为后移算子。 1-3.选取非最小状态空间变量Λ X (k)，形式如下： 其中，Λ x(k)的维数是 m = pX (n_l)+qXn。 将步骤1-2中的模型经过转换后，可得状态空间模型：CN 105159095 A 说明书 6/7 页 B = C = Ax(k+l)、Ay(k+l)分别是第k+Ι时刻当前第1页1 2 本文档来自技高网...

【技术保护点】
多变量过程的蒸馏塔模型预测控制优化PID控制方法，其特征在于包括如下步骤：步骤(1).建立被控对象的扩展非最小状态空间模型；步骤(2).设计多变量过程的PID控制器；步骤1所述的建立被控对象的扩展非最小状态空间模型，具体步骤如下：1‑1.通过采集被控对象的实时数据，用最小二乘法建立多输入多输出系统的模型Y(k)+H1Y(k‑1)+…HnY(k‑n)＝L1U(k)+L2U(k‑1)+……+LnU(k‑n+1)其中，Y(k)=y1(k)y2(k)···yp(k),U(k)=u1(k)u2(k)···uq(k);]]>Y(k)为p维输出，U(k)为q维输入，H1…Hn,L1,L2…Ln是系统需要辨识的系数；通过最小二乘法辨识系统模型结果如下：其中，Yj=yj(n+1)yj(n+2)···yj(n+N)Hj=-yj(n)...-yj(1)UT(n+1)...UT(1)-yj(n+1)...-yj(2)UT(n+2)...UT(2)··················-yj(n+N-1)...-yj(N)UT(n+N)...UT(N)]]>1‑2.将步骤1‑1中辨识得到的模型进一步处理为如下形式：△y(k+1)+H1△y(k)+H2△y(k‑1)+…+Hn△y(k‑n+1)＝L1△u(k)+L2△u(k‑1)+…+Ln△u(k‑n+1)其中，y(k)和u(k)分别是k时刻输出和输入，△为后移算子；1‑3.选取非最小状态空间变量△x(k)，形式如下：△x(k)T＝[△y(k)T,△y(k‑1)T,…,△y(k‑n+1)T,△u(k‑1)T,△u(k‑2)T,…,△u(k‑n+1)T]其中，△x(k)的维数是m＝p×(n‑1)+q×n；将步骤1‑2中的模型经过转换后，可得状态空间模型：△x(k+1)＝A△x(k)+B△u(k)△y(k+1)＝C△x(k+1)其中，A=-H1-H1...-H1-H1-H1...-H1-H1Iq0...000...000Iq...000...00······...·········...······00...Iq00...0000...000...0000...00Iq...00······...·········...······00...000...Iq0]]>B＝[L1 0 0 … 0 Ip 0 0]C＝[Iq 0 0 … 0 0 0 0]△x(k+1)、△y(k+1)分别是第k+1时刻的状态和输出，△u(k)是第k时刻的输入变量增量值，A、B、C分别对应的是状态矩阵、输入矩阵和输出矩阵；Ip为p维的单位矩阵，Iq为q维的单位矩阵；1‑4定义预期的输出r(k)，则跟踪误差表示为：e(k)＝y(k)‑r(k)结合步骤1‑3中的状态空间模型和定义的跟踪误差得到：e(k+1)＝e(k)+CA△x(k)+CB△u(k)‑△r(k+1)其中△u(k)，△r(k+1)分别是经过后移算子后的输入和预期的输出；为了获得扩展非最小状态空间模型，构造一个新的状态变量如下：z(k)=Δx(k)e(k)]]>进一步扩展状态空间模型，形式如下：z(k+1)＝Amz(...

【技术特征摘要】

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