一种频偏条件下新的Chirp-BOK分数域解调方法

宁晓燕， 聂毅帆， 孙志国， 王秋滢

(信息与通信工程学院 哈尔滨工程大学，黑龙江 哈尔滨 150001)

扩频通信是一种发射信号占用带宽远远大于发送信息所需带宽最大值的传输方式，由于其具有抗干扰能力强、隐蔽性好、易于实现码分多址等诸多优点被广泛应用于定位、通信、测距等领域。线性调频信号扩频(chirp spread spectrum, CSS)作为多种扩频技术的一种，在一定的符号间隔内，Chirp信号频率扫过一定的范围，因此占用的频带宽度比信息带宽大得多，实现了对信息的频谱扩展，因此可将其视为一种扩频信号[1]。

Tsai等[2]在匹配滤波相干解调基础上进一步分析了在Rayleigh和Rician衰落信道下采用Chirp信号抗多径干扰的优势。基于Chirp信号的多用户通信、混合扩频通信等的提出，进一步提升通信系统容量和性能[3-5]。基于Chirp信号的Rake接收的提出，进一步提升系统的抗多径性能[6]。近年来，多种技术与扩频通信技术相结合也逐渐成为研究热点。分数域傅里叶变换(fractional Fourier transform, FrFT)作为傅里叶变换的推广，其以Chirp信号为基函数，因此Chirp信号在适当的分数域出现聚敛现象[7-8]。在文献[9-10]中，利用了Chirp信号分数域特性进行参数估计；在文献[11]中提出将Chirp信号用作OFDM系统的同步信号；在文献[12]中提出了一种基于FrFT的Chirp信号扩频调制算法，以二进制数字调制解调为例，探讨了基于分数域傅里叶变换的Chirp非相干解调的相关理论问题。通过与匹配滤波解调法的比较可知，FrFT非相干解调的优势是不受到相位误差的影响[11]，且在对抗多普勒频移和码同步误差的影响时与时域相干解调是相同的[12]；分数域非相干解调的性能劣势是没有利用到相位信息导致的。

为了保留FrFT非相干解调对抗相位偏差的优势，并且提升对抗多普勒频移和码同步误差的能力，本文在FrFT非相干解调的基础上进行了改进，提出了峰值搜索-FrFT非相干解调，与传统方法进行了对比分析，证明了该方法在频偏条件下具有更好的误码率性能。

1 Chirp-BOK调制解调

1.1 Chirp-BOK调制方案

基于二进制正交键控(binary orthogonal keying，BOK)，提出了Chirp-BOK调制方式。它利用了Chirp信号的扫频特性，对Chirp信号的正、反调频率进行信息携带处理，达到信号调制解调的目的，这种方式具有结构简单、低复杂度、低功耗等特点，是Chirp信号的常用调制方式[13-14]。设k为正整数，本文使用调频率为k的升频信号(Up-Chirp)表示“1”，使用调频率为-k的降频信号(Down-Chirp)表示“0”。Chirp信号表达式为：

(1)

Chirp-BOK调制框图如图1所示。

图1 Chirp-BOK调制Fig.1 The bock diagram of Chirp-BOK modulation

1.2 解调方案

1.2.1 匹配滤波相干解调

Up-Chirp与调频率相反的Down-Chirp信号互为匹配滤波函数[11]，因此接收端解调框图如图2所示。

图2 匹配滤波解调Fig.2 The bock diagram of MF demodulation

以升频信号为例，匹配滤波输出如图3(a)所示，非匹配滤波输出如图3(b)所示。

图3 滤波器输出波形Fig.3 Output waveform of filter

1.2.2 分数域非相干解调

分数域傅里叶变换可以看作是传统傅里叶变换的一种推广。正余弦函数是傅里叶变换的基函数，而Chirp函数正是其调频率对应阶次的分数域傅里叶变换的基函数[15]。分数域傅里叶变换是公认的适合处理Chirp信号的工具，分数域傅里叶变换的定义式[16]为：

(2)

式中：下标p为分数域傅里叶变换的阶次；α=pπ/2；u为分数域坐标；t为时域坐标；算子核K(α;u,t)定义为:

2utcscα+t2cotα)]

(3)

调频率为k的Chirp信号经过p=-2arccot(k/π)阶FrFT后出现聚敛峰值，该阶次p被称为Chirp信号的最优阶次[7]；由高斯分布线性变换不变性，高斯白噪声经过p阶分数域傅里叶变换后仍然是高斯白噪声。

根据欧拉公式可以将调频率为k的实Chirp信号拆分为2个调频互为相反数的复Chirp信号：

(4)

同理，调频率为-k的实Chirp信号拆分为2个调频互为相反数的复Chirp信号：

(5)

图4 p阶分数域谱Fig.4 Impact of p order fractional domains

1) 定点采样解调。

根据Chirp信号的分数域性质，对接收到的Chirp信号进行分数域傅里叶变换，然后采样判决[17-18]。FrFT非相干解调框图如图5所示。因为该方法是在固定频点进行采样，为了与本文提出的方法进行区分，称之为定点采样-FrFT非相干解调。

2) 峰值搜索解调。

针对存在多普勒频移和码同步误差时定点采样法性能下降严重的问题，本文提出了峰值搜索-FrFT非相干解调。峰值搜索示意图如图6所示，图中N为奇数。当分数域峰值位置存在偏移，如果还在原来的位置进行采样判决，系统的性能将严重下降，甚至无法解调。但是峰值位置的大致区域是可以确定的，如果将定点采样扩展为区域峰值采样，只要该区域的范围超过最大偏移的范围，即可以保证峰值一定能落在该区域内。然后在采样得到的2N个值中搜索峰值，由峰值所在区域得到判决结果。

图5 定点采样-FrFT非相干解调框图Fig.5 The bock diagram of Fixed point sampling-FrFT non-coherent demodulation

峰值搜索法的解调框图如图7所示。基带信号先经过p阶FrFT，取模值后得到Chirp信号的分数域谱，以cup(t)和cdown(t)对应位置uup和udown为参考点将定点采样扩展为区域峰值搜索。

图6 分数域峰值搜索示意Fig.6 Impact of peak search

图7 峰值搜索-FrFT非相干解调框图Fig.7 The bock diagram of peak search-FrFT non-coherent demodulation

2 性能分析

2.1 AWGN信道下性能分析

2.1.1 匹配滤波法的性能分析

对于匹配滤波相干解调，根据二进制最佳接收公式可以得出其误码率公式：

(6)

式中：ρ为升频信号和降频信号的相关系数；Eb为码元能量；n0为带限高斯白噪声的单边功率谱密度。

2.1.2 定点采样法的性能分析

以调频率为k的Chirp信号为例，经过AWGN信道后接收到的信号为：

r(t)=Ae2πf0t+πkt2+n(t)ejθ

(7)

接收信号r(t)经过p阶FrFT后，在u1处采样可以得到：

Rp(uup)=Aejθ0+n(uup)ejθ1

(8)

而在u2处采样可以得到：

Rp(udown)=ξAejθ0+n(udown)ejθ2

(9)

式中：A为经过最优阶次变换的分数域采样值；ξA是经过非最优阶次变换的分数域采样值，ξ和θ0是与Chirp信号具体参数有关的值；θ1和θ2是由噪声引入的随机相位。

调整相位可以得到：

(10)

两路信号取模值可以得到:

(11)

式中abs(·)表示取模值。

由随机信号分析可知，abs(Rp(uup))的抽样值Rup和abs(Rp(udown))的抽样值Rdown服从莱斯分布。它们的概率密度分布分别为：

(12)

式中I0(·)为第1类零阶修正贝塞尔函数。

以升频信号为例，定点采样法判决正确的概率为：

(13)

当发送信号为降频信号时，判决正确的概率Pdown=Pup。

当发送升频信号和降频信号的概率相等时，定点采样-FrFT非相干解调的误码率为:

(14)

2.1.3 峰值搜索法的性能分析

对于峰值搜索法，由2点比较判决扩展为2N个点互相比较判决。首先对存在峰值区域的N个点进行搜索，寻找最大值。在AWGN信道且同步准确的情况下，第(N+1)/2个点即峰值所在位置。然后对不存在峰值区域的N个点进行搜索，寻找最大值。将这2N个点记为x1,x2,…,xi,…,x2N，xi的采样值为：

xi=ξiA+niejθi

(15)

式中：ξiA为每个采样点的幅值，ξi和θi是与Chirp信号具体参数有关的值，0<ξi≤1。abs(xi)服从莱斯分布，它们的概率密度函数为:

(16)

当xi,i=1,2,…,N为最大值时可以正确解调，则峰值搜索法判决正确的概率为：

xi>x2N|xi)f(xi)d(xi)

(17)

式中2N个采样点之间相互独立，互不影响，式(17)可以化简为：

(18)

式中:

(19)

经过化简，可得峰值搜索法误码率为：

(20)

2.2 频偏对分数域解调法的影响

以调频率为k的Chirp信号为例，分析频率偏移对Chirp信号FrFT非相干解调的影响，以证明将定点采样改进为区域峰值搜索的必要性和优越性。

首先，频偏主要来源于位同步误差、多普勒频移和收发端接收机参数差异。当存在同步误差时，本地同步时间与信号实际到达时间存在时刻τ的延时或超前，信号表达式为：

c(t)=cos(2πf0(t+τ)+πk(t+τ)2)

(21)

把τ看做定值，则上式可以化简为:

c(t)=cos(πfmt+πkt2+Δ)

(22)

式中：fm=f0+kτ/2，Δθ=2πf0τ+πkτ2。可以看到，Chirp信号的时间τ的同步误差可以转化为频率的频移kτ/2和相位的偏移Δθ；同时，由于存在时间偏差，因此只有部分信号被利用到，造成了能量的损失。

当存在多普频移时， Chirp信号中心频率发生偏移，因为入射角的不同多普勒频移的值将会出现变化，取最大值时Chirp信号表达式为：

c(t)=cos(2π(f0+fd)t+πkt2+φ)

(23)

可以看到，多普勒频移造成的影响首先是中心频率偏移fd，与同步误差造成的影响是类似的。

下面分析频偏对采样点的具体影响。首先将同步误差和多普勒频移造成的频偏统一记为fd。根据1.2.2节中的分析，以升频信号为例，在p阶分数域的峰值聚敛位置为uup=f0sinα，能量集中在|u-uup|

c(t)=cos(2πfmt+πkt2)

(24)

新的峰值聚敛位置为um=(f0+fd)sinα，偏移距离fdsinα。当频偏较大时，能量聚敛位置甚至会完全偏离该区域，如果仍以uup为采样点，则无法正确解调出信息。

当采样点移动的区域仍在|u-uup|

c(t)=exp(j2πfmt+jπkt2)

(25)

将c(t)代入分数域变换公式得：

(26)

当k=-cotα，u=f0sinα时，可得：

(27)

化简得：

(28)

取包络得：

(29)

可以看出abs(Cp(u))的值与T和fd有关，当T固定时，fd越大，采样值越小。

综上，同步误差、多普勒频移和接收机收发3端参数差异可以统一转换为中心频率的偏移和能量的损失。映射到分数域，对应着峰值位置的偏移和采样值的降低，严重影响了定点采样-FrFT非相干解调的性能。因此本文提出了峰值搜索-FrFT非相干解调来改善频偏下的性能。

引入频率偏移后，峰值搜索-FrFT非相干解调仍然为2N个点互相比较判决，只是峰值所在位置可能不再是第(N+1)/2个点。推导过程与AWGN信道下相同，故误码公式也与式(20)相同。系数ξi与引入频偏的Chirp信号分数域采样值有关。

3 仿真验证

3.1 AWGN信道下性能仿真

图8对比了Chirp-BOK调制的3种解调方式在AWGN信道下的误码率曲线。仿真参数：带宽20 MHz、码长10 μs、搜索点数为3个采样点。从图8中可以看出相干解调的性能最好；而峰值搜索-FrFT方法的性能相比于定点采样-FrFT性能下降不超过0.5 dB。

图8 AWGN信道下3种解调方式性能对比Fig.8 Performance comparison of three demodulation methods in AWGN channel

图9给出了AWGN信道下峰值搜索-FrFT的理论曲线和仿真曲线。搜索区域长度分别设置为3、21个采样点。从图9中可以看出，当搜索区域的长度越小时，本方法在AWGN信道下的性能越好。因此搜索区域范围的设置要根据频偏动态地调整，使得搜索范围尽可能小。

图9 不同长度搜索区域的误码率曲线Fig.9 Performance comparison of different length search ranges

3.2 抗频偏性能分析

图10和图11给出了存在频偏时3种解调方式的性能对比。通过仿真可以看出，当归一化频偏比为0.001 5时，峰值搜索-FrFT的误码率性能优于定点采样-FrFT与相干解调。当频偏比增大至0.002 5时，定点采样-FrFT与相干解调的性能迅速下降，但是峰值搜索-FrFT的性能却几乎不受影响。

图10 归一化频偏0.001 5时3种解调方式性能对比Fig.10 Performance comparison of three demodulation methods under normalized frequency offset 0.001 5

图11 归一化频偏0.002 5时3种解调方式性能对比Fig.11 Performance comparison of three demodulation methods under normalized frequency offset 0.002 5

4 结论

1)该方法将传统的定点采样改进为采样点周围区域的峰值搜索，以提升对位同步误差、多普勒频移和收发端参数差异导致的信号频率偏移的容忍度。

2)推导了该解调方法下AWGN信道的误码率公式，在AWGN信道下的性能主要受搜索长度的影响，当搜索长度越短时，性能越接近定点采样-FrFT；当存在频偏时，峰值搜索-FrFT可以对峰值位置进行搜索，只要满足搜索范围大于偏移范围，则可以极大程度地改善误码率性能，而匹配滤波相干解调和定点采样-FrFT的性能严重恶化。