光学学报, 2019, 39 (2): 0206003, 网络出版: 2019-05-10   

基于贝克曼分布的星间激光通信链路性能优化 下载: 1066次

Link Performance Optimization for Inter-Satellite Laser Communications Based on Beckmann Distribution
作者单位
国防科技大学电子对抗学院脉冲功率激光技术国家重点实验室, 安徽 合肥 230037
摘要
针对星间激光通信系统易受瞄准误差影响的问题,在瞄准误差服从贝克曼分布的条件下,利用三阶中心矩法将其等效为修正后的Rayleigh分布,分别在中断概率一定和发射功率一定的假设下,建立了发射功率优化模型和中断概率优化模型,并计算出了两种假设条件下高斯光束经天线发射时的最优均方根宽度。通过数值仿真,给出了考虑瞄准误差时的系统中断概率、最小发射功率与最优均方根宽度之间的关系。根据该数值分析结果,在明确瞄准误差分布规律的基础上,选择最优的高斯光束均方根宽度,以实现最优星间激光通信链路性能。
Abstract
Aim

ing at the problem that the inter-satellite laser communication system is easily affected by the pointing errors, the method of third-order central moment is used to make the pointing error equivalent to the modified Rayleigh distribution under the condition that the pointing errors obey the Beckmann distribution. Under the assumptions that the outage probability and the transmission power are fixed, the optimal model of transmitter power and the optimal model of outage probability are established, respectively. The optimal root mean square width of the Gaussian beam emitted by the antenna under the above two assumptions is calculated. Through the numerical simulation, the numerical relationship of the outage probability and the minimum transmitting power with the optimal root mean square width is given. According to the numerical analysis results, the optimal root mean square width of a Gaussian beam can be selected based on the definite pointing error distributions, so that the optimal link performance for the inter-satellite laser communications is obtained.

1 引言

在星间激光通信过程中,必须确保发射端与接收端的光路对准,才能进行有效的数据传输。然而,受星上平台振动、星间相对运动和系统内部噪声等因素的影响,对准的光路发生了偏移,使得发射机和接收机之间存在一定的瞄准误差[1]。在实际应用中,由于光束发散角极窄且通信距离通常非常远,因此造成系统性能恶化的主要因素是瞄准误差[2]

为了抑制瞄准误差对星间激光通信系统性能的影响,当前热门的解决方案是采用捕获、瞄准、跟踪(APT)技术。比如,郑燕红等[3]提出采用控制混合灵敏度设计方法对控制器进行设计,以抑制干扰并解决受控对象不确定性问题;倪小龙等[4]设计并构建了基于反射式液晶光调制器的星载平台振动补偿系统,满足了星际激光通信系统的高动态跟踪精度的要求。但以上方法所采用的技术复杂、成本过高,实现起来比较困难。为此,学者们希望通过优化设计系统参数以降低瞄准误差的影响。丁涛等[5]提出通过控制接收端光斑的形状和大小来抑制振动的影响,但是该方法仅适用于短距离的通信环境;Toyoshima等[6]在瞄准误差服从Rayleigh分布时,通过数值计算给出了系统最优光束发散角,但并没有给出闭合的表达式;Song等[7]提出了采用动态束腰控制技术来抑制瞄准误差的影响,但需要在系统发射端增加高速束腰控制模块。在已有的研究中,均假设瞄准误差在方位向和俯仰向是理想同分布的,没有考虑静态偏置误差的影响,无法精确地描述实际卫星平台的随机振动特性,且没有给出解析形式的参数优化方法,系统实现比较困难。

本文在瞄准误差服从贝克曼分布的情况下,结合星间激光通信链路传输方程,给出了系统中断概率的闭合表达式,并分别在中断概率一定和发射功率一定的条件下,建立了发射功率优化模型和中断概率优化模型,通过模型求解,计算出了这两种条件下的高斯光束最优均方根宽度。最后通过数值仿真,给出了不同瞄准误差角分布情况下的最优高斯光束均方根宽度。

2 星间激光通信基本原理

2.1 星间激光通信信号模型

假定星间激光通信系统采用强度调制直接探测,调制样式为开关键控(OOK)调制,发射光束为高斯光束,此时的星间链路传输方程可以表示为[8]

PR=PTGTGRλ4πd2ηTηRLT(GT,θ),(1)

式中:PR为接收到的光功率;PT为发射功率;GT为发射天线增益,当发射光束为高斯光束时,GT=(2πW/λ)2,W为信号高斯光束经天线往外发射时的均方根宽度,λ为激光波长;GR=(πDR)2为接收天线增益,DR为接收天线的孔径直径,d为链路传输距离;ηTηR分别为发射天线效率和接收天线的效率;LT(GT,θ)=exp(-GTθ2)为瞄准误差损耗因子,θ为瞄准误差偏离角。

在系统接收端,探测器会将接收到的光功率转化为信号电流I,其表达式为

I=RPR,(2)

式中:R为光电响应度, R=ηqλ/(hc),其中,h为普朗克常量,q为电子电荷,η为量子效率,c为光波频率。

综上,光电探测器输出的信号电流为

I=RPTGTGRλ4πd2ηTηRexp(-GTθ2)(3)

2.2 基于贝克曼分布的瞄准误差模型

卫星平台振动等因素将会导致星间激光通信发射机与接收机之间存在瞄准误差。假定指向光束的径向瞄准误差角向量θ=(θx,θy)T,其中θxθy分别表示方位向和俯仰向的瞄准误差角分量,分别用均值非零的高斯分布变量来描述,θx~N(μx, σx2),θy~N(μy, σy2),此时瞄准误差下探测区域与接收光束截面之间的关系如图1所示。图中μ= μx2+μy2表示光通信终端的光学系统轴线和天线法向之间存在的一个夹角,称为静态偏差角。则此时径向的瞄准误差角θ=|θ|= θx2+θy2服从贝克曼分布[9],即

fθ(θ)=θ2πσxσy·02πexp-(θcosϕ-μx)22σx2-(θsinϕ-μy)22σy2dϕ,(4)

式中:μxμy分别表示方位向和俯仰向的静态偏移误差;σxσy分别表示方位向和俯仰向的瞄准误差抖动标准差。

图 1. 瞄准误差下探测区域与接收光束截面之间关系示意图

Fig. 1. Schematic of relationship between detection area and beam cross section of receiving beam under pointing errors

下载图片 查看所有图片

假设将贝克曼分布近似等效成参数为σmod的修正Rayleigh分布,σmod的表达式可以通过三阶中心矩法给出。通过数值比较可以看出,利用三阶中心矩法可以比较准确地计算出σmod的值,使得修正后的Rayleigh分布与贝克曼分布基本保持一致[10]

随机变量χm阶中心矩可以表示为

Ωmχ=E{[χ-E(χ)]m},(5)

式中E(•)表示期望。

为了方便起见,令u=θ2,此时,当瞄准误差角 θ服从修正后的Rayleigh分布时, u的概率密度函数为

f(u)=12σmod2exp-u2σmod2(6)

u的三阶中心矩 Ω3u可以表示为

Ω3u=16σmod6(7)

此时,θ2的三阶中心矩 Ω3θ2同样可以用两个正态随机变量平方和的三阶中心矩来表示,即

Ω3θ2=8σx4(3μx2+σx2)+8σy4(3μy2+σy2)(8)

利用(7)、(8)式,则可以给出 σmod2的表达式,即

σmod2=3μx2σx4+3μy2σy4+σx6+σy621/3(9)

最终,可以将贝克曼分布等效为修正后的Rayleigh分布,θ的概率密度函数(PDF)可近似表示为

fθ(θ)θσmod2exp-θ22σmod2,θ0,(10)

其中,

σmod=3μx2σx4+3μy2σy4+σx6+σy621/6(11)

通过计算,θ的累计分布函数(CDF)为

Fθ(θ)1-exp-θ22σmod2(12)

3 星间激光通信系统优化模型

3.1 中断概率分析

根据信息理论的相关知识,存在一定的概率使信道容量C不足以满足传输速率R0的需求,通常称此时的情况为中断,可用中断概率[11]Pout来衡量,即

Pout=P[C(S)R0],(13)

式中S为信噪比(SNR)。

由于C(·)是随SNR单调递增的,故中断概率可表示为

Pout=P(Sa),(14)

式中:a=C-1(R0)为预先设定好的阈值,C-1(·)为C(·)的反函数。

在星间激光通信中,系统接收端经光电探测器后得到的信号w可以表示为w=RPRX+n,其中X表示发送的信号,n表示加性高斯白噪声,噪声均值为0,方差为 σN2,结合(3)式可得系统接收到的电信号的瞬时信噪比S

S=I22σN2=Bexp(-2GTθ2),0<SB,(15)

其中

B=12σN2λ4πd4(RPTGTGRηTηR)2(16)

为了求出系统的中断概率,需要先求出瞬时信噪比S的CDF,其表达式为

FSs=P(Ss)=PBexp-2GTθ2s=Pθ12GTlnBs=1-Fθ12GTlnBs=sB14GTσmod2(17)

在给定阈值a时,可得系统的中断概率为

Pout=P(Sa)=aB14GTσmod2(18)

3.2 星间激光通信链路优化模型

从(16)式中可以看出,在星间激光通信链路性能计算中涉及到的系统参数比较多,其中,发射功率和发射天线增益是其中相对比较重要的系统参数,为了综合考虑所有涉及到的系统参数,在此定义了两个新的变量[12]:归一化发射功率(y=αPT/σmod2)和归一化发射天线增益(x=GTσmod2),其中α可以表示为

α=122σNλ4πd2RGRηTηR,(19)

B=4x2y2代入 (18)式中可得

Pout=a414x1xy12x(20)

1) 发射功率优化模型。

为了使系统中断概率一定时所需的发射功率最小,建立了发射功率优化模型,可以表示为

miny=f(x),(21)

可以通过选取适当的x值使得y最小。由于y=αPT/σmod2中,ασmod都可以看作是定值,此时所求的最小发射光功率即可转化为使中断概率满足Pout=b时最小的y值。

此时,将Pout=b代入(20)式中可得

a414x1xy12x=b,(22)

利用yx进行求导可得

dydx=a42xln(1/b)-1x2b2x,(23)

dydx=0,可得此时的极值点为

x=12ln(1/b)(24)

为了确定这个极值点是否为极小值点,可以通过求yx的二阶导数来确认,yx的二阶导数为

d2ydx2=a44x2ln(1/b)xln1/b-1+2xx4b2x(25)

将(24)式代入到(25)式中可得

d2ydx2=a4[2ln(1/b)]3b-1ln(1/b),0<b<1(26)

因此,(24)式中的极值点即为所求的极小值点。当中断概率满足Pout=b时,最小的y值为

ymin=aln(1/b)b(1/lnb)(27)

y=αPT/σmod2可得,此时最小的发射功率为

PTmin=22aσNln(1/b)b(1/lnb)σmod2RGRηTηR4πdλ2(28)

由于x=GTσmod2=(2πW/λ)2σmod2,此时高斯光束发射时的最优均方根宽度W

W=λ2πσmod12ln(1/b)(29)

2) 中断概率优化模型

为使在给定发射功率条件下系统的中断概率最小,建立了系统中断概率优化模型,可以表示为

minPout=g(x),(30)

可以通过选取适当的x值使得Pout最小。

在发射功率一定的情况下,令y=k,将其代入(20)式中可得

Pout=a414x1xk12x,(31)

Poutx进行求导可得

dPoutdx=-a4x2k214x2+lna4x2k24x2,(32)

dPoutdx=0,可得此时的极值点为

x=ae2k,(33)

式中e为自然常数,约为2.71828。

为了确定所求的极值点是否为所需的极小值点,求Poutx二阶导数为

d2Poutdx2=a4x2k214x2+lna4x2k22+8x2+lna4x2k2+8x16x4(34)

从(34)式中可以看出,当x= ae/2k时,{2+ln[a/(4x2k2)]}=0,此时有d2Pout/dx2>0,x= ae/(2k)即为所求的极小值点。

此时的最小中断概率为

PO,min=(1/e)kae(35)

同样,由x=GTσmod2=(2πW/λ)2σmod2,可得此时高斯光束经天线往外发射时的最优均方根宽度W

W=λ2πae2kσmod2(36)

4 数值仿真

为了分析星间激光通信系统中各个参数对链路性能的影响,利用(28)、(29)、(35)、(36)式进行了数值仿真。仿真参数[13]设置如下:传输距离d=20000 km,接收天线的直径DR=0.3 m,激光波长λ=1064 nm,发射天线效率ηT=0.5,接收天线效率ηR=0.35,量子效率η=0.7,电子电荷q=1.602176×10-19 C,普朗克常量h=6.626069×10-34 J·s,光波频率c=3×108 m/s,高斯白噪声的方差 σN2=4.27×10-12 A2,设定阈值a=1×10-7。瞄准误差角θ服从贝克曼分布,分别在μx=μy=0和μx=1 μrad, μy=0.5 μrad这两种情况下,当σx=1 μrad,σyx分别取1, 2, 3/2 μrad时进行了仿真分析实验。

图2所示为所需的最小发射功率随中断概率的变化情况,从图中可以看出,当系统需要满足的中断概率越大时,所需的发射功率越小,且随着中断概率的减小,所需的发射功率将迅速增大。从瞄准误差角的分布情况不同时的对比可以看出,静态偏差角为0 (μx=μy=0)时所需的最小发射功率要明显小于静态偏差角不为0(μx=1 μrad, μy=0.5 μrad)时的,并且当σxy的取值不同时,一定的中断概率下,系统所需的发射功率相差较大,如当系统中断概率为0.01,μx=μy=0,σyx=1时,所需的发射功率要比σyx=3/2和σyx=2时分别小12.13 mW和34.38 mW。

图3所示为系统所需发射功率最小时的高斯光束最优均方根宽度随中断概率的变化规律。从图中可以看出,随着中断概率的增大,高斯光束的最优均方根宽度也随着增大,且当瞄准误差角的分布情况不同时,最优均方根宽度的取值差异较大;从图中还可以看出,当静态偏差角为0时的最优均方根宽度要明显大于静态偏差角不为0时的,并且当系统中断概率为0.01,μx=μy=0,σyx分别为1, 3/2,2时的最优均方根宽度依次为55.8,41.1,31.2 mm。结果表明,瞄准误差角的分布情况对星间激光通信系统参数的选取影响很大,在中断概率一定时,为了使系统所需的发射功率最小,必须要在明确瞄准误差角的分布规律下,选择合适的高斯光束均方根宽度。

图 2. 最小发射功率与中断概率间的关系

Fig. 2. Minimum transmitter power versus outage probability

下载图片 查看所有图片

图 3. 最优均方根宽度与中断概率间的关系

Fig. 3. Optimal RMS width versus outage probability

下载图片 查看所有图片

图4、5所示为当发射功率固定时,系统可以达到的最小中断概率以及此时高斯光束的最优均方根宽度随发射功率的变化情况。从图4、5可以看出,随着发射功率的增大,系统的中断概率性能逐渐提升,此时最优均方根宽度也随之减小。从图5可以看出,当瞄准误差角的分布情况不同时,最优均方根宽度的取值分布产生差别,此时静态偏差角为0时的最优均方根宽度要明显大于静态偏差角不为0时的,并且当系统发射功率为10 mW,μx=μy=0,σyx分别为1, 3/2,2时的最优均方根宽度依次为35.1,25.9,19.7 mm。结果表明发射功率固定时,为达到最小的中断概率,必须要结合瞄准误差角分布情况,选择最优的高斯光束均方根宽度。

图 4. 最小中断概率与发射功率间的关系

Fig. 4. Minimum outage probability versus transmitter power

下载图片 查看所有图片

图 5. 最优均方根宽度与发射功率间的关系

Fig. 5. Optimal RMS width versus transmitter power

下载图片 查看所有图片

5 结论

针对瞄准误差影响下的星间激光通信系统,结合星间链路传输方程,给出了当瞄准误差服从贝克曼分布时的中断概率表达式,并分别在固定中断概率和发射功率的条件下建立了发射功率优化模型和中断概率优化模型,通过模型求解,分别给出了这两种假设条件下高斯光束发射时的最优均方根宽度。通过数值仿真,给出了当中断概率一定时,系统所需的最小发射功率和最优均方根宽度随中断概率的变化规律,以及当发射功率一定时,系统所能达到的最小中断概率和最优均方根宽度随发射功率的变化规律。从数值仿真结果可以看出,当瞄准误差角的分布情况不同时,最优均方根宽度的取值相差较大,比如当系统中断概率为0.01,μx=μy=0,σyx分别为1,3/2,2时的最优均方根宽度依次为55.8,41.1,31.2 mm;当系统发射功率为10 mW,μx=μy=0时,σyx分别为1,3/2,2时的最优均方根宽度分别为35.1,25.9,19.7 mm。因此,为了使系统性能最优,必须要结合瞄准误差角分布,选择最优的高斯光束均方根宽度。本研究可为星间激光通信系统的研究和设计提供理论依据,具有一定的参考价值。

参考文献

[1] 周凌林, 郧建平, 杨海峰, 等. 指数韦伯信道下带瞄准误差的相干光通信系统中断概率分析[J]. 激光与光电子学进展, 2017, 54(11): 110602.

    Zhou L L, Yun J P, Yang H F, et al. Outage probability analysis of coherent optical communication system with pointing error in exponentiated Weibull channel[J]. Laser & Optoelectronics Progress, 2017, 54(11): 110602.

[2] 张慧颖, 李洪祚, 肖冬亚, 等. 大气湍流综合效应下空间分集接收性能研究[J]. 中国激光, 2016, 43(4): 0405002.

    Zhang H Y, Li H Z, Xiao D Y, et al. Performance analysis of spatial-diversity reception over combined effects of atmospheric turbulence[J]. Chinese Journal of Lasers, 2016, 43(4): 0405002.

[3] 郑燕红, 王岩, 陈兴林. 卫星光通信APT控制系统H∞设计[J]. 航空学报, 2008, 30(6): 1619-1625.

    Zheng Y H, Wang Y, Chen X L. H∞ control applied for APT system of inter-satellite laser communications[J]. Acta Aeronautica ET Astronautica Sinica, 2008, 30(6): 1619-1625.

[4] 倪小龙, 刘智, 姜会林, 等. 基于液晶光调制器的星际激光通信平台振动补偿系统[J]. 中国激光, 2013, 40(s1): s105007.

    Ni X L, Liu Z, Jiang H L, et al. Space laser communication vibration compensation system based on liquid crystal light modulator[J]. Chinese Journal of Lasers, 2013, 40(s1): s105007.

[5] 丁涛, 许国良, 张旭苹, 等. 空间光通信中平台振动对误码率影响的抑制[J]. 中国激光, 2007, 34(4): 499-502.

    Ding T, Xu G L, Zhang X P, et al. Control of bit error rate introduced by platform vibration for free space optical communication[J]. Chinese Journal of Lasers, 2007, 34(4): 499-502.

[6] Toyoshima M, Jono T, Nakagawa K, et al. Optimum divergence angle of a Gaussian beam wave in the presence of random jitter in free-space laser communication systems[J]. Journal of the Optical Society of America A, 2002, 19(3): 567-571.

[7] Song T Y, Wang Q, Wu M W, et al. Performance of laser inter-satellite links with dynamic beam waist adjustment[J]. Optics Express, 2016, 24(11): 11950-11960.

[8] 刘宏展, 纪越峰, 刘立人. 像差对星间相干光通信接收系统误码性能的影响[J]. 光学学报, 2012, 32(1): 0106002.

    Liu H Z, Ji Y F, Liu L R. Effect of aberration on performance of the bit error rate in an inter- satellite coherent optical communication receiving system[J]. Acta Optica Sinica, 2012, 32(1): 0106002.

[9] Yang F, Cheng J L, Tsiftsis T A. Free-space optical communication with nonzero boresight pointing errors[J]. IEEE Transactions on Communications, 2014, 62(2): 713-725.

[10] Boluda-Ruiz R, García-Zambrana A, Castillo-Vázquez C, et al. Novel approximation of misalignment fading modeled by Beckmann distribution on free-space optical links[J]. Optics Express, 2016, 24(20): 22635-22649.

[11] 韩立强, 江红兵. 一种混合认知RF和MIMO FSO系统的中断概率分析[J]. 中国激光, 2018, 45(4): 0406001.

    Han L Q, Jiang H B. Outage probability analysis of a mixed cognitive RF and MIMO FSO system[J]. Chinese Journal of Lasers, 2018, 45(4): 0406001.

[12] Chen C C, Gardner C S. Impact of random pointing and tracking errors on the design of coherent and incoherent optical intersatellite communication links[J]. IEEE Transactions on Communications, 1989, 37(3): 252-260.

[13] Song T Y, Wang Q, Wu M W, et al. Impact of pointing errors on the error performance of intersatellite laser communications[J]. Journal of Lightwave Technology, 2017, 35(14): 3082-3091.

宛雄丰, 郝士琦, 张岱, 赵青松, 徐晨露, 唐进迎. 基于贝克曼分布的星间激光通信链路性能优化[J]. 光学学报, 2019, 39(2): 0206003. Xiongfeng Wan, Shiqi Hao, Dai Zhang, Qingsong Zhao, Chenlu Xu, Jinying Tang. Link Performance Optimization for Inter-Satellite Laser Communications Based on Beckmann Distribution[J]. Acta Optica Sinica, 2019, 39(2): 0206003.

本文已被 3 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!