主权项 |
一种基于LSSVM及在线更新的电站锅炉烟气软测量方法,其特征在于,所述方法采用LSSVM构建初始烟气排放模型,然后利用样本追加和样本替换来实现模型的增量更新,能够对锅炉烟气各成分进行精确地预测,所述方法包括以下步骤:步骤1:收集初始训练样本来构建LSSVM模型;步骤2:利用此模型对烟气成分含量进行预测;步骤3:当实际的烟气成分含量的传感器测量值y<sub>q</sub>采集到后,计算样本(x<sub>q</sub>,y<sub>q</sub>)的预测误差Er;步骤4:判断预测误差:若Er>Δ,Δ为误差阈值,则执行步骤5,否则执行步骤9;步骤5:从历史运行数据中选取距新采样样本(x<sub>q</sub>,y<sub>q</sub>)最近的样本点(x<sub>k</sub>,y<sub>k</sub>),其中<img file="FDA0000992235660000011.GIF" wi="527" he="107" />步骤6:根据新采样样本(x<sub>q</sub>,y<sub>q</sub>)与其最近的样本(x<sub>k</sub>,y<sub>k</sub>)之间的距离进行判断并确定更新类型;步骤7:根据确定的更新类型,对步骤1获得的初始LSSVM模型进行增量更新,根据步骤6进行样本增加或样本替换来更新;步骤8:所述步骤1中LSSVM模型可描述为以下优化问题,<maths num="0001"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mi>w</mi><mo>,</mo><mi>b</mi><mo>,</mo><mi>ξ</mi></mrow></munder><mi>J</mi><mrow><mo>(</mo><mrow><mi>w</mi><mo>,</mo><mi>ξ</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>γ</mi><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><msub><mi>ξ</mi><mi>i</mi></msub><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000992235660000012.GIF" wi="1825" he="162" /></maths><img file="FDA0000992235660000013.GIF" wi="1201" he="113" />其中,J为目标函数,<img file="FDA0000992235660000015.GIF" wi="92" he="55" />是核空间映射函数,w为权重向量,γ为惩罚系数,ξ<sub>i</sub>为误差变量,b为偏差,上标T表示矩阵的转置;为解此优化问题,定义Lagrange函数如下:<img file="FDA0000992235660000014.GIF" wi="1847" he="154" />其中,α=[α<sub>1</sub>,…,α<sub>n</sub>]<sup>T</sup>为Lagrange乘子;利用Lagrange函数对各变量求偏导,并令导数值为零可得到:<img file="FDA0000992235660000021.GIF" wi="948" he="617" />消去中间变量w和ξ<sub>i</sub>,将其转化为求解线性方程组:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd><mtd><msup><mover><mi>I</mi><mo>→</mo></mover><mi>T</mi></msup></mtd></mtr><mtr><mtd><mover><mi>I</mi><mo>→</mo></mover></mtd><mtd><mrow><mi>Ω</mi><mo>+</mo><mn>1</mn><mo>/</mo><mi>γ</mi><mi>I</mi></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>b</mi></mtd></mtr><mtr><mtd><mi>a</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>y</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000992235660000022.GIF" wi="1327" he="166" /></maths>其中y=[y<sub>1</sub>,…,y<sub>n</sub>]<sup>T</sup>,<img file="FDA0000992235660000025.GIF" wi="283" he="73" />I为n×n阶单位矩阵,Ω={Ω<sub>kl</sub>|k,l=1,…,n},且<img file="FDA0000992235660000023.GIF" wi="774" he="98" />定义为核函数;通过求解方程组得到α和b的值为:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>b</mi><mo>=</mo><mfrac><mrow><msup><mover><mi>I</mi><mo>→</mo></mover><mi>T</mi></msup><msup><mi>H</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>y</mi></mrow><mrow><mover><mi>I</mi><mo>→</mo></mover><msup><mi>H</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mover><mi>I</mi><mo>→</mo></mover></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mi>a</mi><mo>=</mo><msup><mi>H</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><mover><mi>I</mi><mo>→</mo></mover><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000992235660000024.GIF" wi="1927" he="322" /></maths>用求得的新的模型的特征矩阵H<sup>‑1</sup>计算得到相应的模型参数α和b,实现对烟气软测量模型h(x)的更新;步骤9:判断测试样本是否结束,若结束则停止程序;否则执行步骤2,对下一样本进行预测。 |