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

本例为一束激光从地基激光器出发传输到近地轨道上的一个转换镜上。转换镜将激光反射到一个聚焦镜上，然后这束光打到大气层内部一个低海拔的目标上。根据Kolmogorov模型，假设目标值的半径为10cm，就可以计算大气像差。本例中包含了激光扩束器像差、转换镜上的像差以及聚焦镜上的像差。 s :BW}PM  
大气模型假设的波阵面光谱功率为（忽略内部和外部的尺寸限制）：
nnN$?'%~6  
frDMFEXXP  
N-&ZaK  
D)DD6  
其中W^2 (f)是波阵面的光谱功率，r0为可视参数，f是空间频率，L0是外部尺寸，Li是内部尺寸，这些参数的单位分别为rad，m，m-1。由于大气像差和光束扩束器的像散，斯特列尔比SR=0.34。经过一个激励器影响半径为4.0cm的自适应镜校正后，斯特列尔比为0.87。经过全程传播到达目标后，光斑直径为50cm，剩下56%的能量，相对于没有自适应镜时能量的22%，有了明显的提高。由于没有考虑大气对光学元件散射效应的衰减效应，所以实际中传输到的光更加少。 ceD6q~)  
2 'D,1F  
s^IC]sW\%  
图1.地对空激光通信系统示意图 XqUQ{^;aI  
0'.z|
Jg=  
.-mIU.Nwi  
表1.关键参数
.boBb<  
:4Nv6X61  
###激光器光束初始化 Y<(7
u`F  
set/alias/off V_D wHq2  
wavelength/set 1 .48               # 设置激光器波长
=EM<LjO  
array/set 1 256                    # 设置计算初始矩阵大小 G3+e5/0  
units/s 1 .1 :A!EjIL`#  
gauss/c/c 1 1 1.25                 # 定义高斯光束 ^ <`SUBI  
clap/cir/con 1 1.25 DR3om;Uk  
energy/norm 1 1.                   # 能量归一化设置 y6-P6T  
set/density 64                     # 设置画图线条密度 *{j;LA.BR#  
title 1: starting laser distribution zTfjuI|R  
plot/watch ex26_1.plt $BgaLJs/O  
plot/liso nsl=64 xr=1.5 yr=1.5       # 绘制激光束初始强度曲线 GD d'{qE6  
####激光器初始光束相位分布如图1所示： LklE,W  
z;2& d<h  
?I?~BWu  
图1激光器光束初始分布 T}1"  
cJ@fJ|  
##光束扩束器（20X）模拟 }vx,i99W?  
mirror 1 20 focallength            # 扩束器透镜1焦距设置为20cm. 5pDE!6gQ  
dist -420                          # 透镜分离 #W|Obc]K  
mirror 1 400 focallength           #扩束器透镜2焦距设置为400cm. =54D#,[B  
abr/ast 1 .2 45                    # 光束峰谷像散设置 .m8l\h^3  
clap/cir/con 1 25                  # 设置孔径光阑直径50 qxAh8RR;/  
phase/random/kolmogorov 1 10. 7    # 光束附加大气像差 EA1&D^nT  
strehl { pQJ.QI  
title 2: phase after beam expander and atmoshpheric turbulence &{glwVKV  
plot/watch ex26_2.plt R@NFpiw  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 "eOl(TSu/  
hEZo{0:b"  
NF4(+E9g  
图2.经过扩束器及大气干扰后光束相位分布 xoj,>[7 D  
eFS$;3FP1  
####光束传输至自适应镜 vO~w~u5  
adapt 1. 4.                        # 自适应镜对光束的曲率半径影响为4cm X"`[&l1  
strehl bXc*d9]  
title 3: phase after adaptive optic correction 0<uLQVoR2n  
plot/watch ex26_3.plt {dxl8~/I  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 2 I0Allw[  
####光束经过自适应镜之后的相位分布如图3 7A'E+>1d  
K;rgLj0m  
x0{B7/FN  
图3.光束经过自适应镜校正之后的相位分布 4:Ju|g]O  
_unoDoB  
####光束传输至转换镜 8d8jUPFQ  
dist 5e7                           #光束传输500km至转换镜 &s}sA+w  
title 4: irradiance at relay mirror before aperture pCo3%(  
plot/watch ex26_4.plt _%Xp2`
m  
plot/liso xr=300 yr=300 nsl=64 >b9J!'G,(  
####光束到达转换镜前的相位分布如图4所示： 6LCR ;~ ]  
1BW9,Xr  
,YBO}l  
图4.光束到达转换镜前的相位分布 Bfo#N31F}  
o@dTiQK_  
####光束传输经过转换镜 J?dz>3Rhx9  
clap/cir/con 1 75                  # aperture of 1.5 m. dia. for relay -~xQ@+./  
title 5: irradiance at relay mirror after aperture )eZ}Kt+  
plot/watch ex26_5.plt 6)ZaK  
plot/liso xr=300 yr=300 nsl=64 \<xo`2b  
mirror 1 3e8 foc L aA<`  
abr/ast 1 .3 90 +_*NY~  
####光束经过转换镜后相位分布如图5所示： <"`f!k#[  
?q Q.Wj6Mj  
)!E:  
图5.光束经过转换镜后相位分布 W)]&G}U<  
[4YRyx&:++  
####光束传输至聚焦镜 T95FoA  
dist 3e8                           #光束传输3000km至聚焦镜 VB4V[jraCF  
clap/cir/con 1 200                 # 聚焦镜孔径直径4m o$%KbfXO]  
energy n{~&^Nby*I  
title 6: irradiance on focusing mirror |&-*&)iD|w  
plot/watch ex26_6.plt R1Q,
m  
plot/liso xr=300 yr=300 nsl=64 rwW"B  
####光束传输至聚焦镜相位分布如图6所示： )G, S7A  
}1V+8'D  
i?@7>Ca  
图6.光束传输至聚焦镜相位分布 QQ{*j7i)  
CH!\uK22  
####光束经过聚焦镜传输至聚焦靶面后相位分布 mAW(j@5sp  
phase/random 1 .1 50               # 聚焦镜附加随机像差 MNTVG&h  
mirror 1 5e7 focallength           # 聚焦镜焦距500km NRP)'E  
dist 5e7                           # 传输至靶面 BdU .;_K  
clap/cir/con 1 25                  # 设置出瞳孔径直径50cm cp~6\F;c  
energy *&]8rm{  
title 7: target irradiance inside 50 cm. $RF.L
Vc  
plot/watch ex26_7.plt f>cUdEPBb  
plot/liso xr=50 yr=50 nsl=64 {c]
dz7'?  
end ;ZcwgsxTM  
####光束传输至聚焦靶面后相位分布如图7所示： s?1Aj<  
IN#/~[W  
5?Q5cD2]\6  
图7.光束传输至聚焦靶面后相位分布

不错，值得学习。