[转载]GLAD：地对空激光通讯系统 [复制链接]

本例为一束激光从地基激光器出发传输到近地轨道上的一个转换镜上。转换镜将激光反射到一个聚焦镜上，然后这束光打到大气层内部一个低海拔的目标上。根据Kolmogorov模型，假设目标值的半径为10cm，就可以计算大气像差。本例中包含了激光扩束器像差、转换镜上的像差以及聚焦镜上的像差。 n Syq}Y3  
大气模型假设的波阵面光谱功率为（忽略内部和外部的尺寸限制）： +:uz=~mo`  
1Vdi5;dn  
Nwz?*~1  
z,4mg6gt  
其中W^2 (f)是波阵面的光谱功率，r0为可视参数，f是空间频率，L0是外部尺寸，Li是内部尺寸，这些参数的单位分别为rad，m，m-1。由于大气像差和光束扩束器的像散，斯特列尔比SR=0.34。经过一个激励器影响半径为4.0cm的自适应镜校正后，斯特列尔比为0.87。经过全程传播到达目标后，光斑直径为50cm，剩下56%的能量，相对于没有自适应镜时能量的22%，有了明显的提高。由于没有考虑大气对光学元件散射效应的衰减效应，所以实际中传输到的光更加少。 l>}f{az-T  
s6@mXO:H^  
xsB0LUt  
图1.地对空激光通信系统示意图 .dk<?BI#H  
q6)fP4MQ]  
jZzTnmm&?  
表1.关键参数 _yv#v_Z  
!*}UP|8  
###激光器光束初始化 <kdlXS>J.  
set/alias/off s@ r{TXEn  
wavelength/set 1 .48               # 设置激光器波长 \O;2^  
array/set 1 256                    # 设置计算初始矩阵大小 R ,qQC<  
units/s 1 .1 3{$7tck,  
gauss/c/c 1 1 1.25                 # 定义高斯光束 GeTk/tU  
clap/cir/con 1 1.25 /7 Tm2Vj8  
energy/norm 1 1.                   # 能量归一化设置 Si>38vCJ*  
set/density 64                     # 设置画图线条密度 R,G*]/r`  
title 1: starting laser distribution lb.Q^TghU  
plot/watch ex26_1.plt OALNZKP  
plot/liso nsl=64 xr=1.5 yr=1.5       # 绘制激光束初始强度曲线 2Q
p}f^  
####激光器初始光束相位分布如图1所示： "0PrdZMx  
= 'NV3by  
B={_}
f  
图1激光器光束初始分布 &\N>N7/1  
HCs^?s8Pp  
##光束扩束器（20X）模拟 CpO_p%P  
mirror 1 20 focallength            # 扩束器透镜1焦距设置为20cm. (c(F1=K  
dist -420                          # 透镜分离 j4?Qd0z  
mirror 1 400 focallength           #扩束器透镜2焦距设置为400cm. 8x U*j  
abr/ast 1 .2 45                    # 光束峰谷像散设置 ]8\I{LR  
clap/cir/con 1 25                  # 设置孔径光阑直径50 J+`aj8_B  
phase/random/kolmogorov 1 10. 7    # 光束附加大气像差 =Z  
strehl !k)6r6  
title 2: phase after beam expander and atmoshpheric turbulence l~rj7f;  
plot/watch ex26_2.plt >#|%'Us  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 ]=00<~ l*q  
$V(]z`b&  
Reikf}9Q  
图2.经过扩束器及大气干扰后光束相位分布 vd4@jZ5  
dtC@cK/,D  
####光束传输至自适应镜 HU $"o6ap  
adapt 1. 4.                        # 自适应镜对光束的曲率半径影响为4cm 2'Dl$DH  
strehl `FMo;
,j  
title 3: phase after adaptive optic correction "l56?@-x  
plot/watch ex26_3.plt %xg+UW }  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 2 z%YNZ^d  
####光束经过自适应镜之后的相位分布如图3 /l*v *tl  
('5?-  
 vmqa_gU\  
图3.光束经过自适应镜校正之后的相位分布 X_TjJmc  
f8]sjeY  
####光束传输至转换镜 8%7%[WC#  
dist 5e7                           #光束传输500km至转换镜 vFVUdxPOw  
title 4: irradiance at relay mirror before aperture u#@/^h;  
plot/watch ex26_4.plt 
Y<xqws  
plot/liso xr=300 yr=300 nsl=64 6h;(b2p{  
####光束到达转换镜前的相位分布如图4所示： R |c=I}@F  
7AV{ h[J  
6{y7e L3!  
图4.光束到达转换镜前的相位分布 Cs6`lX >  
xnY?<?J"!  
####光束传输经过转换镜 '47 b"uV  
clap/cir/con 1 75                  # aperture of 1.5 m. dia. for relay o&?c,FwN  
title 5: irradiance at relay mirror after aperture ,MCTb'=G  
plot/watch ex26_5.plt =Zq6iMD  
plot/liso xr=300 yr=300 nsl=64 |Kd#pYt%O  
mirror 1 3e8 foc $D^\[^S  
abr/ast 1 .3 90 ,n/]ALz>~  
####光束经过转换镜后相位分布如图5所示： 7ftn gBv?  
PftK>,+,  
$o[-xNn1  
图5.光束经过转换镜后相位分布 m-AF&( ;K  
PJ.\)oP  
####光束传输至聚焦镜 ?"x4u#x  
dist 3e8                           #光束传输3000km至聚焦镜 V?u#WJy/  
clap/cir/con 1 200                 # 聚焦镜孔径直径4m kTb$lLG\xk  
energy u&:N`f  
title 6: irradiance on focusing mirror \=ML*Gi*  
plot/watch ex26_6.plt ~:a1ELqVw  
plot/liso xr=300 yr=300 nsl=64 +U oNJ   
####光束传输至聚焦镜相位分布如图6所示： 24jf`1XFW  
g&$=Y7G  
yLa@27T\A  
图6.光束传输至聚焦镜相位分布 9-42A7g^C  
'c3
5%?]  
####光束经过聚焦镜传输至聚焦靶面后相位分布 ],]Rv#`  
phase/random 1 .1 50               # 聚焦镜附加随机像差 rnIjpc F  
mirror 1 5e7 focallength           # 聚焦镜焦距500km *g(d}C!  
dist 5e7                           # 传输至靶面 ppAbG,7  
clap/cir/con 1 25                  # 设置出瞳孔径直径50cm vzIo2,/7  
energy UQl?
_[G  
title 7: target irradiance inside 50 cm. NL9.J@"b  
plot/watch ex26_7.plt $plk>Khg  
plot/liso xr=50 yr=50 nsl=64 ?2<) Jw  
end %ih\|jRt  
####光束传输至聚焦靶面后相位分布如图7所示： r[L.TX3Ah=  
;+hh|NiQ  
P6gkbtg  
图7.光束传输至聚焦靶面后相位分布

不错，值得学习。