多普勒域和延迟域二维捕获方法及装置制造方法及图纸

本发明专利技术涉及一种多普勒域和延迟域二维捕获方法及装置，属于直接序列扩频通信技术领域。该方法首先对直扩导频信号取FFT变换和圆周移位，然后与本地伪码的FFT变换序列相乘，经过IFFT变换，得到载波多普勒补偿的匹配滤波结果，然后使用内插拟合方法调整匹配滤波结果的采样速率，实现码多普勒补偿，内插拟合的输出结果经过非相干积累和门限判决，得到直扩导频信号的载波频偏和伪码相位。对比现有技术，本发明专利技术使用载波多普勒和码多普勒联合补偿的捕获方法，克服了传统捕获方法在高动态、大伪码长度、弱信号环境下的不足，通过内插拟合的方式修正了相关峰走动，从而避免非相干积累中的平台化效应，使捕获性能大幅度提高。

【技术实现步骤摘要】
多普勒域和延迟域二维捕获方法及装置
本专利技术涉及一种多普勒域和延迟域二维捕获方法，尤其适用于高动态、低信噪比环境下直接序列扩频系统导频信号的检测及搜索其多普勒频偏、码偏和伪码相位，属于直接序列扩频通信

技术介绍
直接序列扩频通信系统基带数字信号处理主要用到捕获和跟踪两个模块。捕获模块主要实现对接收信号载波频率和伪码相位的粗略估计，为跟踪模块的载波环和码环提供频率和码相位的初始值。假设直接序列扩频通信的扩频码为c(l)，l＝0,1,2,…，L-1，L为扩频码长度。c(l)取值为0或1，采用无线电技术的惯例，c(l)＝0表示逻辑正，c(l)＝1表示逻辑负，并采用根升余弦脉冲成形，脉冲成形滤波器冲击函数为h(t)，则成形后的扩频码形状为其中Tc为码片周期。为了便于捕获，直接序列扩频通信的发射端往往在通信开始时发射一段导频信号，即只包含“1”信息码的扩频信号，使基带信号具有周期性。其基带信号形式如下所示：其中Ts为码元周期，Ts＝LTc；M为导频信号长度。由上式可知，Sb,s(t)是x(t)以Ts为周期重复M次的有限周期信号。经过调制后，射频所发出的信号为，Ssend(t)＝cos(2πfRFt)Sb,s(t)上式中fRF为射频频率，并且为了表达简洁方便，认为在发射时载波初相位和码初相位都为0。经过无线传输过程的多普勒效应和信道延迟，到达接收端的直接序列扩频信号是：Srec(t)＝cos(2πfRFδt+τ0)Sb,s(δt+τ0)其中δ是由于多普勒效应所产生的时间系数，v是通信的收发双方相对径向速度，相向运动时v为正，向背运动时v为负。c是无线电传播速度，一般认为是光速，且有v＜＜c。τ0是码相位延迟。接收机对射频信号进行I/Q两路正交下变频，得到正交基带信号I(t)＝cos(2πfRF(δ-1)t+τ0)Sb,s(δt+τ0)Q(t)＝sin(2πfRF(δ-1)t+τ0)Sb,s(δt+τ0)且有称r(t)为复基带信号，它包含两部分，载波残余和基带伪码Sb,s(δt+τ0)。载波残余表现为复正弦波，其频率为称为载波多普勒偏移频率，简称频偏；基带伪码与发送端的基带信号相比具有两方面变化:首先其具有码相位延迟τ0；其次伪码的形状也相对时间t发生了δ倍的展缩，从码片层面上看，伪码的码片速率从变成其差值称为伪码多普勒偏移速率，简称为码偏。由于常规的捕获过程包括对接收信号的载波频率和伪码相位这两个相对独立的参数进行二维搜索。一种常用的捕获算法是基于快速傅里叶变换的捕获算法(邢方剑,程兰,谢刚等.一种基于快速傅里叶变换和圆周移位的捕获算法[J].科学技术与工程,2014,14(6):196-198,202.DOI:10.3969/j.issn.1671-1815.2014.06.044.)，该算法利用了FFT的调制特性和圆周相关性质，其原理如下：r(t)经过频率为fs的采样，过采样率fs/Rc为正整数，则每个码元周期得到的采样数为N＝(fs/Rc)*L,则整段导频信号的正交基带采样序列为以其连续N个采样点为一组，共分为M组，得到n＝0,1,2,…N-1,m＝0,1,2,…,M-1分别进行FFT变换，得到Rm(k)＝FFT[rm(n)]，k＝0,1,2,…N-1,m＝0,1,2,…,M-1。接收机生成本地基带采样序列：并对其取FFT变换，得到Sb,l(k)＝FFT[Sb,l(n)],k＝0,1,2,…,N-1利用了FFT的调制特性和圆周相关性质，将Rm(k)圆周移位q位后，和Sb,l(k)的共轭逐点相乘，并取逆快速傅里叶变换，可得到上其中表示圆周卷积运算。根据圆周卷积的性质，只有当即时，圆周卷积的结果才比较大，否则圆周卷积结果近似于0。通过FFT变换的圆周移位等效实现时域调制的过程称为频偏补偿。基于快速傅里叶变换的捕获算法正是利用了这一原理，通过不断改变Rm(k)圆周移位量q，检测逆快速傅里叶变换输出结果中是否出现峰值，来估计Δf。当时，忽略由于分辨率带来的误差，认为由于v＜＜c，忽略采样点之间的微小的相位展缩，认为因此，根据伪码的自相关性质，可得出，|ym(n)|的峰值出现在n＝(<τ0fs+δmN>)N处，基于快速傅里叶变换的捕获算法正是利用了这一原理，通过检测峰值位置，来估计τ0。在长距离、弱信号等诸多情况下，捕获模块往往需要通过非相干积累来抬高相关峰高度，以满足捕获准确度要求，但传统算法往往会忽略码偏带来的相关峰走动，在多普勒较大或伪码较长的情况下，相关峰走动会致使非相干累积出现平台化效应，给捕获精确度带来巨大影响。平台化效应的原理是：如上文所述，|ym(n)|的峰值出现在n＝(<τ0fs+δmN>)N处，随着m的变化，峰值位置也会随之发生走动，当多普勒较小、伪码长度较短或非相干积累码元数较少时，传统算法往往忽略相关峰的走动，认为(<τ0fs+δmN>)N＝(<τ0fs+mN>)N＝(<τ0fs>)N，所有相关峰都出现在n＝(<τ0fs>)N处；但是在卫星通信等应用环境中，由于高动态、远距离、发射信号功率受限、通信保密需求高等原因，使得多普勒效应强，伪码较长且非相干积累码元数多，致使相关峰不再出现在同一个位置，而呈现小范围的平均分散，非相干积累结果呈现平台状，峰值明显削弱，使捕获精度受到严重影响。
技术实现思路
本专利技术的目的是为了克服传统伪码捕获方法的不足，提供一种载波多普勒和码多普勒联合补偿方法，实现多普勒域和延迟域的二维捕获，可以保证在高动态、低信噪比条件下，直接序列扩频系统导频信号快速捕获成功。本专利技术方法的目的是通过以下技术方案实现的：一种多普勒域和延迟域二维捕获方法，包括以下步骤：1、接收机天线所接收到的导频信号，通过射频前端已完成正交解调及AD采样，并可能经过必要的降采样等步骤，以I/Q两路正交解调基带采样序列的形式提供给数字信号处理部分；I/Q两路正交解调基带采样序列所组成的复序列简称为输入信号r(n)，其形式如前文所述；2、将r(n)分为M组，每组长度为N个采样，对每一组输入信号进行FFT变换，得到Rm(k)，m＝0,1,2,…,M-1，k＝1,2,…,N-1,N；3、M组FFT序列被提供给并行搜索通道。共有P条搜索通道，根据一定的信道先验信息估计的频偏动态范围为(-Δfmax,Δfmax),表示向上取整函数；4、在搜索通道p中(p＝1,2,…P)，Rm(k)圆周移位qp位，其中得到Rm((k+qp))N，并且与本地存储的本地基带采样序列FFT变换的共轭逐点相乘，并进行IFFT变换，得到M组圆周相关结果；5、在搜索通道p中，对步骤4所得到的M组圆周相关结果取模值，进行并串转换，使用插值拟合滤波器对其进行采样速率转换，将采样率从fs变换到插值拟合之后的序列经过必要的补零或删减之后，重新进行串并转换，分为M组，每组N个采样点，进行非相干积累；6、所有搜索通道都有一相同的信号检测门限A，将所有搜索通道非相干积累结果与信号检测门限A进行比较，对所有超过检测门限A的相关峰，记录其峰值及其位置；7、比较所有通道中的所有超过检测门限的相关峰，取其最大值，其对应通道的FFT变换圆周移位位数为qmax，其峰值在该通道非相干积本文档来自技高网...

【技术保护点】
一种多普勒域和延迟域二维捕获方法，其特征在于，包括训练步骤以及分割步骤，具体如下：步骤一、接收机天线所接收到的导频信号，经正交解调及采样率为fs的AD采样得到输入信号r(n)；步骤二、将r(n)分为M组，每组长度为N个采样，对每一组输入信号进行FFT变换，得到Rm(k)，m＝0，1，2，...，M‑1，k＝1，2，...，N‑1，N；步骤三、M组FFT序列被提供给P条并行搜索通道，其频偏动态范围为(‑Δfmax，Δfmax)，其中，Δfmax表示最大频偏，fs表示输入信号的采样率，表示向上取整函数；步骤四、在搜索通道p中(p＝1，2，...P)，Rm(k)圆周移位qp位，其中得到Rm((k+qp))N，并且与本地存储的本地基带采样序列FFT变换的共轭逐点相乘，并进行IFFT变换，得到M组圆周相关结果；(·)*表示共轭函数。步骤五、在搜索通道p中，对步骤四得到的M组圆周相关结果取模值，进行并串转换，使用插值拟合滤波器对其进行采样速率转换，将采样率从fs变换到插值拟合之后的序列经过必要的补零或删减之后，重新进行串并转换，分为M组，每组N个采样点，进行非相干积累；其中，fRF表示射频频率；步骤六、所有搜索通道都有一相同的信号检测门限A，将所有搜索通道非相干积累结果与信号检测门限A进行比较，对所有超过检测门限A的相关峰，记录其峰值及其位置；步骤七、比较所有通道中的所有超过检测门限的相关峰，取其最大值，其对应通道的FFT变换圆周移位位数为qmax，其峰值在该通道非相干积累序列的下标为nmax，则频偏估计值为Δf^=fsNqmax,]]>码相位估计值为τ^0=nmaxfs.]]>...

【技术特征摘要】
1.一种多普勒域和延迟域二维捕获方法，其特征在于，包括训练步骤以及分割步骤，具体如下：步骤一、接收机天线所接收到的导频信号，经正交解调及采样率为fs的AD采样得到输入信号r(n)；步骤二、将r(n)分为M组，每组长度为N个采样，对每一组输入信号进行FFT变换，得到Rm(k)，m＝0，1，2，...，M-1，k＝1，2，...，N-1，N；步骤三、M组FFT序列被提供给P条并行搜索通道，其频偏动态范围为(-Δfmax，Δfmax),其中，Δfmax表示最大频偏，fs表示输入信号的采样率，表示向上取整函数；步骤四、在搜索通道p中,p＝1，2，...，P，Rm(k)圆周移位qp位，其中得到Rm((k+qp))N，并且与本地存储的本地基带采样序列FFT变换的共轭逐点相乘，并进行IFFT变换，得到M组圆周相关结果；(·)*表示共轭函数；步骤五、在搜索通道p中，对步骤四得到的M组圆周相关结果取模值，进行并串转换，使用插值拟合滤波器对其进行采样速率转换，将采样率从fs变换到插值拟合之后的序列经过必要的补零或删减之后，重新进行串并转换，分为M组，每组N个采样点，进行非相干积累；其中，fRF表示射频频率；步骤六、所有搜索通道都有一相同的信号检测门限A，将所有搜索通道非相干积累结果与信号检测门限A进行比较，对所有超过检测门限A的相关峰，记录其峰值及其位置；步骤七、比较所有通道中的所有超过检测门限的相关峰，取其最大值，其对应通道的FFT变换圆周移位位数为qmax，其峰值在该通道非相干积累序列的下标为nmax，则频偏估计值为码相位估计值为2.一种多普勒域和延迟域二维捕获装置，其特征在于，包括依次相连的串并转换/分组模块，FFT模块和多个并联的搜索通道模块以及判决模块；所述串并转换/分组模块用于将输入信...

【专利技术属性】
技术研发人员：王帅，罗士荀，安建平，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京;11