发明名称 一种利用太赫兹时域谱精确测定食用油光学参数的方法
摘要 本发明公开了一种利用太赫兹时域谱精确测定食用油光学参数的方法。利用透射式太赫兹时域谱(THz-TDS)装置,测量没有放置样品时样品池的THz时域谱作为参考信号,再测量放置样品时样品池的THz时域谱作为样品信号;通过求样品信号和参考信号傅立叶变换的比值可得到样品在太赫兹波段传输系数的测量值,然后利用含容器在内的三层传输函数模型,使用Nelder-Mead寻优法和信赖域平滑法拟合样品在太赫兹波段传输系数的测量值,从而精确确定相应太赫兹波段食用油的折射率和吸收系数。本发明充分考虑了在多层介质中THz波传播时的多次反射效应,并考虑了噪声和测量误差影响,有效提高了食用油光学参数的测量精度,有利于将THz-TDS技术应用于食用油成份分析等定量检测场合。
申请公布号 CN101419157B 申请公布日期 2011.12.21
申请号 CN200810120801.6 申请日期 2008.09.05
申请人 中国计量学院 发明人 李九生;李向军
分类号 G01N21/17(2006.01)I;G01N33/03(2006.01)I 主分类号 G01N21/17(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 张法高
主权项 1.一种利用太赫兹时域谱精确测定食用油光学参数的方法,其特征在于包括如下步骤:1)在太赫兹时域谱透射测量装置上,测量没有放置样品时样品池的太赫兹时域谱作为参考信号,再测量放置样品时样品池的太赫兹时域谱作为样品信号;2)利用下式计算样品信号和参考信号傅立叶变换的比值,得到样品在太赫兹波段传输系数的测量值:<maths num="0001"><![CDATA[<math><mrow><msub><mi>H</mi><mi>exp</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>E</mi><mi>Sample</mi><mi>exp</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mi>E</mi><mi>Ref</mi><mi>exp</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中<img file="FSB00000558894700012.GIF" wi="192" he="59" />是样品信号的傅立叶变换和<img file="FSB00000558894700013.GIF" wi="153" he="59" />是参考信号的傅立叶变换;3)根据电磁理论,利用材料的光学参数复折射率<img file="FSB00000558894700014.GIF" wi="254" he="63" />n为实折射率,k为消光系数,建立包含容器在内的三层传输函数模型:<maths num="0002"><![CDATA[<math><mrow><msub><mi>H</mi><mi>Theory</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>E</mi><mi>Sample</mi><mi>th</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mi>E</mi><mi>Ref</mi><mi>th</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>这时<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>Ref</mi><mi>th</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&eta;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>E</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>T</mi><mn>10</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>P</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>F</mi><msubsup><mi>P</mi><mrow><mn>0</mn><mo>,</mo><mi>&delta;</mi></mrow><mn>13</mn></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>T</mi><mn>03</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>Sample</mi><mi>th</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&eta;</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>E</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>T</mi><mn>12</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>2</mn></msub><mo>)</mo></mrow><msubsup><mi>FP</mi><mrow><mn>2</mn><mo>,</mo><mi>&delta;</mi></mrow><mn>13</mn></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>2</mn></msub><mo>)</mo></mrow><msub><mi>T</mi><mn>23</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中H<sub>Theory</sub>(ω)为太赫兹波传输系数的理论值,<img file="FSB00000558894700018.GIF" wi="348" he="127" />和<img file="FSB00000558894700019.GIF" wi="338" he="122" />分别为入射介质A和出射介质B的反射和透射系数;<img file="FSB000005588947000110.GIF" wi="580" he="53" />为电磁波在介质A传播距离L<sub>A</sub>的传播因子,c为真空光速;<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>FP</mi><mrow><mi>B</mi><mo>,</mo><mi>&delta;</mi></mrow><mi>AC</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>&delta;</mi></munderover><msup><mrow><mo>{</mo><msub><mi>R</mi><mi>BC</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msup><msub><mi>P</mi><mi>B</mi></msub><mn>2</mn></msup><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mi>B</mi></msub><mo>)</mo></mrow><msub><mi>R</mi><mi>BA</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>}</mo></mrow><mi>k</mi></msup></mrow></math>]]></maths>为多次反射的FP效应,介质0为空气,介质1、3为SiO<sub>2</sub>容器壁,介质2是食用油,η<sub>x</sub>(ω)为THz波在SiO<sub>2</sub>容器壁的传输、传播和多重反射因子:<maths num="0006"><![CDATA[<math><mrow><msub><mi>&eta;</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>T</mi><mn>01</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><msub><mi>P</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>1</mn></msub><mo>)</mo></mrow><msubsup><mi>FP</mi><mrow><mn>1</mn><mo>,</mo><mi>&delta;</mi></mrow><mrow><mn>0</mn><mi>x</mi></mrow></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>P</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>3</mn></msub><mo>)</mo></mrow><msubsup><mi>FP</mi><mrow><mn>3</mn><mo>,</mo><mi>&delta;</mi></mrow><mrow><mi>x</mi><mn>0</mn></mrow></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><msub><mi>L</mi><mn>3</mn></msub><mo>)</mo></mrow><msub><mi>T</mi><mn>30</mn></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>4)利用含容器在内的三层传输函数模型,使用Nelder-Mead寻优法拟合样品在太赫兹波段传输系数的测量值来确定光学参数;5)通过多次测量得到二个以上样品信号和参考信号,然后根据测量误差理论确定噪声门限和测量重复性引起的统计偏差作为平滑算法的信赖域,平滑光学参数。
地址 310018 浙江省杭州市下沙高教园区学源街