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

本例为一束激光从地基激光器出发传输到近地轨道上的一个转换镜上。转换镜将激光反射到一个聚焦镜上，然后这束光打到大气层内部一个低海拔的目标上。根据Kolmogorov模型，假设目标值的半径为10cm，就可以计算大气像差。本例中包含了激光扩束器像差、转换镜上的像差以及聚焦镜上的像差。 s+G(N$0U  
大气模型假设的波阵面光谱功率为（忽略内部和外部的尺寸限制）： r(wtuD23q  
pjSM7PhQ  
cP@H8
|c=  
np}0OX  
其中W^2 (f)是波阵面的光谱功率，r0为可视参数，f是空间频率，L0是外部尺寸，Li是内部尺寸，这些参数的单位分别为rad，m，m-1。由于大气像差和光束扩束器的像散，斯特列尔比SR=0.34。经过一个激励器影响半径为4.0cm的自适应镜校正后，斯特列尔比为0.87。经过全程传播到达目标后，光斑直径为50cm，剩下56%的能量，相对于没有自适应镜时能量的22%，有了明显的提高。由于没有考虑大气对光学元件散射效应的衰减效应，所以实际中传输到的光更加少。 _ |HA\!  
%P#| }  
vQ$"|8,  
图1.地对空激光通信系统示意图 BZXee>3"  
f6Y-ss;'  
ATmyoN2@>  
表1.关键参数 q%/.+g2-\  
AAB_Ytf  
###激光器光束初始化 uOKdb6]r6  
set/alias/off 1UB.2}/:  
wavelength/set 1 .48               # 设置激光器波长 Zx6h%l,%  
array/set 1 256                    # 设置计算初始矩阵大小 "EWq{l_I5$  
units/s 1 .1 9j5Z!Vsy  
gauss/c/c 1 1 1.25                 # 定义高斯光束 jC?l :m?  
clap/cir/con 1 1.25 BuC\Bd^0  
energy/norm 1 1.                   # 能量归一化设置 N"~P$B1X  
set/density 64                     # 设置画图线条密度 n}4q2x"  
title 1: starting laser distribution As tuM]  
plot/watch ex26_1.plt pB%oFWqK  
plot/liso nsl=64 xr=1.5 yr=1.5       # 绘制激光束初始强度曲线 j^f54Ky.  
####激光器初始光束相位分布如图1所示： 37M,Os1(  
X .K*</(g  

]}]+aB  
图1激光器光束初始分布 nI+.De~  
u'Hh||La"  
##光束扩束器（20X）模拟 g^i\7'  
mirror 1 20 focallength            # 扩束器透镜1焦距设置为20cm. {sq:vu@NC  
dist -420                          # 透镜分离 7v
.O Lp  
mirror 1 400 focallength           #扩束器透镜2焦距设置为400cm. x&EMg!  
abr/ast 1 .2 45                    # 光束峰谷像散设置 b1."mT!p  
clap/cir/con 1 25                  # 设置孔径光阑直径50 o{mVXidE  
phase/random/kolmogorov 1 10. 7    # 光束附加大气像差 k@[[vj|W  
strehl #}tdA( -  
title 2: phase after beam expander and atmoshpheric turbulence 65mfq&"P?  
plot/watch ex26_2.plt {)Wa"|+  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 uOd1:\%*  
Zl]@;*u  
x{rjngp2  
图2.经过扩束器及大气干扰后光束相位分布 =i/r:  
-|=)
 
####光束传输至自适应镜 ##1/{9ywy  
adapt 1. 4.                        # 自适应镜对光束的曲率半径影响为4cm n+vv %  
strehl Mdu\ci)lr  
title 3: phase after adaptive optic correction Sj8fo^K50  
plot/watch ex26_3.plt C 8d9(u  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 2 jpMMnEVj6P  
####光束经过自适应镜之后的相位分布如图3 d9T:0A`M  
RyWfoLc  
2Z`
Jr/  
图3.光束经过自适应镜校正之后的相位分布 q+w] Xs;  
Vk76cV D  
####光束传输至转换镜 _C'VC#Sy  
dist 5e7                           #光束传输500km至转换镜 3u;0,:X&  
title 4: irradiance at relay mirror before aperture AthR|I|8  
plot/watch ex26_4.plt #Ny+6XM  
plot/liso xr=300 yr=300 nsl=64 ("<4Ry.u  
####光束到达转换镜前的相位分布如图4所示： D>-Pv-f/  
|Fq\%y#  
RE]*fRe7#  
图4.光束到达转换镜前的相位分布 $)=`Iai  
{OS[0LB  
####光束传输经过转换镜 JX{rum  
clap/cir/con 1 75                  # aperture of 1.5 m. dia. for relay v|3mbApv  
title 5: irradiance at relay mirror after aperture ZA'0q  
plot/watch ex26_5.plt C MGDg}  
plot/liso xr=300 yr=300 nsl=64 NZwi3  
mirror 1 3e8 foc :e}j$vF  
abr/ast 1 .3 90 f\Q_]%^W  
####光束经过转换镜后相位分布如图5所示： v~YGef;D  
d%p{l)Hd  
9h6siK(F  
图5.光束经过转换镜后相位分布 /-=h|A#Kh  
KHeeB`V>J  
####光束传输至聚焦镜 1ZvXRJ)%  
dist 3e8                           #光束传输3000km至聚焦镜 B? XK;*])  
clap/cir/con 1 200                 # 聚焦镜孔径直径4m tC7 4=  
energy zIU6bMMT3u  
title 6: irradiance on focusing mirror Go[anf  
plot/watch ex26_6.plt I.%EYAai  
plot/liso xr=300 yr=300 nsl=64 m\|E
M'@k  
####光束传输至聚焦镜相位分布如图6所示： ~cfvL*~5  
xi)M8\K  
5mm&l+N)  
图6.光束传输至聚焦镜相位分布 }0OQm?xh  
X%`:waR  
####光束经过聚焦镜传输至聚焦靶面后相位分布 QS-X_  
phase/random 1 .1 50               # 聚焦镜附加随机像差 @U =~c9  
mirror 1 5e7 focallength           # 聚焦镜焦距500km $vnx)#r3  
dist 5e7                           # 传输至靶面 Z)}2bJwA  
clap/cir/con 1 25                  # 设置出瞳孔径直径50cm %+C6#cj  
energy ^<j =.E  
title 7: target irradiance inside 50 cm. U.N&~S  
plot/watch ex26_7.plt V7 hO}  
plot/liso xr=50 yr=50 nsl=64 veS) j?4  
end !0v3Lu~j  
####光束传输至聚焦靶面后相位分布如图7所示： 6O*lZNN  
NK%Ok  
]qEg5:yY  
图7.光束传输至聚焦靶面后相位分布

不错，值得学习。