[发明专利]热力站节能改造控制对象建模方法

权利要求书

1、一种热力站节能改造控制对象建模方法，其特征在于所述建模方法是按照以下步骤实现的：

步骤一、取得单位脉冲信号、单位阶跃信号和单位斜坡信号的响应模型：

供热过程采用二阶加滞后来描述，二阶加滞后对象的模型结构为G(s)=K(T1s+1)(T2s+1)e-τs,]]>其被辨识的参数为K，T₁，T₂和τ；其中s为拉普拉斯算子，单位为秒^-1；K为比例系数；T₁，T₂均为时间常数；τ为滞后时间常数；

以下分别计算二阶系统的单位脉冲信号、单位阶跃信号、单位斜坡信号的响应模型，并分别令其为A(t)，B(t)，C(t)；当单位脉冲信号的响应模型A(t)、单位阶跃信号的响应模型B(t)、单位斜坡信号的响应模型C(t)三者之一由实验获取，则另两个响应模型可通过式(1)求得：

B(t)=∫t0tA(ξ)dξC(t)=∫t0tB(ξ)dξ---(1)]]>

热力站供热对象的模型为G(s)=K(T1s+1)(T2s+1)e-τs]]>时，其中令T₂＝βT₁，β为比例系数，则单位脉冲信号的响应模型为：

A(t)=KT1(1-β)[exp(-t-τT1)-exp(-t-τβT1)]+ω(t)---(2)]]>

上式中ω(t)为噪声项，t为时间，以下各式同；

单位阶跃信号的响应模型为：

B(t)=∫t0tA(ξ)dξ=K+K1-β[βexp(-t-τβT1)-exp(-t-τT1)]+∫t0tω(ξ)dξ---(3)]]>

单位斜坡信号的响应模型为：

C(t)=∫t0tB(ξ)dξ]]>

=K(t-τ)-(β+1)KT1+KT11-β[exp(-t-τT1)-β2exp(-t-τβT1)]+∫∫ω(ξ)dξdσ---(4)]]>

步骤二、由典型传递函数取得辨识方程

由步骤一中所述的式(2)、式(3)和式(4)

可求得供热对象的二阶加滞后最小二乘辨识模型方程：

C(t)=t-1-A(t)-B(t)KKτT1T2T1+T2+Ω(t)---(5)]]>

其中Ω(t)为噪声项，是对ω(t)的滤波；

取不同时刻t₁，t₂，…，t_m，可得线性差分方程组

C＝Uθ+Ω                (6)

其中，C=C(1)C(2)···C(m),]]>U=1-1-A(1)-B(1)2-1-A(2)-B(2)············m-1-A(m)-B(m),]]>θ=KKτT1T2T1+T2,]]>Ω=Ω1Ω2···Ωm]]>

差分方程中θ未知，通过测量C和U来估计未知参数θ，即被辨识的参数K，T₁，T₂和τ，为了估计θ的n个参数，必须有m≥n；

步骤三、利用最小二乘法辨识二阶加滞后最小二乘辨识模型方程中的未知参数

由于存在模型误差和观测噪声，未知参数θ的估计要选择最小误差平方最小性能指标，定义误差向量平方J＝Ω^TΩ，由式(6)可得：

J＝(C-Uθ)^T(C-Uθ)

估计，使误差平方J最小，即使∂J∂θ=0,]]>得最小二乘参数估计值：

θ^=(UTU)-1UTC---(7)]]>

当未知参数估计值取得后，可得对象的传递函数模型：

G(s)=K(T1s+1)(T2s+1)e-τs---(8)]]>

步骤四、经离散化后取得供热对象的差分方程

将步骤三中得到的传递函数模型G(s)进行z变换，离散化后得到脉冲传递函数：

G(z)=Y(z)R(z)=bm′zm′+bm′-1zm′-1+···+b0an′zn′+an′-1zn′-1+···+a0z-k0---(9)]]>

其中z为z变换算子；Y(z)和R(z)分别为对象的输出与输入z变换；m’，n’分别为分子、分母的阶数，且m′≥0，n′≥0；a₀，a₁，a₂，…，a_n′与b₀，b₁，b₂，…，b_m′分别是分母、分子多项式系数；k₀为输出纯滞后步数，且k₀≥0；

将所述脉冲传递函数的分子、分母同除以zⁿ′并进行数学变换，再将所得的式子中的Y(z)替换为y(k)、z^-iY(z)替代为y(k-i)、z^-jR(z)替换为r(k-j)，最终得到对象的差分方程：

y(k)=[-Σi=0n′-1aiy(k-n′+i)+Σj=0m′bjr(k-n′+j-k0)]/an′---(10)]]>

其中k为第k个时间点。