[发明专利]印染车间空气温度控制方法

权利要求书

1.印染车间空气温度控制方法，其特征在于该方法的具体步骤是：

步骤1.建立被控对象的多预测模型，具体是：

首先，以冷水二通阀、加热二通阀和蒸汽二通阀的阀门开度值为输入控制量，以温度传感器采集到的染色车间空气温度值为输出量，通过系统辨识方法，建立染色过程在各个典型温度下控制系统的离散时间传递函数模型

A_i(z^-1)y(k)＝B_i(z^-1)u(k)

其中y(k)表示染色车间空气温度值；u(k)＝[u₁(k),u₂(k),u₃(k)]^T为k时刻的控制输入变量，其中u₁(k)表示k时刻冷水二通阀阀门开度值，u₂(k)表示k时刻加热二通阀阀门开度值，u₃(k)表示k时刻蒸汽二通阀阀门开度值；在染色过程中，通过控制冷水二通阀、加热二通阀和蒸汽二通阀的阀门开度值来实现染色车间空气温度控制；A_i(z^-1)和B_i(z^-1)表示第i个典型工况对应的多项式矩阵，可以通过辨识得到，其形式为

Ai(z-1)=I+Σj=1naAijz-i,Bi(z-1)=Σj=0nbBijz-i]]>

其中A_ij,B_ij表示需要辨识的模型参数，I表示具有合适维数的单位矩阵，n_a、n_b表示采样个数；

然后，将前述离散时间传递函数模型通过状态空间实现，建立下述基于状态空间描述的印染车间空气温度控制系统模型

x(k+1)＝A_ix(k)+B_iΔu(k)

y(k)＝C_ix(k)

其中x(k)＝[y^T(k)…y^T(k-n_a)Δu^T(k-1)…Δu^T(k-n_b)]^T为系统k时刻的状态向量，Δu(k)＝u(k)-u(k-1)为k时刻控制输入的增量；上述模型中的系数矩阵A_i,B_i和C_i分别为

Bi=[Bi0T0...0|I0...0]]]>

Ci=[I0...0|00...0]]]>

步骤2.预测输出值，具体是：由第1步建立的被控对象的状态空间模型，通过迭代计算，可得输出值

y(k+j)=CiAijx^(k)+Σl=0j-1CiAij-l-1Δu(k+l)]]>

其中为状态向量x(k)的估计值，y(k+j)表示利用k时刻的数据估计k+j时刻的输出值，j为不小于1的正整数；

令N和N_u分别表示预测时域和控制时域，则可将预测输出的状态空间模型表示为如下形式

Y(k)=Fix^(k)+GiU(k)]]>

其中

Y(k)=y^(k+1|k)y^(k+2|k)...y^(k+N|k),U(k)=Δu(k)Δu(k+1)...Δu(k+Nu-1)]]>

式中为k+s时刻输出量的预测值，且s＝1,2,...,N；

由上式知，对输出量的预测依赖于系统状态向量的估计值状态估计的准确性及收敛速度将直接影响输出值y的预测；

步骤3.切换系统的有限时间域状态估计，具体是：通过设计切换系统的有限时间域收敛状态观测器，以解决多模型切换系统状态估计的精度和收敛快速性问题，可用于印染车间空气温度的优化控制；

对于可观的矩阵对(A_i,C_i)，设计两个观测器

z_i(k+1)＝(A_i-L_jiC_i)z_i(k)+B_iu(k)+L_jiy(k),j＝1,2,i＝1,2,…n

其中zi=z1iz2i,]]>z_1i(k)和z_2i(k)分别为k时刻两个观测器的状态向量，L_1i,L_2i为待设计的观测器增益矩阵；

引入滞后时间d，d为正整数，并记M_i＝A_i-L_jiC_i，则

z_i(k+1)＝M_iz_i(k)+N_iu(k)+L_jiy(k)

x^(k)=K[zi(k)-Midzi(k-d)]]]>

其中

Mi=M1i00M2i=Ai-LjiCi00Ai-LjiCi]]>

Ni=BiBi,Lji=L1iL2i,T=II,Ki=[I0][TMidT]-1]]>

选择合适的矩阵L_1i,L_2i，使得矩阵M_i稳定，则有

若初始状态为x(0)＝x₀，及可得

zi(k)=Tx(k)+Mid[Tx^0-Tx0]x(k)=Kizi(k)-KiMidzi(k-d)^=Ki[Tx(k)+MidT(x^0-x0)]-KiMidTx^0]]>

显然K_iT＝I和所以当k=d时，下式成立

x^(k)-x(k)=KiMikT(x^0-x0)=0]]>

对于染色过程温度控制系统，选择绝对值大于零且小于1的实数标量利用标准的极点配置方法，设计矩阵L_1i，将M_1i的极点配置到单位圆内的则存在非奇异矩阵Γ，将M_1i做相似变换，得

M‾1i=Γ-1(Ai-LjiCi)Γ=diag{λ1,λ2,...,λna+nb+1}]]>

式中diag{λ1,λ2,...,λna+nb+1}]]>表示元素为λ1,λ2,...,λna+nb+1]]>的对角矩阵；

类似的，再设计矩阵L_2i，使A_i-L_2iC_i＝α(A_i-L_1iC_i)，其中标量α满足0<α<1，则可将M_2i的极点配置为单位圆内的且有

M‾2i=Γ-1(Ai-L2iCi)Γ=αdiag{λ1,λ2,...,λna+nb+1}]]>

M2i=ΓM‾2iΓ-1=αdiag{λ1,λ2,...,λna+nb+1}=αM‾1i]]>

上式表明，对任意有限大小的正整数d，矩阵[TMidT]]]>的行列式为det[TMidT]=det[M2id-M1id]≠0,]]>即矩阵[TMidT]]]>可逆，切换系统状态向量x(k)的估计误差将在任意的有限时间步长k=d时收敛至零；

步骤4.优化控制求解；

（1）建立切换系统二次型优化目标函数

Ji=Σj=1NuPiΔu2(k+j-1)+Σj=1NQiy^2(k+j|k)]]>

其中P_i和Q_i分别为多模型切花系统中控制量和预测输出的加权矩阵；

当取Q_i＝I时，前述优化目标函数可写为

Ji=[GiUi(k)-Fix^]T[GiU(k)-Fix^(k)]+UT(k)PiU(k)]]>

（2）切换系统优化控制问题的求解

利用QP算法，当矩阵Q_i＝I时，通过求解可得最优控制输入为

U(k)=-(GiTGi+Pi)-1GiTFix^(k)]]>

其中I表示维数合适的单位矩阵。