发明名称 基于单频激光辐照的双层参与性介质光谱辐射特性测量方法
摘要 基于单频激光辐照的双层参与性介质光谱辐射特性测量方法,属于参与性介质光学参数测量技术领域。本发明为了解决现有双层参与性介质光谱辐射特性的测量成本高及测量结果不准确的问题。它利用单频激光先后从两侧辐照双层参与性介质表面,利用探测器获得样品表面的频域复半球反射信号和复半球透射信号,最后利用反演的方法获得双层参与性介质的光谱吸收系数和光谱散射系数。本发明用于测量双层参与性介质光谱辐射特性。
申请公布号 CN103439283A 申请公布日期 2013.12.11
申请号 CN201310412796.7 申请日期 2013.09.11
申请人 哈尔滨工业大学 发明人 齐宏;任亚涛;张彪;孙双成;阮立明
分类号 G01N21/31(2006.01)I;G01N21/49(2006.01)I 主分类号 G01N21/31(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张利明
主权项 1.一种基于单频激光辐照的双层参与性介质光谱辐射特性测量方法,待测双层参与性介质由第一层介质和第二层介质组成,其特征在于,它包括以下步骤:步骤一:将频率为ω的单频调制激光束,从待测双层参与性介质的第一层介质表面侧垂直入射,然后采用探测器分别采集该第一层介质表面侧的一组频域复半球反射信号<img file="FDA0000380554300000011.GIF" wi="98" he="78" />和第二层介质表面侧的一组频域复半球透射信号<img file="FDA0000380554300000012.GIF" wi="230" he="78" />步骤二:由设定的第一层介质的光谱吸收系数κ<sub>a1</sub>和光谱散射系数κ<sub>s1</sub>的取值,根据辐射传输方程获得单频调制激光束计算域内的辐射强度场;然后执行步骤四;步骤三:由设定的第二层介质的光谱吸收系数κ<sub>a2</sub>和光谱散射系数κ<sub>s2</sub>的取值,根据辐射传输方程获得单频调制激光束计算域内的辐射强度场;然后执行步骤四;步骤四:根据当前计算获得的单频调制激光束计算域内的辐射强度场,当i=1时,通过下述公式:<maths num="0001"><![CDATA[<math><mrow><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>est</mi><mn>1</mn></msubsup><mo>=</mo><mo>-</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>I</mi><mi>in</mi></msub></mfrac><msub><mo>&Integral;</mo><mrow><mi>&mu;</mi><mo>&lt;</mo><mn>0</mn></mrow></msub><msub><mover><mi>I</mi><mo>^</mo></mover><mi>d</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>&mu;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>d&mu;</mi><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mover><mi>T</mi><mo>^</mo></mover><mi>est</mi><mn>1</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><msub><mi>I</mi><mi>in</mi></msub></mfrac><mo>[</mo><mn>2</mn><mi>&pi;</mi><msub><mo>&Integral;</mo><mrow><mi>&mu;</mi><mo>></mo><mn>0</mn></mrow></msub><msub><mover><mi>I</mi><mo>^</mo></mover><mi>d</mi></msub><mrow><mo>(</mo><mi>L</mi><mo>,</mo><mi>&mu;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>d&mu;</mi><mo>+</mo><msub><mover><mi>I</mi><mo>^</mo></mover><mi>c</mi></msub><mrow><mo>(</mo><mi>L</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>]</mo><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>计算获得第一层介质表面侧的一组频域复半球反射信号估计值<img file="FDA0000380554300000015.GIF" wi="70" he="78" />和第二层介质表面侧的一组频域复半球透射信号估计值<img file="FDA0000380554300000016.GIF" wi="98" he="78" />当i=2,通过下述公式:<maths num="0003"><![CDATA[<math><mrow><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>est</mi><mn>2</mn></msubsup><mo>=</mo><mo>-</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>I</mi><mi>in</mi></msub></mfrac><msub><mo>&Integral;</mo><mrow><mi>&mu;</mi><mo>&lt;</mo><mn>0</mn></mrow></msub><msub><mover><mi>I</mi><mo>^</mo></mover><mi>d</mi></msub><mrow><mo>(</mo><mi>L</mi><mo>,</mo><mi>&mu;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>d&mu;</mi><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mover><mi>T</mi><mo>^</mo></mover><mi>est</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><msub><mi>I</mi><mi>in</mi></msub></mfrac><mo>[</mo><mn>2</mn><mi>&pi;</mi><msub><mo>&Integral;</mo><mrow><mi>&mu;</mi><mo>></mo><mn>0</mn></mrow></msub><msub><mover><mi>I</mi><mo>^</mo></mover><mi>d</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>&mu;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>d&mu;</mi><mo>+</mo><msub><mover><mi>I</mi><mo>^</mo></mover><mi>c</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>]</mo><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>计算获得第二层介质表面侧的一组频域复半球反射信号估计值<img file="FDA0000380554300000019.GIF" wi="80" he="78" />和第一层介质表面侧的一组频域复半球透射信号估计值<img file="FDA00003805543000000110.GIF" wi="99" he="78" />式中I<sub>in</sub>为单频调制激光束的辐射强度,μ为待测双层参与性介质壁面外法向量与x轴正方向的夹角余弦,该x轴正方向为步骤一中单频调制激光束的入射方向,<img file="FDA00003805543000000111.GIF" wi="57" he="78" />为待测双层参与性介质的扩散辐射强度,L为待测双层参与性介质的厚度,<img file="FDA00003805543000000112.GIF" wi="49" he="78" />为单频调制激光束平行于入射方向的辐射强度;步骤五:计算获得逆问题算法中的目标函数值<img file="FDA0000380554300000021.GIF" wi="98" he="75" /><maths num="0005"><![CDATA[<math><mrow><msubsup><mi>F</mi><mi>obj</mi><mi>i</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><msup><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>est</mi><mi>i</mi></msubsup><mo>|</mo></mrow><mrow><mo>|</mo><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>mea</mi><mi>i</mi></msubsup><mo>|</mo></mrow></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><msubsup><mover><mi>T</mi><mo>^</mo></mover><mi>est</mi><mi>i</mi></msubsup><mo>|</mo></mrow><mrow><mo>|</mo><msubsup><mover><mi>T</mi><mo>^</mo></mover><mi>mea</mi><mi>i</mi></msubsup><mo>|</mo></mrow></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><msubsup><mi>&theta;</mi><mrow><mi>est</mi><mo>,</mo><mi>r</mi></mrow><mi>i</mi></msubsup><msubsup><mi>&theta;</mi><mrow><mi>mea</mi><mo>,</mo><mi>r</mi></mrow><mi>i</mi></msubsup></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><msubsup><mi>&theta;</mi><mrow><mi>est</mi><mo>,</mo><mi>t</mi></mrow><mi>i</mi></msubsup><msubsup><mi>&theta;</mi><mrow><mi>mea</mi><mo>,</mo><mi>t</mi></mrow><mi>i</mi></msubsup></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA0000380554300000023.GIF" wi="97" he="98" />为频域复半球反射信号估计值的幅值,<img file="FDA0000380554300000024.GIF" wi="94" he="98" />表示频域复半球透射信号估计值的幅值,<img file="FDA0000380554300000025.GIF" wi="112" he="98" />为采用探测器获得的频域复半球反射信号的幅值,<img file="FDA0000380554300000026.GIF" wi="110" he="98" />为采用探测器获得的频域复半球透射信号的幅值,<img file="FDA0000380554300000027.GIF" wi="93" he="78" />为频域复半球反射信号相位角的估计值,<img file="FDA0000380554300000028.GIF" wi="93" he="78" />为频域复半球透射信号相位角的估计值,<img file="FDA0000380554300000029.GIF" wi="109" he="78" />为采用探测器获得的频域复半球反射信号相位角,<img file="FDA00003805543000000210.GIF" wi="112" he="76" />为采用探测器获得的频域复半球透射信号相位角;步骤六:判断步骤五中获得的目标函数值<img file="FDA00003805543000000211.GIF" wi="77" he="78" />是否小于预设置的阈值ξ,若是,当i=1时,执行步骤七,当i=2时,执行步骤九;否则,当i=1时,根据逆问题算法重新设定第一层介质的光谱吸收系数κ<sub>a1</sub>和光谱散射系数κ<sub>s1</sub>的取值,返回步骤二,当i=2时,根据逆问题算法重新设定第二层介质的光谱吸收系数κ<sub>a2</sub>和光谱散射系数κ<sub>s2</sub>的取值,返回步骤三;步骤七:将设定的当前第一层介质的光谱吸收系数κ<sub>a1</sub>和光谱散射系数κ<sub>s1</sub>的取值作为待测双层参与性介质第一层介质的测量光谱吸收系数κ<sub>a1</sub>和光谱散射系数κ<sub>s1</sub>的值;步骤八:再将频率为ω的单频调制激光束,从待测双层参与性介质的第二层介质表面侧垂直入射,再采用探测器分别采集该第二层介质表面侧的一组频域复半球反射信号<img file="FDA00003805543000000212.GIF" wi="92" he="78" />和第一层介质表面侧的一组频域复半球透射信号<img file="FDA00003805543000000213.GIF" wi="227" he="74" />返回执行步骤三;步骤九:将设定的当前第二层介质的光谱吸收系数κ<sub>a2</sub>和光谱散射系数κ<sub>s2</sub>的取值作为待测双层参与性介质第二层介质的测量光谱吸收系数κ<sub>a2</sub>和光谱散射系数κ<sub>s2</sub>的值;由此实现双层参与性介质光谱辐射特性的测量。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号