Performance Analysis of DCO-OFDM and ACO-OFDM Systems in Indoor Visible Light Communications
1 引言
可见光通信(VLC)具有照明和提供通信的双重功能。与传统的射频(RF)通信相比,VLC可提供更高的数据率,通信保密性好,在邻近房间可以实现频率复用,拥有不经许可即可使用的近乎无限的频谱资源,对人体健康几乎没有危害。此外,VLC不会与RF相互干扰,适用于医院、矿井、飞行器等对电磁干扰敏感的区域[1]。
VLC一般设计为强度调制直接检测(IM/DD)系统,即信息只包含在光信号的强度中。由于光电检测器(PD)的尺寸通常比可见光波长大上千倍,光信号在PD表面就形成了类似空间分集的效果,因此不存在多径衰落现象。干扰对每条光信号的影响相当于被求取空间平均,所以VLC信道可以看作时间稳定的、随PD位置变化而缓慢变化的信道。VLC系统常作为线性时不变系统处理,完全可以用冲激响应来描述。在研究和设计VLC时,只有对室内光无线传播信道的特征有充分的了解,才能确保所设计的通信系统有令人满意的性能。到目前为止,尚没有公认的、能够包括复杂反射的高效专用VLC信道表征方案[2]。现有研究大多是基于红外光通信的信道模型,国内外学者提出了多种室内红外光无线信道建模方法[3-9],文献[
1]则给出了VLC中视线传输(LOS)和一次反射时的直流增益计算方法。随着通信符号速率的提高,光信号经过室内墙壁、屋顶等反射物反射后所引起的码间干扰(ISI)将不可忽视。为了准确评价多径效应对系统性能的影响,必须建立室内VLC多径信道模型。
另一方面,IM/DD系统要求传输的信号为非负实信号,单极性的脉冲位置调制(PPM)和脉冲幅度调制(PAM)可以满足该要求,但是当系统符号速率较高时,严重的ISI会使系统的性能降低。此外,室内荧光灯等常用人造光源所产生的窄带干扰会对PAM和PPM产生影响。正交频分复用(OFDM)技术可以有效地解决光信号漫射引起的ISI和窄带干扰等问题,同时可以提高系统的频谱利用率。
本文研究了两种最常用的OFDM技术,即直流偏置光OFDM(DCO-OFDM)和非对称限幅光OFDM(ACO-OFDM)[10-12]。首先在单光源VLC系统中,提出了一种多径信道建模方法。该方法用迭代法[13]计算LOS和反射传输冲激响应,得到了接收端抽样周期和发光二极管(LED)调制带宽之间的关系,并定义了ISI。以LOS信道时延作为建模起点,将接收端相邻抽样间隔之间的所有冲激响应之和作为多径信道路径增益,建立了多径信道模型。其次,考虑非线性限幅噪声时,在多径信道下,推导了DCO-OFDM和ACO-OFDM的理论误码率(BER)。最后,分析PD室内位置、LED发射角、PD的视场角(FOV)和直流偏置等参数变化时,DCO-OFDM和ACO-OFDM系统的性能变化。
2 VLC多径信道模型
采用单光源的室内VLC系统的几何模型如图1所示。屋顶安装用于照明和通信的LED,PD随机分布在室内工作平台上,光信号经过LOS和反射传播入射到PD。
图 1. 室内VLC几何模型
Fig. 1. Geometrical model of indoor VLC
下载图片 查看所有图片
当信源S、接收器R和室内反射环境一定时,包含多次反射的信道冲激响应可表示为
式中h(k)(t;S,R)为光信号经过k次反射的信道冲激响应,k=0表示LOS信道。当LED服从朗伯辐射模式时,LOS信道冲激响应为
式中d为LED到PD的距离,AR为PD的面积,ϕ为LOS光线的出射角,φ为光信号的入射角,ΨFOV为PD的视场角,c为光速,δ(x)为狄拉克函数,κ=-ln2/ln(cosθ1/2)为表征光源辐射方向性的辐射模式指数,θ1/2为光源半功率角,rect(·)为矩形函数,其定义为
用迭代法[13-14]计算经过k次(k>0)反射的信道冲激响应为
式中r为微反射单元的位置矢量,
为r处微反射单元的单位法向矢量,符号⊗代表卷积运算。对反射面
上的所有微反射单元进行积分,实际计算时,将所有反射平面划分为面积为ΔA的小反射单元,积分运算数字化后得到
式中Nref为反射单元的总数,ρi为第i个反射单元的反射率,D为信源到反射点的距离。
从最先到达PD的第1路光信号开始,时间延迟大于符号周期一半的光信号将引起ISI[15],建立多径信道模型原理图如图2所示,图中Tsp=Tsym/2为接收端信号采样周期,Tsym为驱动LED发光的符号周期。根据奈奎斯特定理,Tsym最小值为WLED/2,其中WLED为LED的调制带宽,τ0=d/c为最先到达PD的光信号时延。
图 2. 多径信道建模原理图
Fig. 2. Schematic of the multipath channel modeling
下载图片 查看所有图片
建立多径信道增益矢量为
式中[·]T为矩阵的转置,L为多径信道路径总数,第l路信道路径增益表示为
3 光OFDM系统
3.1 发送端
可见光DCO-OFDM和ACO-OFDM系统原理如图3所示。信源经M阶正交幅度调制(MQAM),调制符号满足E[|X(n)|2]=1,其中E[]为功率归一化函数,对调制符号X(n)进行映射操作,ACO-OFDM调制时输出为[
]T,其中,N为快速傅里叶逆变换(IFFT)的长度,*表示共轭运算,所有的奇数位子载波包含信息,偶数位子载波都为零。DCO-OFDM调制时映射输出为[
]T。显然,映射信号Xmapping满足厄米对称性,即
图 3. 可见光OFDM系统原理框图
Fig. 3. Diagram of visible light OFDM system
下载图片 查看所有图片
频域映射信号输入到IFFT模块,输出时域信号xIFFT表示为
式中H表示矩阵共轭转置,F为N×N的归一化离散傅里叶变换矩阵。
在ACO-OFDM系统中,xIFFT满足反对称性[12],即有
式中k=0,1,…,N/2-1。
在DCO-OFDM系统中,xIFFT(k)可表示为
通常为了避免直流偏移,令Xmapping(0)=Xmapping(N/2)=0,因此xIFFT(k)为双极性的实数信号。
在IM/DD系统中,要求驱动LED的信号为非负实数信号。根据中心极限定理(CLT),当N≥64时,xIFFT服从高斯分布,因此需要对xIFFT进行限幅。限幅信号经并串(P/S)转换后,加上循环前缀(CP)以消除ISI的影响。在ACO-OFDM中,经数/模(D/A)转换的信号xACO(t)直接驱动LED。在DCO-OFDM中,需要加上直流偏置信号BDC,得到LED的前置驱动信号xDCO(t)=xclip(t)+BDC。
限幅对两种调制系统的影响说明如下。 在ACO-OFDM系统中,xIFFT有反对称性,因此将xIFFT中小于零的值删除,可得到单极性实数信号xclip,限幅输出信号可表示为[12]
式中nclip(k)为限幅噪声。由文献[
12]可知,限幅后的信号所对应的频域信号为
由(13)式可以看出,有用的信号在奇数位子载波上,其幅度是未限幅频域信号的一半,而限幅产生的噪声都落在偶数位子载波上。
在DCO-OFDM中,限幅操作的上下门限分别为εbottom=-BDC,εtop=¥。根据Bussgang理论[16],将高斯分布信号限幅表示为
式中nclip为与xIFFT不相关的非高斯分布的限幅噪声,η为信号衰减因子,可表示为
式中
为xIFFT的方差,λbottom=εbottom/σ0为归一化限幅门限,Q(υ)=(1/)
exp(-u2/2)du为互补累积函数。
通常直流偏置信号可以表示为
式中μ为比例常数,用分贝(dB)表示BDC的大小为10lg(μ2+1) dB。
3.2 接收端
光信号经过室内VLC多径信道传输后,被PD检测到并转换为电信号,表示为
式中γ为光电转换因子,不失一般性,设γ=1 A/W;n(t)为背景光与电路热噪声之和,是独立于信号且服从高斯分布的加性高斯白噪声(AWGN)[17]。假设噪声单边功率谱密度为N0,经过A/D转换和S/P转换后删除CP,由于CP的长度大于等于信道冲激响应的长度,(17)式中线性卷积变为循环卷积,可得到
式中xLED和y分别为长度为N的发送和接收符号矢量,n为噪声矢量。
为(6)式中多径矢量h所对应的循环矩阵[14]。
y信号输入到N点快速傅里叶变换(FFT)模块,输出的频域信号分别为
式中H=FFH为频域信道响应,是一个N×N的对角矩阵,用H(n)表示其主对角线上的第n个元素;N为方差
=WN0的AWGN矢量的傅里叶变换,其中W为OFDM信号的带宽;B=[
]T为BDC的傅里叶逆变换;NDCO,clip和NACO,clip分别为非高斯分布的随机限幅噪声的傅里叶变换。根据CLT,非高斯分布的随机限幅噪声转变为高斯分布的信号,其均值为零,方差为
。叠加在DCO-OFDM载荷信息的子载波上的噪声方差[16]为
式中φ(u)=(1/
)exp(-u2/2)。可以看出,限幅噪声方差与IFFT/FFT的长度N和调制阶数M无关,仅与归一化门限λbottom有关。
与发送端相对应,在ACO-OFDM系统中,提取YACO的前一半奇数位子载波信号,得到长度为N/4的提取信号,
式中N(m)为PD在第m个子载波上的噪声。因为限幅噪声项不对奇数位子载波产生影响,所以(23)式去掉了限幅噪声项。
在DCO-OFDM系统中,提取YDCO的1~N/2-1个子载波信号作为解调信息的信号,得到长度为N/2-1的提取矢量,
因为直流偏置信号的傅里叶变换只在第0个子载波非零,对提取的子载波没有影响,因此直流项被删除。
为了均衡多径信道的影响,采用迫零(ZF)检测算法。ACO-OFDM系统中,输入到检测器的信号为
式中Ysd,ACO为(N/4)×(N/4)维等效频域传输矩阵;Xex,ACO为Xmapping的第m行元素组成的(N/4)×1维矩阵,m=1,3,…,N/2-1;
为均值为零、方差为1的高斯噪声。
DCO-OFDM系统中,输入到检测器的信号为
式中Hsd,DCO为(N/2-1)×(N/2-1)维等效频域传输矩阵;Xex,DCO为Xmapping取m=1,2,3,…,N/2-1行元素组成的(N/2-1)×1维的矩阵;
为均值为零、方差为1的高斯噪声。
ZF检测器的加权矩阵为
式中(·)-1表示矩阵求逆。
在ACO-OFDM系统中,信源符号X(l),l=0,1,2,…,N/4-1的估计符号为
式中WZF,ACO[l,:]表示(N/4)×(N/4)的ZF检测器加权矩阵的第l行。
在DCO-OFDM中,WZF,DCO为(N/2-1)×(N/2-1)维的加权矩阵。信源符号X(l),l=0,1,2,…,N/2-1的估计符号为
将估计符号输入到MQAM解调器,解调出二进制序列。采用最大似然检测的MQAM解调BER[18]为
式中erfc(·)为误差函数;ΓSNR为信噪比,表示MQAM解调器输入的符号能量和噪声功率谱密度之比。由(29)、(30)式计算
(l)的信噪比分别为
将ΓACO(l)和ΓDCO(l)分别代入(31)式可以得到
(l)的BER,系统的总BER为所有数据流BER的平均值, 即
4 数值仿真和分析
在长、宽和高分别为6,6,4 m的房间内,安装垂直指向地面的距屋顶中心0.5 m的LED,PD位于高度为0.85 m的工作平台上,垂直指向屋顶。将墙面在三维坐标方向上按间隔0.1 m划分成小反射单元,其他仿真参数如表1所示。
假设LED调制带宽为50 MHz, LED的调制符号周期为10 ns,则相对于最先到达PD的光信号,时间延迟大于接收端采样周期(Tsp=5 ns)的光信号就认为发生了ISI。在接收端接收的光功率中LOS信道和光信号一次反射约占所有接收光功率的90%。为了简单起见,仅考虑LOS和一次反射光功率。图4所示为PD位于室内3个典型位置(3,3,0.85)、(1.5,1.5,0.85)和(0.5,0.5,0.85)时,LED和PD之间的多径信道模型。可以看出,当PD在房间中心时,LOS路径延迟小,多径分量h0(第1径)较大,多径分量衰减快;当PD在房间墙角时,LOS路径延迟变大,多径分量h0变小,多径信道路径分量衰减较慢,且相对于h0较大。
表 1. 仿真参数
Table 1. Simulation parameters
Parameter | Value |
---|
LED modulation bandwidth BLED /MHz | 50 | Semi-angle at half power θ1/2 /(°) | 60 | FOV at PD ΨFOV /(°) | 80 | Detector physical area of PD AR /cm2 | 1 | Reflectivity of wall ρi | 0.8 | Reflectivity area of small region ΔA /m2 | 0.01 |
|
查看所有表
图 4. PD位于不同位置时的多径信道增益。 (a) (3,3,0.85);(b) (1.5,1.5,0.85);(c) (0.5,0.5,0.85)
Fig. 4. Multipath channel gain for different positions of PD. (a) (3, 3, 0.85); (b) (1.5, 1.5,0.85); (c) (0.5, 0.5, 0.85)
下载图片 查看所有图片
ACO-OFDM和DCO-OFDM系统分别采用4QAM、16QAM和64QAM调制,IFFT/FFT长度为N=256。图5为ACO-OFDM系统BER性能,可以看出,随着M的增大,系统BER性能变差。当PD在房间中心位置时,信道增益较强,第1径h0最大,其多径分量相对很小,多径效应对系统性能影响小,系统性能最好。当PD移动到墙角位置(0.5,0.5,0.85)时,与PD在位置(1.5,1.5,0.85)相比较,大量反射路径变短,多径分量中第1径h0相比于其他路径较大,多径效应影响较小,BER性能比PD位于位置(1.5,1.5,0.85)时好。
图 5. ACO-OFDM系统BER性能比较
Fig. 5. Comparison of BER performance of ACO-OFDM system
下载图片 查看所有图片
图6和图7所示为BDC=7 dB和BDC=13 dB时DCO-OFDM系统的BER性能。可以看出,当BDC=7 dB时BER性能比BDC=13 dB时好。这是由于信噪比一定时,直流偏置越大,信号中包含信息的分量越小,虽然大的直流偏置对应的限幅噪声较小,但此时限幅噪声对系统性能的影响相对很小。另一方面,给定接收位置时,ACO-OFDM系统的BER性能比DCO-OFDM系统好。
图 6. BDC=7 dB时DCO-OFDM系统BER性能比较
Fig. 6. Comparison of BER performance of DCO-OFDM system when BDC=7 dB
下载图片 查看所有图片
图 7. BDC=13 dB时DCO-OFDM系统BER性能比较
Fig. 7. Comparison of BER performance of DCO-OFDM system when BDC=13 dB
下载图片 查看所有图片
图 8. FOV变化时BER性能比较
Fig. 8. Comparison of BER performance with different FOV
下载图片 查看所有图片
图 9. LED半功率角变化时BER性能比较
Fig. 9. Comparison of BER performance with different θ1/2 of LED
下载图片 查看所有图片
图8和图9所示为PD在(1.5,1.5,0.85)位置时, ACO-OFDM系统的BER性能和BDC=7 dB时DCO-OFDM系统的BER性能随FOV和LED半功率角的变化。可以看出,随着FOV的增加,BER性能变差,这是由于有更多的反射信号被PD接收,多径增益相对h0变大,多径效应的影响增大。另一方面,随着LED半功率角增大,BER性能也变差,这是因为半功率角变大时,发射光信号越发散,反射信号越强,多径干扰越大。
5 结论
结合LED的调制带宽,介绍了一种单光源室内VLC多径信道建模方法。随着接收端位置在室内的变化,对应的信道模型也相应变化,信道模型的多路径增益反映了反射光信号的情况,以及多径干扰对系统性能的影响。在多径信道下研究了两种常用光OFDM系统的BER性能,建立了Matlab蒙特卡罗仿真模型,分析了PD的FOV、LED半功率角和直流偏置等参数对系统的性能影响,仿真结果和理论分析相吻合,验证了理论的正确性,也验证了信道建模的正确性和实用性。
贾科军, 靳斌, 郝莉, 张守琴. 室内可见光通信中DCO-OFDM和ACO-OFDM系统性能分析[J]. 中国激光, 2017, 44(8): 0806003. Jia Kejun, Jin Bin, Hao Li, Zhang Shouqin. Performance Analysis of DCO-OFDM and ACO-OFDM Systems in Indoor Visible Light Communications[J]. Chinese Journal of Lasers, 2017, 44(8): 0806003.