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

本例为一束激光从地基激光器出发传输到近地轨道上的一个转换镜上。转换镜将激光反射到一个聚焦镜上，然后这束光打到大气层内部一个低海拔的目标上。根据Kolmogorov模型，假设目标值的半径为10cm，就可以计算大气像差。本例中包含了激光扩束器像差、转换镜上的像差以及聚焦镜上的像差。 #p]V?  
大气模型假设的波阵面光谱功率为（忽略内部和外部的尺寸限制）： `
Q#)N0  
Y3s8@0b3  
Zsmv{p  
Gm\)1b  
其中W^2 (f)是波阵面的光谱功率，r0为可视参数，f是空间频率，L0是外部尺寸，Li是内部尺寸，这些参数的单位分别为rad，m，m-1。由于大气像差和光束扩束器的像散，斯特列尔比SR=0.34。经过一个激励器影响半径为4.0cm的自适应镜校正后，斯特列尔比为0.87。经过全程传播到达目标后，光斑直径为50cm，剩下56%的能量，相对于没有自适应镜时能量的22%，有了明显的提高。由于没有考虑大气对光学元件散射效应的衰减效应，所以实际中传输到的光更加少。 ecO$L<9>  
XsDZ<j%x89  
q&_\A0  
图1.地对空激光通信系统示意图 f9 :=6  
~b0l?P*Ff  
vK+!m~kDu  
表1.关键参数 }2:q#}"  
7FD,TJs  
###激光器光束初始化 G l2WbY  
set/alias/off SuNc&e#(  
wavelength/set 1 .48               # 设置激光器波长 !&3"($-U3G  
array/set 1 256                    # 设置计算初始矩阵大小 b\zq,0%  
units/s 1 .1 J!G92A~*]  
gauss/c/c 1 1 1.25                 # 定义高斯光束 Fy!s$!\C0  
clap/cir/con 1 1.25 +nim47  
energy/norm 1 1.                   # 能量归一化设置 |T/s>OW  
set/density 64                     # 设置画图线条密度 {K}Dpy  
title 1: starting laser distribution qh&q<
M  
plot/watch ex26_1.plt [*I7^h%  
plot/liso nsl=64 xr=1.5 yr=1.5       # 绘制激光束初始强度曲线 jTZi< Y:bB  
####激光器初始光束相位分布如图1所示： #W l^!)#j?  
,fN <I  
?<Hgq8J  
图1激光器光束初始分布 aY;34SF  
O1\25D  
##光束扩束器（20X）模拟 )8'v@8;-  
mirror 1 20 focallength            # 扩束器透镜1焦距设置为20cm. o?~27   
dist -420                          # 透镜分离 9}Za_ZgG  
mirror 1 400 focallength           #扩束器透镜2焦距设置为400cm. -[pfLo  
abr/ast 1 .2 45                    # 光束峰谷像散设置 ~l.C
-  
clap/cir/con 1 25                  # 设置孔径光阑直径50 %cDDu$9;  
phase/random/kolmogorov 1 10. 7    # 光束附加大气像差 +2}Ar<elP  
strehl :#_k`{WG  
title 2: phase after beam expander and atmoshpheric turbulence cxvO,8NiB  
plot/watch ex26_2.plt vjh'<5w9Wi  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 ^k{/Yl  
+lZ-xU1  
c* ~0R?  
图2.经过扩束器及大气干扰后光束相位分布 $:1/`m19  
lw?C:-m  
####光束传输至自适应镜 HZS.%+2  
adapt 1. 4.                        # 自适应镜对光束的曲率半径影响为4cm $G9E=wn  
strehl k>8OxpaWv?  
title 3: phase after adaptive optic correction + 65<|0  
plot/watch ex26_3.plt ,
KF>@3f  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 2 2n2,MB  
####光束经过自适应镜之后的相位分布如图3 ZCb@!V}=  
| z_av  

B#.xs>{N  
图3.光束经过自适应镜校正之后的相位分布 -]Mk} z$  
&e#pL`N  
####光束传输至转换镜 +ut%C.1  
dist 5e7                           #光束传输500km至转换镜 g2*}XS3  
title 4: irradiance at relay mirror before aperture 3-n&&<  
plot/watch ex26_4.plt \IzZJ
Gi  
plot/liso xr=300 yr=300 nsl=64 Z=%u:K}[  
####光束到达转换镜前的相位分布如图4所示： v&%W*M0q@  
s>WqVuXmn  
bNtOqhi  
图4.光束到达转换镜前的相位分布 .L^;aL  
iEy2z+/"^  
####光束传输经过转换镜 #)#'^MZX  
clap/cir/con 1 75                  # aperture of 1.5 m. dia. for relay IM[=]j.?  
title 5: irradiance at relay mirror after aperture I\rjw$V#  
plot/watch ex26_5.plt `/wXx5n5<  
plot/liso xr=300 yr=300 nsl=64 A)d0Z6G`  
mirror 1 3e8 foc glKPjL*  
abr/ast 1 .3 90 E%stFyr9`/  
####光束经过转换镜后相位分布如图5所示： @y0kX<M  
8JW0;H<  
7)Tix7:9S;  
图5.光束经过转换镜后相位分布 d~QJ}a  
wLt0Fq6QG  
####光束传输至聚焦镜 %$Q!'+YW  
dist 3e8                           #光束传输3000km至聚焦镜 Y2N$&]O{  
clap/cir/con 1 200                 # 聚焦镜孔径直径4m L;b-=mF  
energy 98_os2`  
title 6: irradiance on focusing mirror YY'[PXP$Y  
plot/watch ex26_6.plt NA/Sv"7om  
plot/liso xr=300 yr=300 nsl=64 /^&$ma\  
####光束传输至聚焦镜相位分布如图6所示： _n4`mL8>kH  
!ueh%V Ky  
M$f_I +  
图6.光束传输至聚焦镜相位分布 I>-}ys`[  
a"1LF`  
####光束经过聚焦镜传输至聚焦靶面后相位分布 AP=SCq;  
phase/random 1 .1 50               # 聚焦镜附加随机像差 \S~<C[P  
mirror 1 5e7 focallength           # 聚焦镜焦距500km &qa16
bz  
dist 5e7                           # 传输至靶面 &;GoCU Le  
clap/cir/con 1 25                  # 设置出瞳孔径直径50cm ;[Esop  
energy +l7Bu}_?  
title 7: target irradiance inside 50 cm. /\1Q :B3W  
plot/watch ex26_7.plt uHDUuK:Ur  
plot/liso xr=50 yr=50 nsl=64 )coA30YR  
end S%7bM~J@  
####光束传输至聚焦靶面后相位分布如图7所示： v yP_qG  
4L,&a+)  
{z'Gg  
图7.光束传输至聚焦靶面后相位分布

不错，值得学习。