CN103610227A - 一种烘丝机头尾段工艺变量优化控制方法 - Google Patents

Abstract

本发明公开了一种烘丝机头尾段工艺变量优化控制方法，依据烘丝过程头尾段筒温、风温、排潮风门等工艺变量的历史数据，采用三次函数作为径向基函数的Cubic-RBF-ARX模型对烘丝动态特性进行建模；所建模型具有自调节能力，能反映不同模式下的入口流量以及入口水分的变化对出口水分的影响，可根据头尾段不同模式的入口流量及入口水分的变化来预测未来出口水分的变化情况；根据所建模型对各工艺变量进行优化设定，可使头尾段叶丝出口水分的控制达到较好的效果。本发明方法综合考虑了来料量与各输入变量间的动态特性，可以更有效地克服来料流量和水分变化对烘丝过程头尾段的影响，适用于不同模式下叶丝入口流量与入口水分时的头尾段控制。

Description

一种烘丝机头尾段工艺变量优化控制方法

技术领域

本发明涉及烘丝机头尾段工艺变量优化控制方法。

背景技术

烘丝过程是香烟制丝生产中最重要的一道加工工序，它主要是通过对叶丝进行加热干燥，降低叶丝的含水率，使烘烤后叶丝的含水率、温度均匀一致，并控制在一定的数值范围内，以满足生产工艺要求。烘丝的工艺流程主要分为预热、干头、中间以及干尾过程四个部分。在干头阶段，叶丝入口流量不断增加，但无叶丝出口水分的检测值，难以进行反馈控制，容易造成干头阶段出口水分控制品质差、干料多；在干尾阶段，由于叶丝入口流量骤然减少，而烘丝筒具有较大热容，筒壁内部温度难以按规定的速率下降等问题，也容易造成干尾阶段出口水分控制性能低且干料多。因此，“干头干尾”是目前烘丝过程出口水分控制的难点所在。

现有的干头干尾过程控制方法主要有：

（1）利用进入和输出烘丝机的物料和干燥介质作为热质平衡对象建立数学模型，结合前馈PID调节筒温的控制方式。但前馈数学模型仅考虑了进料的含水率和流量，并没有考虑热风温度等其他对出口水分有重要影响的因素，不能完全反应真实过程，造成头尾段烘丝机出口水分波动大，需要操作人员进行人工干预，对于头尾段不同模式下的、具有不同入口流量和入口水分的来料难以获得满意的控制效果。

（2）在上述前馈控制的基础上，在头尾段增加蒸喷加湿装置对头尾料施加蒸汽水来提高头尾料的含水率，以降低干料量。但此方法仅对叶丝表层进行加湿，仅提高了叶丝表层湿度，仍然会造成烟丝内在质量的降低，且增加了出口水分控制的难度与稳定性。

（3）通过多次试验、寻求最佳头尾阶段热风温度值和调整排潮阀门开度等工艺参数来减少干料量。此方法缺乏自调节能力，无法保证对于不同模式下具有不同入口流量和入口水分的来料时，该组工艺参数均为最优设定值；

（4）在PID控制策略的基础上，将模糊控制的思想应用到烘丝机水分控制中。

仅仅用单纯的二维模糊控制器来解决烘丝过程出口水分的控制问题仍然无法获得最优的工艺参数设定值，而且对于不同模式下的入口流量与入口水分的变化，还需对模糊控制规则表进行调整，这对工业生产带来不便。

发明内容

本发明所要解决的技术问题是，针对现有技术不足，提供一种烘丝机头尾段工艺变量优化控制方法，使干头阶段叶丝出口水分尽可能快地上升、并快速到达稳定状态，使干尾阶段叶丝出口水分尽可能缓慢地下降，从而有效地减少头尾段的干料量，提高烘丝过程的控制性能；更有效地克服来料流量和水分变化对烘丝过程头尾段的影响，避免人工整定输入工艺变量参数的不便。

为解决上述技术问题，本发明所采用的技术方案是：一种烘丝机头尾段工艺变量优化控制方法，该方法为：

1）根据烘丝机的运行流程，建立烘丝过程中叶丝入口流量、入口水分、筒温、风温、排潮风门、出口水分的时序关系，同时根据烘丝过程干头阶段无叶丝出口水分检测值、干尾阶段无叶丝入口流量与入口水分检测值的特点，采用三次函数作为径向基函数的Cubic-RBF-ARX模型，分别建立烘丝过程干头阶段与干尾阶段的Cubic-RBF-ARX模型；

2）根据烘丝机头尾段的历史运行数据，采用结构化非线性参数优化方法分别优化烘丝过程干头阶段与干尾阶段的Cubic-RBF-ARX模型；

3）依据优化的烘丝过程干头阶段与干尾阶段的Cubic-RBF-ARX模型，采用双S型函数描述干头阶段的排潮风门、风温、筒温的最优输入曲线；采用阶跃函数描述干头阶段的入口流量的最优输入曲线；采用指数函数描述干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线；

4）采用列维布格奈奎尔特方法，通过使优化的干头阶段与干尾阶段的Cubic-RBF-ARX模型计算出的出口水分预测值与出口水分设定值的误差最小，寻找出烘丝过程干头阶段与干尾阶段的最优输入曲线的参数，以适应来料情况的变化，减少干尾阶段的干料量。

所述步骤1）中，烘丝机干头阶段Cubic-RBF-ARX模型为：

其中：

其中，y^H(t^H)表示烘丝机干头阶段Cubic-RBF-ARX模型的出口水分；

分别表示干头阶段Cubic-RBF-ARX模型的排潮风门开度、风温、筒温、入口流量及入口水分；X^H(t^H-1)为入口流量和入口水分的状态变量；np^H，nq^H，d^H和m^H均表示干头阶段Cubic-RBF-ARX模型的阶次；

分别为干头阶段Cubic-RBF-ARX模型输出项与输入项的RBF神经网络的中心；

为干头阶段Cubic-RBF-ARX模型的标量权系数；||·||_F表示矩阵的Frobenius范数；ξ^H(t^H)是干头阶段Cubic-RBF-ARX模型的建模误差，为高斯白噪声；T₀ ^H为烘丝机干头阶段Cubic-RBF-ARX模型建模采样时间，T₁为从有入口流量检测值到有入口水分检测值的时间，T₂为从有入口水分检测值到有出口水分检测值的时间，T₃为从有入口水分检测值到烘丝筒入口的时间，T₄为叶丝在烘丝筒烘干的时间。

所述步骤1）中，烘丝机干尾阶段Cubic-RBF-ARX模型为：

其中：

其中，y^T(t^T)表示烘丝机干尾阶段Cubic-RBF-ARX模型的出口水分；

分别表示干尾阶段Cubic-RBF-ARX模型的筒温、热风风温、排潮风门开度、入口流量、入口水分及筒体电机频率；X^T(t^T-1)为热风风温和筒体电机频率的状态变量；np^T,nq^T,d^T和m^T均表示干尾阶段Cubic-RBF-ARX模型的阶次；

分别为干尾阶段Cubic-RBF-ARX模型输出项与输入项的RBF神经网络的中心；

为干尾阶段Cubic-RBF-ARX模型的标量权系数；ξ^T(t^T)是干尾阶段Cubic-RBF-ARX模型建模误差，为高斯白噪声；T₀ ^T为烘丝机干尾阶段Cubic-RBF-ARX模型建模采样时间。

所述步骤2）中，烘丝机干头阶段Cubic-RBF-ARX模型优化如下：

$({\widehat{\mathrm{\θ}}}_N^H,{\widehat{\mathrm{\θ}}}_L^H)=\arg \underset{{\mathrm{\θ}}_N^H,{\mathrm{\θ}}_L^H}{\min }\underset{t^{\mathrm{oh}}=1}{\overset{N^H}{\mathrm{\Σ}}}{({\stackrel{\‾}{y}}^H\left(t^{\mathrm{oh}}\right)-{\hat{y}}^H\left(t^{\mathrm{oh}}\right))}^2$

其中，是烘丝机干头阶段出口水分的实际值，

是在实际输入作用下，由烘丝机干头阶段Cubic-RBF-ARX模型计算出的出口水分的预测值； ${\widehat{\mathrm{\θ}}}_L^H=\{{\mathrm{\ω}}_0^{H,0},{\mathrm{\ω}}_{i^H,0}^{y^H},{\mathrm{\ω}}_{n,j^H,0}^{u^H},{\mathrm{\ω}}_{k^H}^{H,0},{\mathrm{\ω}}_{i^H,k^H}^{y^H},{\mathrm{\ω}}_{j^H,k^H}^{u^H}|i^H=1,...,{\mathrm{np}}^H;j^H=1,...,{\mathrm{nq}}^H;k^H=1,...,m^H\}$ 为烘丝机干头阶段Cubic-RBF-ARX模型的线性参数；

为烘丝机干头阶段Cubic-RBF-ARX模型的非线性参数；N^H为烘丝机干头阶段Cubic-RBF-ARX模型建模数据长度。

烘丝机干尾阶段Cubic-RBF-ARX模型优化如下：

$({\widehat{\mathrm{\θ}}}_N^T,{\widehat{\mathrm{\θ}}}_L^T)=\arg \underset{{\mathrm{\θ}}_N^T,{\mathrm{\θ}}_L^T}{\min }\underset{t^{\mathrm{ot}}=1}{\overset{N^T}{\mathrm{\Σ}}}{({\stackrel{\‾}{y}}^T\left(t^{\mathrm{ot}}\right)-{\hat{y}}^T\left(t^{\mathrm{ot}}\right))}^2$

其中，

是烘丝机干尾过程中出口水分的实际值；是在实际输入作用下，由烘丝机干尾阶段Cubic-RBF-ARX模型计算出的出口水分的预测值； ${\mathrm{\θ}}_L^T=\{{\mathrm{\ω}}_0^{T,0},{\mathrm{\ω}}_{i^T,0}^{y^T},{\mathrm{\ω}}_{n,j^T,0}^{u^T},{\mathrm{\ω}}_{k^T}^{T,0},{\mathrm{\ω}}_{i^T,k^T}^{y^T},{\mathrm{\ω}}_{j^T,k^T}^{u^T}|i^T=1,...,{\mathrm{np}}^T;j^T=1,...,{\mathrm{nq}}^T;k^T=1,...,m^T\}$ 为烘丝机干尾阶段Cubic-RBF-ARX模型的线性参数，

为烘丝机干尾阶段Cubic-RBF-ARX模型的非线性参数；N^T为烘丝机干头阶段Cubic-RBF-ARX模型建模数据长度。

所述步骤3）中：

用于描述烘丝机干头阶段排潮风门、风温、筒温的最优输入曲线的双S型函数表达式为：

$U_{\mathrm{sc}}\left(t^s\right)=\frac{{\mathrm{\λ}}_1}{1+e^{\frac{t^s-{\mathrm{\λ}}_2}{{\mathrm{\λ}}_3}}}+{\mathrm{\λ}}_4+\frac{{\mathrm{\λ}}_5}{1+e^{\frac{t^s-{\mathrm{\λ}}_6}{{\mathrm{\λ}}_7}}}$

其中，t^s为输入的时间，单位为s；λ₁,λ₄,λ₅分别为双S型函数的起点、转折点及终点值；λ₂,λ₆分别为双S型函数的两条对称轴中心位置；λ₃,λ₇分别为双S型函数上升或下降的速度；λ₃,λ₇大于0时表示S型函数上升，λ₃,λ₇小于0时表示S型函数下降；c=1,2,3，U_s1(t^s)是排潮风门的设定值；U_s2(t^s)是风温的设定值；U_s3(t^s)是筒温的设定值。

用于描述烘丝机干头阶段入口流量的最优输入曲线的阶跃函数表达式为：

$U_T\left(t^T\right)=\left\{\begin{array}{cc}\frac{{\mathrm{\κ}}_1{t}^T}{{\mathrm{\κ}}_2}& t^T\∈[1,{\mathrm{\κ}}_2]\\ {\mathrm{\κ}}_1& t^T\∈[{\mathrm{\κ}}_2+1,{\mathrm{\κ}}_3]\end{array}\right.;$

其中，t^T为输入的时间，单位为s；κ₁,κ₂,κ₃分别为阶跃函数的上升速度、上升时间与终值。

所述步骤4）中，烘丝机干头阶段Cubic-RBF-ARX模型计算出的出口水分预测值

为： ${\stackrel{\‾}{y}}^H\left(t^a\right)=f(U_{s1}\left(t^a\right),U_{s2}\left(t^a\right),U_{s3}\left(t^a\right),U_T\left(t^a\right)),$

通过将烘丝机干头阶段各工艺变量的优化设定曲线代入所构建的干头阶段Cubic-RBF-ARX的输入变量

中得到。通过使干头阶段Cubic-RBF-ARX模型计算出的出口水分预测值

与出口水分设定值y_set(t^a)的误差e_H(t^a)最小，即采用列维布格奈奎尔特方法求解优化问题

寻找出干头阶段排潮风门、风温、筒温的输入曲线的参数λ_x和入口流量输入曲线的参数κ₁,κ₂,κ₃；其中，x=1,2,…,7；g=1,2,3；M是干头阶段持续的时间。

所述步骤3）中，用于描述干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线的指数函数的表达式为：

$\begin{array}{c}{U}_{\mathrm{zp}}\left(t^z\right)={\mathrm{\α}}_{p1}\×{\left({\mathrm{\α}}_{p2}\right)}^{t^z}+{\mathrm{\α}}_{p3}\\ p=\mathrm{1,2,3,4}\end{array};$

式中U_z1(t^z)、U_z2(t^z)、U_z3(t^z)、U_z4(t^z)分别表示干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线。

所述步骤4）中，干尾阶段的Cubic-RBF-ARX模型计算出的出口水分预测值

为： ${\stackrel{\‾}{y}}^T\left(t^b\right)=f(U_{z1}\left(t^b\right),U_{z2}\left(t^b\right),U_{z3}\left(t^b\right),U_{z4}\left(t^b\right)),$

通过将烘丝机干尾阶段各工艺变量的优化设定曲线代入所构建的干尾阶段Cubic-RBF-ARX模型的输入变量中得到；通过使干尾阶段的Cubic-RBF-ARX模型计算出的出口水分预测值

与出口水分设定值y'_set(t^b)的误差e_T(t^b)最小，即采用列维布格奈奎尔特方法求解优化问题寻找出干尾阶段排潮风门、风温、筒温和筒体电机频率最优输入曲线的参数α_pg；其中，g=1,2,3；M'是干尾阶段持续时间。

与现有技术相比，本发明所具有的有益效果为：本发明方法可使干头阶段叶丝出口水分尽可能快地上升、并快速到达稳定状态，可使干尾阶段叶丝出口水分尽可能缓慢地下降，从而有效地减少头尾段的干料量，提高烘丝过程的控制性能，具有较大的经济价值；本发明方法综合考虑了来料量与各输入变量间的动态特性，可以更有效地克服来料流量和水分变化对烘丝过程头尾段的影响，适用于不同模式下叶丝入口流量与入口水分时的头尾段控制；本发明方法基于辨识的模型优化出最优的输入设定曲线，避免了人工整定输入工艺变量参数的不便。

附图说明

图1为烘丝机工艺过程示意图。

具体实施方式

烘丝机工艺过程如图1所示。叶丝进入烘丝工序之前，首先检测叶丝的入口流量u₄和入口水分u₅。经过T₃时间，叶丝到达烘丝机入口处。叶丝在烘丝筒烘干时，***会定时采样筒体的排潮风门开度u₁、风温u₂、筒温u₃等工艺变量参数值。烘干过程持续T₄时间，烘干后的叶丝从烘丝筒出口倒出，并在出口处测量叶丝出口水分值y。从有入口流量检测值到有出口水分检测值需经历一段较长时间，例如某烘丝生产线大约需340s。另外，烘丝机的输入/输出变量间也具有较大的时滞。

当检测到有入口流量时，表明烘丝过程开始运行。运行初期烘丝过程有叶丝入口流量与入口水分检测值，没有叶丝出口水分检测值，此时烘丝过程干头阶段开始。根据烘丝过程干头阶段的特性，建立Cubic-RBF-ARX模型结构：

其中：

其中，y^H(t^H)表示烘丝机干头阶段Cubic-RBF-ARX模型的出口水分；

分别表示干头阶段Cubic-RBF-ARX模型的排潮风门开度、风温、筒温、入口流量及入口水分；X^H(t^H-1)为入口流量和入口水分的状态变量；np^H，nq^H，d^H和m^H均表示干头阶段Cubic-RBF-ARX模型的阶次；

分别为干头阶段Cubic-RBF-ARX模型输出项与输入项的RBF神经网络的中心；

为干头阶段Cubic-RBF-ARX模型的标量权系数；||·||_F表示矩阵的Frobenius范数；ξ^H(t^H)是干头阶段Cubic-RBF-ARX模型的建模误差，为高斯白噪声；T₀ ^H为烘丝机干头阶段Cubic-RBF-ARX模型建模采样时间，T₁为从有入口流量检测值到有入口水分检测值的时间，T₂为从有入口水分检测值到有出口水分检测值的时间，T₃为从有入口水分检测值到烘丝筒入口的时间，T₄为叶丝在烘丝筒烘干的时间。

当入口流量由正常值变为0时，标志着干尾过程的开始，当出口水分下降到3%时，标志着烘丝机整个烘丝过程的结束。干尾过程中无入口流量检测值，但有出口水分检测值。根据烘丝机干尾过程段的特性，建立如下的Cubic-RBF-ARX模型：

其中：

其中，y^T(t^T)表示烘丝机干尾阶段Cubic-RBF-ARX模型的出口水分；

分别表示干尾阶段Cubic-RBF-ARX模型的筒温、热风风温、排潮风门开度、入口流量、入口水分及筒体电机频率；X^T(t^T-1)为热风风温和筒体电机频率的状态变量；np^T,nq^T,d^T和m^T均表示干尾阶段Cubic-RBF-ARX模型的阶次；

分别为干尾阶段Cubic-RBF-ARX模型输出项与输入项的RBF神经网络的中心；

为干尾阶段Cubic-RBF-ARX模型的标量权系数；ξ^T(t^T)是干尾阶段Cubic-RBF-ARX模型建模误差，为高斯白噪声；T₀ ^T为烘丝机干尾阶段Cubic-RBF-ARX模型建模采样时间。

本发明采用结构化非线性参数优化方法（SNPOM）方法对模型进行估计。为了使得上面所构造的Cubic-RBF-ARX模型能够描述烘丝过程头尾段的全局动态特性，我们首先采用SNPOM方法来优化模型的、一步预测误差最小情形下的参数，并以此参数作为长期预测优化目标下的模型参数初始值。然后，采用列维布格奈奎尔特方法（LMM）来进行长期预测性能最优的模型参数的优化。

烘丝机干头阶段Cubic-RBF-ARX模型（1）的参数优化问题如下：

$({\widehat{\mathrm{\θ}}}_N^H,{\widehat{\mathrm{\θ}}}_L^H)=\arg \underset{{\mathrm{\θ}}_N^H,{\mathrm{\θ}}_L^H}{\min }\underset{t^{\mathrm{oh}}=1}{\overset{N^H}{\mathrm{\Σ}}}{({\stackrel{\‾}{y}}^H\left(t^{\mathrm{oh}}\right)-{\hat{y}}^H\left(t^{\mathrm{oh}}\right))}^2---\left(5\right)$

其中，

是烘丝机干头阶段出口水分的实际值，

是在实际输入作用下，由烘丝机干头阶段Cubic-RBF-ARX模型计算出的出口水分的预测值； ${\widehat{\mathrm{\θ}}}_L^H=\{{\mathrm{\ω}}_0^{H,0},{\mathrm{\ω}}_{i^H,0}^{y^H},{\mathrm{\ω}}_{n,j^H,0}^{u^H},{\mathrm{\ω}}_{k^H}^{H,0},{\mathrm{\ω}}_{i^H,k^H}^{y^H},{\mathrm{\ω}}_{j^H,k^H}^{u^H}|i^H=1,...,{\mathrm{np}}^H;j^H=1,...,{\mathrm{nq}}^H;k^H=1,...,m^H\}$ 为烘丝机干头阶段Cubic-RBF-ARX模型的线性参数；

为烘丝机干头阶段Cubic-RBF-ARX模型的非线性参数；N^H为烘丝机干头阶段Cubic-RBF-ARX模型建模数据长度。

烘丝机干尾阶段Cubic-RBF-ARX模型（3）的参数优化问题如下：

$({\widehat{\mathrm{\θ}}}_N^T,{\widehat{\mathrm{\θ}}}_L^T)=\arg \underset{{\mathrm{\θ}}_N^T,{\mathrm{\θ}}_L^T}{\min }\underset{t^{\mathrm{ot}}=1}{\overset{N^T}{\mathrm{\Σ}}}{({\stackrel{\‾}{y}}^T\left(t^{\mathrm{ot}}\right)-{\hat{y}}^T\left(t^{\mathrm{ot}}\right))}^2---\left(6\right)$

其中，

是烘丝机干尾过程中出口水分的实际值；

是在实际输入作用下，由烘丝机干尾阶段Cubic-RBF-ARX模型计算出的出口水分的预测值； ${\mathrm{\θ}}_L^T=\{{\mathrm{\ω}}_0^{T,0},{\mathrm{\ω}}_{i^T,0}^{y^T},{\mathrm{\ω}}_{n,j^T,0}^{u^T},{\mathrm{\ω}}_{k^T}^{T,0},{\mathrm{\ω}}_{i^T,k^T}^{y^T},{\mathrm{\ω}}_{j^T,k^T}^{u^T}|i^T=1,...,{\mathrm{np}}^T;j^T=1,...,{\mathrm{nq}}^T;k^T=1,...,m^T\}$ 为烘丝机干尾阶段Cubic-RBF-ARX模型的线性参数，

为烘丝机干尾阶段Cubic-RBF-ARX模型的非线性参数；N^T为烘丝机干头阶段Cubic-RBF-ARX模型建模数据长度。

依据估计出的烘丝过程干头阶段Cubic-RBF-ARX模型来设计各工艺变量的最优输入曲线，以适应来料情况的变化，尽量减少干头阶段的干料量。本发明采用双S型函数来描述干头阶段排潮风门、风温、筒温的最优输入曲线，采用阶跃型函数来描述入口流量的最优输入曲线。

双S型曲线公式如下：

$U_{\mathrm{sc}}\left(t^s\right)=\frac{{\mathrm{\λ}}_1}{1+e^{\frac{t^s-{\mathrm{\λ}}_2}{{\mathrm{\λ}}_3}}}+{\mathrm{\λ}}_4+\frac{{\mathrm{\λ}}_5}{1+e^{\frac{t^s-{\mathrm{\λ}}_6}{{\mathrm{\λ}}_7}}}---\left(7\right)$

其中，t^s为输入的时间，单位为s；λ₁,λ₄,λ₅分别为双S型函数的起点、转折点及终点值；λ₂,λ₆分别为双S型函数的两条对称轴中心位置；λ₃,λ₇分别为双S型函数上升或下降的速度；λ₃,λ₇大于0时表示S型函数上升，λ₃,λ₇小于0时表示S型函数下降；c=1,2,3，U_s1(t^s)是排潮风门的设定值；U_s2(t^s)是风温的设定值；U_s3(t^s)是筒温的设定值。

描述入口流量输入曲线的阶跃型函数公式如下：

$U_T\left(t^T\right)=\left\{\begin{array}{cc}\frac{{\mathrm{\κ}}_1{t}^T}{{\mathrm{\κ}}_2}& t^T\∈[1,{\mathrm{\κ}}_2]\\ {\mathrm{\κ}}_1& t^T\∈[{\mathrm{\κ}}_2+1,{\mathrm{\κ}}_3]\end{array}\right.---\left(8\right)$

其中，t^T为输入的时间，单位为s；κ₁,κ₂,κ₃分别为阶跃函数的上升速度、上升时间与终值。

将各工艺变量的优化设定曲线（7-8）代入所构建的Cubic-RBF-ARX模型（1）的输入变量

中，可得到干头阶段出口水分的预测值

${\stackrel{\‾}{y}}^H\left(t^a\right)=f(U_{s1}\left(t^a\right),U_{s2}\left(t^a\right),U_{s3}\left(t^a\right),U_T\left(t^a\right))---\left(9\right)$

采用列维布格奈奎尔特（Levenberg-Marquardt Method，LMM）方法，通过使模型计算出的出口水分预测值与出口水分设定值的误差最小，寻找出干头阶段排潮风门、风温、筒温最优输入曲线的λ_i(i=1,2,…,7)参数和入口流量最优输入曲线的κ_j(j=1,2,3)参数。干头阶段出口水分设定值与基于干头动态模型预测值（9）之间的误差为：

$e_H\left(t^a\right)=y_{\mathrm{set}}\left(t^a\right)-{\stackrel{\‾}{y}}^H\left(t^a\right)---\left(10\right)$

y_set(t^a)是出口水分设定值。

干头阶段工艺变量最优设定的优化问题如下：

$\underset{{\mathrm{\λ}}_x,{\mathrm{\κ}}_g}{\min }J=\underset{t^a=1}{\overset{M}{\mathrm{\Σ}}}{e}_H^2\left(t^a\right)---\left(11\right)$

M是干头阶段持续时间。通过求解上述优化问题可得到最优设定曲线的参数值，从而设计出烘丝机干头阶段各个工艺变量的最优输入曲线。

依据估计出的烘丝过程干尾阶段Cubic-RBF-ARX模型来设计各工艺变量的最优输入曲线，以适应来料情况的变化，尽量减少干尾阶段的干料量。采用指数型函数来描述干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线，该指数型曲线公式如下所示：

$\begin{array}{c}{U}_{\mathrm{zp}}\left(t^z\right)={\mathrm{\α}}_{p1}\×{\left({\mathrm{\α}}_{p2}\right)}^{t^z}+{\mathrm{\α}}_{p3}\\ p=\mathrm{1,2,3,4}\end{array}---\left(12\right)$

式中U_z1(t^z)、U_z2(t^z)、U_z3(t^z)、U_z4(t^z)分别表示干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线。将各工艺变量的优化设定曲线（12）代入所构建的Cubic-RBF-ARX模型（3）的输入变量

中，可得到干尾阶段出口水分的预测值：

${\stackrel{\‾}{y}}^T\left(t^b\right)=f(U_{z1}\left(t^b\right),U_{z2}\left(t^b\right),U_{z3}\left(t^b\right),U_{z4}\left(t^b\right))---\left(13\right)$

采用列维布格奈奎尔特（LMM）方法，通过使模型计算出的出口水分预测值与出口水分设定值的误差最小，寻找出干尾阶段排潮风门、风温、筒温和筒体电机频率最优输入曲线的α_pg；其中，g=1,2,3。干尾阶段出口水分设定值与基于干尾动态模型预测值（13）之间的误差为：

$e_T\left(t^b\right)={y^{\′}}_{\mathrm{set}}\left(t^b\right)-{\stackrel{\‾}{y}}^T\left(t^b\right)---\left(14\right)$

y_set(t)是出口水分设定值。

干尾阶段工艺变量最优设定的优化问题如下：

$\underset{{\mathrm{\α}}_{\mathrm{pg}}}{\min }{J}^{\′}=\underset{k=1}{\overset{M^{\′}}{\mathrm{\Σ}}}{e}_T^2\left(t^b\right)---\left(15\right)$

M是干尾阶段持续时间。通过求解上述优化问题可得到最优设定曲线的参数值，从而设计出烘丝机干尾阶段各个工艺变量的最优输入曲线。

Claims (8)

1.一种烘丝机头尾段工艺变量优化控制方法，其特征在于，该方法为：

1）根据烘丝机的运行流程，建立烘丝过程中叶丝入口流量、入口水分、筒温、风温、排潮风门、出口水分的时序关系，同时根据烘丝过程干头阶段无叶丝出口水分检测值、干尾阶段无叶丝入口流量与入口水分检测值的特点，采用三次函数作为径向基函数的Cubic-RBF-ARX模型，分别建立烘丝过程干头阶段与干尾阶段的Cubic-RBF-ARX模型；

2）根据烘丝机头尾段的历史运行数据，采用结构化非线性参数优化方法分别优化烘丝过程干头阶段与干尾阶段的Cubic-RBF-ARX模型；

3）依据优化的烘丝过程干头阶段与干尾阶段的Cubic-RBF-ARX模型，采用双S型函数描述干头阶段的排潮风门、风温、筒温的最优输入曲线；采用阶跃函数描述干头阶段的入口流量的最优输入曲线；采用指数函数描述干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线；

4）采用列维布格奈奎尔特方法，通过使优化的干头阶段与干尾阶段的Cubic-RBF-ARX模型计算出的出口水分预测值与出口水分设定值的误差最小，寻找出烘丝过程干头阶段与干尾阶段的最优输入曲线的参数，以适应来料情况的变化，减少干尾阶段的干料量。

2.根据权利要求1所述的烘丝机头尾段工艺变量优化控制方法，其特征在于，所述步骤1）中，烘丝机干头阶段Cubic-RBF-ARX模型为：

其中：

其中，y^H(t^H)表示烘丝机干头阶段Cubic-RBF-ARX模型的出口水分；

分别表示干头阶段Cubic-RBF-ARX模型的排潮风门开度、风温、筒温、入口流量及入口水分；X^H(t^H-1)为入口流量和入口水分的状态变量；np^H，nq^H，d^H和m^H均表示干头阶段Cubic-RBF-ARX模型的阶次；

分别为干头阶段Cubic-RBF-ARX模型输出项与输入项的RBF神经网络的中心；

为干头阶段Cubic-RBF-ARX模型的标量权系数；||·||_F表示矩阵的Frobenius范数；ξ^H(t^H)是干头阶段Cubic-RBF-ARX模型的建模误差，为高斯白噪声；T₀ ^H为烘丝机干头阶段Cubic-RBF-ARX模型建模采样时间，T₁为从有入口流量检测值到有入口水分检测值的时间，T₂为从有入口水分检测值到有出口水分检测值的时间，T₃为从有入口水分检测值到烘丝筒入口的时间，T₄为叶丝在烘丝筒烘干的时间。

3.根据权利要求1或2所述的烘丝机头尾段工艺变量优化控制方法，其特征在于，所述步骤1）中，烘丝机干尾阶段Cubic-RBF-ARX模型为：

其中：

其中，y^T(t^T)表示烘丝机干尾阶段Cubic-RBF-ARX模型的出口水分；

分别表示干尾阶段Cubic-RBF-ARX模型的筒温、热风风温、排潮风门开度、入口流量、入口水分及筒体电机频率；X^T(t^T-1)为热风风温和筒体电机频率的状态变量；np^T,nq^T,d^T和m^T均表示干尾阶段Cubic-RBF-ARX模型的阶次； 分别为干尾阶段Cubic-RBF-ARX模型输出项与输入项的RBF神经网络的中心；

为干尾阶段Cubic-RBF-ARX模型的标量权系数；ξ^T(t^T)是干尾阶段Cubic-RBF-ARX模型建模误差，为高斯白噪声；T₀ ^T为烘丝机干尾阶段Cubic-RBF-ARX模型建模采样时间。

4.根据权利要求3所述的烘丝机头尾段工艺变量优化控制方法，其特征在于，所述步骤2）中，烘丝机干头阶段Cubic-RBF-ARX模型优化如下：

$({\widehat{\mathrm{\θ}}}_N^H,{\widehat{\mathrm{\θ}}}_L^H)=\arg \underset{{\mathrm{\θ}}_N^H,{\mathrm{\θ}}_L^H}{\min }\underset{t^{\mathrm{oh}}=1}{\overset{N^H}{\mathrm{\Σ}}}{({\stackrel{\‾}{y}}^H\left(t^{\mathrm{oh}}\right)-{\hat{y}}^H\left(t^{\mathrm{oh}}\right))}^2$

其中，

是烘丝机干头阶段出口水分的实际值，

是在实际输入作用下，由烘丝机干头阶段Cubic-RBF-ARX模型计算出的出口水分的预测值； ${\widehat{\mathrm{\θ}}}_L^H=\{{\mathrm{\ω}}_0^{H,0},{\mathrm{\ω}}_{i^H,0}^{y^H},{\mathrm{\ω}}_{n,j^H,0}^{u^H},{\mathrm{\ω}}_{k^H}^{H,0},{\mathrm{\ω}}_{i^H,k^H}^{y^H},{\mathrm{\ω}}_{j^H,k^H}^{u^H}|i^H=1,...,{\mathrm{np}}^H;j^H=1,...,{\mathrm{nq}}^H;k^H=1,...,m^H\}$ 为烘丝机干头阶段Cubic-RBF-ARX模型的线性参数；

为烘丝机干头阶段Cubic-RBF-ARX模型的非线性参数；N^H为烘丝机干头阶段Cubic-RBF-ARX模型建模数据长度。

烘丝机干尾阶段Cubic-RBF-ARX模型优化如下：

$({\widehat{\mathrm{\θ}}}_N^T,{\widehat{\mathrm{\θ}}}_L^T)=\arg \underset{{\mathrm{\θ}}_N^T,{\mathrm{\θ}}_L^T}{\min }\underset{t^{\mathrm{ot}}=1}{\overset{N^T}{\mathrm{\Σ}}}{({\stackrel{\‾}{y}}^T\left(t^{\mathrm{ot}}\right)-{\hat{y}}^T\left(t^{\mathrm{ot}}\right))}^2$

其中，

是烘丝机干尾过程中出口水分的实际值；是在实际输入作用下，由烘丝机干尾阶段Cubic-RBF-ARX模型计算出的出口水分的预测值； ${\mathrm{\θ}}_L^T=\{{\mathrm{\ω}}_0^{T,0},{\mathrm{\ω}}_{i^T,0}^{y^T},{\mathrm{\ω}}_{n,j^T,0}^{u^T},{\mathrm{\ω}}_{k^T}^{T,0},{\mathrm{\ω}}_{i^T,k^T}^{y^T},{\mathrm{\ω}}_{j^T,k^T}^{u^T}|i^T=1,...,{\mathrm{np}}^T;j^T=1,...,{\mathrm{nq}}^T;k^T=1,...,m^T\}$ 为烘丝机干尾阶段Cubic-RBF-ARX模型的线性参数，

为烘丝机干尾阶段Cubic-RBF-ARX模型的非线性参数。N^T为烘丝机干头阶段Cubic-RBF-ARX模型建模数据长度。

5.根据权利要求4所述的烘丝机头尾段工艺变量优化控制方法，其特征在于，所述步骤3）中：

用于描述烘丝机干头阶段排潮风门、风温、筒温的最优输入曲线的双S型函数表达式为：

$U_{\mathrm{sc}}\left(t^s\right)=\frac{{\mathrm{\λ}}_1}{1+e^{\frac{t^s-{\mathrm{\λ}}_2}{{\mathrm{\λ}}_3}}}+{\mathrm{\λ}}_4+\frac{{\mathrm{\λ}}_5}{1+e^{\frac{t^s-{\mathrm{\λ}}_6}{{\mathrm{\λ}}_7}}}$

其中，t^s为输入的时间，单位为s；λ₁,λ₄,λ₅分别为双S型函数的起点、转折点及终点值；λ₂,λ₆分别为双S型函数的两条对称轴中心位置；λ₃,λ₇分别为双S型函数上升或下降的速度；λ₃,λ₇大于0时表示S型函数上升，λ₃,λ₇小于0时表示S型函数下降；c=1,2,3，U_s1(t^s)是排潮风门的设定值；U_s2(t^s)是风温的设定值；U_s3(t^s)是筒温的设定值；

用于描述烘丝机干头阶段入口流量的最优输入曲线的阶跃函数表达式为：

$U_T\left(t^T\right)=\left\{\begin{array}{cc}\frac{{\mathrm{\κ}}_1{t}^T}{{\mathrm{\κ}}_2}& t^T\∈[1,{\mathrm{\κ}}_2]\\ {\mathrm{\κ}}_1& t^T\∈[{\mathrm{\κ}}_2+1,{\mathrm{\κ}}_3]\end{array}\right.;$

其中，t^T为输入的时间，单位为s；κ₁,κ₂,κ₃分别为阶跃函数的上升速度、上升时间与终值。

6.根据权利要求5所述的烘丝机头尾段工艺变量优化控制方法，其特征在于，所述步骤4）中，将烘丝机干头阶段优化设定曲线U_sc(t^s)、U_T(t^T)代入所述烘丝机干头阶段Cubic-RBF-ARX模型的输入变量

中，得到烘丝机干头阶段Cubic-RBF-ARX模型计算出的出口水分预测值

通过使干头阶段Cubic-RBF-ARX模型计算出的出口水分预测值

与出口水分设定值y_set(t^a)的误差e_H(t^a)最小，即采用列维布格奈奎尔特方法求解优化问题

寻找出干头阶段排潮风门、风温、筒温的输入曲线的参数λ_x和入口流量输入曲线的参数κ₁,κ₂,κ₃；其中，x=1,2,…,7；g=1,2,3；M是干头阶段持续的时间。

7.根据权利要求6所述的烘丝机头尾段工艺变量优化控制方法，其特征在于，所述步骤3）中，用于描述干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线的指数函数的表达式为：

$\begin{array}{c}{U}_{\mathrm{zp}}\left(t^z\right)={\mathrm{\α}}_{p1}\×{\left({\mathrm{\α}}_{p2}\right)}^{t^z}+{\mathrm{\α}}_{p3}\\ p=\mathrm{1,2,3,4}\end{array};$

式中U_z1(t^z)、U_z2(t^z)、U_z3(t^z)、U_z4(t^z)分别表示干尾阶段排潮风门、风温、筒温和筒体电机频率的最优输入曲线。

8.根据权利要求7所述的烘丝机头尾段工艺变量优化控制方法，其特征在于，所述步骤4）中，将烘丝机干尾阶段优化设定曲线U_zp(t^z)代入所述烘丝机干尾阶段Cubic-RBF-ARX模型的输入变量

中，得到烘丝机干尾阶段的Cubic-RBF-ARX模型计算出的出口水分预测值

通过使干尾阶段的Cubic-RBF-ARX模型计算出的出口水分预测值与出口水分设定值y'_set(t^b)的误差e_T(t^b)最小，即采用列维布格奈奎尔特方法求解优化问题

寻找出干尾阶段排潮风门、风温、筒温和筒体电机频率最优输入曲线的参数α_pg；其中，g=1,2,3；M'是干尾阶段持续时间。

Images