发明名称 线性调频双光束激光外差测量入射角度的方法
摘要 线性调频双光束激光外差测量入射角度的方法,涉及工程领域。本发明是为了解决现有的角度测量方法需要动镜的方向信息,使得信号的解调和调制复杂,并且测量精度低的问题。本发明所述的该方法是基于下列装置实现的,所述装置包括线性调频激光器、第一反射镜、第二反射镜、会聚透镜、平面标准镜、光电探测器和信号处理系统,在测量过程中,信号处理系统接收光电探测器发出的光电流I,并对光电流I进行处理,采用傅里叶变换获得该光电流I中的外差信号频率f<sub>IF</sub>,由折射定律获得折射角θ和待测激光入射角θ0之间的关系,进而获得待测的激光入射角θ<sub>0</sub>。它可用于激光雷达、机械、仪器仪表和电子产品制造业中。
申请公布号 CN103954266A 申请公布日期 2014.07.30
申请号 CN201410205959.9 申请日期 2014.05.15
申请人 黑龙江大学 发明人 李彦超;冉玲苓;杨九如;高扬;柳春郁;杨瑞海;杜军;丁群;王春晖;马立峰;于伟波
分类号 G01C1/00(2006.01)I 主分类号 G01C1/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 岳泉清
主权项 线性调频双光束激光外差测量入射角度的方法,该方法是基于下列装置实现的,所述装置包括线性调频激光器(1)、第一反射镜(2)、第二反射镜(3)、会聚透镜(5)、平面标准镜(4)、光电探测器(6)和信号处理系统(7),所述线性调频激光器(1)发出的线性调频偏振光经第一反射镜(2)反射后入射到第二反射镜(3),经第二反射镜(3)反射后的光束以入射角θ<sub>0</sub>斜入射到平面标准镜(4)的前表面,平面标准镜(4)前表面将该激光分束成一号反射光和折射光,平面标准镜(4)前表面将该折射光折射至平面标准镜(4)的后表面,该折射光经该平面标准镜(4)的后表面与前表面多次反射后获得多束二号反射光,该多束二号反射光经该平面标准镜(4)的前表面透射之后,与一号反射光通过会聚透镜(5)汇聚至光电探测器(6)的光敏面上,所述光电探测器(6)输出电信号给信号处理系统(7),其特征在于,基于上述装置的线性调频双光束激光外差测量角度的方法为:线性调频激光器(1)以入射角θ<sub>0</sub>斜入射至平面标准镜(5)前表面时的入射光场为:E(t)=E<sub>0</sub>exp{i(ω<sub>0</sub>t+kt<sup>2</sup>)}       (公式1)其中,<img file="FDA0000505782170000011.GIF" wi="165" he="129" />为调频带宽的变化率,T为调频周期,△F为调频带宽,E<sub>0</sub>为入射光场振幅,t为时间,ω<sub>0</sub>为光场角频率,i表示虚数,线性调频激光器发出的线性调频偏振光到达平面标准镜前表面的光程为l,则t‑l/c时刻到达平面标准镜前表面的反射光场为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>&alpha;</mi><msub><mi>E</mi><mn>0</mn></msub><mi>exp</mi><mo>{</mo><mi>i</mi><mo>[</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>l</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>l</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>}</mo></mrow>]]></math><img file="FDA0000505782170000012.GIF" wi="813" he="125" /></maths>       (公式2)同一时刻经过平面标准镜后表面反射后透射出前表面的光场表达式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&alpha;&alpha;</mi><mn>1</mn><mn>2</mn></msubsup><msub><mi>E</mi><mn>0</mn></msub><mi>exp</mi><mo>{</mo><mi>i</mi><mo>[</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>l</mi><mo>+</mo><mn>2</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>l</mi><mo>+</mo><mn>2</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>&omega;</mi><mn>0</mn></msub><mi>nd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>]</mo><mo>}</mo></mrow>]]></math><img file="FDA0000505782170000013.GIF" wi="1610" he="130" /></maths> (公式3)其中,α和α<sub>1</sub>分别为平面标准镜的反射系数和透射系数,d为平面标准镜厚度,θ为入射光折射角,n为平面标准镜折射率,c为光速,光电探测器接收到的总光场表示为:E′(t)=E<sub>1</sub>(t)+E<sub>2</sub>(t)    (公式4)则光电探测器输出的光电流表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>I</mi><mo>=</mo><mfrac><mi>&eta;e</mi><mi>hv</mi></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>D</mi></munder><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><msup><mrow><mo>[</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mrow><mo>*</mo></msup><mi>ds</mi></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mi>&eta;e</mi><mrow><mn>2</mn><mi>hv</mi></mrow></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>D</mi></munder><mo>[</mo><msubsup><mi>E</mi><mn>1</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>E</mi><mn>2</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><msub><mi>E</mi><mn>2</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><msub><mi>E</mi><mn>1</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo><mi>ds</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000505782170000014.GIF" wi="1239" he="291" /></maths>   (公式5)其中,e为电子电量,Z为探测器表面介质的本征阻抗,η为量子效率,D为探测器光敏面的面积,h为普朗克常数,v为激光频率,*号表示复数共轭,直流项经过低通滤波器滤除后,对公式5进行整理获得只含交流项的中频电流为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>I</mi><mi>IF</mi></msub><mo>=</mo><mfrac><mi>&eta;e</mi><mrow><mn>2</mn><mi>hv</mi></mrow></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>D</mi></munder><mo>[</mo><mrow><mo>(</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><msub><mi>E</mi><mn>2</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><msub><mi>E</mi><mn>1</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo><mi>ds</mi></mrow>]]></math><img file="FDA0000505782170000021.GIF" wi="950" he="141" /></maths>     (公式6)将(公式2)和(公式3)代入(公式6),获得计算结果为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='' close='' separators=' '><mtable><mtr><mtd><msub><mi>I</mi><mi>IF</mi></msub><mo>=</mo><mfrac><mi>&eta;e</mi><mi>hv</mi></mfrac><mfrac><mi>&pi;</mi><mi>Z</mi></mfrac><msubsup><mi>E</mi><mn>0</mn><mn>2</mn></msubsup><msup><mi>&alpha;</mi><mn>2</mn></msup><msubsup><mi>&alpha;</mi><mn>1</mn><mn>2</mn></msubsup><mi>cos</mi><mo>{</mo><mo>[</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>l</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mi>l</mi><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>-</mo></mtd></mtr><mtr><mtd><mo>[</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>l</mi><mo>+</mo><mn>2</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mo>+</mo><mi>k</mi><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mfrac><mrow><mi>l</mi><mo>+</mo><mn>2</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>&omega;</mi><mn>0</mn></msub><mi>nd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mi>&eta;e</mi><mi>hv</mi></mfrac><mfrac><mi>&pi;</mi><mi>Z</mi></mfrac><msubsup><mi>E</mi><mn>0</mn><mn>2</mn></msubsup><msup><mi>&alpha;</mi><mn>2</mn></msup><msubsup><mi>&alpha;</mi><mn>1</mn><mn>2</mn></msubsup><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>4</mn><mi>knd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mi>t</mi><mo>-</mo><mfrac><mrow><mn>4</mn><mi>knd</mi><mi>cos</mi><mi>&theta;</mi><mrow><mo>(</mo><mi>l</mi><mo>+</mo><mi>nd</mi><mi>cos</mi><mi>&theta;</mi><mo>)</mo></mrow></mrow><msup><mi>c</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow></mtd></mtr></mtable><mrow><mo>]</mo><mo>}</mo></mrow></mfenced>]]></math><img file="FDA0000505782170000022.GIF" wi="1665" he="397" /></maths> (公式7)信号处理系统(7)接收光电探测器(6)发出的光电流I,并对光电流I进行处理,采用傅里叶变换获得该光电流I中的外差信号频率f<sub>IF</sub>,所述外差信号频率f<sub>IF</sub>与线性调频激光器(1)入射到平面标准镜(4)前表面后的折射角θ之间的关系式为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>f</mi><mi>IF</mi></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mi>knd</mi><mi>cos</mi><mi>&theta;</mi></mrow><mi>&pi;c</mi></mfrac><mo>=</mo><mi>K</mi><mi>cos</mi><mi>&theta;</mi></mrow>]]></math><img file="FDA0000505782170000023.GIF" wi="549" he="130" /></maths>      (公式8)式中,K为比例系数,该比例系数根据公式:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>knd</mi></mrow><mi>&pi;c</mi></mfrac></mrow>]]></math><img file="FDA0000505782170000024.GIF" wi="224" he="121" /></maths>  (公式10)计算获得,由折射定律获得折射角θ和待测激光入射角θ<sub>0</sub>之间的关系为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>cos</mi><mi>&theta;</mi><mo>=</mo><mi>cos</mi><mo>[</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><mi>sin</mi><msub><mi>&theta;</mi><mn>0</mn></msub></mrow><mi>n</mi></mfrac><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000505782170000025.GIF" wi="537" he="130" /></maths>   (公式9)进而获得待测的激光入射角θ<sub>0</sub>。
地址 150080 黑龙江省哈尔滨市南岗区学府路74号