-
文中激光通信终端粗瞄机构(Coarse scanning mechanism,CSM)选用高精度力矩电机,精瞄机构(Refined scanning mechanism,RSM)采用压电陶瓷偏摆镜(Piezoelectric ceramic mirror, PZT)。总体技术参数如表1所示,指标满足研究要求。根据卫星平台的姿态控制与安装定位指标要求,开展扫描不确定域[13](Field of uncertainty, FOU)的计算。
表 1 总体仿真参数
Table 1. Overall simulation parameters
Parameter Value 1 Communication distance/km ≥3000 2 Satellite attitude stability/(°)·s-1 ≤0.02(3σ) 3 Attitude measurement accuracy/(°) 0.01 4 Satellite positioning error/m ≤50(3σ) 5 Signal light wavelength λ/nm 1550 6 Field of vision/mrad ±2.5 7 Signal light divergence angle α/μrad 50-70 8 CSM control frequency fmotor/ Hz ≤10 9 CSM control error δmotor/ μrad ≤30(3σ) 10 RSM control frequency fPZT/ Hz ≤800 11 RSM control error δPZT/ μrad ≤8(3σ) 12 Camera refresh rate fCCD/ Hz ≥120 卫星姿态数据由星上综合电子通过总线进行分发,姿态更新周期为1 Hz,根据表1数据显示,由测量精度和姿态变化引入的误差为524 μrad(3σ)。考虑到双星链路的距离,由卫星定位误差造成的最大偏差约为50 μrad(3σ)。此外,终端装星误差、平台轴系误差、回转误差带来的偏差为365 μrad(3σ)。综上,由卫星姿态稳定度、测量精度和终端安装定位引入的总偏差约为940 μrad,近似为1 mrad。因此,计算得到FOU范围为±1 mrad,大小为2 mrad。
根据表1所示,终端相机刷新频率高,并且终端探测视域(Field of vision, FOV)大于FOU,终端FOV可覆盖FOU。因此,文中扫描方式选用连续扫描[14],即终端在扫描过程中无需在每个扫描点驻留足够长的时间来等待相机探测与处理对方卫星终端发来的光信号。进一步地,双星无信标终端的捕跟方式采用扫描-凝视方式,如图1所示,终端A和B分别为扫描端和凝视端,捕跟流程具体如下:
(1)根据预定控制参数,终端A与终端B分别发出信号光并指向对方卫星终端。受姿轨预测精度和控制误差等因素影响,终端瞄准指向方向在FOU内随机分布,双方连线在FOU内;
(2)终端A按照预定扫描方法在不确定区域内由内至外扫描,终端B保持凝视状态,一旦终端A发来的信号光扫描到终端B相机视域内,则立刻被终端B捕获;
(3)终端B根据捕获到的光斑在其相机视域中心的偏移量,迅速调整位置指向终端A光束发来的方向。由于光束是双向可逆的,此时,终端A相机也能探测到终端B调整后的回光;
(4)终端A一旦探测到终端B的回光信号,则完成扫描捕获过程,双方进入到链路跟踪阶段。
-
根据图1可知,扫描捕获过程实际从捕跟步骤(2)开始,无信标扫描捕获总时长Ttotal由三部分组成,表达式如下:
$$ {T_{total}} = {T_{S tep 2}} + {T_{S tep 3}} + {T_{S tep 4}} $$ (1) 式中:TStep2、TStep3、TStep4分别为步骤(2)~(4)的时长。
与步骤(2)相比,步骤(3)、(4)的耗时时长较短。步骤(3)时长主要由凝视端B转台偏转调整耗时、光束空间传输时延和扫描端A相机探测处理时延组成。转台偏转调整时间由电机控制带宽决定,探测处理时延由相机刷新频率决定,假设2~3个电机控制带宽内转台完成偏转调整,步骤(3)用时如公式(2)所示:
$$ {T_{S tep3}} = n \cdot {({f_{motor}})^{ - 1}} + \frac{D}{c} + {({f_{CCD}})^{ - 1}} $$ (2) 式中:fmotor为电机控制频率;D为双星空间距离;c为光束;fCCD为相机刷新频率;n为电机控制带宽数量,取值2~3。
步骤(4)的时长主要由终端B机构微调耗时组成,假设一个电机和PZT控制带宽内可完成微调,步骤(4)用时表达式如下:
$$ {T_{S tep4}}{\text{ = }}{({f_{PZT}})^{ - 1}} + {({f_{motor}})^{ - 1}} $$ (3) 式中:fPZT为PZT控制频率,由于PZT需要和相机进行匹配控制,因此,PZT控制频率与相机刷新频率保持相同。
根据公式(2)、(3)计算,步骤(3)和步骤(4)总耗时约为0.4 s,说明无信标扫描捕获时间主要由步骤(2)决定。
扫描端A的扫描方法是影响步骤(2)耗时的主要因素。首先,采用在星间激光通信领域被广泛研究应用的扫描方法——等步长螺旋扫描方法[15],以几十微弧度束散角的信号光在FOU内进行全覆盖扫描仿真分析。考虑星上微振动环境,单步扫描步长Iθ与信号光束散角α关系如图2所示,其表达式如公式(4)所示:
$$ {I_\theta } = \frac{{\sqrt 2 }}{2}\alpha - x $$ (4) 式中:x为补偿星上微振动设置的重叠量,大小与星上微振动幅值有关。
基于等步长螺旋扫描方法,束散角假设为60 μrad,全场螺旋扫描点超过2300个,如图3所示。
单步扫描时长由执行机构控制带宽决定,若按照传统方法,将电机作为唯一的扫描执行机构,则单步扫描时间约为0.12 s,全场扫描时间超过276 s (4.6 min),无法对光斑进行快速捕获。因此,实现无信标快速扫描捕获主要有以下两种技术方案途径:
(1)利用PZT控制频率高的特点,采用PZT作为终端唯一扫描执行机构,在FOU内完成高频快速全覆盖扫描;
(2)利用电机扫描范围大和PZT控制频率高的特点,采用粗精机构复合扫描的方法在FOU内完成快速全覆盖扫描。
将方案(1)应用到文中,考虑终端光学压缩率,PZT镜面与光路安装夹角、PZT镜面在光路中的安装位置等因素。PZT最大扫描走点距离lPZT,如公式(5)所示:
$$ {l_{PZT}} = \frac{{{\varphi _{PZT}} \cdot \varepsilon }}{\eta } $$ (5) 式中:φPZT为PZT最大偏摆角度,取值1~1.5 mrad;η为光路压缩率,取值8~10;ε为光学反射传递系数,数值与镜面与光路夹角和镜面安装位置有关,文中取值为2。
根据公式(5),文中终端PZT最大走点距离不超过300 μrad,不能在FOU内单独完成全场扫描。基于方案(1),PZT最大偏摆角度要超过4 mrad才能实现全场快速扫描功能,然而该类PZT尺寸较大,需要更大体积的光学回转体进行内嵌安装,从而增加了激光通信终端的体积重量,这对终端轻小型化设计产生不利影响。因此,方案(1)无法适用于文中无信标扫描方法的应用研究。
根据上述分析,方案(2)同样可实现无信标快速扫描捕获功能,并且能适应终端低功耗、轻小型的设计要求。因此,文中以上述等步长螺旋扫描为基础,开展无信标粗精复合螺旋扫描方法设计工作。
-
根据方案(2),文中粗精复合螺旋扫描方法的设计原则为:在FOU内规划大步长/大螺距螺旋线,利用电机转动范围大的优势,使电机按照螺旋点进行大角度走点;以螺旋点为中心,利用PZT控制频率高的特点,使PZT快速扫描遍历各螺旋点周围区域。所提方法主流程如图4所示。
具体包括如下步骤:
(1)以卫星终端初始指向角度为扫描原点,终端方位角度和俯仰角度为坐标,建立二维直角坐标系,在二维坐标系内设置大螺距等步长螺旋扫描线,各螺旋点为电机目标走点,以各螺旋点为中心将FOU划分为多个扫描子区域;
(2)在扫描子区域内以螺旋点为原点建立子区域极坐标系;
(3)规划PZT在子区域内的扫描路径,将PZT扫描路径坐标点参数从子区域极坐标系依次转换到终端极坐标系、终端二维坐标系,最终得到PZT俯仰和方位角度实际控制量;
(4)控制PZT在子区域内以实际角度扫描,如扫描到对方卫星终端并捕获到对方回光,则停止;如没有扫描到对方卫星终端并捕获到对方回光,则继续控制PZT在子区域内完成全覆盖扫描;
(5)待PZT完成子区域扫描,则控制电机按等步长螺旋线方式走点至下一个螺旋点,如电机走点后扫描到对方卫星终端视域并捕获到对方回光,则停止;如电机没有扫描到对方卫星终端视域并捕获到对方回光,则重复步骤(3)、(4),即控制PZT在下一个子区域内进行扫描;
(6)当电机完成全部的螺旋线走点,即完全遍历覆盖扫描FOU,仍未扫描到对方卫星终端视域并捕获到对方回光,则停止,重复步骤(1),重新进行粗精复合扫描。
螺旋线步长与子区域扫描路径是影响所提复合螺旋扫描方法扫描捕获性能的重要参数。针对上述所提方法,对目标参数进行数值计算与优选。
FOU内的等步长螺旋线总长度l表达式如下:
$$ l = \frac{L}{{4\pi }}\left[ {\theta \sqrt {1 + {\theta ^2}} + \ln \left( {\theta + \sqrt {1 + {\theta ^2}} } \right)} \right] $$ (6) 式中:L为螺旋点步长;θ螺旋角,大小为2π;l为螺旋线总长度。当θ大于π时,可近似简化为:
$$ l \approx \frac{{\pi r_{FOU}^2}}{L} = n \cdot L $$ (7) 式中:rFOU为FOU半角;n为螺旋点总数。
螺旋点在FOU的极坐标表达式为:
$$ \left\{ {\begin{array}{*{20}{c}} {{\rho _i} = \sqrt {\dfrac{i}{\pi }} \cdot L} \\ {{\theta _i} = 2\sqrt {i \cdot \pi } } \end{array}} \right.{\text{ }}i \gt 2 $$ (8) 式中:ρi为螺旋点极径;θi为螺旋点极角;i为第i个螺旋点。将螺旋点极坐标表达式转换为终端二维坐标系表达式,如公式(9)所示:
$$ \left\{ {\begin{array}{*{20}{c}} {{x_i} = \sqrt {\dfrac{i}{\pi }} \cdot L\cos \left( {2\sqrt {i \cdot \pi } } \right)} \\ {{y_i} = \sqrt {\dfrac{i}{\pi }} \cdot L\sin \left( {2\sqrt {i \cdot \pi } } \right)} \end{array}} \right.{\text{ }}i \gt 2 $$ (9) 式中:xi为第i个螺旋点方位轴角度;yi为第i个螺旋点俯仰轴角度。
由于采用电机进行螺旋线走点,因此,扩大螺旋点步长,减少电机螺旋线走点数量,可降低整体扫描时间。然而,在扩大螺旋点步长的同时,每个螺旋点子区域的面积也在增加。子区域面积和受PZT最大走点距离和子区域扫描路径的制约,并且PZT在扫描阶段结束后的跟踪阶段还担负着微调终端对准光斑的功能。因此,在考虑PZT最大走点距离和最大偏摆角度范围的前提下需对PZT在扫描和跟踪阶段的偏摆角度进行分配,如公式(10)所示:
$$ {\varphi _{PZT}} = {\varphi _{PZT\_scan}} + {\varphi _{PZT\_tracking}} $$ (10) 式中:φPZT_scan为扫描阶段PZT最大偏摆角度;φPZT_tracking跟踪阶段PZT最大偏摆角度。
文中假设φPZT_scan与φPZT_tracking平均分配PZT最大偏摆角度量程φPZT,将φPZT_scan代入公式(5)进行计算,得到所提复合扫描方法的PZT实际最大走点距离lPZT小于150 μrad。因此,在考虑lPZT不超过150 μrad的前提下,PZT的走点区域被约束在以螺旋点为中心,半径小于lPZT的圆内。同样地,螺旋点子区域则是以螺旋点为中心的中心对称区域。
光栅-螺旋方法是一种高效遍历中心对称图形区域的规划方法[16],文中采用光栅-螺旋方法在子区域内进行扫描路径规划。基于该方法,以PZT最大走点距离lPZT为约束条件,PZT在子区域内最多可扫描两圈,共24步,如图5所示。根据上述分析可知,螺旋点步长越大,则电机走点数量越少,整体扫描时间也越少。根据这一论述,在满足PZT实际走点距离小于lPZT的前提下,文中选取PZT在子区域内的光栅-螺旋扫描走点数量为24,保证电机螺旋点数量最少,整体扫描耗时最少。
子区域极坐标系下所规划的24个光栅-螺旋扫描路径坐标点矩阵Nmatrix如下:
$$ N_{\text {matrix }}(:, 1)=\rho_{N_{\text {marrix }}}^{\prime}, N_{\text {matrix }}(:, 2)=\theta_{N_{\text {marix }}}^{\prime} $$ (11) 式中:ρ′Nmatirx为子区域极坐标系下PZT所有扫描点的极径;θ′Nmatirx为子区域极坐标系下PZT所有扫描点的极角,如图5所示。
考虑电机控制误差δmotor和大螺距造成的相邻子区域径向漏扫,螺旋点步长L需设置重叠量β,如图6所示,其表达式为:
$$ \left\{ {\begin{array}{*{20}{c}} {L = 5{I_\theta } - \beta } \\ {\beta {\text{ = }}{\delta _{motor}}{\text{ + }}\varepsilon } \end{array}} \right. $$ (12) 式中:δmotor为电机控制误差;ε为径向漏扫补偿量,根据螺旋线步长、螺距和子区域面积等合理取值。
所提复合扫描方法在FOU内的全场最大扫描时间Tstep2_max计算如下:
$$ {T_{step2\_\max }} = n \cdot {({F_{motor}})^{ - 1}} + 24 \cdot n \cdot {({F_{PZT}})^{ - 1}}{\text{ ,}} \frac{{\pi \cdot {r_{FOU}}^2}}{{{L^2}}} = n $$ (13) 式中:rFOU为FOU半角。
将子区域极坐标系下的PZT扫描路径点坐标转换到终端二维坐标系,可转化为PZT的角度实际控制量,如图7所示,转换关系式如下:
$$ \left\{\begin{array}{c} \rho_{N_{\text {marix }}} \cos \left(\theta_{N_{\text {mari }}}\right)=\rho_{i} \cdot \cos \left(\theta_{i}\right)+\rho_{N_{\text {maric }}}^{\prime} \cdot \cos \left(\theta_{i}+\theta_{N_{\text {marix }}}^{\prime}\right) \\ \rho_{N_{\text {marix }}} \sin \left(\theta_{N_{\text {mari }}}\right)=\rho_{i} \cdot \sin \left(\theta_{i}\right)+\rho_{N_{\text {marix }}}^{\prime} \cdot \sin \left(\theta_{i}+\theta_{N_{\text {marix }}}^{\prime}\right) \end{array}\right. $$ (14) $$ \left\{ {\begin{array}{*{20}{c}} {{x_{{N_{matrix}}}} = {\rho _{{N_{matrix}}}}\cos \left( {{\theta _{{N_{matrix}}}}} \right)} \\ {{y_{{N_{matrix}}}} = {\rho _{{N_{matrix}}}}\sin \left( {{\theta _{{N_{matrix}}}}} \right)} \end{array}} \right. $$ (15) 式中:
$\rho_{N_{matirx} }$ 为子区域内PZT所有扫描点在终端极坐标系下的极径;$\theta _{N_{matirx} }$ 为子区域内PZT所有扫描点在终端极坐标系下的极角;$x _{N_{matirx} } $ 为子区域内PZT所有扫描点在终端二维坐标系下的方位角度控制量;$y _{N_{matirx} } $ 为子区域内PZT所有扫描点在终端二维坐标系下的俯仰角度控制量。 -
由于受到卫星姿态抖动等因素的影响,终端初始指向方向在FOU内呈正态随机分布[13],相同地,对方卫星终端位置在FOU内也呈正态随机分布。以扫描端终端为例,当扫描端信号光扫到对方凝视端视域内,凝视端则根据信号光在其视域内的位置进行角度调整,并沿着信号光反方向发射回光,由于这一过程时间很短(Tstep3<0.4 s),因此,可假设在FOU内,扫描端扫到凝视端视域的位置与凝视端回光位置近似相同。
文中仿真参数如表1所示,典型仿真环境为终端通信距离3000 km,星上微振动平均振幅15 μrad,凝视端视域位置与回光光斑位置相同,并且在扫描端FOU内呈正态随机分布。基于典型环境,开展所提复合扫描方法1000次,10000次,100000次重复性捕获回光光斑的仿真测试,光斑在FOU的正态随机分布情况,如图8所示。
文中定义:捕获成功率=成功捕获次数/总扫描场次。1000次,10000次,100000次重复性扫描捕获仿真测试结果,如表2所示。图9(a1)~(a2)、(b1)~(b2)为所提方法的部分仿真示意结果。由表2结果可知,所提方法平均捕获时长不高于12.2 s,全场最大捕获时长不超过36.5 s,最小捕获时长仅为0.0083 s,捕获成功率大于98.7%。根据第2节分析可知,步骤(3)和步骤(4)共耗时约0.4 s,与本节仿真得到的平均时长结果,即步骤(2)的平均时长累加,所提方法的整体平均耗时低于12.6 s,极大地降低了无信标扫描捕获时长,可有效提升无信标星间激光通信终端的在轨建链效率。
表 2 文中所提方法仿真结果
Table 2. Simulation results of the proposed method
No. Simulation item Average acquisition time/s Minimum acquisition time/s Maximum acquisition time/s Acquisition probability 1 1000 simulations 12.1421 0.0165 36.3100 99.1% 2 10000 simulations 12.1239 0.0083 36.4800 98.7% 3 100000 simulations 12.1241 0.0083 36.4800 98.7% 图 9 (a)~(b)所提无信标方法的部分仿真结果;(c)有信标方法的部分仿真结果
Figure 9. (a)-(b) Partial simulation results of the proposed beaconless method; (c) Partial simulation result of the beacon method
为了进一步证明文中所提方法的工程实用性,所提方法与经过在轨验证的有信标螺旋扫描方法进行仿真对比。信标光束散角为230 μrad,其余仿真参数与所提方法一致。仿真表明,有信标方法的最大捕获时间为27.8333 s,最小时间为0 s,捕获概率大于99.0%,图9(c1)~(c2)为有信标螺旋扫描方法的部分仿真示意结果。
与有信标扫描方法结果相比,在最大捕获时长这项指标上,所提方法比有信标方法平均多耗时近9 s,这是由于信号光比信标光小,导致全场走点数量比信标光多十几倍,造成了这项指标不能媲美有信标扫描方法。而在平均捕获时长、最小捕获时长和捕获成功率指标上,所提方法与有信标方法差距不大,进一步证明了所提方法的实用性和有效性。
Coarse-refined combined scanning method of beaconless inter-satellite laser communication system
-
摘要: 无信标星间激光通信系统无需额外设置信标光组件,直接采用小束散角的信号光进行瞄准捕获跟踪,有利于降低终端的质量、体积、功耗和制造成本等重要指标,满足商业航天和低轨卫星空间激光网络的发展要求。为了克服现有无信标激光通信系统扫描时间长、捕获难度大的问题,综合考虑星上姿态变化、测量精度和终端安装定位误差等因素,对扫描不确定域进行了数值计算,设计了无信标扫描捕获流程,定量分析了粗精扫描机构控制带宽、扫描步长、规划路径与光束覆盖面积的耦合关系,提出了一种简单可靠的无信标粗精复合螺旋扫描方法。典型环境仿真结果表明,采用所提方法平均扫描捕获时长小于20 s,捕获概率大于95%,有效提升了无信标系统的扫描效率,可满足未来低轨卫星空间激光网络的快速建链要求。Abstract: The beaconless inter-satellite laser communication system (BILCS) does not to configure beacon light components, and directly applies signal light with small divergence angle to pointing, acquisition and tracking (PAT), which is beneficial to reduce the weight, volume, power consumption and manufacturing cost of the terminal, and meet the development requirements of commercial spaceflight and LEO satellite space laser network. To overcome the existing shortcomings of long scanning time and difficult acquisition of BILCS, satellite attitude change, measurement accuracy and terminal positioning error were taken into consideration, the scanning field of uncertainty was numerically computed, the beaconless scanning-acquisition process was designed, the coupling relationship of the control bandwidth of the coarse and refined scanning mechanism, scanning step size, planning path and the light coverage area was quantitatively analyzed, a simple and reliable coarse-refined combined spiral scanning method based on beaconless was proposed. Simulation results of typical environment demonstrate that the average acquisition time is less than 20 s, and the acquisition probability is more than 95% by this method, which effectively improves the BILCS scanning efficiency, and can satisfy the requirements of the future LEO satellite space laser network to rapidly establish laser communication link.
-
Key words:
- inter-satellite laser communication /
- LEO satellite /
- beaconless /
- scanning method
-
表 1 总体仿真参数
Table 1. Overall simulation parameters
Parameter Value 1 Communication distance/km ≥3000 2 Satellite attitude stability/(°)·s-1 ≤0.02(3σ) 3 Attitude measurement accuracy/(°) 0.01 4 Satellite positioning error/m ≤50(3σ) 5 Signal light wavelength λ/nm 1550 6 Field of vision/mrad ±2.5 7 Signal light divergence angle α/μrad 50-70 8 CSM control frequency fmotor/ Hz ≤10 9 CSM control error δmotor/ μrad ≤30(3σ) 10 RSM control frequency fPZT/ Hz ≤800 11 RSM control error δPZT/ μrad ≤8(3σ) 12 Camera refresh rate fCCD/ Hz ≥120 表 2 文中所提方法仿真结果
Table 2. Simulation results of the proposed method
No. Simulation item Average acquisition time/s Minimum acquisition time/s Maximum acquisition time/s Acquisition probability 1 1000 simulations 12.1421 0.0165 36.3100 99.1% 2 10000 simulations 12.1239 0.0083 36.4800 98.7% 3 100000 simulations 12.1241 0.0083 36.4800 98.7% -
[1] 高世杰, 吴佳彬, 刘永凯等. 微小卫星激光通信系统发展现状与趋势[J]. 中国光学, 2020, v. 13;No. 73(06): 7-17. Gao Shijie, Wu Jiabin, Liu Yongkai, et al. Development status and trend of micro-satellite laser communication systems [J]. Chinese Optics, 2020, 13(73): 7-17. (in Chinese) [2] 顾健, 艾勇, 陈晶等. 扰动观测器在空间光通信PAT系统中的应用[J]. 红外与激光工程, 2016, 45(1): 0122003. doi: 10.3788/irla201645.0122003 Gu Jian, Ai Yong, Chen Jing, et al. Application of distribution observer for space optical communication PAT system [J]. Infrared and Laser Engineering, 2016, 45(1): 0122003. (in Chinese) doi: 10.3788/irla201645.0122003 [3] Zhang M, Tong S F, Song Y S, et al. Modeling and reality analysis of beaconless spatial acquisition for intersatellite laser communication system [J]. Optical Engineering, 2017, 56(10): 106108. [4] Yu Huilong, Bao Zhikang, Wang Xuan, et al. XY-2 setallite laser communication equipment PAT test in orbit [J]. Infrared and Laser Engineering, 2021, 50(5): 20200327. (in Chinese) [5] 王俊尧, 宋延嵩, 佟首峰等. 空间激光通信组网反射镜联动跟踪控制技术[J]. 中国光学, 2020, 13(3): 537-546. Wang Junxiao, Song Yansong, Tong Shoufeng, et al. Linkage tracking control technology of space laser communication network mirror [J]. Chinese Optics, 2020, 13(3): 537-546. (in Chinese) [6] 郭超勇, 刘继奎, 程俊波, 等. 卫星激光通信粗跟踪系统复合控制策略[J]. 光学精密工程, 2020, 28(04): 195-202. Guo Chaoyong, Liu Jikui, Chen Junbo, et al. Compound control strategy of satellite laser communication coarse tracking system [J]. Optics and Precision Engineering, 2020, 28(4): 195-202. (in Chinese) [7] 温涛, 魏急波, 马东堂. 无线光通信中的步进扫描捕获方法研究[J]. 红外与激光工程, 2006, A-0502-05: 1007-2276. Wen Tao, Wei Jibo, Ma Dongtang. A step method for laser acquisition in atmosphere laser communication [J]. Infrared and Laser Engineering, 2006, 35(S1): 1007-2276. (in Chinese) [8] Yu S Y, Wu F, Tan L Y, et al. Static position errors correction on the satellite optical communication terminal [J]. Optical Engineering, 2017, 56(2): 026112. doi: 10.1117/1.OE.56.2.026112 [9] Teng Y J, Zhang M, Tong S F. A novel approach for optimization control of dynamic tracking with platform vibration in the ground test of satellite laser communication systems [J]. IEEE Access, 2019, 7: 35049-35059. doi: 10.1109/ACCESS.2019.2903527 [10] 武凤. 窄信标卫星光通信链路快速捕获和稳定跟踪方法研究[D]. 哈尔滨工业大学, 2018. Wu Feng. Research on fast acquisition and stable tracking method with narrow beacon in satellite optical communication link[D]. Harbin: Harbin Institute of Technology, 2018. (in Chinese) [11] 国爱燕, 高文军, 周傲松等. 星间光通信无信标捕跟瞄技术[J]. 红外与激光工程, 2017, 46(10): 293-299. Guo Aiyan, Gao wenjun, Zhou Aosong, et al. Beaconless acquisition tracking and pointing of inter-satellite optical communication [J]. Infrared and Laser Engineering, 2017, 46(10): 1022002. (in Chinese) [12] Teng Y J, Zhang M, Tong S F. The optimization design of sub-regions scanning and vibration analysis for beaconless spatial acquisition in the inter-satellite laser communication system [J]. IEEE Photonics Journal, 2018, 10(6): 1-11. [13] 于思源. 卫星光通信瞄准捕获跟踪技术[M]. 科学出版社. 2016. Yu Siyuan. Satellite Optical Communication Pointing, Acquisition and Tracking Technology[M]. Beijing: Science Press, 2016. (in Chinese) [14] Xue Zhengyan, Qi Bo, Ren Ge. Vibration-induced jitter control in satellite optical communication[C]//Proc of SPIE, 2013, 8906: 89061. [15] Li Xin, Song Quanyou, Ma Jing, et al. Spatial acquisition optimization based on average acquisition time for intersatellite optical communications[C]//IEEE 10 th Russian-Chinese Symposium on Laser Physics and Laser Technologies (RCSLPLT) and 2010 Academic Symposium on Optoelectronics Technology (ASOT), 2010: 244-248. [16] Tu Bo, Liu Lu, Liu Yihui. Acquisition probability analysis of ultra-wide FOV acquisition scheme in optical links under impact of atmospheric turbulence [J]. Applied Optics, 2013, 52(14): 3147-3155. doi: 10.1364/AO.52.003147