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

本例为一束激光从地基激光器出发传输到近地轨道上的一个转换镜上。转换镜将激光反射到一个聚焦镜上，然后这束光打到大气层内部一个低海拔的目标上。根据Kolmogorov模型，假设目标值的半径为10cm，就可以计算大气像差。本例中包含了激光扩束器像差、转换镜上的像差以及聚焦镜上的像差。 jWK>=|)=c  
大气模型假设的波阵面光谱功率为（忽略内部和外部的尺寸限制）： ;oOv~YB7H  
N`LY$U+N|  
~sTn?~  
]:Pkh./  
其中W^2 (f)是波阵面的光谱功率，r0为可视参数，f是空间频率，L0是外部尺寸，Li是内部尺寸，这些参数的单位分别为rad，m，m-1。由于大气像差和光束扩束器的像散，斯特列尔比SR=0.34。经过一个激励器影响半径为4.0cm的自适应镜校正后，斯特列尔比为0.87。经过全程传播到达目标后，光斑直径为50cm，剩下56%的能量，相对于没有自适应镜时能量的22%，有了明显的提高。由于没有考虑大气对光学元件散射效应的衰减效应，所以实际中传输到的光更加少。 !*}E  
CIQ9dx7>  
<o@__l.  
图1.地对空激光通信系统示意图 ?}No'E1!I  
x) R4_3  
_3{8Zg  
表1.关键参数 A s8IjGNs{  
9L>ep&u)^  
###激光器光束初始化 BLb'7`t  
set/alias/off c1 1?Kq  
wavelength/set 1 .48               # 设置激光器波长 jsq|K=x,  
array/set 1 256                    # 设置计算初始矩阵大小 8< z  
units/s 1 .1
|al'_s}I  
gauss/c/c 1 1 1.25                 # 定义高斯光束 ! ~+mf^D  
clap/cir/con 1 1.25 FB O_B  
energy/norm 1 1.                   # 能量归一化设置 bK|nxL  
set/density 64                     # 设置画图线条密度 _!k\~4U  
title 1: starting laser distribution e*39/B0S  
plot/watch ex26_1.plt 1r<'&f5  
plot/liso nsl=64 xr=1.5 yr=1.5       # 绘制激光束初始强度曲线 R".*dC,0'B  
####激光器初始光束相位分布如图1所示： &_3o1<  
)SfM`W)Y  
=!=DISPo  
图1激光器光束初始分布 *s!T$oc  
+Rq]_sDu  
##光束扩束器（20X）模拟 BWzo|isv  
mirror 1 20 focallength            # 扩束器透镜1焦距设置为20cm. 2`
V(w[zTr  
dist -420                          # 透镜分离 B";Dj~y  
mirror 1 400 focallength           #扩束器透镜2焦距设置为400cm. 1(/rg  
abr/ast 1 .2 45                    # 光束峰谷像散设置 I}\`l+  
clap/cir/con 1 25                  # 设置孔径光阑直径50 FwDEYG  
phase/random/kolmogorov 1 10. 7    # 光束附加大气像差 (!T\[6  
strehl z[0t%]7l  
title 2: phase after beam expander and atmoshpheric turbulence ;RW5XnVx  
plot/watch ex26_2.plt `'ak/%Krh  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 7qg. :h  
$~T|v7Y%  
ORt)sn&~d  
图2.经过扩束器及大气干扰后光束相位分布 Cp#}x1{  
T>m|C}yy  
####光束传输至自适应镜 Jy0(g T  
adapt 1. 4.                        # 自适应镜对光束的曲率半径影响为4cm D`PA
@t  
strehl D/tFN+|P  
title 3: phase after adaptive optic correction '^BTa6W}m  
plot/watch ex26_3.plt <KZ J  
plot/liso/phase xr=30 yr=30 nsl=64 min=-6 max=6 2 \m+;^_;5GW  
####光束经过自适应镜之后的相位分布如图3 4x" je  
=Ct$!uun  
_x<7^^VT  
图3.光束经过自适应镜校正之后的相位分布 CMYkxU  
*m&(h@l  
####光束传输至转换镜 IKnXtydeI}  
dist 5e7                           #光束传输500km至转换镜 l6T5]$  
title 4: irradiance at relay mirror before aperture ,sn 9&E  
plot/watch ex26_4.plt m"vWu0
/#  
plot/liso xr=300 yr=300 nsl=64 n ZzGak  
####光束到达转换镜前的相位分布如图4所示： eit%U  
JjHQn=3AJ  
xUDXg*  
图4.光束到达转换镜前的相位分布 3NrWt2?  
:qvaI,  
####光束传输经过转换镜 hFDo{yI  
clap/cir/con 1 75                  # aperture of 1.5 m. dia. for relay vVH*\&H\T  
title 5: irradiance at relay mirror after aperture *"j3x} U<  
plot/watch ex26_5.plt II(P  
plot/liso xr=300 yr=300 nsl=64 {y
|j**NZ  
mirror 1 3e8 foc Kk/cI6`W  
abr/ast 1 .3 90 <is%lx(GDX  
####光束经过转换镜后相位分布如图5所示： 8-q4'@(  
7G23D  
|$c~Jq  
图5.光束经过转换镜后相位分布 L_fiE3G|>  
9TV1[+JWe  
####光束传输至聚焦镜 j.UO>1{7  
dist 3e8                           #光束传输3000km至聚焦镜 k E-+#p  
clap/cir/con 1 200                 # 聚焦镜孔径直径4m EV N:3  
energy .Yxf0y?uv  
title 6: irradiance on focusing mirror 3mXRLx=0>  
plot/watch ex26_6.plt tnC,1HV0[  
plot/liso xr=300 yr=300 nsl=64 oqy}?<SQ  
####光束传输至聚焦镜相位分布如图6所示： "@@Z{  
e",0Er FT  
=kn-F T  
图6.光束传输至聚焦镜相位分布 7 JVonruaR  
hi9@U]H#  
####光束经过聚焦镜传输至聚焦靶面后相位分布 p}h9>R  
phase/random 1 .1 50               # 聚焦镜附加随机像差 O-]^_LV`  
mirror 1 5e7 focallength           # 聚焦镜焦距500km %s[ n2w  
dist 5e7                           # 传输至靶面 \_3#%%z  
clap/cir/con 1 25                  # 设置出瞳孔径直径50cm TRQH{O\O  
energy x%,!px3s  
title 7: target irradiance inside 50 cm. 1'9YY")#  
plot/watch ex26_7.plt i UW.$1l  
plot/liso xr=50 yr=50 nsl=64 JAI;7  
end D,, x<JG|  
####光束传输至聚焦靶面后相位分布如图7所示： R5MY\^H/A  
Z;J{&OJ3qM  
"9!CsloWhz  
图7.光束传输至聚焦靶面后相位分布

不错，值得学习。