CN110730149B - 一种平移式特殊qpsk系统联合捕获方法 - Google Patents

一种平移式特殊qpsk系统联合捕获方法 Download PDF

Info

Publication number
CN110730149B
CN110730149B CN201910947509.XA CN201910947509A CN110730149B CN 110730149 B CN110730149 B CN 110730149B CN 201910947509 A CN201910947509 A CN 201910947509A CN 110730149 B CN110730149 B CN 110730149B
Authority
CN
China
Prior art keywords
signal
sequence
frequency offset
power detection
received
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910947509.XA
Other languages
English (en)
Other versions
CN110730149A (zh
Inventor
姚如贵
王圣尧
吴啟洪
左晓亚
张远
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201910947509.XA priority Critical patent/CN110730149B/zh
Publication of CN110730149A publication Critical patent/CN110730149A/zh
Application granted granted Critical
Publication of CN110730149B publication Critical patent/CN110730149B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/32Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/3405Modifications of the signal space to increase the efficiency of transmission, e.g. reduction of the bit error rate, bandwidth, or average power
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/32Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/3405Modifications of the signal space to increase the efficiency of transmission, e.g. reduction of the bit error rate, bandwidth, or average power
    • H04L27/3411Modifications of the signal space to increase the efficiency of transmission, e.g. reduction of the bit error rate, bandwidth, or average power reducing the peak to average power ratio or the mean power of the constellation; Arrangements for increasing the shape gain of a signal set
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/32Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/38Demodulator circuits; Receiver circuits
    • H04L27/3809Amplitude regulation arrangements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/32Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/38Demodulator circuits; Receiver circuits
    • H04L27/3818Demodulator circuits; Receiver circuits using coherent demodulation, i.e. using one or more nominally phase synchronous carriers
    • H04L27/3836Demodulator circuits; Receiver circuits using coherent demodulation, i.e. using one or more nominally phase synchronous carriers in which the carrier is recovered using the received modulated signal or the received IF signal, e.g. by detecting a pilot or by frequency multiplication
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/32Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/38Demodulator circuits; Receiver circuits
    • H04L27/3845Demodulator circuits; Receiver circuits using non - coherent demodulation, i.e. not using a phase synchronous carrier
    • H04L27/3854Demodulator circuits; Receiver circuits using non - coherent demodulation, i.e. not using a phase synchronous carrier using a non - coherent carrier, including systems with baseband correction for phase or frequency offset

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明提供了一种平移式特殊QPSK系统联合捕获方法,在发送端,采用帧同步序列对信号的幅度进行调制,将帧同步序列隐藏在信息序列中;在接收端,首先对捕获信号进行功率检测,当进行了设定次数的功率检测或者功率检测捕获成功后进行频偏估计,对经过频偏校准后的信号作相关检测,当相关检测结果满足设定条件时,接收端成功捕获到接收信号的帧起始位置,对接收信号进行解调,解调过程中,在抽样判决的时候加入m序列进行辅助判决。本发明使系统在存在多普勒频移的情况下仍然具有良好的捕获性能和系误比特率。

Description

一种平移式特殊QPSK系统联合捕获方法
技术领域
本发明属于无线通信技术领域,针对将帧同步序列隐藏于信息序列中的平移式特殊QPSK调制解调系统,提出了一种功率检测、频偏估计和相关检测相结合的联合捕获方案,从而使系统具有良好的捕获性能。
背景技术
在数字通信系统中,尤其是突发通信系统,为了在接收端能够正确识别数据帧的起始,通常需要在发送端插入每帧数据的起始标记,在接收端检测并获取这一标志的过程称为帧同步。帧同步的主要作用是提取数据流中每一帧的帧头位置,并将成帧的数据按照要求作适当处理后送给后端模块作数据处理。因此,同步帧的捕获是通信系统中一项关键技术。
在无线通信系统中,由于发送机和接收机之间的长距离和高速移动会导致动态多普勒频移,严重影响捕获性能和误比特率性能,因此多普勒频移效应是导致系统性能退化的主要因素之一。在实际通信系统中,如何在存在多普勒频移的情况下,对接收信号进行频偏估计和同步帧捕获是信号正确恢复的前提条件。
文献1“崔冕,王宇.一种中高速突发QPSK数字解调系统捕获与载波同步算法的研究[J].空间电子技术,2010,7(1):30-33.”针对一种中高速突发解调处理系统,研究了突发捕获与载波同步算法实现,采用相位域突发捕获与同步算法,解决了卫星上行信号幅度变化对突发捕获解调性能的影响,具有突发捕获概率高、载波同步快和解调动态范围大等优点。
文献2“Zeng R,Huang H,Yang L,et al.Frequency Offset and Doppler ShiftJointly Estimation Algorithm in High Mobility Environment Based on OrthogonalAngle Domain Subspace Projection[J].IEEE Transactions on VehicularTechnology,2018,67(3):2254-2266.”提出了一种基于正交角度域子空间投影的高移动环境下联合频偏和多普勒频移估计算法,该算法能够跟踪快速时变多普勒频移,并将多普勒频移与频偏相区分。
文献3“Zhang Z,Liu H,Zhang H.Joint-Detection for Doppler-ShiftAcquisition in Space Communications[J].IEEE Communications Letters,2015,19(2):215-218.”研究了一种存在多普勒频移的通信系统中的帧捕获算法,该算法使用多个组合周期来缩小多普勒频移捕获的搜索范围,以提高捕获概率。
文献4“张毅,欧阳志新,邓云凯.高性能无数据辅助QPSK频偏估计新算法[J].西安电子科技大学学报,2013,40(2):187-193.”提出了一种新的QPSK频偏估计算法,并结合离散傅里叶变换粗频偏估计方法对大频偏进行估计,提高了对频偏变化的适应能力,使之在频偏缓变的情况下依然可以准确估计频偏。
由于特殊QPSK的特殊性,当存在多普勒频移时,上述文献所提同步帧捕获方法将不再适用,多普勒频移将会对捕获性能和误比特率性能造成严重影响。在先的专利申请“一种帧同步码隐含于信息码中传输的方法,201910174035.X,2019-03-08.”中针对辅助的特殊同步码序列或者是循环前缀序列占用信道资源的问题,提出了一种对称式特殊QPSK映射方法。通过改变信息码元的能量,将帧同步码隐藏于信息码中,释放了帧同步码序列所占用的信道资源,实现了同步序列对信道资源的零占用,有效提高了信道资源的利用率。但由于特殊QPSK的特殊性,当接收信号存在频率偏移时,隐藏于符号中的m序列将会受到严重的影响。较小频偏导致的符号旋转都将会对m序列提取的造成很大干扰,频偏稍大时完全无法进行正常捕获。
发明内容
为了克服现有技术的不足,本发明提供一种平移式特殊QPSK系统联合捕获方法,能够使平移式特殊QPSK系统在存在频率偏移时仍具有良好的捕获性能。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤一,在发送端,采用帧同步序列对信号的幅度进行调制,将帧同步序列隐藏在信息序列中;
步骤二,在接收端,首先对捕获信号进行功率检测,当进行了设定次数的功率检测或者功率检测捕获成功后进行频偏估计,对经过频偏校准后的信号作相关检测,当相关检测结果满足设定条件时,接收端成功捕获到接收信号的帧起始位置;
步骤三,对接收信号进行解调,解调过程中,在抽样判决的时候加入m序列进行辅助判决。
所述的步骤一中,设信息序列为{b0,b1,b2,…,bN-1},帧同步序列m为{m0,m1,…,mN-1},其中,N为数据帧长度且为偶数,bi,mi∈{0,1},i=0,1,…,N-1,b′i=2bi-1,m′i=2mi-1∈{-1,1};发射端采用平移式特殊QPSK调制,第k个特殊QPSK调制符号的同向分量为sI(k)=a×b′2k+m′2k×d,正交分量为sQ(k)=a×b′2k+1+m′2k+1×d,其中k=0,1,…,N/2-1,d代表幅度偏移量,a代表QPSK的星座点参数;发送的调制符号
Figure BDA0002224573220000031
其中,
Figure BDA0002224573220000032
Figure BDA0002224573220000033
Figure BDA0002224573220000034
为由信息序列决定的调制相位,
Figure BDA0002224573220000035
θk是由同步序列决定的相位,
Figure BDA0002224573220000036
且θk∈{-θ0,0,θ0},其中θ0=arctan(a2/a1)-π/4。
所述的步骤二中,在接收端,假设接收信号经过下变频后存在Δf的频率偏移,则第k个接收符号
Figure BDA0002224573220000037
其中频率偏移引起的角频率为Δω=2πΔf;θk是由同步序列决定的相位;n0(k)为叠加在第k个接收符号上的复高斯白噪声;
功率检测模块采用一个FIFO结构,取信号的绝对值代替信号的平方,检测一个帧长度内的信号的功率;功率检测模块使用一个定时计数变量n,在功率检测中用于判断是否进行频偏估计;当
Figure BDA0002224573220000038
或者功率检测捕获成功时,令n=0,输出频偏估计的启动信号,开始执行频偏估计算法;否则,令定时计数变量n值加1,继续执行功率检测;
频偏
Figure BDA0002224573220000039
其中fs是符号速率;
设经过频偏校准后的信号为y′(k),y′I(k)和y′Q(k)分别为其实部和虚部;在接收端接收信号时,提取出隐藏在信号中的m序列,其第2k个和第2k+1个比特对应的符号分别估计为
Figure BDA00022245732200000310
当相关检测结果出现相关峰且相关峰值大于判决门限时,接收端成功捕获到接收信号的帧起始位置。
所述的步骤二中,在理想信道下接收的信息序列为{s(0),s(1),s(2),…,s(k),…},实际信道下接收的信息序列为{y(0),y(1),y(2),…,y(k),…};根据理想信道下的一帧数据的接收功率估算为
Figure BDA0002224573220000041
对于正常接收,实际接收信息序列在FIFO中每隔一个采样时刻,按照先入先出的顺序逐一移位,进行一次功率检测;当处于检测时刻k0时,实际信道下的一帧数据的接收功率
Figure BDA0002224573220000042
根据功率检测的阈值Pth,当Pr(k′)>Pth时,认为数据帧从第k′个采样开始,功率检测模块输出功率捕获成功信号。
所述的步骤二中,功率检测的阈值设为Pth=0.8Pp
所述的步骤二中,设实际信道中接收信号中的m序列与本地m序列完全对齐与未对齐时的归一化相关峰值分别为Rc和Re,将提取的m序列
Figure BDA0002224573220000043
与本地m序列{m′0,m′1,…,m′N-1}进行相关检测并判决,数据帧捕获的判决门限设置为Rc≥max(0.6d,1.6max(Re))。
所述的步骤三中,设经过频偏校正后的第k个接收符号的同相分量和正交分量分别为y′I(k)、y′Q(k),则第2k个信息比特估计为
Figure BDA0002224573220000044
即当
Figure BDA0002224573220000045
判决为1,否则,判决为0;同理,第2k+1个信息比特估计为
Figure BDA0002224573220000046
本发明的有益效果是:针对平移式特殊QPSK调制解调系统,提出了一种将功率检测、频偏估计和相关检测相结合的联合捕获方案,从而使系统在存在多普勒频移的情况下仍然具有良好的捕获性能和系误比特率。提出的联合捕获方法可以按照实例中的步骤扩展到对称式特殊QPSK的应用中,不失一般性。
附图说明
图1是平移式特殊QPSK星座图;
图2是平移式特殊QPSK信号调制原理框图;
图3是平移式特殊QPSK的联合捕获方案流程图;
图4是平移式特殊QPSK信号解调原理框图;
图5是m序列提取和相关捕获原理框图;
图6是SNR=4dB时,联合捕获的相关效果图;
图7是SNR=4dB时,频偏估计的误差图;
图8是不同信噪比下的BER性能图;
图9是不同信噪比下的捕获概率图。
具体实施方式
本发明所提的方法适用于正交振幅调制(Quadrature Amplitude Modulation,QAM)、多相制(multiple phase shift keying,MPSK)等调制方式,这里仅以正交相移键控(Quadrature Phase Shift Keying,QPSK)为例对该方法进行详细介绍。
QPSK是无线通信中一种广泛应用的数字调制方式。针对同步码元序列占用信道资源的问题,本发明结合一种平移式特殊QPSK星座图,提出了一种联合功率检测、频偏估计和相关检测三部分的捕获方法,从而使平移式特殊QPSK系统在存在频率偏移时仍具有良好的捕获性能。
为了更加清楚的描述本发明提出的联合捕获方案,下面先简单介绍下平移式特殊QPSK的映射方案。
本发明研究平移式特殊QPSK在多普勒频移环境下的联合捕获。平移式特殊QPSK基于QPSK调制方式,使用m序列作为帧同步序列,通过m序列控制I、Q两路信号在初始幅度a的基础上加上或者减去幅度为d的偏移量,完成将帧同步序列隐藏于信息序列的功能。如图1所示,a1、a2表示单路信号发生偏移后的幅度,d代表幅度偏移量,a代表传统QPSK的星座点参数,其中d=(a2-a1)/2,a=(a2+a1)/2。平移式特殊QPSK调制信号的星座图以QPSK的星座图坐标E为基准,根据m序列将符号映射到E0、E1、E2或E3的位置上。其中E0、E1、E2和E3关于虚线圈E中心对称,每个象限中的映射坐标均可以由其他象限通过平移得到,而在先专利中的对称式特殊QPSK中的映射坐标是通过轴对称而来的,这也是两者的区别所在。在发射端进行信号调制时,信息序列决定了星座图上的具体象限,m序列决定了某一象限内的具体位置。在接收端进行信号解调时,星座图中的象限解调为发送的信息序列,而象限内具体位置则对应m序列。
本发明的技术方案包括以下步骤:
步骤一:平移式特殊QPSK调制
设信息序列为{b0,b1,b2,…,bN-1},m序列为{m0,m1,…,mN-1},其中,N为数据帧长度且为偶数,bi,mi∈{0,1},i=0,1,…,N-1。本发明中采用双极性码,所以b′i=2bi-1,m′i=2mi-1∈{-1,1}。发射端采用的平移式特殊QPSK调制方案如图2所示,第k个特殊QPSK调制符号的同向分量为sI(k)=a×b′2k+m′2k×d,正交分量为sQ(k)=a×b′2k+1+m′2k+1×d,其中k=0,1,…,N/2-1。因此,发送的调制符号可以表示为
Figure BDA0002224573220000061
其中,
Figure BDA0002224573220000062
Figure BDA0002224573220000063
为由信息序列决定的调制相位,可表示为
Figure BDA0002224573220000064
θk是由同步序列决定的相位,可计算为
Figure BDA0002224573220000065
且θk∈{-θ0,0,θ0},其中θ0=arctan(a2/a1)-π/4。
步骤二:联合捕获方案
在接收端采用本发明提出的联合捕获方法。如图3所示,主要包含功率检测,频偏估计和相关检测三部分。
在接收端,假设接收信号经过下变频后存在Δf的频率偏移,则第k个接收符号可表示为
Figure BDA0002224573220000066
其中频率偏移引起的角频率为Δω=2πΔf;θk是由同步序列决定的相位;n0(k)为叠加在第k个接收符号上的复高斯白噪声。
1)功率检测
用功率检测作为第一级的粗捕获,可以有效的降低特殊QPSK的捕获复杂度。功率检测模块采用一个FIFO结构,取信号的绝对值代替信号的平方,检测一个帧长度内的信号的功率。功率检测模块使用一个定时计数变量n,在功率检测中用于判断是否进行频偏估计。当
Figure BDA0002224573220000067
或者功率检测捕获成功时,令n=0,输出频偏估计的启动信号,开始执行频偏估计算法;否则,令定时计数变量n值加1,继续执行功率检测。
在理想信道下,接收的信息序列为{s(0),s(1),s(2),…,s(k),…},与发送的信息序列一致。实际信道下接收的信息序列为{y(0),y(1),y(2),…,y(k),…},其中s(k)和y(k)表示在第k个采样时刻接收信号的采样值。根据理想信道下的一帧数据的接收功率估算为
Figure BDA0002224573220000068
对于正常接收,实际接收信息序列在FIFO中每隔一个采样时刻,按照先入先出的顺序逐一移位,进行一次功率检测。当处于检测时刻k0时,实际信道下的一帧数据的接收功率可以表示为
Figure BDA0002224573220000071
根据工程经验,将功率检测的阈值设为Pth=0.8Pp,当Pr(k′)>Pth时,可认为数据帧从第k′个采样开始,功率检测模块输出功率捕获成功信号。
2)频偏估计
在高信噪比下,忽略n0(k),对y(k)进行四次方后,
Figure BDA0002224573220000072
由于
Figure BDA0002224573220000073
所以化简可得
Figure BDA0002224573220000074
直接对y4(k)进行FFT可以估算出
Figure BDA0002224573220000075
可以表示为
Figure BDA0002224573220000076
其中fs是符号速率。
3)相关检测
设经过频偏校准后的信号为y′(k),y′I(k)和y′Q(k)分别为其实部和虚部。在接收端接收信号时,需要提取出隐藏在信号中的m序列,其第2k个和第2k+1个比特对应的符号可估计为
Figure BDA0002224573220000077
其中sign[·]为取符号函数,即
Figure BDA0002224573220000078
设实际信道中接收信号中的m序列与本地m序列完全对齐与未对齐时的归一化相关峰值分别为Rc和Re,这里的归一化定义为相关峰值与N的比值。将提取的m序列
Figure BDA0002224573220000079
与本地m序列{m′0,m′1,…,m′N-1}进行相关检测并判决,根据数据帧捕获的工程经验,当出现一个相关峰值大于0.6d且大于一帧数据长度内其他相关峰值的1.6倍时,即可判定为成功捕获。因此,本发明将数据帧捕获的判决门限设置为Rc≥max(0.6d,1.6max(Re)),当相关检测结果出现尖锐的相关峰且相关峰值大于判决门限时,接收端成功捕获到接收信号的帧起始位置。
步骤三:平移式特殊QPSK解调
接收端采用如图4所示的解调方式。与传统的QPSK解调相比,平移式特殊QPSK仅需在抽样判决的时候,加入m序列进行辅助判决即可。设经过频偏校正后的第k个接收符号的同相分量和正交分量分别为y′I(k)、y′Q(k),则第2k个信息比特可估计为
Figure BDA0002224573220000081
即当y′I(k)>m′2kd,
Figure BDA0002224573220000082
判决为1,否则,判决为0。同理,第2k+1个信息比特可估计为
Figure BDA0002224573220000083
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
在发射端,对信号进行平移式特殊QPSK调制以发送。经信道后,在接收端,对接收信号先进行联合捕获同步,然后再对平移式特殊QPSK进行解调。实际调制解调将在第一部分描述,联合捕获是本专利的核心,将在第二部分进行详细描述,第二部分具体阐述了一种平移式特殊QPSK系统的联合捕获方案,包括功率检测,频偏估计和相关检测三部分。第三部分总结了联合捕获方案的具体实施过程。
I.平移式特殊QPSK的调制解调模型
本发明研究的平移式特殊QPSK星座图如图1所示,a1、a2表示单路信号发生偏移后的幅度,幅度偏移量d=(a2-a1)/2,传统QPSK的星座点参数a=(a1+a2)/2。设信息序列为{b0,b1,b2,…,bN-1},m序列为{m0,m1,…,mN-1},其中bi,mi∈{0,1},i=0,1,…,N-1;{b′0,b′1,…,b′N-1}和{m′0,m′1,…,m′N-1}分别为信息序列和m序列的双极性形式,即b′i=2bi-1,m′i=2mi-1∈{-1,1}。
1)平移式特殊QPSK调制
在进行调制时,第k个特殊QPSK调制符号的同向分量为sI(k)=a×b′2k+m′2k×d,正交分量为sQ(k)=a×b′2k+1+m′2k+1×d,其中
Figure BDA0002224573220000084
因此,发送的调制符号可以表示为
Figure BDA0002224573220000085
其中
Figure BDA0002224573220000086
Figure BDA0002224573220000087
Figure BDA0002224573220000088
为由信息序列决定的调制相位且
Figure BDA0002224573220000089
θk为由同步序列决定的相位,
Figure BDA00022245732200000810
且θk∈{-θ0,0,θ0},其中θ0=arctan(a2/a1)-π/4。
可以看到,Ak和θk的变化均由m序列和信息序列共同作用,且满足
Figure BDA0002224573220000091
其中
Figure BDA0002224573220000092
表示事件发生的概率。
2)平移式特殊QPSK解调
在经过接收端联合捕获之后,第k个接收符号的同相分量和正交分量可以分别表示为
Figure BDA0002224573220000093
其中Δf和
Figure BDA0002224573220000094
分别表示实际的多普勒频移及其估计值,nI(k)和nQ(k)是第k个接收符号叠加的高斯白噪声且服从于均值为0、方差为σ2/2的正态分布。平移式特殊QPSK信号的解调原理框图如图4所示,第2k个信息比特可估计为
Figure BDA0002224573220000095
即当y′I(k)>m′2kd,
Figure BDA0002224573220000096
判决为1,否则,判决为0。同理,第2k+1个信息比特可估计为
Figure BDA0002224573220000097
II.联合捕获方案
在实际的通信系统中,频率的偏移总是存在的。由于特殊QPSK的特殊性,当接收信号存在频率偏移时,隐藏于符号中的m序列将会受到严重的影响,较小频偏导致的符号旋转都将会对m序列提取的造成很大干扰,频偏稍大时完全无法进行正常捕获。为克服这一问题,本发明提出一种平移式特殊QPSK系统联合捕获方案,主要包含功率检测、频偏估计和相关检测三部分,具体流程如图3所示。功率检测中,我们采用绝对值的方式去表示信号功率,当功率检测器中达到判决门限时,输出一个功率捕获的信号,启动频偏估计;对输入信号进行四次方,去除符号的调制信息,然后采用基于FFT的频偏估计算法;然后,利用估计的频率值,对输入数据进行频偏补偿并输出到相关检测中;最后,在相关检测中提取出m序列并将其与本地m序列进行相关运算,对相关结果进行判决,完成帧同步的功能。
1)功率检测
功率检测模块主要采用一个FIFO结构,取信号的绝对值代替信号的平方,检测一个帧长度内信号的功率。令Pp为理想信道下一帧接收数据的功率,Pr(k0)为实际信道中k0时刻一帧接收数据的功率。定义n为一个定时计数变量,在功率检测中用于辅助进行频偏估计,防止由于功率检测的漏捕获而导致数据帧丢失的现象。当功率检测捕获成功或者
Figure BDA0002224573220000098
时,令n=0,输出频偏估计的启动信号,开始执行频偏估计算法;否则,定时计数变量n=n+1,继续执行功率检测。
在理想信道下,接收的信息序列为{s(0),s(1),s(2),…,s(k),…},与发送的信息序列一致,实际信道下接收的信息序列为{y(0),y(1),y(2),…,y(k),…},其中s(k)和y(k)表示在第k个采样时刻接收信号的采样值。根据理想信道下的一帧数据的接收功率估算为
Figure BDA0002224573220000101
对于正常接收,实际接收信息序列在FIFO中每隔一个采样时刻,按照先入先出的顺序逐一移位,进行一次功率检测。当处于检测时刻k0时,实际信道下的一帧数据的接收功率可以表示为
Figure BDA0002224573220000102
根据工程经验,将功率检测的阈值设为Pth=0.8Pp,当Pr(k′)>Pth时,可认为数据帧从第k′个采样开始,功率检测模块输出功率捕获成功信号。
功率检测具有低时间复杂度、低实现复杂度等优点,是平移式特殊QPSK信号捕获的第一步,可以为有效的降低信号捕获过程中的复杂度。
2)频偏估计
为了找到合适的频偏估计算法,本发明分析了特殊QPSK信号的特性。假设信道为理想信道,接收端接收到的采样信号
Figure BDA0002224573220000103
可以表示为
Figure BDA0002224573220000104
其中
Figure BDA0002224573220000105
θk为由同步序列决定的相位,且θk∈{-θ0,0,θ0},其中θ0=arctan(a2/a1)-π/4。对
Figure BDA0002224573220000106
进行四次方可得
Figure BDA0002224573220000107
由于
Figure BDA0002224573220000108
所以有
Figure BDA0002224573220000109
可以发现,通过对
Figure BDA00022245732200001010
进行四次方,得到的
Figure BDA00022245732200001011
中不再含有调制信息对应的相位,实现了去除调制符号的功能。
针对
Figure BDA00022245732200001012
信号的频谱进行分析,单个
Figure BDA00022245732200001013
符号中的平均能量为
Figure BDA00022245732200001014
其直流分量的能量为
Figure BDA0002224573220000111
所有交流信号的总能量为
Figure BDA0002224573220000112
设交流信号能量与直流信号能量的比例为η,则
Figure BDA0002224573220000113
根据在先专利,
Figure BDA0002224573220000114
可以得到η≈0.69。由此可见,
Figure BDA0002224573220000115
信号中的直流能量大于交流能量。此外,交流信号是由相位和幅度的随机性偏移的产生的,没有固定频率,所以,其频谱将会分散到整个带宽中。因此,在
Figure BDA0002224573220000116
频谱图上,将只会在零频处出现一个极大的峰值,直流分量的幅度大于其他交流信号的幅度值。
假设信道是高斯白噪声信道,并且接收信号具有Δf的频率偏移,则第k个接收符号可以表示为
Figure BDA0002224573220000117
其中频率偏移引起的角频率Δω=2πΔf,n0(k)为叠加在第k个接收符号上的复高斯白噪声。在高信噪比下,忽略n0(k),对y(k)进行四次方后化简可得
Figure BDA0002224573220000118
采用最大似然估计方法,建立Δω的最大似然函数如下
Figure BDA0002224573220000119
其中
Figure BDA00022245732200001110
yI(k)和yQ(k)分别为y4(k)的实部和虚部。将yI(k)、yQ(k)、YI(k)和YQ(k)带入化简可得
Figure BDA00022245732200001111
对于已接收数据来说,
Figure BDA0002224573220000121
为定值,
Figure BDA0002224573220000122
恒成立,所以最大似然函数可以简化为
Figure BDA0002224573220000123
因为
Figure BDA0002224573220000124
Figure BDA0002224573220000125
Figure BDA0002224573220000126
具有相同的单调性,式(12)可以简化为为
Figure BDA00022245732200001215
其中Re[x]代表x的实部。因此,Δω的最大似然估计为
Figure BDA0002224573220000127
其中
Figure BDA0002224573220000128
从前面分析的结果可以得出,当复信号y4(k)中相位和幅度的发生随机性偏移时,不会影响频谱中的最大值。因此,式(14)可以修改为
Figure BDA0002224573220000129
基于离散傅立叶变换(DFT)特性,可以发现,
Figure BDA00022245732200001210
对应于y4(k)的DFT运算的最大值。由于Δω∈[-π/4,π/4],
Figure BDA00022245732200001211
可以重新表示为
Figure BDA00022245732200001212
Δω接近于Cramer-Rao边界。
最后,我们得到
Figure BDA00022245732200001213
的估计值
Figure BDA00022245732200001214
其中fs是符号速率。
3)相关检测
特殊QPSK调制将同步码序列隐藏在每一个符号的能量中。在接收端接收信号时,需要提取出隐藏在信号中的m序列,将提取的m序列
Figure BDA0002224573220000131
与本地m序列{m′0,m′1,…,m′N-1}进行相关检测并判决,从而完成特殊QPSK调制解调系统的帧同步功能。本发明在m序列末尾补0,使得m序列中0和1的个数相等,即长度N为偶数。因此,理想信道下,应用于平移式特殊QPSK系统的m序列的归一化自相关函数可定义为
Figure BDA0002224573220000132
本发明设计的m序列提取和相关捕获的原理框图如图5所示。接收端m序列的提取是发送端m序列调制的逆过程,接收信号中隐藏的信号可以用来表征m序列。所以,m序列中的第2k个和第2k+1个比特对应的符号可估计为
Figure BDA0002224573220000133
设实际信道中接收信号中的m序列与本地m序列完全对齐与未对齐时的归一化相关峰值分别为Rc和Re。根据数据帧捕获的工程经验,当出现一个相关峰值大于0.6d且大于一帧数据长度内其他相关峰值的1.6倍时,即可判定为成功捕获数据帧的起始位置。因此,本发明将数据帧捕获的判决门限设置为Rc≥max(0.6d,1.6max(Re))。将接收端提取的m序列
Figure BDA0002224573220000134
与本地m序列{m′0,m′1,…,m′N-1}做相关并进行判决,当相关结果出现尖锐的相关峰且相关峰值Rc大于等于判决门限时,接收端成功捕获到接收信号的帧起始位置。
III.联合捕获算法
综上所述,算法1总结了联合捕获方案的具体实施过程,其中步骤1-2是功率检测,4-5是频偏补偿,6-7是相关检测。此外,我们可以看到频偏估计独立于时间轴上的其他模块。
算法1参数说明:
Figure BDA0002224573220000135
是频偏估计的值,
Figure BDA0002224573220000136
是上一次的频偏估计值。
Figure BDA0002224573220000137
和y′(k)是第k个接收符号及其经过频偏补偿后的符号。N为数据帧长度且为偶数。
算法1流程:
(1)由功率检测模块计算出Pr(k′)。
(2)如果Pr(k′)≥Pth或者
Figure BDA0002224573220000141
则执行步骤3;否则n=n+1,执行步骤5。
(3)令n=0,并根据式(17)计算频偏估计值
Figure BDA0002224573220000142
(4)如果
Figure BDA0002224573220000143
不等于
Figure BDA0002224573220000144
Figure BDA0002224573220000145
(5)令y′(k)=y(k)exp(-j2πΔf′);
(6)计算相关峰值Rc
Figure BDA0002224573220000146
(7)当Rc≥max(0.6d,1.6max(Re)),则说明捕获成功,跳出循环;否则继续循环步骤1-7;
在实施例中,为了验证本发明所提方法的可行性,对所设计的调制解调系统进行了性能仿真。后续仿真中,若为重复实验,均为100,000次独立实验的结果。
给定数据帧长度N=1024,帧捕获的噪声环境为SNR=4dB,本发明将数据帧捕获的判决门限设置为Rc≥max(1.6max(Re),0.6d),在此基础上通过MATLAB进行仿真。
将多普勒频移对信息速率进行归一化。当归一化多普勒频率偏移为0.1时,平移式特殊QPSK在SNR=4dB下的捕获效果如图6所示。从图中可以看出,基于上述联合捕获方法,接收信号中提取出的m序列与本地m序列具有非常好的相关特性。当提取的信号与本地的m序列完全对齐时,其归一化相关峰值会远远大于其他峰值。
图7展示了基于Nf点FFT算法的频偏估计性能,其中Nf=2048。FFT的最小归一化频率间隔是
Figure BDA0002224573220000147
理论上的最大归一化频偏估计误差是2.44×10e-4。从图中可以看出,平移式特殊QPSK在SNR=4dB时的最大归一化频率估计误差为1.6×10e-4,小于理论推导的2.44×10e-4,具有较高的频偏估计精度。
图8对平移式特殊QPSK、对称式特殊QPSK和QPSK调制解调系统在不同SNR下的误比特性能进行了仿真。在加性高斯白噪声信道下,对称式特殊QPSK性能较差,而平移式特殊QPSK性能与传统QPSK的性能非常接近,仅有微小的BER性能损失。从图中可以看出,即便是在信噪比为10dB的时候,平移式特殊QPSK的BER性能损失也非常的微小。
图9对平移式特殊QPSK调制解调系统在不同SNR下的联合捕获性能进行了仿真。当归一化多普勒频率偏移为0.1时,从图中可以看出,当SNR=1dB时,捕获概率也能有50%左右,而随着SNR的增加,捕获概率急剧增大,当SNR≥4dB时,捕获概率达到了100%,说明了本发明提出的联合捕获方案具有良好的捕获性能和可行性。
结论:本发明以平移式特殊QPSK系统为例,给出了平移式特殊QPSK系统的调制、解调以及联合捕获的方法,考虑存在多普勒频移的情况下,提出了一种将功率检测,频偏估计和相关检测相结合的联合捕获方案。给定数据帧长度为512和具体捕获率要求,通过MATLAB仿真验证了所提方法的可行性。本发明所提出的方法,具有良好的捕获性能和误比特性能,可在实际工程中应用。本发明中的方案设计与优化方法不失一般性,同样适用于对称式特殊QPSK中。

Claims (7)

1.一种平移式特殊QPSK系统联合捕获方法,其特征在于包括以下步骤:步骤一,在发送端,采用帧同步序列对信号的幅度进行调制,将帧同步序列隐藏在信息序列中;步骤二,在接收端,首先对捕获信号进行功率检测,当进行了设定次数的功率检测或者功率检测捕获成功后进行频偏估计,对经过频偏校准后的信号作相关检测,当相关检测结果满足设定条件时,接收端成功捕获到接收信号的帧起始位置;步骤三,对接收信号进行解调,解调过程中,在抽样判决的时候加入m序列进行辅助判决。
2.根据权利要求1所述的平移式特殊QPSK系统联合捕获方法,其特征在于:所述的步骤一中,设信息序列为{b0,b1,b2,…,bN-1},帧同步序列m为{m0,m1,…,mN-1},其中,N为数据帧长度且为偶数,bi,mi∈{0,1},i=0,1,…,N-1,b′i=2bi-1,m′i=2mi-1∈{-1,1};发射端采用平移式特殊QPSK调制,第k个特殊QPSK调制符号的同向分量为sI(k)=a×b′2k+m′2k×d,正交分量为sQ(k)=a×b′2k+1+m′2k+1×d,其中k=0,1,…,N/2-1,d代表幅度偏移量,a代表QPSK的星座点参数;发送的调制符号
Figure FDA0003067062410000011
其中,
Figure FDA0003067062410000012
Figure FDA0003067062410000013
a1、a2表示单路信号发生偏移后的幅度;
Figure FDA0003067062410000014
为由信息序列决定的调制相位,
Figure FDA0003067062410000015
θk是由同步序列决定的相位,
Figure FDA0003067062410000016
且θk∈{-θ0,0,θ0},其中θ0=arctan(a2/a1)-π/4。
3.根据权利要求2所述的平移式特殊QPSK系统联合捕获方法,其特征在于:所述的步骤二中,在接收端,假设接收信号经过下变频后存在Δf的频率偏移,则第k个接收符号
Figure FDA0003067062410000017
其中频率偏移引起的角频率为Δω=2πΔf;θk是由同步序列决定的相位;n0(k)为叠加在第k个接收符号上的复高斯白噪声;功率检测模块采用一个FIFO结构,取信号的绝对值代替信号的平方,检测一个帧长度内的信号的功率;功率检测模块使用一个定时计数变量n,在功率检测中用于判断是否进行频偏估计;当
Figure FDA0003067062410000021
或者功率检测捕获成功时,令n=0,输出频偏估计的启动信号,开始执行频偏估计算法;否则,令定时计数变量n值加1,继续执行功率检测;频偏
Figure FDA0003067062410000022
其中fs是符号速率;设经过频偏校准后的信号为y′(k),y′I(k)和y′Q(k)分别为其实部和虚部;在接收端接收信号时,提取出隐藏在信号中的m序列,其第2k个和第2k+1个比特对应的符号分别估计为
Figure FDA0003067062410000023
当相关检测结果出现相关峰且相关峰值大于判决门限时,接收端成功捕获到接收信号的帧起始位置。
4.根据权利要求3所述的平移式特殊QPSK系统联合捕获方法,其特征在于:所述的步骤二中,在理想信道下接收的信息序列为{s(0),s(1),s(2),…,s(k),…},实际信道下接收的信息序列为{y(0),y(1),y(2),…,y(k),…};根据理想信道下的一帧数据的接收功率估算为
Figure FDA0003067062410000024
对于正常接收,实际接收信息序列在FIFO中每隔一个采样时刻,按照先入先出的顺序逐一移位,进行一次功率检测;当处于检测时刻k0时,实际信道下的一帧数据的接收功率
Figure FDA0003067062410000025
根据功率检测的阈值Pth,当Pr(k′)>Pth时,认为数据帧从第k′个采样开始,功率检测模块输出功率捕获成功信号。
5.根据权利要求4所述的平移式特殊QPSK系统联合捕获方法,其特征在于:所述的步骤二中,功率检测的阈值设为Pth=0.8Pp
6.根据权利要求3所述的平移式特殊QPSK系统联合捕获方法,其特征在于:所述的步骤二中,设实际信道中接收信号中的m序列与本地m序列完全对齐与未对齐时的归一化相关峰值分别为Rc和Re,将提取的m序列
Figure FDA0003067062410000026
与本地m序列{m′0,m′1,…,m′N-1}进行相关检测并判决,数据帧捕获的判决门限设置为Rc≥max(0.6d,1.6max(Re))。
7.根据权利要求2所述的平移式特殊QPSK系统联合捕获方法,其特征在于:所述的步骤三中,设经过频偏校正后的第k个接收符号的同相分量和正交分量分别为y′I(k)、y′Q(k),则第2k个信息比特估计为
Figure FDA0003067062410000031
即当y′I(k)>m′2kd,
Figure FDA0003067062410000032
判决为1,否则,判决为0;同理,第2k+1个信息比特估计为
Figure FDA0003067062410000033
CN201910947509.XA 2019-10-08 2019-10-08 一种平移式特殊qpsk系统联合捕获方法 Active CN110730149B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910947509.XA CN110730149B (zh) 2019-10-08 2019-10-08 一种平移式特殊qpsk系统联合捕获方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910947509.XA CN110730149B (zh) 2019-10-08 2019-10-08 一种平移式特殊qpsk系统联合捕获方法

Publications (2)

Publication Number Publication Date
CN110730149A CN110730149A (zh) 2020-01-24
CN110730149B true CN110730149B (zh) 2021-07-09

Family

ID=69219691

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910947509.XA Active CN110730149B (zh) 2019-10-08 2019-10-08 一种平移式特殊qpsk系统联合捕获方法

Country Status (1)

Country Link
CN (1) CN110730149B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112019233B (zh) * 2020-08-20 2022-05-24 西安烽火电子科技有限责任公司 一种短波通信的多路信号同步头快速捕获方法
CN115037327B (zh) * 2022-05-31 2023-12-22 江苏屹信航天科技有限公司 初始捕获跟踪的接收方法、接收装置、终端

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103812813A (zh) * 2012-11-12 2014-05-21 德尔福电子(苏州)有限公司 基于软件定义无线电cmmb解调器的时域同步方法
US9999009B1 (en) * 2009-07-21 2018-06-12 Marvell International Ltd. Simultaneous uplink transmission in a wireless network
CN108512791A (zh) * 2018-03-29 2018-09-07 中国人民解放军国防科技大学 基于定时频偏补偿的星载ais解调方法
CN110099024A (zh) * 2019-03-08 2019-08-06 西北工业大学 一种帧同步码隐含于信息码中传输的方法
CN110278168A (zh) * 2018-03-13 2019-09-24 意法半导体国际有限公司 用于通信接收器的载频偏移补偿电路和处理

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5220584A (en) * 1990-12-21 1993-06-15 Mikros Systems Corp. System for demodulation and synchronizing multiple tone waveforms
CN101291311B (zh) * 2007-04-18 2013-01-16 中兴通讯股份有限公司 多输入多输出正交频分复用系统的同步实现方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9999009B1 (en) * 2009-07-21 2018-06-12 Marvell International Ltd. Simultaneous uplink transmission in a wireless network
CN103812813A (zh) * 2012-11-12 2014-05-21 德尔福电子(苏州)有限公司 基于软件定义无线电cmmb解调器的时域同步方法
CN110278168A (zh) * 2018-03-13 2019-09-24 意法半导体国际有限公司 用于通信接收器的载频偏移补偿电路和处理
CN108512791A (zh) * 2018-03-29 2018-09-07 中国人民解放军国防科技大学 基于定时频偏补偿的星载ais解调方法
CN110099024A (zh) * 2019-03-08 2019-08-06 西北工业大学 一种帧同步码隐含于信息码中传输的方法

Also Published As

Publication number Publication date
CN110730149A (zh) 2020-01-24

Similar Documents

Publication Publication Date Title
US11108612B2 (en) Anti-interference signal detection and synchronization method for wireless broadband communication system
US9954701B2 (en) Bluetooth signal receiving method and device using improved packet detection and symbol timing acquisition
CN104852876B (zh) 一种航空无线突发通信系统
CN109379314B (zh) 高速突发数字解调方法和设备
CN110113276B (zh) 基于ieee802.11的ofdm频偏估计方法、系统及装置
CN110099024B (zh) 一种帧同步码隐含于信息码中传输的方法
WO2005074223A1 (en) Frame synchronization and initial symbol timing acquisition system and method
CN111884685A (zh) 数字通信信号同步解调方法及其装置
CN103023831B (zh) 一种适用于突发波形的载波频偏估计方法
CN110730149B (zh) 一种平移式特殊qpsk系统联合捕获方法
CN107257324A (zh) 一种ofdm系统中的时频联合同步方法及装置
CN105871765A (zh) 一种基于fft辅助s-pll的无线通信载波跟踪方法
CN112399551B (zh) 一种面向短时突发信号的高精度同步方法
CN107342960B (zh) 一种适合幅度相移键控的非数据辅助频偏估计方法
CN107682294B (zh) 一种基于FPGA的高速16apsk信号的相位模糊校正方法
CN111147406B (zh) 基于恒虚警信号检测及频偏校正的asm体制解调方法
CN105610755B (zh) 一种突发信号的频偏估计方法及装置
CN114205200A (zh) 一种实现vdes系统帧头捕获和载波同步的方法
CN107318157B (zh) 基于前导序列的定时误差估计方法
CN110266625B (zh) 一种基于帧同步码隐含于信息码的qpsk调制的动态门限判决的解调方法及系统
JP5624527B2 (ja) シングルキャリア受信装置
CN109889461B (zh) 一种低复杂度并行的载波恢复系统及其方法
CN108650203B (zh) 一种基于侦察接收机的调制模式识别方法
CN113115430A (zh) 一种高速突发数字解调系统
CN114338334B (zh) 一种基于伪码导频的相位噪声估计与补偿方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant