Performance Analysis for Mixed RF/FSO Airborne Communication Systems over Atmospheric Turbulence and Pointing Error
1 引言
自由空间光(FSO)通信具有高速率传输、安全抗干扰及架设方便快捷等优势,与传统射频无线通信相比,具有广泛的民用和**应用前景,近年来受到越来越多的关注。FSO系统具有高性价比和大吞吐量特点,可将其应用于航空通信系统中与射频(RF)通信互补,以提高系统容量及通信性能[1-2]。然而FSO链路易受大气湍流及航空平台运动的影响,尤其是在长距离传输条件下,大气湍流引起的光强闪烁及航空平台抖动导致的指向误差(PE)都会影响航空通信系统性能[3-4],因此采用混合RF/FSO通信系统可实现两者优势互补。在混合RF/FSO航空通信系统中,RF链路用于子节点接入骨干链路,而宽带FSO链路用于航空骨干链路的建立。为了实现RF子节点对FSO骨干链路的无缝接入,同时增强混合RF/FSO通信系统的可靠性,可采用双跳中继技术来补偿由大气湍流及指向误差导致的性能衰落。近年来,许多专家学者基于双跳中继的混合RF/FSO系统开展研究[5-8]。
当前大气激光通信中的大气湍流模型有对数正态(LN)分布、Gamma-Gamma(GG)分布及 M分布等。Samimi等[9]基于放大-转发中继方式,研究了混合RF/FSO系统的中断性能,其中FSO信道服从M分布,RF通信链路中信号衰落服从瑞利分布。Ansari等[10]分析了放大转发中继混合RF/FSO系统的误码率性能,RF链路及FSO链路分别基于Nakagami-m和GG分布。Zedini等[11]针对大气湍流及指向误差影响下的双跳中继混合系统性能进行研究,分析了采用放大转发固定增益中继和可变增益中继方式的中继系统性能。LN分布模型适用于弱湍流条件下;M分布及GG分布模型可有效描述基于点接收机的弱-强湍流条件下的光强起伏,但不适用于孔径平均接收条件。为了准确描述从弱到强各种湍流强度下基于有效接收孔径的光强起伏,Barrios等[12]提出了适用于孔径平均条件下的指数型Weibull (EW)衰落模型。而基于指数型Weibull湍流分布模型及指向误差影响的混合RF/FSO航空通信系统中断及误码性能鲜有报道。
在无线光通信系统中,与一般开关键控(OOK)调制方案相比,副载波强度调制(SIM)方式在大气湍流信道中具有优势[13-15]。本文采用副载波多进制相移键控调制构建混合RF/FSO航空通信系统,RF及FSO信道分别采用Nakagami-m及指数型Weibull衰落分布模型,综合考虑大气湍流、指向误差及衰落指数对混合传输系统性能的影响,通过Meijer's G函数推导出混合传输系统的平均误码率和中断概率的闭合表达式,通过数值分析研究了不同湍流强度和接收孔径尺寸条件下航空通信参数与系统中断性能的关系,并通过蒙特卡罗仿真验证了推导公式的正确性。
2 系统模型
考虑双跳中继混合RF/FSO航空通信系统,系统采用放大-转发中继方式,通信系统如图1所示,其中源节点S到中继节点R间采用RF链路,中继节点R与目的节点间经过FSO链路。在源节点S发送端,信源数据被调制为电信号,经过RF链路传输在中继节点R对光源进行强度调制,将光信号经过大气湍流FSO信道发送至接收端。
图 1. 双跳中继混合RF/FSO航空通信系统
Fig. 1. Dual-hop relaying mixed RF/FSO aeronautical communication system
下载图片 查看所有图片
副载波相移键控调制(M-PSK)强度调制/直接检测的混合RF/FSO航空通信系统框图如图2所示。在源节点S发射端(Tx),对信源产生的信息序列进行信道编码和电调制;调制后的电信号经过RF信道传输,信道衰落服从Nakagami-m分布;在中继节点R信号经过带通滤波器(BPF),再进行直流(DC)偏压保证光信号输出,通过前置放大器(PA)放大后对光信号进行强度调制;调制后的光信号经过FSO大气湍流信道传输到达目的节点;在接收端(Rx)上信号先经过光学带通滤波器(OBPF)滤波,再通过光电探测器将光信号变成电信号,最后经过电解调和信道解码得到原信息。
图 2. 副载波强度调制系统示意图
Fig. 2. Diagram of SIM system
下载图片 查看所有图片
2.1 RF链路
在中继节点接收到的信号可表示为
式中x为源节点S发送的信号,平均功率为ES,R,hS,R为RF信道增益,nS,R为方差为
的加性高斯白噪声。RF链路服从Nakagami-m衰落分布,根据(1)式中定义,系统信噪比可表示为γS,R=ES,R|hS,R|2/
,RF链路信噪比γS,R的概率密度分布(PDF)公式为[16]
式中m为衰落指数(m≥1/2),Γ(m)为伽玛函数,γ为信噪比,平均信噪比
=ES,R/。将(2)式代入FγS,R(γ)=fγS,R(x)dx中,并利用低阶不完全Gamma函数定义γ(s,x)=ts-1exp(-t)dt进行计算(s,x为伽玛函数定义中的变量),可得γS,R的累积分布函数(CDF)为
2.2 FSO 链路
目的节点D的接收信号可表示为
式中
为中继节点处接收信号的估计值,ηR,D为光电转换效率,nR,D为均值为0、方差为
的加性高斯白噪声,IR,D为FSO信道衰落增益。FSO信道衰落包括路径损耗、指向误差及大气湍流等导致的衰落。
2.2.1 仅考虑大气湍流影响
考虑孔径平均效应,假设FSO链路服从指数型Weibull分布模型,则根据文献[
12]可得信噪比γR,D的PDF表达式为
相应的CDF表达式为
式中α,β和η分别为与闪烁因子及光强均值相关的参数,满足条件α>0,β>0及η>0,文献[
17]中给出了通过曲线拟合得到的参数经验公式,具体参数值可通过公式计算得到。
2.2.2 考虑大气湍流及指向误差影响
采用指数型Weibull分布模型描述大气湍流衰落,考虑FSO链路指向误差服从瑞利分布,根据牛顿一般二项式定理(1+z)δ=zi(z,δ,i为牛顿一般二项式定理公式中的变量),推导可得信噪比γR,D的PDF表达式为
式中ρ=ωeq/(2σs)表示接收端等效光束半径与瞄准误差偏移量标准差之比;ωeq为等效光束半径;
=ω2erf(v)/[2vexp(-v2)],其中erf()为误差函数,v=a/与接收端探测器半径相关;σs为接收端指向误差的标准差;A0=[erf(v)]2表示径向距离r=0处的部分收集信号功率[18];
为Merjer G函数。
将(7)式代入FγR,D(γ)=fγR,D(x)dx中可得相应的CDF表达式为
3 信噪比模型
分别考虑大气湍流影响下及大气湍流与指向误差衰落共同作用影响下的FSO信道,基于解码-转发中继方式,给出了混合RF/FSO通信系统中端到端信噪比(SNR)的累积分布函数。系统端到端信噪比可表示为
3.1 仅考虑大气湍流影响
假设FSO信道服从EW大气湍流分布模型,则变量γz的累积分布函数可写为[7]
将(3)式及(6)式代入(10)式,可得其累积分布函数CDF为
3.2 考虑大气湍流及指向误差影响
综合考虑大气湍流及指向误差概率密度分布函数下的FSO信道,将(3)式及(8)式代入(10)式,可得相应的CDF表达式为
4 系统性能分析
4.1 仅考虑大气湍流影响
4.1.1 中断概率
在混合RF/FSO通信系统中,中断概率描述系统端到端信噪比低于某一目标门限值γth的概率[3],根据(11)式可将解码-转发中继系统的中断概率表示为
4.1.2 平均误码率
系统采用M-PSK调制方式时,系统误码率(BER)可表示为[19]
式中D(m)=max(lb M,2),T=max(M/4,1),aτ=sin[(2τ-1)π/M],erfc(·)为误差函数,则系统平均误码率为
对(14)式求导可得系统误码率的一阶导数为
根据牛顿二项式定理,可将(11)式中第二项包含的多项式
展开得到
(16)式中的指数项exp(x)可写作Meijer's G函数[20]。根据函数定义,可将不完全Gamma函数展开为无穷级数形式[21]:
为了对(15)式中的积分形式进行求解,将(11)、(16)及(17)式代入(15)式中,利用上述无穷级数展开式和Meijer's G函数相关运算性质[文献[
22]中(5.6.3.1)及文献[
20]中(2.24.3.1)式],可得到系统平均误码率的具体表达式为
根据Meijer's G函数性质[20],(19)式中l和k为满足l/k=β/2的正整数,有k∈Ν+∧l∈Ν+∧gcd(k,l)=1(N+为正整数集合,∧为逻辑与,gcd为最大公约数);Δ(K,A)=,
,…,
。
4.2 考虑大气湍流及指向误差影响
4.2.1 中断概率
考虑大气湍流及指向误差影响,系统中断概率可由(12)式得到:
4.2.2 平均误码率
联合(12)、 (15)~(17)式,利用Meijer's G函数相关运算性质[文献[
22]中(5.6.3.1)式及文献[
20]中(2.24.3.1)式],可推导得到混合RF/FSO通信系统的平均误码率表达式为
5 仿真及结果分析
为了分析孔径平均效应影响下混合RF/FSO通信系统传输性能,根据以上推导的中断概率及误码率表达式(13)、(19)~(21)式进行了仿真分析。(19)~(21)式中,中断概率及误码率的表达式为无穷级数形式,在计算时选择i=30,则公式基本收敛,此时公式截断误差近似值分别为2×10-7,6×10-8,8×10-7。同时进行蒙特卡罗仿真来验证公式的正确性,蒙特卡罗仿真中,通过混合RF/FSO系统累积分布函数的反函数,对自变量在[0,1]内随机取值,计算得到关于自变量的符合CDF分布的样本,为了降低仿真中的统计不确定性,在每次仿真中产生108个随机数进行计算。图3和图4所示为混合微波/激光传输系统的中断概率随平均信噪比的变化规律,图5~8所示为系统平均误码率与平均信噪比之间的关系。通信系统中链路具体参数如表1 所示,假设RF链路和FSO链路的平均信噪比相同。
表 1. 系统参数
Table 1. System parameters
Parameter | Symbol | Value |
---|
Link distance /km | L | 100 | Wavelength /nm | λ | 1550 | Optical-to-electrical conversion coefficient | η | 0.9 | Noise standard deviation /(A·Hz-1) | σn | 10-7 | Corresponding jitter standard deviation /cm | σs | 10 | Corresponding beam radius /cm | ω | 10 |
|
查看所有表
图3所示为弱湍流大气信道影响下接收机孔径分别为15 cm和20 cm时系统中断概率与平均信噪比的关系曲线。设大气结构常数
的数值为2.7×10-18 m-2/3,则相应的弱湍流条件下Rytov指数为
=0.24;设信噪比门限值为10 dB。由图3可知,在两种接收孔径条件下系统中断概率均受指向误差的影响而增大。如在平均信噪比30 dB处,接收孔径为15 cm时,无指向误差影响条件下系统的中断概率为5.2×10-6,而大气湍流和指向误差联合作用下系统中断概率数值增加了4个数量级,约为5.1×10-2。同时可以看到系统中断概率随接收孔径增大而降低,在大气湍流及指向误差联合作用下,当接收孔径尺寸从D=15 cm增加至D=20 cm时,中断概率降低了1个数量级。
中等强度湍流条件下系统中断性能如图4所示,相应的大气结构常数为
=1.9×10-17 m-2/3(
=1.79)。从图4中数值可知数值计算结果与蒙特卡罗仿真结果基本符合,验证了(13)式和(20)式的理论正确性。同时可从图4中看到,与弱湍流条件相比,中等湍流强度条件下指向误差对系统中断性能影响较小。这是因为在强大气湍流影响下,系统性能衰落严重,此时指向误差导致的信道衰落对系统中断性能影响较小。由图3中数值结果分析可知,在弱湍流条件下,指向误差对系统中断性能的影响效果明显,当接收孔径为15 cm、平均信噪比为30 dB时,指向误差影响下中断概率从5.2×10-6增大到5.1×10-2,约提高了4个数量级,显著高于中等湍流强度条件下的数值结果。
图 3. 弱湍流条件下中断概率随平均信噪比变化规律
Fig. 3. Outage probability versus average SNR for weak turbulence conditions
下载图片 查看所有图片
图 4. 中等强度湍流条件下中断概率随平均信噪比变化规律
Fig. 4. Outage probability versus average SNR for moderate turbulence conditions
下载图片 查看所有图片
图5所示为混合中继通信系统中不同调制方式下平均误码率的变化规律。该仿真在中等湍流强度(
=1.79)条件下进行,分别选取接收孔径为点接收机和D=20 cm,且不考虑指向误差影响,图5所示仿真结果由(19)式计算得出。由图5可知,BPSK调制方式下系统误码性能优于高阶PSK调制方式下的误码率结果,同时可以看出孔径平均效应可以明显改善链路误码性能,例如,当系统平均信噪比为30 dB时,BPSK调制方式下接收孔径D=20 cm时的系统误码率为1.8×10-5,相应的点接收机条件下的误码率为8.5×10-5。
考虑指向误差影响,不同接收孔径尺寸对应的平均误码率结果如图6所示,该仿真采用BPSK调制方式,考虑弱湍流大气信道条件,相关信道参数与图3中所用参数相同。由图6可以看出,随着接收孔径尺寸的增大,混合RF/FSO航空通信系统误码性能逐渐被改善。当平均信噪比取14 dB时,考虑大气湍流及指向误差影响,接收孔径D=20 cm时系统误码率为9×10-3,而D=15 cm时该数值增大为4.2×10-2。这是因为在孔径平均效应的作用下,接收端孔径尺寸大于光强起伏的相干长度,则通信系统接收的光强为若干个散斑场的平均值,因此与点接收孔径相比,系统误码性能随接收孔径尺寸的增大而明显改善[23]。
图 5. 不同调制方式下平均误码率随平均信噪比变化规律
Fig. 5. Average BER versus average SNR for different modulation techniques
下载图片 查看所有图片
图 6. 弱湍流条件下平均误码率随平均信噪比变化规律
Fig. 6. Average BER versus average SNR for weak turbulence conditions
下载图片 查看所有图片
分别考虑有无指向误差影响,混合系统平均误码率在中等强度湍流条件下随平均信噪比的变化规律如图7 所示。该仿真分别选取D=15 cm和D=20 cm的接收孔径,系统相关参数与图4 中给出的数值相同。其中在有无指向误差影响条件下,分别得到系统误码性能,由图7可知,数值计算结果与相应的蒙特卡罗仿真结果在各种条件下均可匹配,验证了(19)式和(21)式的准确性。
图8所示为平均误码率与平均信噪比在不同大气湍流强度及衰落指数影响下的变化关系。接收孔径取D=20 cm且不考虑指向误差影响,数值仿真结果由(19)式计算得到。由图8可知,在强RF链路衰落(m=1)及中等湍流强度条件影响下的系统误码性能优于相应的弱RF链路衰减(m=3)及强湍流影响下得到的数值结果,由此可以看出,FSO信道对混合系统性能起主要作用。当平均信噪比取30 dB时,在中等湍流强度 (α=3.57, β=2, η=0.74) 及强RF信道衰落(m=1)条件下,系统误码率为2.8×10-4;而当信道条件为强湍流(α=5.97, β=0.45, η=0.1)及弱RF信道衰落(m=3)时,系统误码率增长为2×10-3。
图 7. 中等强度湍流条件下平均误码率随平均信噪比变化规律
Fig. 7. Average BER versus average SNR for moderate turbulence conditions
下载图片 查看所有图片
图 8. 不考虑指向误差时平均误码率随平均信噪比变化规律
Fig. 8. Average BER versus average SNR with pointing errors not considered
下载图片 查看所有图片
6 结论
分析了基于解码-转发中继方式的混合RF/FSO航空通信系统中断及误码性能,其中FSO信道服从指数型Weibull衰落分布模型,RF信道服从Nakagami-m衰落模型。基于Meijer's G函数推导得到混合系统端到端信噪比的累积分布函数表达式,并得到了中断概率及误码率的具体表达式。通过仿真验证了数值计算方法的准确性,分析了有无指向误差影响下,不同湍流强度条件下孔径平均效应对混合RF/FSO通信系统性能的影响。仿真结果表明,指向误差严重影响系统中断特性,而孔径平均效应可明显改善混合系统的通信性能,另外大气湍流FSO信道在混合RF/FSO航空通信系统中起主要作用。
赵静, 赵尚弘, 赵卫虎, 王翔, 陈柯帆. 大气湍流和指向误差下混合RF/FSO航空通信系统性能分析[J]. 中国激光, 2017, 44(9): 0906001. Zhao Jing, Zhao Shanghong, Zhao Weihu, Wang Xiang, Chen Kefan. Performance Analysis for Mixed RF/FSO Airborne Communication Systems over Atmospheric Turbulence and Pointing Error[J]. Chinese Journal of Lasers, 2017, 44(9): 0906001.