-
Golay3阵列结构如图1所示。
Golay3光学系统的点扩散函数(point spread function, PSF)为:
$$ \begin{split} PSF( {x,y} ){\rm{ = }}&PS{F_{{\rm{one}}}}( {x,y} ) \cdot \Bigg\{ {3 + } \Bigg. 2{\rm{cos}}\Bigg[ \frac{{2\pi }}{{\lambda f}}\Bigg( \frac{\sqrt 3 ( {D - d} )}{2}x + \\ & \frac{{D - d}}{2}y \Bigg) \Bigg] + {2{\rm{cos}}\Bigg[ {\frac{{2\pi }}{{\lambda f}}\Bigg( { \pm \frac{{\sqrt 3 ( {D - d} )}}{{\rm{4}}}x + \frac{{{\rm{3}}( {D - d} )}}{{\rm{4}}}y} \Bigg)} \Bigg]} \Bigg\} \\ \end{split} $$ (1) 式中:
$PS{F_{{\rm{one}}}}\left( {x,y} \right)$ 为单孔径PSF;$d$ 为子孔径直径;$D$ 为外接圆直径;$\lambda $ 为入射光波长;$f$ 为系统焦距。Golay3光学系统的调制传递函数(modulation transfer function, MTF) [11]为:
$$ \begin{split} MTF\left( {{f_x},{f_y}} \right){\rm{ = }}&MT{F_{{\rm{one}}}}\left( {{f_x},{f_y}} \right) * \left[ {\delta \left( {{f_x},{f_y}} \right) + } \right. \\ &\dfrac{1}{3}\delta \left( {{f_x} \pm \dfrac{{\sqrt 3 \left( {D - d} \right)}}{{2\lambda f}},{f_y} \pm \dfrac{{D - d}}{{2\lambda f}}} \right) + \\ &\dfrac{1}{3}\delta \left( {{f_x} \pm \dfrac{{\sqrt 3 \left( {D - d} \right)}}{{4\lambda f}},{f_y} \pm \dfrac{{3\left( {D - d} \right)}}{{4\lambda f}}} \right) + \\ &\left. {\dfrac{1}{3}\delta \left( {{f_x} \mp \dfrac{{\sqrt 3 \left( {D - d} \right)}}{{4\lambda f}},{f_y} \pm \dfrac{{3\left( {D - d} \right)}}{{4\lambda f}}} \right)} \right] \\ \end{split} $$ (2) 式中:
$MT{F_{{\rm{one}}}}\left( {{f_x},{f_y}} \right)$ 表示单孔径MTF;$({f_x},{f_y})$ 表示空间频率。由公式(2)可知,当子孔径间存在相位误差时,孔径间的相对位置关系不再对称,导致系统的MTF发生改变,成像质量下降。
子孔径间的相位误差对光学系统成像质量影响非常大,在实验调试阶段必须对其进行严格校正。图2给出了Golay3结构中子孔径3存在活塞误差或倾斜误差时的MTF分布图以及相应的面物体成像。
图 2 MTF分布((a)无误差;(b)活塞误差π;(c)倾斜误差0.012 μrad)和面物体成像((d)无误差;(e)活塞误差π;(f)倾斜误差0.012 μrad)
Figure 2. MTF distribution( (a) Error free; (b) Piston error π; (c) Tilt error 0.012 μrad) and surface object imaging ((d) Error free; (e) Piston error π; (f) Tilt error 0.012 μrad)
由图2可以看出,无误差的多孔径MTF分布呈中心对称分布,中心能量最大,周围对称分布六个旁瓣,这与子孔径的排列方式有关。当最大活塞误差为π时,系统MTF偏移量最大。倾斜误差不改变MTF中心位置,但使得MTF分布地更加杂乱无章。多孔径子光束间存在活塞误差和倾斜误差都会导致成像结果模糊。所以,使用多孔径光学系统观测前必须对子孔径进行调整,以排除系统自身的活塞误差和倾斜误差。
-
湍流是大气的一种运动形式,是一种复杂的现象,在湍流大气中,折射率随时间、空间随机变化。大多数情况下,湍流会充满整个视场空间,但由于大气环境复杂,湍流的尺度变化很大,因此还会存在仅影响视场部分区域的小尺度湍流。局部小尺度湍流虽然没有充满整个视场,但是仍会对目标观察产生干扰。文中重点针对这种外尺度较小的湍流情况进行分析。
实际观测中若多孔径望远系统全视场角为0.2°,观测距离为10 km目标,则子孔径可观测到的视场范围约为350 m。当成像光路中出现近似为团状湍流,湍流外尺度小于子孔径的视场时,则望远系统视场中有一部分近似圆形的区域入射光相位值与其他部分不同发生了畸变。为了分析近似圆形湍流对成像系统的影响,对其进行了简化模型构建,如图3(a)中阴影部分所示。
图 3 (a)近似圆形湍流在视场中的示意图及(b)湍流区域相位分布
Figure 3. (a) Schematic diagram of similar circular turbulence in the field of view and (b) phase distribution in turbulent region
近似圆形湍流引起的相位差由以下方法计算。根据参考文献[12]使用各向异性大气湍流功率谱
${\varPhi _n}({\kappa _x},{\kappa _y})$ :$${\varPhi _n}({\kappa _x},{\kappa _y}) = \dfrac{{A(\alpha )\overline C_n^2{\mu _x}{\mu _y}}}{{{{(\mu _x^2\kappa _x^2 + \mu _y^2\kappa _y^2)}^{{\alpha / 2}}}}}$$ (3) 式中:
${\mu _x}$ 和${\mu _y}$ 分别表示$x$ 和$y$ 方向各向异性分子;$({\kappa _x},{\kappa _y})$ 表示空间频率;$\bar C_n^2$ 表示折射率广义结构常数,单位为${\rm{m} ^{3 - \alpha }}$ ,是表示湍流强度的一个重要参数;$A(\alpha )$ 表示保持其指数结构函数与其功率谱一致性的常数。$$A(\alpha ) = \frac{1}{{4{\pi ^2}}}\varGamma (\alpha - 1)\cos ({{\alpha \pi } / 2})\quad (3 < \alpha < 4)$$ (4) 式中:
$\varGamma (x)$ 表示$\gamma $ 函数;$\alpha $ 表示与高度有关的光谱指数。通过对大气湍流的功率谱进行滤波,再做逆傅里叶变换可以计算出湍流影响的相位值[13-14]分布如图3(b)所示。
随着湍流随机运动,如果空间存在一个分布方向正好垂直于视场的近似长条状湍流,例如对流湍流中山波湍流[15-16],望远系统视场中有一部分近似长条状的区域入射光相位值与其他部分不同,发生了畸变,对其进行简化模型构建,如图4(a)阴影部分所示。
图 4 (a)近似长条状湍流在视场中的示意图;(b)湍流气体密度分布及(c)相位分布
Figure 4. (a) Schematic diagram of similar elongated turbulence in the field of view; (b) Turbulent gas density distribution and (c) phase distribution
文中近似长条状湍流引起的相位差由以下方法计算。对流湍流的分析模型需要结合空气的流体动力学情况来建立,利用Fluent软件计算强风通过山脉产生的山波湍流,得到湍流区域的密度分布,如图4(b)所示,再利用Gladstone-Dale关系:
$$n = 1 + K\rho $$ (5) 式中:
$n$ 为折射率;$\;\rho $ 为密度;$K$ 为Gladstone-Dale系数[17]。$$K = 2.23 \times {10^{ - 4}}\left( {1 + \frac{{7.52 \times {{10}^{ - 3}}}}{{{\lambda ^2}}}} \right)$$ (6) 式中:
$\lambda $ 为入射光波长。计算出湍流区域的折射率场分布,进而可求出光场经过近似长条状湍流的相位分布,如图4(c)所示。
-
近似圆形湍流影响下的系统PSF可表示为:
$$ \begin{split} PSF\left( {x,y} \right) =& \left| {\sum\limits_{n = 1}^3 {\left\{ {\left( {\frac{{\pi {{{d}}^2}}}{{4\lambda f}}} \right)\left[ {\frac{{2{J_1}\left( {{{\pi {{d}}{r_n}} / {\lambda f}}} \right)}}{{{{\pi {{d}}{r_n}} / {\lambda f}}}}} \right]} \right.} } \right. + \\ &{\left. {\left( {\frac{{\pi {l^2}}}{{4\lambda f}}} \right)\left[ {\frac{{2{J_1}\left( {{{\pi l{r_{ln}}} / {\lambda f}}} \right)}}{{{{\pi l{r_{ln}}} / {\lambda f}}}}} \right] \cdot \left( {{{\rm{e}}^{i\varphi }} - 1} \right)} \right|^2} \\ \end{split} $$ (7) 式中:
${J_1}$ 表示一阶贝塞尔函数;${r_n} = \sqrt {x_n^2 + y_n^2} $ ,$({x_n},{y_n})$ 表示子孔径中心坐标;$d$ 表示子孔径直径;${r_{ln}} = \sqrt {x_{ln}^2 + y_{ln}^2} $ ,$({x_{{l_n}}},{y_l}_{_n})$ 表示近似圆形湍流中心坐标;$\varphi $ 表示湍流影响附加相位差。后文中湍流影响导致的附加相位差都用相位差表示。公式(7)等号右边绝对值内第一、二项分别表示没有湍流影响和有湍流影响的振幅扩散函数,振幅扩散函数模平方为PSF。湍流变化引起的大气折射率起伏导致光波波前畸变,改变望远镜系统视场内相位均匀性,造成光学图像的模糊,影响光学系统PSF,进而恶化系统MTF。
系统的MTF为:
$$MTF\left( {{f_x},{f_y}} \right) = \int {\int_{ - \infty }^\infty {PSF\left( {x,y} \right)} } {{\rm{e}}^{ - i2\pi \left( {{f_x}x + {f_y}y} \right)}}{\rm{d}}x{\rm{d}}y$$ (8) 通过改变近似圆形湍流在视场中的不同位置,计算得出湍流位于子孔径边缘对系统成像性能影响最大。并且经过计算发现,当湍流区域直径小于d/6时,系统MTF分布基本不变,即尺度较小的湍流对系统成像影响不大。
经过计算发现,湍流在视场中的位置和尺度大小会影响系统MTF,但其主要的影响规律相同。如图3(a)所示,近似圆形湍流直径为d/3,湍流尺度约为100 m,取湍流层的传输距离为512 m,分析湍流变化对系统成像的影响。
图5给出了近似圆形湍流区域相位变化
$0 \sim 2\pi $ 的MTF分布以及相应的面物体成像和不同相位差的MTF曲线。图 5 近似圆形湍流对系统的MTF影响(相位差(a) π/2 (3π/2);(b) π)及面物体成像(相位差(c) π/2 (3π/2);(d) π;(e)不同相位差MTF曲线)
Figure 5. Influence of similar circular turbulence on MTF of system ((a) π/2 (3π/2); (b) π) and plane object imaging phase difference ((c) π/2 (3π/2); (d) π; (e) MTF curves with different phase differences)
由图5可知,当近似圆形湍流影响改变,系统MTF中心位置并没有发生移动,随着相位差由0增加到π,MTF畸变程度越来越严重,中心由圆形变成椭圆,椭圆两尖端逐渐分裂开,周围的六个旁瓣范围逐渐缩小。从图5(e)定量分析,以相位差为
$0\pi (2\pi )$ 的MTF曲线为基准,相位差为${\pi / 2}({{3\pi } / 2})$ 和$\pi $ 的MTF分别在归一化空间频率0.16处下降了0.06和0.19,在频率0.45处下降了0.04和0.09,在频率0.69下降了0.034和0.08。另外,当湍流强度为
$\overline C_n^2 = {10^{ - 8}}\;{{\rm{m}}^{3 - \alpha }}$ ,距离入瞳距离分别为h=8、10、15、20 km时,由参考文献[12]可知当湍流强度一定时,湍流导致的相位差值与湍流层高度有关,所以在近似圆形湍流区域构建不同高度的随机相位屏,相位分布的值放入512×512的矩阵中,矩阵中每点的相位差值可表示对应空间位置的相位起伏,以此来研究近似圆形湍流影响下的MTF。便于分析计算出对应随机相位屏的相位均值分别为70.82 π、8.09 π、0.61 π、0.027 π。由图6可知,当高度为8、10、15 km时,MTF的畸变比较明显,中央类似一个椭圆形,周围六个旁瓣位置并没有发生明显变化。当高度为20 km时,圆形湍流对MTF分布的影响不大。
图 6 不同高度湍流影响的MTF ((a) 8 km; (b) 10 km;(c) 15 km; (d) 20 km)
Figure 6. MTF affected by turbulence at different heights ((a) 8 km; (b) 10 km; (c) 15 km; (d) 20 km)
当湍流高度为8 km时改变湍流强度,令
$ {\overline{C}}_{n}^{2}= {\rm{10}}^{\rm{-8}}\rm{ }{\text{、}} $ $ {\rm{ 10}}^{\rm{-13}}\rm{ }{\text{、}}{\rm{ 10}}^{\rm{-15}}\;{\rm{m}}^{3-\alpha }$ ,在近似圆形湍流区域构建不同湍流强度的随机相位屏,并计算出对应随机相位屏的相位均值分别为70.82 π、1.11 π、0.14 π。图7(d)所示MTF曲线为过MTF分布中心且平行于
${f_x}$ 方向的数据,结合随机湍流相位屏均值,可知同一高度湍流强度越小对系统的MTF影响越小。以理想情况的MTF曲线为基准,湍流强度最小的MTF曲线和理想情况相似,而另外两条强度大的MTF曲线,在归一化空间频率0.16、0.45和0.69处都有所下降。强度为${\rm{1}}{{\rm{0}}^{{\rm{ - 13}}}}\;{{\rm{m}}^{3 - \alpha }}$ 和${\rm{1}}{{\rm{0}}^{{\rm{ - 8}}}}\;{{\rm{m}}^{{\rm{3 - }}\alpha }}$ MTF曲线分别在频率0.16处下降0.06和0.09,0.45处下降0.04和0.05,0.69处下降0.03和0.04。 -
近似长条状湍流影响下的系统PSF可表示为:
$$ \begin{split} PSF\left( {x,y} \right) =& \left| {\sum\limits_{n = 1}^3 {\left\{ {\left( {\frac{{\pi {d^2}}}{{\lambda f}}} \right)\left[ {\frac{{2{J_1}\left( {{{\pi d{r_n}} / {\lambda f}}} \right)}}{{{{\pi d{r_n}} / {\lambda f}}}}} \right]} \right.} } \right. + b\sqrt {{d^2} - {b^2}} {\rm{sin}}c\cdot \\ & {\left. {\left( {\frac{{x - {x_{ln}}}}{{\lambda f}}b,\frac{{y - {y_{ln}}}}{{\lambda f}}\sqrt {{d^2} - {b^2}} } \right) \cdot \left( {{{\rm{e}}^{i\varphi }} - 1} \right)} \right|^2} \\ \end{split} $$ (9) 式中:
$({x_{{l_n}}},{y_l}_{_n})$ 表示近似长条状湍流中心坐标;$b$ 表示近似长条状湍流宽度。公式(9)等号右边绝对值内第一、二项分别表示没有湍流影响和有湍流影响的振幅扩散函数。通过改变近似长条状湍流在视场中的位置,计算得出位于子孔径中心的湍流,能更好地突出对光学系统性能的影响。并且计算发现,湍流区域宽度越窄,对系统MTF分布变化影响越大。
如图4(a)所示,近似长条状湍流宽度为
${d / {10}}$ ,取湍流层的传输距离为512 m,分析湍流变化对系统成像的影响。图8给出了近似长条状湍流区域相位变化$0 \sim 2\pi $ 的MTF分布以及相应的面物体成像和不同相位差的MTF曲线。图 8 近似长条状湍流对系统MTF影响(相位差(a) π/2 (3π/2);(b) π面物体成像相位差;(c) π/2 (3π/2);(d) π;(e)不同相位差的MTF曲线
Figure 8. Influence of similar elongated turbulence on MTF of system ((a) π/2 (3π/2); (b) π plane object imaging phase difference; (c) π/2 (3π/2); (d) π; (e) MTF curves with different phase differences
由图8可知,当近似长条状湍流引起相位差逐渐改变时,MTF的中心区域没有移动,但中心区域范围减小,中心亮斑和旁瓣产生线条状的畸变,除中心处外,线条处的能量最大。图8(c)、(d)为面物体成像,相比近似圆形湍流,近似长条状湍流对系统的成像质量影响更加严重。从图8(e)定量分析,以相位差
$0\pi $ $ (2\pi )$ 为基准,相位差${\pi / 2}({{3\pi } / 2})$ 的MTF在归一化空间频率0.17处下降0.537,0.39处(MTF为0)下降0.28,0.69处下降0.144;相位差π的MTF畸变程度非常大,MTF值在频率0.078、0.21、0.34、0.48、0.63和0.75处下降为0,同时MTF除了在中央零频有极大值,分别在频率0.14、0.27、0.42、0.54、0.68和0.81处有极大值,分别为0.21、0.40、0.14、0.33、0.028和0.044。由此可知,近似长条状湍流对系统的MTF影响非常明显。另外,利用Fluent软件计算光波长550 nm、风速分别为10、20、30 m/s的强风通过山脉产生的山波湍流,再通过计算得到近似长条状湍流区域的相位差分布,以此来研究近似长条状湍流对系统MTF的影响。便于分析计算出不同风速湍流对应相位差的相位均值分别为9.39 π、9.54 π、9.63 π。图9给出不同风速影响下的MTF分布和MTF曲线。
图 9 不同风速影响的MTF ((a) 10 m/s;(b) 20 m/s;(c) 30 m/s)及(d)不同风速影响的MTF曲线
Figure 9. MTF affected by different wind speeds ((a) 10 m/s; (b) 20 m/s; (c) 30 m/s) and (d) MTF curves of different wind speed influence
由图9结合Fluent仿真,对流湍流的密度分布在两种气流的交界处变化梯度非常大,进而不同风速下近似长条状湍流,除了交界处相位差有明显差异,其余湍流域相位差相差不大,因此MTF只在较小空间频率处有明显的差异。由图9(d)MTF曲线定量分析,以理想情况的MTF为基准,风速为10、20、30 m/s的MTF曲线在归一化空间频率0.16处分别下降0.23、0.24、0.25,在频率0.42处均下降约0.089,在频率0.69处均下降0.046。同时三条MTF曲线除在中央零频有极大值,分别在频率0.25和0.50处出现极大值,分别为0.37和0.32。
由图3和图4的近似湍流示意图分析,成像系统实际上是一个低通滤波器,根据空域与频域的对应关系可知,通光孔径不同空间位置将对应不同的空间频率范围。在近似圆形和近似长条状湍流面积占成像系统视场比相差不大情况下,近似圆形湍流只对目标图像某一频段产生影响,近似长条状湍流由于贯穿了整个孔径,因此对目标图像低、中、高频整个频段产生影响。同时由图5和图8分析,附加相位差相同时,近似长条状湍流影响的MTF变化更加明显。由图7和图9MTF曲线分析,除去多孔径系统附加相位差2 π周期效应,近似圆形湍流强度为
${\rm{1}}{{\rm{0}}^{{\rm{ - 8}}}}\; {{\rm{m}}^{{\rm{3 - }}\alpha }}$ 的有效均值相位差为0.82 π,近似长条状湍流风速为20 m/s的有效相位差为0.54 π,虽然近似圆形湍流的有效相位差较大,但是近似长条状湍流对MTF的影响更显著。
Influence of turbulence on imaging characteristics of multi-aperture optical system
-
摘要: 为了研究大气湍流变化对多孔径光学系统成像质量的影响,针对Golay3结构建立了一个理论分析模型,推导了湍流影响下望远系统的点扩散函数(point spread function, PSF)的表达式。分别针对近似圆形和近似长条状分布的湍流结构对成像特性的影响进行了讨论,具体构建了两种湍流影响下的光场相位计算模型,并计算比较了两种湍流的调制传递函数(modulation transfer function, MTF)。结果表明,光路中存在湍流将导致系统成像质量下降。在近似圆形湍流中,湍流强度越小,系统MTF影响越小。不同强度湍流影响的MTF在归一化空间频率0.16、0.45和0.69处均下降0.05左右。在近似长条状湍流中,对流湍流风速越小,系统MTF影响越小。不同风速的湍流影响的MTF在空间频率0.16、0.42和0.69处分别下降0.25、0.09和0.05左右。比较两种湍流表明,近似长条状湍流对成像系统MTF影响更为明显。Abstract: In order to study the influence of atmospheric turbulence on the imaging quality of multi-aperture optical system, a theoretical analysis model was established for the Golay3 structure. The expression of point spread function (PSF) and modulation transfer function (MTF) of this structure were derived. To discuss the influence of similar circular and elongated turbulence structures on the imaging characteristics, the MTF under two turbulence conditions were calculated and compared, and the phase calculation models were constructed. The results show that when there is presence of turbulence in the optical path, the image quality of the system will decrease. As for similar circular turbulence, the smaller the intensity of turbulence is, the less the impact on the MTF is. Specially, the MTF affected by turbulence of different intensities decreases by about 0.05 at the normalized spatial frequencies of 0.16, 0.45, and 0.69. As for similar elongated turbulence, the smaller turbulent wind speed is, the less the impact on the MTF is. The MTF affected by turbulence of different wind speeds decreases by about 0.25, 0.09, and 0.05 at the spatial frequencies of 0.16, 0.42, and 0.69, respectively. The comparison of the influence of the two turbulences on the system shows that the similar elongated turbulence has a more obvious influence on the MTF.
-
-
[1] Golay M J E. Point arrays having compact, nonredundant autocorrelations [J]. JOSA, 1971, 61(2): 272-273. doi: 10.1364/JOSA.61.000272 [2] Lampkin C M, Flint G W, MacFarlane M J. Optical and mechanical design aspects of a four-telescope array for combined imaging [J]. Opt Eng, 1988, 27(9): 749-754. [3] Christopher R, De Hainaut K P, Hentz L D, et al. Design of a wide field of view phased array telescope [J]. Opt Eng, 1988, 27(9): 736-739. [4] Flores Jorge L, Strojnik Marija, Páez Gonzalo, et al. Effects of misalignment errors on the optical transfer functions of synthetic aperture telescopes [J]. Applied Optics, 2004, 43(32): 5926-5932. doi: 10.1364/AO.43.005926 [5] Kendrick R L, Aubrun J N, Bell R, et al. Wide-field fizeau imaging telescope: Experimental results [J]. Appl Opt, 2006, 45(18): 4235-4240. doi: 10.1364/AO.45.004235 [6] Ni M, Benson L, Camp J, et al. Autonomous tip/tilt alignment and phasing of adistributed aperture imaging testbed [J]. Opt Express, 2010, 18(12): 13051-13056. doi: 10.1364/OE.18.013051 [7] 吴泉英. 稀疏孔径光学系统成像研究[D]. 苏州大学, 2006. Wu Quanying. Research on imaging of sparse aperture optical system[D]. Suzhou: Soochow University, 2006. (in Chinese) [8] Zhang Wei, Deng Jian, Long Funian, et al. Research on image quality evaluation of sub-aperture synthetic optical imaging system [J]. Acta Optica Sinica, 2008, 28(4): 687-691. (in Chinese) doi: 10.3321/j.issn:0253-2239.2008.04.014 [9] Sun Jiangqin, Qian Lin, Wu Quanying. Modulation transfer function and imaging properties of axisymmetric synthetic aperture optical system [J]. Applied Optics, 2008, 29(4): 548-552, 564. (in Chinese) [10] Xie Zongliang, Ma Haotong, Li Bo. Experimental demonstration of enhanced resolution of a Golay3 sparse-aperture telescope [J]. Chinese Optics Letters, 2017, 15(4): 34-37. [11] Fiete R D, Tantalo T A, Calus J R, et al. Image quality of sparse-aperture designs for remote sensing [J]. Opt Eng, 2002, 41(8): 1957-1969. doi: 10.1117/1.1490555 [12] Tang Shanfa, Yan Jiawei, Yong Kangle, et al. Propagation characteristics of vortex beams in anisotropic atmospheric turbulence [J]. Opt Soc Am, 2020, 37(1): 133-137. doi: 10.1364/JOSAB.37.000133 [13] Zhang Huimin, Li Xinyang. Research on numerical simulation method of atmospheric turbulence distortion phase screen [J]. Optoelectronic Engineering, 2006, 33(1): 14-19. (in Chinese) doi: 10.3969/j.issn.1003-501X.2006.01.004 [14] 张滔. 大气湍流中涡旋光束轨道角动量谱的扩散及修正的研究[D]. 电子科技大学, 2017. Zhang Tao. Diffusion and correction of orbital angular momentum spectrum of vortex beams in atmospheric turbulence[D]. Chengdu: University of Electronic Science and Technology of China, 2017. (in Chinese) [15] Victor Klimenko, Jimmy Krozel. Clear-air turbulence impact modeling based on flight route analysis [C]//AIAA Guidance, Navigation, and Control Conference, 2011, 6: 5199-5211. [16] 范源丹. 基于垂直载荷因子的晴空湍流检测方法研究[D]. 中国民航大学, 2019. Fan Yuandan. Study on the clear-air turbulence detection method based on vertical load factor[D]. Tianjin: Civil Aviation University of China, 2019. (in Chinese) [17] 曹原. 喷流冷却条件下高速飞行器光学窗口气动特性研究[D]. 哈尔滨工业大学, 2020. Cao Yuan. Aerodynamic characteristics of high speed aircraft optical window under jet cooling condition[D]. Harbin: Harbin Institute of Technology, 2020. (in Chinese)