高炉顶压控制
正常工况下,高炉煤气经过通过透平机静叶产生能量供发电机组发电;TRT紧急停机时,通过快开、旁开阀进入煤气总管。此处仅对正常工况下高炉顶压控制展开研究。
设高炉顶压为P1,静叶流通面积为F,静叶开度为L,由工艺可设: F=f(L) , P1=f(F)
根据静叶阀门的直线结构特性,可得F与L的直线关系:
dF?Kf 积分得:F?KfL?k0 dL设高炉顶压煤气流量为Q1,通过静叶的煤气流量为Q2,由阀门的气体压差关系可得: Q2?
dp1RT?K(Q1?Q2) 由气体方程, ? ? n 则有 P ?dtV2AF?p??
dp12AF??Kdt?p1?p2??KQ1
所以
dp1p?p22A??K(KfL?k0)1?KQ1 dt??系统输入输出之间的关系:
随机模型,为SISO系统,将模型改写为:
z(k)?h(k)??v(k)
z(k)?h(k)?
z(k)为系统的输出向量,即高炉顶压,z(k)为输出估计向量,?为系统待辨识的参数向量,v(k)是均值为零的随机噪声。
选取1000组连续输入向量Li=[Li], i=1,2….1000, 代入判别方程,
?
?1000??110001J(?)?{?[z(k)?z(k)]2}?{?[z(k)?f(?)]2}2i?12i?1?最小二乘估计,求得使J(?)最小的参数估计向量?
Matlab仿真程序如下
?