发明名称 结合卡尔曼滤波器的腔衰荡光谱技术气体浓度监测方法
摘要 本发明属于气体检测领域,具体是一种结合卡尔曼滤波器的腔衰荡光谱技术气体浓度监测方法。一种结合卡尔曼滤波器的腔衰荡光谱技术气体浓度监测方法,由采用基于腔衰荡光谱技术的气体浓度监测装置采集数据的方法以及对采集到的数据采用卡尔曼滤波器进行分析的方法组成;所述基于腔衰荡光谱技术的气体浓度监测装置包括一个激光发射装置,激光发射装置的出射端通过光纤连接有一个声光调制器,声光调制器的出射光路上顺次设有耦合透镜以及腔体上设有进气口与出气口的光学腔,光学腔的出射光路上设有第一光电探测器。本发明结合了卡尔曼滤波器,对腔衰荡信号进行后续处理,从而减小了系统误差,提高了腔衰荡光谱的准确度。
申请公布号 CN103884679B 申请公布日期 2016.01.20
申请号 CN201410156271.6 申请日期 2014.04.18
申请人 山西大学 发明人 赵刚;马维光;李志新;谭巍;尹王保;贾锁堂
分类号 G01N21/39(2006.01)I 主分类号 G01N21/39(2006.01)I
代理机构 太原科卫专利事务所(普通合伙) 14100 代理人 朱源
主权项 一种结合卡尔曼滤波器的腔衰荡光谱技术气体浓度监测方法,其特征在于,由采用基于腔衰荡光谱技术的气体浓度监测装置采集数据的方法以及对采集到的数据采用卡尔曼滤波器进行分析的方法组成;所述基于腔衰荡光谱技术的气体浓度监测装置包括一个激光发射装置,激光发射装置的出射端通过光纤连接有一个声光调制器(3),声光调制器(3)的出射光路上顺次设有耦合透镜(4)以及腔体上设有进气口与出气口的光学腔(5),所述光学腔(5)位于光路上的两端均为高反镜且两个高反镜的反射面相对,其中一个高反镜固定在压电陶瓷上,压电陶瓷顺次连接有一个高压放大器(8)和函数发生器(9);光学腔(5)的出射光路上设有第一光电探测器(6),第一光电探测器(6)的信号输出端分为两路,一路连接有计算机(15),另一路通过一个阈值电路(7)与声光调制器(3)的信号输入端相连接;所述阈值电路(7)包括一个定时器以及一个比较器;所述比较器的同向输入端与第一光电探测器(6)的信号输出端相连接,比较器的反向输入端连接有一个电位器R5,比较器的信号输出端与定时器的第六管脚相连接;所述定时器的第三管脚与声光调制器(3)的信号输入端相连接;比较器的第七管脚以及定时器的第八管脚均与外部电源相连接;定时器的第三管脚还通过电阻R3连接有电容C2,定时器的第二管脚与电容C2相连接;电容C2的另一端接地;定时器的第四管脚与外部电源相连接,定时器的第七管脚通过电阻R1与外部电源相连接;定时器的第一管脚接地;所述比较器的第四管脚接地,定时器的第六管脚通过电阻R2接地;电位器R5的一端接地,另一端通过电阻R4与外部电源相连接;电位器R5两端并联有稳压管D1;所述定时器采用LM555CN,比较器采用TLV3501;所述计算机(15)在内设的基于labview平台的CRDS采集及拟合程序软件的支持下,采用卡尔曼滤波器分析方法对由上述装置采集并得到的气体浓度数据进行分析,进而得出优化后待测气体的浓度;所述卡尔曼滤波器分析方法包括以下步骤:设k‑1时刻真实的气体浓度为C<sub>k‑1</sub>,浓度的变化为w<sub>k‑1</sub>,则k时刻真实的浓度值可表示为:C<sub>k</sub>=C<sub>k‑1</sub>+w<sub>k‑1</sub>通过实验系统测量得到浓度值为z<sub>k</sub>,测量误差为v<sub>k</sub>,则有:z<sub>k</sub>=C<sub>k</sub>+v<sub>k</sub>其中,w<sub>k</sub>和v<sub>k</sub>为高斯白噪声,二者方差分别为σ<sub>w</sub><sup>2</sup>和σ<sub>v</sub><sup>2</sup>;卡尔曼滤波过程是一个递归过程,包括5个公式:预测过程:在k‑1时刻估计k时刻的值,浓度不会发生突变,则假设k时刻的浓度值等于k‑1时刻得到的最佳估计值:C<sub>k|k‑1</sub>=C<sub>k‑1|k‑1</sub>计算先验误差,先验误差用来表征使用k‑1时刻值估算k时刻值的估计误差,它等于k‑1时刻的后验误差加浓度C的变化方差:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>w</mi><mn>2</mn></msubsup></mrow>]]></math><img file="FDA0000800900270000021.GIF" wi="516" he="124" /></maths>更新过程:更新k时刻的卡尔曼增益K<sub>g</sub>,K<sub>g</sub>用来权衡估计值与测量值的可信程度<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>K</mi><mi>g</mi></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>v</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow>]]></math><img file="FDA0000800900270000031.GIF" wi="654" he="133" /></maths>更新k时刻的最佳估计值,Z<sub>k</sub>为k时刻实际测量值C<sub>k|k</sub>=C<sub>k|k‑1</sub>+K<sub>g</sub>(z<sub>k</sub>‑C<sub>k|k‑1</sub>)计算k时刻的后验误差,后验误差用来表征最佳估计值的估计误差P<sub>k|k</sub>=(1‑K<sub>g</sub>)P<sub>k|k‑1</sub>最终得到的C<sub>k|k</sub>就是经过卡尔曼滤波方法滤波后的待测气体浓度数据;其中浓度变化误差方差与测量误差方差σ<sub>w</sub><sup>2</sup>、σ<sub>v</sub><sup>2</sup>可以根据下式进行估计:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&sigma;</mi><mi>&nu;</mi><mn>2</mn></msubsup><mo>=</mo><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mi>z</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>*</mo><mi>z</mi><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000800900270000032.GIF" wi="807" he="188" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>&sigma;</mi><mi>w</mi><mn>2</mn></msubsup><mo>=</mo><mi>R</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>-</mo><msubsup><mi>&sigma;</mi><mi>v</mi><mn>2</mn></msubsup><mo>.</mo></mrow>]]></math><img file="FDA0000800900270000033.GIF" wi="409" he="94" /></maths>
地址 030006 山西省太原市坞城路92号