CN107835035B - 低信噪比短帧突发通信开环解调方法及装置 - Google Patents
低信噪比短帧突发通信开环解调方法及装置 Download PDFInfo
- Publication number
- CN107835035B CN107835035B CN201711041223.2A CN201711041223A CN107835035B CN 107835035 B CN107835035 B CN 107835035B CN 201711041223 A CN201711041223 A CN 201711041223A CN 107835035 B CN107835035 B CN 107835035B
- Authority
- CN
- China
- Prior art keywords
- data
- frequency offset
- phase
- section
- estimation
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004891 communication Methods 0.000 title claims abstract description 40
- 230000008859 change Effects 0.000 claims abstract description 78
- 238000012545 processing Methods 0.000 claims abstract description 20
- 230000001427 coherent effect Effects 0.000 claims abstract description 18
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims description 8
- 230000000875 corresponding effect Effects 0.000 claims description 5
- 101100313471 Streptomyces sp getA gene Proteins 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 3
- 239000010931 gold Substances 0.000 claims description 3
- 229910052737 gold Inorganic materials 0.000 claims description 3
- 238000007476 Maximum Likelihood Methods 0.000 claims description 2
- 230000002596 correlated effect Effects 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 230000009191 jumping Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 238000012937 correction Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000007480 spreading Effects 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000007493 shaping process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details 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/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
- H04L27/22—Demodulator circuits; Receiver circuits
- H04L27/227—Demodulator circuits; Receiver circuits using coherent demodulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明低信噪比短帧突发通信开环解调方法及装置,属于直接序列扩频通信解调技术领域。本发明采用分段数据处理的方式,先对解扩数据解线调,之后做补零FFT,估计出频偏变化率和频偏并进行补偿。利用穷举法对每段频偏进行精确估计并且补偿,利用帧头极性对相位进行解模糊,完成数据的正确解调。对比现有技术,本发明克服传统锁相环相干解调装置具有相位模糊以及收敛速度慢等不足,能够实现对频偏和频偏变化率快速搜索并补偿,完成相位鉴别,相位纠正,保证短时间,低信噪比条件下,直接序列扩频信号能够正常解调。
Description
技术领域
本发明涉及一种载波相位跟踪与相干解调方法(下文简称解调装置),尤其涉及低信噪比短帧突发通信开环解调方法及装置,适用于高动态,低信噪比等环境下的短帧突发直接序列扩频***载波频偏、频偏变化率的开环搜索、跟踪和符号相干解调,属于直接序列扩频通信中的解调技术领域。
背景技术
高速飞行器间与地面设备或高速飞行器之间的通信具有通信距离远、链路衰减大、信噪比低的特点,因此常采用直接序列扩频通信体制;且该种通信往往要求保证通信隐蔽性、具备防侦测能力,因而采用短帧突发体制。由于飞行器的运动速度快,运动特性复杂,通信链路上的发射机和接收端之间的相对径向运动存在一阶甚至高阶加速度,一方面由于是短帧突发的通信体制,帧通信的时间很短,所以高阶运动加速度可忽略,但另一方面接收机接收信号存在很大的多普勒频偏和多普勒一阶频偏变化率(下文简称频偏变化率)不可忽略;由于接收信号的信噪比很低,所以采用相干解调的方式,如果采用锁相环的方式去实现对载波多普勒频偏和频偏变化率的跟踪,在短帧突发通信体制下,会出现环路无法及时收敛,导致数据不能正常解调,因此应采用搜索速度更快,时间更短的开环搜索方式,完成相干解调。
通常的直接序列扩频接收机主要包括捕获模块、解扩模块、解调模块,译码模块等。捕获模块的作用是获得接收信号的频偏和相位信息,送入解扩模块;解扩模块利用捕获信息,恢复扩频前的数据;解调模块利用解扩输出数据,恢复调制前的基带信息;译码模块是利用解调模块的输出完成译码。
假设直接序列扩频通信接收机扩频码为c[l],l=0,1,2,…,L-1,L为扩频码长度。c(l)取值为0或者1,采用无线电技术惯例,c[l]=0表示逻辑正,c[l]=1表示逻辑负,并且采用根升余弦脉冲成形,脉冲成形滤波器冲击函数为h(t),则成形之后的基带波形表达式为:
其中,b[m]为基带数据符号,其周期Ts=LTc,Tc为码片周期,短帧突发通信体制中,一般对突发帧帧结构有一定要求;为了便于捕获,帧通信开始先要发射导频信号,即b[m]恒为0。导频信号结束后,b[m]为长度固定的M序列或者是Gold码的帧头,帧头结束后是数据段。经过调制后,射频所发出的信号为:
Ssrnd(t)=cos(2πfRFt)Sb,s(t)
上式中,fRF为射频频率,并且为了表达方便,假设发射信号载波初相位和码初相位均为0。
经过无线传输过程的多普勒效应和信号延迟,到达接收机,对射频信号进行I/Q两路正交下变频,得到复基带信号为:
上式中fβ和fσ分别为,通信链路双方相对径向加速运动产生的多普勒频偏变化率和多普勒频偏;τ0为码相位延迟。通过接收机的捕获模块完成多普勒频偏和伪码相位二维搜索。其中,一旦考虑频偏变化率,“频偏”和“初始频偏”就是两个不同的概念。默认情况下“频偏”指的是瞬时频偏;但是对于存在频偏变化率的场合,瞬时频偏是时变函数,只有“初始频偏”(帧突发通信开始时的瞬时频偏)才是恒定不变的。捕获装置得到频偏是瞬时频偏。
利用扩频码周期性特点对输入数据做相关操作,即解扩。只有当rde_In[n]和y[n]为同一种序列并且他们的相位又相互一致的时候,此时相关值最大,此时d[n]为解扩输出结果,按照符号速率更新。考虑到频偏变化率和残余频偏的影响,所以解扩数据输出可以写成:
上式中,β是对符号速率归一化的数字化频偏变化率;ωσ为捕获模块得到频偏粗估计相对于码片速率归一化数字频率,Rc为码片速率,Rc=Rs·L。代表残余载波初始(数字)频偏,假定初始相位为bn∈0/1代表发送的二进制符号(BPSK调制);解调模块的作用是完成频偏变化率,频偏的搜索、补偿以及相位鉴别,解调后符号记为z(n),其表达式为:
令
传统BPSK相干解调模块通常是利用锁相环,比如平方环、科斯塔斯(Costas)环等,实现载波跟踪的,但都不能有效解决“倒π”问题,即恢复的本地载波与所需的相干载波可能同相,也可能反相,这种相位不确定性将会造成解调出来的数字基带信号正好相反,判决出来的数字符号全部出错。
发明内容
本装置的目的是为了克服传统锁相环相干解调装置具有相位模糊以及收敛速度慢为主的技术缺陷,提出了低信噪比短帧突发通信开环解调方法及装置,实现对频偏和频偏变化率搜索并补偿,完成相位鉴别,相位纠正,保证短时间,低信噪比条件下,直接序列扩频信号能够被正常解调。
低信噪比短帧突发通信开环解调方法及装置包括低信噪比短帧突发通信开环解调装置,简称本装置;以及低信噪比短帧突发通信开环解调方法,简称本方法。
本装置包括二维搜索模块、精细搜索模块1、精细搜索模块2、精细搜索模块3和帧同步模块;
其中,精细搜索模块1、精细搜索模块2与精细搜索模块3是复用的,本装置中复用的3个精细搜索模块完全一样;
本装置的输入数据为解扩模块的输出,本装置的输出数据为解调结果,送入译码模块;
本装置各模块的连接关系如下:
二维搜索模块连接精细搜索模块1;精细搜索模块1与精细搜索模块2相连,精细搜索模块2与精细搜索模块3相连,精细搜索模块3和帧同步模块相连;
本装置中各模块的功能如下:
二维搜索模块的功能是进行二维搜索,完成对接收信号初始频偏和频偏变化率进行联合估计,得到初始频偏和频偏变化率的估计,并且对多普勒频偏变化率和多普勒频偏进行补偿;精细搜索模块的功能是完成频偏的估计并补偿、相位估计以及相位补偿;帧同步模块的功能是解相位模糊。
本方法是通过以下技术方案实现的:
本装置对输入数据采用分段数据处理的方法,包括如下步骤:
步骤一、基于本装置的输入数据运行多普勒频偏变化率和多普勒频偏变化率二维搜索模块,完成对接收信号初始频偏和频偏变化率进行联合估计,得到初始频偏和频偏变化率的估计;
其中,本装置的输入数据表示为如下公式(1),也称为第一段数据;
能够通过二维网格搜索进行近似求解,将频偏和频偏变化率全部不确定范围按照一定精度划分成二维平面网格,得到二维平面的方式是:对于频偏变化率,在其变化范围内,取Q个即按照对x[n]进行Q次“解线性调频调制”,简称“解线调”,对应的符号xq[n]表达为如下公式(3):
对xq[n]进行补零至长度为K点,再进行K点FFT运算,得到xq[k];对xq[k]取模,结果为|xq[k]|;再从全部Q×K个取模结果|xq[k]|找到最大值,根据这个最大值在二维网格上的位置,得到β和Δω估计量和其中,表示步骤一对频偏变化率的估计结果,表示步骤一对初始频偏的估计结果;
完成对接收信号初始频偏和频偏变化率进行联合估计,得到初始频偏和频偏变化率的估计,从第二段数据开始,频偏变化率的补偿均使用第一段数据得到的频偏变化率的估计结果,并且默认频偏变化率已经补偿完毕。即使步骤一得出的频偏变化率估计结果存在误差,即存在残余的多普勒频偏变化率,只要第二段数据处理长度和Q取值合适,完全可以忽略残余的多普勒频偏变化率对频偏估计造成的影响;
各段数据只需要对频偏进行估计并补偿即可,完成数据的解调;
步骤二、基于第一段数据得到的初始频偏的估计结果,运行精细搜索模块,基于第二段数据得到的频偏的精确估计结果并进行补偿,从而完成数据解调;
步骤二、具体为:
步骤2.1基于步骤一中对初始频偏的估计值进行补偿,得出公式(4)的第二段数据:
步骤2.2运行精细搜索模块,得到频偏的精确估计;
运行本装置精细搜索模块,精细搜索模块采用穷举法得到频偏的精确估计,进而得出第二段数据的相干解调结果,具体为:
步骤2.21先对第二段数据(4)进行平方去调制,再与本地内积模板的平方做内积,本地内积模板表达式如公式(5):
与本地内积模板的平方做内积的具体方法是将Δω2等间隔分成P个小区间,用P个数字频率不同的复正弦波的平方与x2[n]内积,再取内积模值最大者为复正弦波数字频率为-2Δω2的一个估计,基于此估计得到用于补偿频偏;相位按照下式(6)估计:
上式(6)中的“除以2”操作会导致相位模糊,即和都满足公式(6)的要求,由于初始值未知,在这里取或者都可以;当数据处理到帧头段时,利用帧头的相关结果极性,来确定绝对相位,不妨先取得到第二段数据相干解调结果为下式(7):
步骤三、对本段数据频偏进行精确估计并且补偿,对本段数据相位进行解模糊,即需要与前一段数据相位保持一致;
步骤三、具体为:
第三段数据为:
本段数据为:
其中,第一次运行步骤三时,本段数据指第三段数据,其后依次类推;
步骤3.1:运行精细搜索模块,得到频偏的精确估计;精细搜索模块采用穷举法得到频偏的精确估计,进而得出本段数据的相干解调结果,具体为:
先对本段数据进行平方去调制,再与本地内积模板的平方做内积,本地内积模板表达式如公式(5);与本地内积模板的平方做内积的具体方法是将Δω2等间隔分成P个小区间,用P个数字频率不同的复正弦波的平方与x2[n]内积,再取内积模值最大者为复正弦波数字频率为-2Δω3的一个估计,基于此估计得到用于补偿频偏;相位按照下式(9)估计:
步骤3.2:对步骤3.1输出的相位进行解模糊,解模糊的参考相位根据前一段频偏和初相位的结果来构造,因为前一段的末相位就是本段的初相位;
这因为末相位和初相位之间隔了N-1个符号;如果考虑扩频***的解扩过程,则第三段的初相位和第二段的末相位还稍有区别,应作(10):
由于第二段本身就存在相位模糊,所以第三段和以后各段只需要严格跟前一段相位保持一段相位连续即可;具体的处理过程如下(11):
步骤四、基于后续数据段,重复运行步骤三,直到导频段数据处理结束;所得的相位估计结果能和前一段数据相位保持一致;
步骤五、步骤一到步骤四将导频段数据处理结束后,进入帧头段的数据处理,具体的处理方式是是对帧头数据取相关,得出帧头相关结果;
其中,帧头通常为M序列或者Gold码序列,帧头是数据段的开始;
帧头通过解扩、频率估计、补偿和相位估计,输出数据z[n],与本地的帧头模板做如下公式(14)的相关运算:
此公式(14)的相关过程与解扩过程的运算一样,区别是解扩中是在码片级别上做的相关操作,解扩是相关长度为一个符号的时间长度;而帧同步是在符号级别上做的相关运算,相关运算的长度取决于帧头的长度;rheader[n]为帧头相关结果;sheader_tep[n]为帧头本地模板;M为帧头长度;
步骤六、将帧头相关结果与本地模板做相关,依据相关结果绝对值是否超过预定门限来确定帧同步是否成功,并进行相应操作,具体为:对帧头相关运算结果rheader[n]取模操作,如果取模结果超过门限,则认为帧同步成功,跳至步骤七,启动译码装置;否则认为本帧数据接收失败,返回至步骤一等待下一帧数据到来;
步骤七、帧同步成功时,根据相关结果极性,对相位进行解模糊;具体为:根据为解决“倒π”问题,当超过门限的时刻,需要确认超过门限时相关结果的正负;如果极性为正,说明此时保持的相位为0,直接输出z[n]进译码装置;如果相关结果极性为负,说明此时保持的相位为π,需要将z[n]取反,之后给译码模块;译码结束后,本帧数据接收成功,本装置返回至步骤一等待下一帧数据到来;
至此,从步骤一到步骤七,完成了低信噪比短帧突发通信开环解调方法。
有益效果
本发明提出的低信噪比短帧突发通信开环解调方法及装置,与现有技术相比,具有如下有益效果:
1.利用解线调和二维搜索等方法能够在很短时间内完成对多普勒频偏和频偏变化率的精确估计;
2.通过分段处理方式,利用穷举法将每段的频偏的精细估计并进行补偿,实现对载波频率跟踪,完成相位鉴别;
3、利用帧头极性对相位进行纠正,保证数据信息能够正确解调,适用于高动态,低信噪比等环境下的短帧突发直接序列扩频***载波频偏、频偏变化率的开环搜索、跟踪和符号相干解调。
附图说明
图1是本发明低信噪比短帧突发通信开环解调方法及装置及实施例1中的本装置的原理框图;
图2是本发明低信噪比短帧突发通信开环解调方法及装置中实施例步骤一中频偏和频偏变化率二维搜索平面;
图3是本发明低信噪比短帧突发通信开环解调方法及装置中实施例步骤一中解线调原理;
图4是本发明低信噪比短帧突发通信开环解调方法及装置中穷举法频率精细搜索原理框图。
具体实施方式
下面结合实例和附图对本装置做进一步详细说明。
本实例以直接序列扩频BPSK通信接收机为例,其基带数字信号处理部分如图1所示。
从图1可以看出进入本装置的是解扩数据,依次运行本装置的多普勒频偏变化率和多普勒频偏二维搜索模块、精细搜索模块以及帧同步模块,完成解调功能,将数据输出至接收机的译码模块,完成译码。
将解扩后的数据输入至本装置,数据频偏和频偏变化率的先验信息为Δf=-130.9140Hz;ΔΔf=-3253.3Hz/s。第一段数据长度为128个符号,且全部为导频(基带数据全“0”)。进入多普勒频偏变化率和多普勒频偏二维搜索模块后,将初始的频偏和频偏变化率的全部不确定范围各自按照一定精度划分为二维平面网格,先按照对x[n]进行解线调,如图2所示,频偏变化率ΔΔf搜素范围±6kHz/s,搜索精度250Hz/s(相邻两个q值之间相差250Hz/s),Q值为 解线调表达式为
上式中对49组解线调结果xq[n]补128个零,做256点FFT,得到49组xq[k],每组xq[k]长度256。将全部的49×256的取模结果|xq[k]|画在网格上,根据最大值对应的坐标,得出此段信号的频偏和频偏变化率。如附图3,得到的频偏为-136.7188Hz,频偏变化率-3000Hz/s,即转化为数字域上得到的频偏估计量频偏变化率分别对频偏变化率和频偏进行补偿。
从第二段数据开始,数据处理的长度都为64,并且默认为频偏变化率已经利用第一段频偏变化率估计结果进行补偿,即使在第一段估计出来的有误差,误差最大为500Hz/s,在64个符号的时间长度,频率积累为如果本段频偏估计得出的频率超过12.8Hz,认为本段频率估计有错,维持前一次频率估计结果,只对每段数据估计出的频偏进行更新并补偿。由数字信号处理知识可知,第一段数据处理后参与频偏补偿的精度为所以从第二段数据开始采用穷举法来得到频偏更加精确的估计。
第二段数据进入精细搜索模块的数据为:
从第三段数据开始,需要对数据进行解模糊,保证与前一段相位保持一致。
第三段的初相位为
第四段以及以后各段的处理流程与第三段完全一样,当帧同步成功时,此时相关结果极性为负,帧头结束后,数据段解调结果需要将z[n]符号取反,输出至译码装置,完成译码。
试验结果
本发明提出的解调装置,能够在较低信噪比下,短时间内,对频偏和频偏变化率完成估计并且实时补偿,实现接收机对信号的频偏和频偏变化率准确跟踪,完成相位鉴别和相位纠正,保证解调数据误码率与理论值相符。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.低信噪比短帧突发通信开环解调装置,其特征在于:包括多普勒频偏变化率和多普勒频偏变化率二维搜索模块、精细搜索模块和帧同步模块;
其中,精细搜索模块是复用的,本装置中复用的精细搜索模块完全一样;
本装置的输入数据为解扩模块的输出,本装置的输出数据为解调结果,送入译码模块;
本装置各模块的连接关系如下:
二维搜索模块连接精细搜索模块1;精细搜索模块1与精细搜索模块2相连,精细搜索模块2与精细搜索模块3相连,精细搜索模块3和帧同步模块相连;
本装置中各模块的功能如下:
多普勒频偏变化率和多普勒频偏变化率二维搜索模块的功能是进行二维搜索,完成对接收信号初始频偏和频偏变化率进行联合估计,得到初始频偏和频偏变化率的估计,并且对多普勒频偏变化率和多普勒频偏进行补偿;精细搜索模块的功能是完成频偏的估计并补偿、相位估计以及相位补偿;帧同步模块的功能是解相位模糊。
2.如权利要求1所述的低信噪比短帧突发通信开环解调装置中采用低信噪比短帧突发通信开环解调方法,其特征在于:包括如下步骤:
步骤一、基于本装置的输入数据运行多普勒频偏变化率和多普勒频偏变化率二维搜索模块,完成对接收信号初始频偏和频偏变化率进行联合估计,得到初始频偏和频偏变化率的估计;
其中,本装置的输入数据也称为第一段数据;
步骤二、基于第一段数据得到的初始频偏的估计结果,运行精细搜索模块1、精细搜索模块2以及精细搜索模块3得到第二段数据,基于第二段数据得到的频偏的精确估计结果并进行补偿,从而完成数据解调;
步骤三、对本段数据频偏进行精确估计并且补偿,对本段数据相位进行解模糊,即需要与前一段数据相位保持一致;
步骤四、基于后续数据段,重复运行步骤三,直到导频段数据处理结束;所得的相位估计结果能和前一段数据相位保持一致;
步骤五、步骤一到步骤四将导频段数据处理结束后,进入帧头段的数据处理,具体的处理方式是是对帧头数据取相关,得出帧头相关结果;
其中,帧头通常为M序列或者Gold码序列,帧头是数据段的开始;
帧头通过解扩、频率估计、补偿和相位估计,输出数据z[n],与本地的帧头模板做如下公式(14)的相关运算:
此公式(14)的相关过程与解扩过程的运算一样,区别是解扩中是在码片级别上做的相关操作,解扩是相关长度为一个符号的时间长度;而帧同步是在符号级别上做的相关运算,相关运算的长度取决于帧头的长度;rheader[n]为帧头相关结果;sheader_tep[n]为帧头本地模板;M为帧头长度;
步骤六、将帧头相关结果与本地模板做相关,依据相关结果绝对值是否超过预定门限来确定帧同步是否成功,并进行相应操作,具体为:对帧头相关运算结果rheader[n]取模操作,如果取模结果超过门限,则认为帧同步成功,跳至步骤七,启动译码装置;否则认为本帧数据接收失败,返回至步骤一等待下一帧数据到来;
步骤七、帧同步成功时,根据相关结果极性,对相位进行解模糊;具体为:根据为解决“倒π”问题,当超过门限的时刻,需要确认超过门限时相关结果的正负;如果极性为正,说明此时保持的相位为0,直接输出z[n]进译码装置;如果相关结果极性为负,说明此时保持的相位为π,需要将z[n]取反,之后给译码模块;译码结束后,本帧数据接收成功,本装置返回至步骤一等待下一帧数据到来;
至此,从步骤一到步骤七,完成了低信噪比短帧突发通信开环解调方法。
3.根据权利要求2所述的低信噪比短帧突发通信开环解调方法,其特征在于:步骤一中,本装置的输入数据表示为如下公式(1):
具体通过二维网格搜索进行近似求解,将频偏和频偏变化率全部不确定范围按照一定精度划分成二维平面网格,得到二维平面的方式是:对于频偏变化率,在其变化范围内,取Q个即按照对x[n]进行Q次“解线性调频调制”,简称“解线调”,对应的符号xq[n]表达为如下公式(3):
对xq[n]进行补零至长度为K点,再进行K点FFT运算,得到xq[k];对xq[k]取模,结果为|xq[k]|;再从全部Q×K个取模结果|xq[k]|找到最大值,根据这个最大值在二维网格上的位置,得到β和Δω估计量和其中,表示步骤一对频偏变化率的估计结果,表示步骤一对初始频偏的估计结果;
完成对接收信号初始频偏和频偏变化率进行联合估计,得到初始频偏和频偏变化率的估计,从第二段数据开始,频偏变化率的补偿均使用第一段数据得到的频偏变化率的估计结果,并且默认频偏变化率已经补偿完毕;即使步骤一得出的频偏变化率估计结果存在误差,即存在残余的多普勒频偏变化率,只要第二段数据处理长度和Q取值合适,完全可以忽略残余的多普勒频偏变化率对频偏估计造成的影响;各段数据只需要对频偏进行估计并补偿即可,完成数据的解调。
4.根据权利要求3所述的低信噪比短帧突发通信开环解调方法,其特征在于:步骤二、具体为:
步骤2.1基于步骤一中对初始频偏的估计值进行补偿,得出公式(4)的第二段数据:
步骤2.2运行精细搜索模块,得到频偏的精确估计;
运行本装置精细搜索模块,精细搜索模块采用穷举法得到频偏的精确估计,进而得出第二段数据的相干解调结果,具体为:
步骤2.21先对第二段数据(4)进行平方去调制,再与本地内积模板的平方做内积,本地内积模板表达式如公式(5):
与本地内积模板的平方做内积的具体方法是将Δω2等间隔分成P个小区间,用P个数字频率不同的复正弦波的平方与x2[n]内积,再取内积模值最大者为复正弦波数字频率为-2Δω2的一个估计,基于此估计得到用于补偿频偏;相位按照下式(6)估计:
上式(6)中的“除以2”操作会导致相位模糊,即和都满足公式(6)的要求,由于初始值未知,在这里取或者都可以;当数据处理到帧头段时,利用帧头的相关结果极性,来确定绝对相位,不妨先取得到第二段数据相干解调结果为下式(7):
5.根据权利要求4所述的低信噪比短帧突发通信开环解调方法,其特征在于:步骤三具体为:
本段数据为:
其中,第一次运行步骤三时,本段数据指第三段数据,其后依次类推;
步骤3.1:运行精细搜索模块,得到频偏的精确估计;精细搜索模块采用穷举法得到频偏的精确估计,进而得出本段数据的相干解调结果,具体为:
先对本段数据进行平方去调制,再与本地内积模板的平方做内积,本地内积模板表达式如公式(5);与本地内积模板的平方做内积的具体方法是将Δω2等间隔分成P个小区间,用P个数字频率不同的复正弦波的平方与x2[n]内积,再取内积模值最大者为复正弦波数字频率为-2Δω3的一个估计,基于此估计得到用于补偿频偏;相位按照下式(9)估计:
步骤3.2:对步骤3.1输出的相位进行解模糊,解模糊的参考相位根据前一段频偏和初相位的结果来构造,因为前一段的末相位就是本段的初相位;
这因为末相位和初相位之间隔了N-1个符号;如果考虑扩频***的解扩过程,则第三段的初相位和第二段的末相位还稍有区别,应作(10):
由于第二段本身就存在相位模糊,所以第三段和以后各段只需要严格跟前一段相位保持一段相位连续即可;具体的处理过程如下(11):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711041223.2A CN107835035B (zh) | 2017-10-30 | 2017-10-30 | 低信噪比短帧突发通信开环解调方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711041223.2A CN107835035B (zh) | 2017-10-30 | 2017-10-30 | 低信噪比短帧突发通信开环解调方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107835035A CN107835035A (zh) | 2018-03-23 |
CN107835035B true CN107835035B (zh) | 2020-06-16 |
Family
ID=61651133
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711041223.2A Expired - Fee Related CN107835035B (zh) | 2017-10-30 | 2017-10-30 | 低信噪比短帧突发通信开环解调方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107835035B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108900452A (zh) * | 2018-05-25 | 2018-11-27 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 减少频率窗的同步检测方法 |
CN108964783B (zh) * | 2018-07-30 | 2020-03-10 | 中国电子科技集团公司第五十四研究所 | 大频偏条件下相干光接收机载波同步方法 |
CN109728827B (zh) * | 2018-12-13 | 2022-11-11 | 航天恒星科技有限公司 | 一种序列辅助捕获低信噪比tpc编译码*** |
CN111884700B (zh) * | 2020-06-03 | 2021-11-23 | 中国人民解放军军事科学院国防科技创新研究院 | 基于低轨卫星的飞行员定位报文处理装置及处理方法 |
CN111884984B (zh) * | 2020-06-29 | 2022-09-02 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 快变载波多普勒频移捕获*** |
CN111865865B (zh) * | 2020-08-04 | 2021-06-15 | 北京空天智数科技有限公司 | 适用于高灵敏度星载ads-b接收机的频偏及相偏估计方法 |
CN112751797B (zh) * | 2020-12-29 | 2023-11-03 | 厦门城市职业学院(厦门开放大学) | Ofdma上行载波频偏盲蔽式估测法 |
CN114285709B (zh) * | 2021-12-31 | 2023-04-25 | 北京中科晶上科技股份有限公司 | 接收信号的相位跟踪方法、装置和信号处理*** |
CN115209519B (zh) * | 2022-06-02 | 2024-01-26 | 四川大学 | 短帧开环结构的无线时间同步设备 |
CN115379550B (zh) * | 2022-07-29 | 2024-05-31 | 西安空间无线电技术研究所 | 一种基于离散导频的突发帧同步方法、装置及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6125136A (en) * | 1997-12-31 | 2000-09-26 | Sony Corporation | Method and apparatus for demodulating trellis coded direct sequence spread spectrum communication signals |
CN102316058A (zh) * | 2011-03-18 | 2012-01-11 | 中国科学院上海微***与信息技术研究所 | 非静止轨道卫星dqpsk通信的相干解调装置 |
CN102546500A (zh) * | 2012-03-20 | 2012-07-04 | 西安电子科技大学 | 基于导频和软信息联合辅助的soqpsk载波同步方法 |
CN104158582A (zh) * | 2014-07-04 | 2014-11-19 | 航天恒星科技有限公司 | 一种用于高速飞行器天基测控的数据处理器*** |
CN105721375A (zh) * | 2016-03-28 | 2016-06-29 | 电子科技大学 | 一种低信噪比短前导突发信号的解调***及方法 |
-
2017
- 2017-10-30 CN CN201711041223.2A patent/CN107835035B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6125136A (en) * | 1997-12-31 | 2000-09-26 | Sony Corporation | Method and apparatus for demodulating trellis coded direct sequence spread spectrum communication signals |
CN102316058A (zh) * | 2011-03-18 | 2012-01-11 | 中国科学院上海微***与信息技术研究所 | 非静止轨道卫星dqpsk通信的相干解调装置 |
CN102546500A (zh) * | 2012-03-20 | 2012-07-04 | 西安电子科技大学 | 基于导频和软信息联合辅助的soqpsk载波同步方法 |
CN104158582A (zh) * | 2014-07-04 | 2014-11-19 | 航天恒星科技有限公司 | 一种用于高速飞行器天基测控的数据处理器*** |
CN105721375A (zh) * | 2016-03-28 | 2016-06-29 | 电子科技大学 | 一种低信噪比短前导突发信号的解调***及方法 |
Non-Patent Citations (2)
Title |
---|
一种直接序列扩频***的大频偏二次捕获算法;薛斌等;《北京理工大学学报》;20111130;第31卷(第11期);1351-1354 * |
高动态短时突发扩频信号的快速捕获;贾鹏等;《通信技术》;20150630;第48卷(第6期);657-661 * |
Also Published As
Publication number | Publication date |
---|---|
CN107835035A (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107835035B (zh) | 低信噪比短帧突发通信开环解调方法及装置 | |
CN113238261B (zh) | 低轨卫星扩频通信体制信号捕获跟踪*** | |
CN111884685B (zh) | 数字通信信号同步解调方法及其装置 | |
CN109104390B (zh) | 一种高速信号的捕获和跟踪方法及装置 | |
CN104253774B (zh) | 一种高动态环境下多普勒频偏估计***和方法 | |
CN107342960B (zh) | 一种适合幅度相移键控的非数据辅助频偏估计方法 | |
CN106603451B (zh) | 一种基于延时自相关的高动态多普勒频偏及频偏变化率估计方法 | |
CN108011653B (zh) | 基于自适应的快速捕获跟踪***及方法 | |
CN101969321B (zh) | 基于fft的直接序列扩频***的大频偏二次捕获方法 | |
CN107026810A (zh) | 突发直扩***及其突发直扩波形的伪码同步方法 | |
CN109088838B (zh) | 一种高动态下直扩dpsk信号的伪码—多普勒快速捕获方法 | |
CN110518936B (zh) | 高超声速飞行器大动态多普勒快速捕获方法及通信*** | |
CN106998586A (zh) | 一种高动态环境中无线通信***的同步捕获方法 | |
CN109586761B (zh) | 一种高动态扩频信号的跟踪解调方法 | |
CN106302296A (zh) | 高动态窄带信号频率跟踪方法 | |
CN105763500A (zh) | 一种连续相位调制信号的频偏、时延和相偏的联合同步方法 | |
CN109660309A (zh) | 利用双向通信传输帧同步码测量距离速度钟差频差的方法 | |
CN109655847B (zh) | 一种适于动态信号的快速捕获方法 | |
CN105610755B (zh) | 一种突发信号的频偏估计方法及装置 | |
CN114205200A (zh) | 一种实现vdes***帧头捕获和载波同步的方法 | |
CN106230473B (zh) | 一种dsss_qpsk载波相位正交误差接收补偿***及方法 | |
CN107820212B (zh) | 一种基于移动多媒体广播***实现的定位方法及定位终端 | |
CN110943956B (zh) | 一种星载自动识别***ais的信号解调方法及*** | |
CN115865127B (zh) | 一种直扩信号的参数估计和解调方法 | |
CN117040564A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200616 |
|
CF01 | Termination of patent right due to non-payment of annual fee |