发明名称 一种远场光场全矢量计算的方法
摘要 本发明涉及一种远场光场全矢量计算的方法,通过其他全矢量计算方法得到光学器件的近场分布后,通过全矢量积分的方式,得到光场远场任意点的分布。既保留了全矢量计算的精度,又兼具了衍射算法的速度和远场可计算性。计算结果比标量近似方法大为精确,不仅可以准确地计算出远场光场的强度,也可以精确计算远场光场的位相。同时,由于积分公式占用内存很小,且与远场距离无关,可以计算其他全矢量方法无法计算的远场,兼具快速高效的优点。可广泛应用于光学领域。
申请公布号 CN104090364B 申请公布日期 2016.03.02
申请号 CN201410209050.0 申请日期 2014.05.16
申请人 北京理工大学 发明人 胡滨;刘娟;王涌天;韩剑
分类号 G02B27/00(2006.01)I 主分类号 G02B27/00(2006.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 李迪
主权项 一种远场光场全矢量计算的方法,其特征在于,首先通过全矢量计算方法得到光学器件的某一近场平面Γ的近场光场分布;然后通过基于格林函数的全矢量积分方法得到远场任意一点P<sub>0</sub>的光场分布,计算方法如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mo>&Integral;</mo><mi>s</mi></msub><mi>i</mi><mi>&omega;</mi><mi>G</mi><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>,</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>&times;</mo><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>+</mo><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>&times;</mo><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>&times;</mo><mo>&dtri;</mo><mi>G</mi><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>,</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>+</mo><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>&dtri;</mo><mi>G</mi><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>,</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mi>d</mi><mi>&Gamma;</mi></mrow>]]></math><img file="FDA0000847065960000011.GIF" wi="1525" he="100" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>=</mo><msub><mo>&Integral;</mo><mi>s</mi></msub><mi>i</mi><mi>&omega;</mi><mi>G</mi><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>,</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>&times;</mo><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>+</mo><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>&times;</mo><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>&times;</mo><mo>&dtri;</mo><mi>G</mi><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>,</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>+</mo><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>&dtri;</mo><mi>G</mi><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>,</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>)</mo></mrow><mi>d</mi><mi>&Gamma;</mi></mrow>]]></math><img file="FDA0000847065960000012.GIF" wi="1477" he="101" /></maths>其中<img file="FDA0000847065960000014.GIF" wi="51" he="69" />为电场分量,<img file="FDA0000847065960000015.GIF" wi="60" he="70" />为磁场分量,ω为入射光的角频率,<img file="FDA0000847065960000016.GIF" wi="43" he="66" />为法向方向,<img file="FDA0000847065960000017.GIF" wi="41" he="69" />为远场某一点P<sub>0</sub>点位置,<img file="FDA0000847065960000018.GIF" wi="52" he="68" />为近场平面Γ上某一点的位置,S为近场平面Γ的面积,G为格林函数,可以表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>G</mi><mo>=</mo><mfrac><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>ik</mi><mn>0</mn></msub><mo>|</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>-</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>|</mo></mrow></msup><mrow><mn>4</mn><mi>&pi;</mi><mo>|</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>-</mo><msup><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&prime;</mo></msup><mo>|</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000847065960000013.GIF" wi="294" he="206" /></maths>其中k<sub>0</sub>为光在空间中传播的波矢e为自然底数,约为2.718281828。
地址 100081 北京市海淀区中关村南大街5号北京理工大学