一种聚烯烃过程的非线性预测函数控制方法技术

本发明专利技术属于工业自动化领域，公开了一种聚烯烃过程的非线性预测函数控制方法，首先根据采集到的非线性化工过程的阶跃响应数据对过程模型的时间常数以及滞后时间进行辨识，然后结合非线性过程的机理模型公式对增益进行公式推导，之后设计非线性化工过程的模型预测控制器，最后将得到的精确最优控制律实施于被控的非线性化工过程。本发明专利技术的技术方案是通过非线性化工过程阶跃响应数据采集、部分模型参数辨识、增益公式推导、非线性模型预测控制器设计等手段，建立基于机理模型的化工过程非线性模型预测控制方法，避免传统模型预测控制中由于模型线性化带来的预测误差，最终实现非线性化工过程的精准控制，进一步提升非线性化工过程生产的稳定性。程生产的稳定性。

【技术实现步骤摘要】
一种聚烯烃过程的非线性预测函数控制方法


[0001]本专利技术属于工业自动化领域，尤其涉及一种聚烯烃过程的非线性预测函数控制方法。

技术介绍

[0002]对工业中的聚烯烃过程而言，其分为动态模型与稳态模型，对动态模型而言，其时间常数以及滞后时间可以通过过程数据简单辨识得到，其稳态模型一般通过机理模型来进行表示。对这些过程来说，由稳态模型计算得到的过程变量间的模型增益大部分情况下是时变的，因此在预测函数控制技术的应用过程中，会由于增益参数的时变导致预测的不准确，最终影响聚烯烃过程的整体控制效果。如果能将这些时变的增益参数引入到预测函数控制中进行精确的预测推导，最终将从根本上解决传统预测函数控制在聚烯烃控制过程中存在预测偏差的问题，最终推动非线性预测函数控制技术在聚烯烃生产过程中的应用发展。

技术实现思路

[0003]本专利技术目的在于提供一种聚烯烃过程的非线性预测函数控制方法，以解决传统预测函数控制在无法有效解决聚烯烃生产过程中由于增益参数时变带来的预测偏差的技术问题。
[0004]为解决上述技术问题，本专利技术的一种聚烯烃过程的非线性预测函数控制方法的具体技术方案如下：
[0005]一种聚烯烃过程的非线性预测函数控制方法，包括如下步骤：步骤1：采集聚烯烃过程的阶跃响应数据，拟合出过程的动态模型参数；步骤1.1：给聚烯烃过程的输入端一个阶跃信号，并开始记录聚烯烃过程对应的阶跃响应数据；
[0006]步骤1.2：将获得的阶跃响应数据y
a
(k)转换成无纲量形式y<br/>a*
(k)；
[0007]步骤1.3：选取满足y
a*
(k1)＝0.39和y
a*
(k1)＝0.63的两个点，计算聚烯烃过程的时间常数和滞后时间；
[0008]步骤1.4：根据具体聚烯烃过程的工艺机理列出对应的机理模型，该机理模型也就被称为稳态模型，聚烯烃过程的实时增益从该机理模型中获得；步骤1.5：最终得到的聚烯烃过程传递函数模型；
[0009]步骤2：设计聚烯烃过程的非线性预测函数控制器；
[0010]步骤2.1：在采样时间T
s
以及零阶保持器下，聚烯烃过程的传递函数模型转换为离散方程；
[0011]步骤2.2：引入smith预估对离散方程中的时滞进行补偿，修正的无时滞的模型；
[0012]步骤2.3：选择预测函数控制的基函数为阶跃函数，得到基于无时滞模型的过程预测输出；
[0013]步骤2.4：预测函数控制器的设计选取目标函数；
[0014]步骤2.5：将预测输出式子进行变换；
[0015]步骤2.6：将得到的预测函数控制的最优控制律u(k)实施于聚烯烃过程，在下个采样周期按照步骤2.2～2.5中的步骤依次循环求解最新的最优控制律。
[0016]进一步地，所述步骤1.2将获得的阶跃响应数据y
a
(k)转换成无纲量形式y
a*
(k)：
[0017]y
a*
(k)＝y
a
(k)/y
s
[0018]其中，y
s
为阶跃测试中y
a
(k)的稳态值。
[0019]进一步地，所述步骤1.3依据下面的公式计算聚烯烃过程的时间常数和滞后时间：
[0020]T＝2(k2‑
k1)
[0021]τ＝2k1‑
k2[0022]其中，T和τ分别为聚烯烃过程模型的时间常数和滞后时间。
[0023]进一步地，所述步骤1.4的具体计算式子如下：
[0024]y(k)＝f(u(k))
[0025]K(k)＝f
′
(u(k))
[0026]其中，y(k)是聚烯烃过程的模型输出，u(k)是聚烯烃过程的输入，K(k)是聚烯烃过程的实时增益，f,f
′
分别为聚烯烃过程的机理模型以及对应的一阶导数。
[0027]进一步地，所述步骤1.5最终得到的聚烯烃过程传递函数模型为
[0028][0029]其中，G(k),s分别为聚烯烃过程的实时传递函数以及拉普拉斯算子。
[0030]进一步地，所述步骤2.1聚烯烃过程的传递函数模型转换为如下的离散方程：
[0031]y(k)＝αy(k
‑
1)+K(k
‑
1)(1
‑
α)u(k
‑1‑
d)
[0032]其中，进一步地，所述步骤2.2修正的无时滞的模型如下：
[0033]y
c
(k)＝αy
c
(k
‑
1)+K(k
‑
1)(1
‑
α)u(k
‑
1)
[0034]其中，y
c
(k)为修正后的无时滞的模型输出；
[0035]实际输出经过校正后为
[0036]y
ac
(k)＝y
a
(k)+y
c
(k)
‑
y
c
(k
‑
d)
[0037]其中，y
ac
(k)为修正后的实际输出。
[0038]进一步地，所述步骤2.3得到基于无时滞模型的过程预测输出如下：
[0039]y
c
(k+P)＝α
P
y
c
(k)+K(k)Bu(k)
[0040]其中，P为预测时域，
[0041][0042]进一步地，所述步骤2.4预测函数控制器的设计选取如下的目标函数：
[0043][0044]其中，y
r
(k)为对应的参考轨迹点，Q为跟踪误差的加权矩阵，e(k)为补偿的预测误差，e(k)＝y
ac
(k)
‑
y
c
(k)；
[0045]参考轨迹一般取如下式子：
[0046]y
r
(k+i)＝β
i
y
a
(k)+(1
‑
β
i
)c(k)
[0047]其中，β为参考轨迹柔化系数，c(k)为对应的设定值。
[0048]进一步地，所述步骤2.5将预测输出式子进行变换，得到如下形式：
[0049]y
c
(k+P)＝α
P
y
c
(k)+Bθ(k)
[0050]其中，θ(k)＝K(k)u(k)；
[0051]对上述目标函数求导，可得最优解为
[0052]θ(k)＝((β
i
‑
1)y
a
(k)+(1
‑
β
i
)c(k)
‑
α
P
y
...

【技术保护点】

【技术特征摘要】
1.一种聚烯烃过程的非线性预测函数控制方法，其特征在于，包括如下步骤：步骤1：采集聚烯烃过程的阶跃响应数据，拟合出过程的动态模型参数；步骤1.1：给聚烯烃过程的输入端一个阶跃信号，并开始记录聚烯烃过程对应的阶跃响应数据；步骤1.2：将获得的阶跃响应数据y
a
(k)转换成无纲量形式y
a*
(k)；步骤1.3：选取满足y
a*
(k1)＝0.39和y
a*
(k1)＝0.63的两个点，计算聚烯烃过程的时间常数和滞后时间；步骤1.4：根据具体聚烯烃过程的工艺机理列出对应的机理模型，该机理模型也就被称为稳态模型，聚烯烃过程的实时增益从该机理模型中获得；步骤1.5：最终得到的聚烯烃过程传递函数模型；步骤2：设计聚烯烃过程的非线性预测函数控制器；步骤2.1：在采样时间T
s
以及零阶保持器下，聚烯烃过程的传递函数模型转换为离散方程；步骤2.2：引入smith预估对离散方程中的时滞进行补偿，修正的无时滞的模型；步骤2.3：选择预测函数控制的基函数为阶跃函数，得到基于无时滞模型的过程预测输出；步骤2.4：预测函数控制器的设计选取目标函数；步骤2.5：将预测输出式子进行变换；步骤2.6：将得到的预测函数控制的最优控制律u(k)实施于聚烯烃过程，在下个采样周期按照步骤2.2～2.5中的步骤依次循环求解最新的最优控制律。2.根据权利要求1所述的聚烯烃过程的非线性预测函数控制方法，其特征在于，所述步骤1.2将获得的阶跃响应数据y
a
(k)转换成无纲量形式y
a*
(k)：y
a*
(k)＝y
a
(k)/y
s
其中，y
s
为阶跃测试中y
a
(k)的稳态值。3.根据权利要求1所述的聚烯烃过程的非线性预测函数控制方法，其特征在于，所述步骤1.3依据下面的公式计算聚烯烃过程的时间常数和滞后时间：T＝2(k2‑
k1)τ＝2k1‑
k2其中，T和τ分别为聚烯烃过程模型的时间常数和滞后时间。4.根据权利要求1所述的聚烯烃过程的非线性预测函数控制方法，其特征在于，所述步骤1.4的具体计算式子如下：y(k)＝f(u(k))K(k)＝f
′
(u(k))其中，y(k)是聚烯烃过程的模型输出，u(k)是聚烯烃过程的输入，K(k)是聚烯烃过程的实时增益，f,f
′
分别为聚烯烃过程的机理模型以及对应的一阶导数。5.根据权利要求1所述的聚烯烃过程的非线性预测函数控制方法，其特征在于，所述步骤1.5最终得到的聚烯烃过程传递函数模型为
其中，G(k),s分别为聚烯烃过程的实时传递函数以及拉普拉斯算子。6.根据权利要求1所述的聚烯烃过程的非线性预测函数控制方法，其特征在于，所述步骤2.1聚烯烃过程的传递函数模型转换为如下的离散方程：y(...

【专利技术属性】
技术研发人员：吴胜，柴俊沙，王元华，周鹏，
申请(专利权)人：杭州司南智能技术有限公司，
类型：发明
国别省市：