发明名称 基于M<sup>2</sup>因子和光闪烁指数确定大气湍流参数方法及装置
摘要 本发明涉及一种基于测量M<sup>2</sup>因子和光闪烁指数<img file="dda0000031114770000011.GIF" wi="119" he="46" />来确定大气湍流参数的方法及其实验装置,属于空间激光通信领域。该方法是先将激光器发出的高斯光束经过光学发射天线准直扩束后发射,使发射出的高斯准直光束通过大气湍流再被光学接收天线接收;然后将接收到的光束分成两束,一束通过M<sup>2</sup>因子测量仪,另一束通过光闪烁测量仪;最后利用测得的M<sup>2</sup>因子和光闪烁指数来确定大气湍流参数。实现上述方法的实验装置包括激光器,光学发射天线,光学接收天线,分束器,M<sup>2</sup>因子测量仪,光闪烁测量仪,数据处理器。本发明的方法及其实验装置不仅能确定大气湍流折射率结构常数和大气湍流内尺度参数,且实验装置结构简单,操作测量方便。
申请公布号 CN101980056B 申请公布日期 2012.10.03
申请号 CN201010533777.6 申请日期 2010.11.06
申请人 四川大学 发明人 张彬;但有全;潘平平;齐娜
分类号 G01W1/00(2006.01)I;G01N21/17(2006.01)I;G01N21/41(2006.01)I;G01B11/00(2006.01)I;H04B10/10(2006.01)I 主分类号 G01W1/00(2006.01)I
代理机构 成都科海专利事务有限责任公司 51202 代理人 刘双兰
主权项 1.一种基于M<sup>2</sup>因子和光闪烁指数确定大气湍流参数的方法,包括以下步骤:(1)将激光器发出的高斯光束经过光学发射天线准直扩束后发射,使发射出的高斯准直光束通过大气湍流后被光学接收天线接收;(2)然后用分束器将光学接收天线接收到的高斯准直光束分成两束,一束通过M<sup>2</sup>因子测量仪,另一束通过光闪烁测量仪;(3)最后根据M<sup>2</sup>因子测量仪和光闪烁测量仪所测得的数值来确定大气湍流参数中大气湍流折射率结构常数和大气湍流内尺度值;所述高斯光束均为完全相干高斯光束,且M<sup>2</sup>因子和光闪烁指数满足以下两个函数关系式:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>n</mi><mn>2</mn></msubsup><msubsup><mi>l</mi><mn>0</mn><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>M</mi><mn>2</mn></msup><mo>)</mo></mrow><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><mi>C</mi><mi>n</mi><mn>2</mn></msubsup><mi>g</mi><mrow><mo>(</mo><msub><mi>l</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&sigma;</mi><mi>I</mi><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>公式(1)中,<img file="FDA00001856314400013.GIF" wi="45" he="48" />为大气湍流折射率结构常数,l<sub>0</sub>为大气湍流内尺度,f(M<sup>2</sup>)表示M<sup>2</sup>因子的函数关系式,所述公式(1)在弱起伏条件和强起伏条件下均成立;公式(2)中,<img file="FDA00001856314400014.GIF" wi="58" he="48" />为光闪烁指数,g(l<sub>0</sub>)表示大气湍流内尺度l<sub>0</sub>的函数关系式;所述弱起伏条件和强起伏条件的公式表示为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>R</mi><mn>2</mn></msubsup><mo>=</mo><mn>1.23</mn><msubsup><mi>C</mi><mi>n</mi><mn>2</mn></msubsup><msup><mi>k</mi><mrow><mn>7</mn><mo>/</mo><mn>6</mn></mrow></msup><msup><mi>L</mi><mrow><mn>11</mn><mo>/</mo><mn>6</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>公式(3)中,σ<sub>R</sub><sup>2</sup>为Rytov方差,k为光波数,L为光束传输距离,σ<sub>R</sub><sup>2</sup>≥1时,表示强起伏条件成立;σ<sub>R</sub><sup>2</sup>&lt;1时,表示弱起伏条件成立;通过公式(1)和公式(2)计算得出大气湍流折射率结构常数<img file="FDA00001856314400016.GIF" wi="44" he="48" />和大气湍流内尺度l<sub>0</sub>,其相应的基本步骤如下:1)为公式(1)和公式(2)的参数赋值:高斯光束波长λ,束腰宽度w<sub>0</sub>,及传输距离L,并输入M<sup>2</sup>因子和光闪烁指数测量值;2)判断光闪烁指数测量值是否小于1,如果是则执行弱起伏条件下的计算步骤A);否则执行强起伏条件下的计算步骤B);A)弱起伏条件下的计算步骤:A1)利用公式(1)计算出f(M<sup>2</sup>)的值;A2)利用公式(1)和公式(2)得方程:<img file="FDA00001856314400021.GIF" wi="425" he="57" />并用数值方法求解此方程,得大气湍流内尺度l<sub>0</sub>;A3)将计算出的大气湍流内尺度l<sub>0</sub>代入公式(1),即可算出大气湍流折射率结构常数<img file="FDA00001856314400022.GIF" wi="72" he="57" />A4)进一步将计算出的大气湍流内尺度l<sub>0</sub>和大气湍流折射率结构常数<img file="FDA00001856314400023.GIF" wi="74" he="57" />代入公式(3),计算Rytov方差<img file="FDA00001856314400024.GIF" wi="74" he="57" />并判断它是否小于1,如果小于1,则将计算出的<img file="FDA00001856314400025.GIF" wi="50" he="57" />和l<sub>0</sub>作为最终结果直接输出,同时输出<img file="FDA00001856314400026.GIF" wi="51" he="57" />的值;否则,执行强起伏条件下的计算步骤B);B)强起伏条件下的计算步骤:B1)将大气湍流内尺度l<sub>0</sub>取为通常值,即l<sub>0</sub>=1mm;B2)将大气湍流内尺度l<sub>0</sub>=1mm及M<sup>2</sup>因子测量值代入公式(1),即可计算出大气湍流折射率结构常数<img file="FDA00001856314400027.GIF" wi="74" he="57" />并输出其值作为强起伏条件下测量<img file="FDA00001856314400028.GIF" wi="50" he="57" />的结果,同时输出大气湍流内尺度l<sub>0</sub>=1mm,以及由公式(3)计算出的Rytov方差<img file="FDA00001856314400029.GIF" wi="50" he="57" />的值。
地址 610207 四川省成都市双流县川大路二段2号