发明名称 一种分布反馈式激光器时频响应曲线的确定方法
摘要 本发明公开了一种分布反馈式激光器时频响应曲线的确定方法,相比于现有技术在采用吸收光谱法测量气体浓度时需要预先使用光学干涉仪来获取DFB激光器在扫描电压调制下的时频响应曲线,本发明激光器时频响应曲线的确定方法在不使用干涉仪的条件下,对峰值归一化的吸收线型进行最小二乘拟合,也能够快速准确获得激光器扫描过程的时频响应曲线,并利用该时频响应曲线实现气体浓度的准确测量,从而有效降低了测量成本。
申请公布号 CN106053021A 申请公布日期 2016.10.26
申请号 CN201610328995.3 申请日期 2016.05.18
申请人 东南大学 发明人 周宾;李可;许康;程禾尧;陆勇
分类号 G01M11/02(2006.01)I;G01N21/31(2006.01)I;G01N21/39(2006.01)I 主分类号 G01M11/02(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 李倩
主权项 一种分布反馈式激光器时频响应曲线的确定方法,其特征在于,包括如下步骤:步骤1,搭建气体浓度测量系统;步骤2,向气体吸收池内通入一定浓度的待测气体,使用DFB激光器作为光源,调节DFB激光器的工作温度,使加载扫描电压调制后的DFB激光器出射光频率覆盖所需目标气体吸收峰,出射光经准直后穿过气体吸收池并由探测器接收,对测量到的原始信号进行基线拟合,获得时域上的吸光度曲线测量值α<sub>m</sub>;步骤3,在设定工作温度条件下,DFB激光器在给定的调制电压范围内时频响应曲线基本保持不变,由经验公式(1)表示:v(t)=a+b·t+c·t<sup>2</sup>+d·t<sup>3</sup>   (1);式中,v(t)表示激光器出射光在时刻t时的频率;步骤4,不同频率v处的吸光度α(v)用式(2)表示:<maths num="0001"><math><![CDATA[<mrow><mi>&alpha;</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>A&phi;</mi><mi>C</mi></msub><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>A</mi><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mfrac><mrow><msub><mi>&Delta;v</mi><mi>C</mi></msub></mrow><mrow><msup><mrow><mo>(</mo><mi>v</mi><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&Delta;v</mi><mi>C</mi></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000993108820000011.GIF" wi="1140" he="183" /></maths>式中,A为待测气体的积分吸光度,v<sub>0</sub>为吸收峰的中心频率,Δv<sub>C</sub>为碰撞加宽;其中,积分吸光度A用式(3)表示:A=PXLS(T)   (3);式中,P为气体压力,X为气体浓度,L为光程长度,S(T)为谱线在温度T下的线强;当v=v<sub>0</sub>时,可得峰值处吸光度α(v<sub>0</sub>)=2A/(π·Δv<sub>C</sub>);此时峰值归一化吸光度α′可表示为:<maths num="0002"><math><![CDATA[<mrow><msup><mi>&alpha;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>(</mo><mfrac><mrow><mi>v</mi><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub></mrow><mrow><msub><mi>&Delta;v</mi><mi>C</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000993108820000012.GIF" wi="974" he="191" /></maths>对于混合气体,吸收线型的碰撞加宽Δv<sub>C</sub>由式(5)计算得到:<maths num="0003"><math><![CDATA[<mrow><msub><mi>&Delta;v</mi><mi>C</mi></msub><mo>=</mo><mi>P</mi><munder><mo>&Sigma;</mo><mi>j</mi></munder><msub><mi>X</mi><mi>j</mi></msub><mn>2</mn><msub><mi>&gamma;</mi><mi>j</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000993108820000013.GIF" wi="942" he="116" /></maths>式中,X<sub>j</sub>为组分j的浓度,2γ<sub>j</sub>为组分j的碰撞加宽系数;步骤5,将式(1)代入式(4),可得:<maths num="0004"><math><![CDATA[<mrow><msup><mi>&alpha;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>v</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>(</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub></mrow><mrow><msub><mi>&Delta;v</mi><mi>C</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>(</mo><mfrac><mrow><mi>a</mi><mo>+</mo><mi>b</mi><mo>&CenterDot;</mo><mi>t</mi><mo>+</mo><mi>c</mi><mo>&CenterDot;</mo><msup><mi>t</mi><mn>2</mn></msup><mo>+</mo><mi>d</mi><mo>&CenterDot;</mo><msup><mi>t</mi><mn>3</mn></msup><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub></mrow><mrow><msub><mi>&Delta;v</mi><mi>C</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000993108820000021.GIF" wi="1381" he="199" /></maths>对测量得到的低浓度气体的吸光度曲线α<sub>m</sub>进行峰值归一化处理,得到时域上的测量值α’<sub>m</sub>,设置目标函数:<maths num="0005"><math><![CDATA[<mrow><mi>O</mi><mi>B</mi><mi>J</mi><mo>=</mo><mo>|</mo><mo>|</mo><mfrac><mrow><msup><mi>&alpha;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>v</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msub><msup><mi>&alpha;</mi><mo>&prime;</mo></msup><mi>m</mi></msub></mrow><mrow><msup><mi>&alpha;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>v</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000993108820000022.GIF" wi="997" he="147" /></maths>利用最小二乘算法拟合多项式系数a、b、c、d,当OBJ<ε时,ε取10<sup>‑5</sup>,即可获得激光器在该扫描电压下的时频响应曲线。
地址 210096 江苏省南京市玄武区四牌楼2号