CN108809354B - 一种高动态多普勒捕获的长时间累积方法 - Google Patents
一种高动态多普勒捕获的长时间累积方法 Download PDFInfo
- Publication number
- CN108809354B CN108809354B CN201810563389.9A CN201810563389A CN108809354B CN 108809354 B CN108809354 B CN 108809354B CN 201810563389 A CN201810563389 A CN 201810563389A CN 108809354 B CN108809354 B CN 108809354B
- Authority
- CN
- China
- Prior art keywords
- doppler
- accumulation
- long
- signal
- representing
- 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
Links
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
- H04B1/7075—Synchronisation aspects with code phase acquisition
- H04B1/70758—Multimode search, i.e. using multiple search strategies
Abstract
本发明公开了一种高动态多普勒捕获的长时间累积方法,包括将接收到的MLopt个码元信号划分到连续的M个短时间累积周期内,并进行非相干累积;根据M个累积结果进行最大似然搜索,得到M个多普勒捕获结果;对M个多普勒捕获结果进行最小二乘拟合,得到多普勒加速度的估计值;根据多普勒加速度的估计值,计算相邻码元信号的频率漂移;补偿码元信号的频率漂移;将补偿后的码元信号合并到一个长时间累积周期,并进行非相干累积;根据长时间累积周期的累积结果进行最大似然搜索,得到最终的多普勒捕获结果。本发明抑制了长时间累积导致的能量扩散问题,有效地提高多普勒捕获性能,更适合低信噪比、高动态环境下的多普勒捕获,为后续载波同步和解调奠定基础。
Description
技术领域
本发明涉及一种高动态多普勒捕获的长时间累积方法,属于空间通信领域。
背景技术
随着航天领域的蓬勃发展,通信载体的活动范围已经从地面扩大到空中乃至外太空。一方面,通信距离的增大会造成很大的路径损耗,从而严重降低接收信号的信噪比。另一方面,通信载体的高速(超过第一宇宙速度7.9km/s)和高加速度(可达几十g,g=9.8m/s2)运动使得信号载波产生很大的多普勒偏移。在接收端,信号载波解调需要补偿多普勒,而多普勒补偿的核心之一就是多普勒捕获技术。因此,低信噪比和高动态环境下信号载波的多普勒捕获就成为一个亟待解决的问题。
当前,航天领域的信号主要是扩频体制且多普勒捕获主要分为以下两类:滑动相关和快速傅里叶变换(fast Fourier transform,FFT)。
滑动相关方法首先将频率搜索范围划分为若干个小间隔,然后通过不断地调整本地频率并与接收信号进行相关处理,直至相关值超过捕获门限。但是,在多普勒变化范围很大的情况下,滑动相关捕获方法需要搜索的次数过多,从而导致捕获时间很长。
FFT方法首先将接收到的码元信号从时域转换到频域,然后在频域上进行最大似然搜索得到频域谱峰,是当前应用最为广泛的方法。
低信噪比和高动态环境会严重影响FFT方法的多普勒捕获性能。首先,在低信噪比环境下,对单个码元信号进行快速傅里叶变换后的频域谱峰淹没在强噪声中从而无法凸显出来,无法完成多普勒捕获。针对这种情况,D.Borio,C.O'Driscoll和G.Lachapell等在“Coherent,noncoherent and differentially coherent combining techniques foracquisition of new composite GNSS signals”(IEEE Transactions on Aerospace andElectronic Systems,vol.45,no.3,pp.1227-1240,2009.)中总结了三种信号累积方式(相干、非相干和差分)来提高检测信噪比。通过逐渐积累码元信号能量,多个信号累积的方法使频域谱峰从噪声中凸显出来,从而提高多普勒的捕获概率。然而,当通信载体在高速运动,尤其是存在很大加速度的情形下,不同码元信号上的多普勒并不是恒定的而是随时间快速变化(即,高动态多普勒)。在长时间信号累积过程中,高动态多普勒会导致严重的能量扩散问题,从而降低检测信噪比。针对载体运动加速度对信号累积效果的影响,A.Yasotharan和T.Thayaparan在“Strengths and limitations of the Fourier methodfor detecting accelerating targets by pulse Doppler radar”(IEE Proceedings ofRadar,Sonar and Navigation,vol.149,no.2,pp.83-88,2002)中发现检测信噪比是信号累积个数的凸函数并且存在唯一的最优累积周期长度。另外,通过大量数据分析,该项研究还给出了最优累积周期长度的经验表达式(即,最优累积周期长度式中,fa和fs分别表示多普勒加速度和采样频率)。由此可见,运动加速度越大,最优累积周期长度越小,检测信噪比越小,捕获概率也会随之越低。
通过以上的分析,我们可以发现:一方面,低信噪比要求长时间信号累积以提高检测信噪比;另一方面,高加速度引起的能量扩散问题使得最优信号累积长度较小。因此,在低信噪比和高加速度环境下,如何提高多普勒的捕获性能就成为一个极富挑战性的技术难题。
发明内容
本发明提供了一种高动态多普勒捕获的长时间累积方法,解决了低信噪比、高动态(高速、高加速)环境下信号载波的多普勒捕获问题。
为了解决上述技术问题,本发明所采用的技术方案是:
一种高动态多普勒捕获的长时间累积方法,包括以下步骤,
将接收到的MLopt个码元信号划分到连续的M个短时间累积周期内,并进行非相干累积;
根据M个累积结果进行最大似然搜索,得到M个多普勒捕获结果;
对M个多普勒捕获结果进行最小二乘拟合,得到多普勒加速度的估计值;
根据多普勒加速度的估计值,计算相邻码元信号的频率漂移;
补偿码元信号的频率漂移;
将补偿后的码元信号合并到一个长时间累积周期,并进行非相干累积;
根据长时间累积周期的累积结果进行最大似然搜索,得到最终的多普勒捕获结果。
接收的码元信号为扩频信号,先对其进行解扩处理,然后对解扩后的码元信号进行快速傅里叶变换。
第i个解扩后的码元信号为x[i]=[x0,...,xN-1],
其中,i∈[1,MLopt],a表示振幅,b表示比特速率为rb的二进制码元数据信息,fa和fs分别表示多普勒加速度和采样频率,fd表示码元上的多普勒初始频偏,表示初始相位,wn表示高斯白噪声,N=fs/rb表示码元信号长度。
对每个长度为N的码元信号进行Nf点的快速傅里叶变换,Nf≥N;
x[i]对应的Nf点频域向量I[i]为,
I[i]=fft(x[i])
其中,fft(·)表示Nf点的快速傅里叶变换。
根据M个累积结果进行谱峰搜索,得到M个多普勒捕获结果;
频率漂移公式为,
第i′个码元信号的频率漂移补偿量为(i′-1)Δf,i′∈[2,MLopt],Δf表示频率漂移。
本发明所达到的有益效果:本发明抑制了长时间累积导致的能量扩散问题,有效地提高多普勒捕获性能,更适合低信噪比、高动态环境下的多普勒捕获,为后续载波同步和解调奠定基础。
附图说明
图1为本发明的流程图;
图2为多普勒捕获的长时间累积原理图;
图3为多普勒捕获概率提升曲线仿真图;
图4为传统傅里叶变换方法下多普勒捕获误差仿真图;
图5为本发明方法下多普勒捕获误差仿真图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1和2所示,一种高动态多普勒捕获的长时间累积方法,包括以下步骤:
步骤1,通信载体将接收到的MLopt个码元信号进行解扩处理。
由于接收的码元信号为扩频信号,为了便于后续处理,需先对其进行解扩,第i个解扩后的码元信号可表示为x[i]=[x0,...,xN-1],i∈[1,MLopt],M≥3,其中,
式中,a表示振幅,b表示比特速率为rb的二进制码元数据信息,fa和fs分别表示多普勒加速度和采样频率,fd表示码元上的多普勒初始频偏,表示初始相位,wn表示高斯白噪声,N=fs/rb表示码元信号长度。
步骤2,对每个长度为N的码元信号进行Nf点的快速傅里叶变换,Nf≥N。
x[i]对应的Nf点频域向量I[i]为,
I[i]=fft(x[i])
其中,fft(·)表示Nf点的快速傅里叶变换。
步骤3,将接收到的MLopt个码元信号划分到连续的M个短时间累积周期内,并进行非相干累积。
步骤4,根据M个累积结果进行谱峰搜索,得到M个多普勒捕获结果。
频率漂移公式为:
式中,Δf表示频率漂移。
步骤7,补偿码元信号的频率漂移,具体如下:
以第1个码元信号I[1]为基准,补偿第2到MLopt个码元信号的频率漂移,其中,第i′个码元信号的频率漂移补偿量为(i′-1)Δf,i′∈[2,MLopt],得到MLopt个补偿后的信号。
步骤8,将补偿后的码元信号合并到一个长时间累积周期,并进行非相干累积。
步骤9,根据长时间累积周期的累积结果进行最大似然搜索,得到最终的多普勒捕获结果。
为了进一步验证上述方法,以空间飞行器测控为例,信号载波频率30GHz,飞行器运动速度和加速度的最大值分别为7.9×103m/s和50g(g=9.8m/s2),根据运动速度和加速度,可以得到多普勒及其加速度的最大值分别为790kHz和49kHz/s。扩频C/A码长度为1023且码率为10.23Mcp/s,码元二进制数据速率为10kb/s,采样频率20.46MHz,码元信号初始相位0,多普勒捕获可容忍误差门限1.0×103Hz,FFT点数2048,信噪比范围-33dB~-44dB。
从附图3可以看出,时频分析方法通常要求较高的信噪比条件,其捕获概率随信噪比的下降而急剧下降。传统FFT在长时间累积过程中导致严重的能量扩散问题,从而造成捕获概率下降。本发明方法利用多个短时间累积周期来构建一个长时间累积周期,可以抑制加速度引起的能量扩散问题,从而明显地提高多普勒捕获概率。从附图4和附图5可以看出,与传统傅里叶方法相比,本发明方法下多普勒捕获误差也会极大地减小,其幅度大约减小为原来的5%。因此,本发明提出的多普勒捕获的长时间累积方法能够提高捕获性能,即,不仅提高捕获概率而且减小捕获误差,适用于低信噪比和高动态环境。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (4)
1.一种高动态多普勒捕获的长时间累积方法,其特征在于:包括以下步骤,
将接收到的MLopt个码元信号划分到连续的M个短时间累积周期内,并进行非相干累积;
根据M个累积结果进行最大似然搜索,得到M个多普勒捕获结果;
对M个多普勒捕获结果进行最小二乘拟合,得到多普勒加速度的估计值;
根据多普勒加速度的估计值,计算相邻码元信号的频率漂移;
频率漂移公式为:
补偿码元信号的频率漂移;
以第1个码元信号I[1]为基准,补偿第2到MLopt个码元信号的频率漂移,其中,第i′个码元信号的频率漂移补偿量为(i′-1)Δf,i′∈[2,MLopt],得到MLopt个补偿后的信号;
将补偿后的码元信号合并到一个长时间累积周期,并进行非相干累积;
根据长时间累积周期的累积结果进行最大似然搜索,得到最终的多普勒捕获结果;
2.根据权利要求1所述的一种高动态多普勒捕获的长时间累积方法,其特征在于:接收的码元信号为扩频信号,先对其进行解扩处理,然后对解扩后的码元信号进行快速傅里叶变换。
4.根据权利要求3所述的一种高动态多普勒捕获的长时间累积方法,其特征在于:对每个长度为N的码元信号进行Nf点的快速傅里叶变换,Nf≥N;
x[i]对应的Nf点频域向量I[i]为,
I[i]=fft(x[i])
其中,fft(·)表示Nf点的快速傅里叶变换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810563389.9A CN108809354B (zh) | 2018-06-04 | 2018-06-04 | 一种高动态多普勒捕获的长时间累积方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810563389.9A CN108809354B (zh) | 2018-06-04 | 2018-06-04 | 一种高动态多普勒捕获的长时间累积方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108809354A CN108809354A (zh) | 2018-11-13 |
CN108809354B true CN108809354B (zh) | 2020-05-26 |
Family
ID=64090365
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810563389.9A Active CN108809354B (zh) | 2018-06-04 | 2018-06-04 | 一种高动态多普勒捕获的长时间累积方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108809354B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111884984B (zh) * | 2020-06-29 | 2022-09-02 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 快变载波多普勒频移捕获系统 |
CN111901266B (zh) * | 2020-07-14 | 2023-02-14 | 南京邮电大学 | 面向低信噪比和高动态场景的载波多普勒频偏捕获方法 |
CN113972951B (zh) * | 2021-10-22 | 2022-05-17 | 金陵科技学院 | 一种低信噪比高动态信号载波多普勒频偏的梳状搜索方法 |
CN114039823B (zh) * | 2021-11-05 | 2022-10-28 | 金陵科技学院 | 一种低信噪比下高动态多普勒频偏及其加速度估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104253774A (zh) * | 2014-09-16 | 2014-12-31 | 重庆邮电大学 | 一种高动态环境下多普勒频偏估计系统和方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104378316A (zh) * | 2014-10-30 | 2015-02-25 | 深圳市国创新能源研究院 | 一种多普勒频偏估计方法和装置 |
CN104880718A (zh) * | 2015-05-31 | 2015-09-02 | 西安电子科技大学 | 基于迭代搜索范围校正的多普勒捕获方法 |
CN106291615B (zh) * | 2016-07-28 | 2019-03-29 | 西安空间无线电技术研究所 | 一种高动态多普勒频偏的两阶段捕获方法 |
-
2018
- 2018-06-04 CN CN201810563389.9A patent/CN108809354B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104253774A (zh) * | 2014-09-16 | 2014-12-31 | 重庆邮电大学 | 一种高动态环境下多普勒频偏估计系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108809354A (zh) | 2018-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108809354B (zh) | 一种高动态多普勒捕获的长时间累积方法 | |
CN103616679B (zh) | 基于差波束调制和波形分析的pd雷达测距测角方法 | |
CN106291615B (zh) | 一种高动态多普勒频偏的两阶段捕获方法 | |
CN107132534B (zh) | 一种高速雷达目标频域检测的优化方法 | |
CN111736128B (zh) | 基于skt-siaf-mscft的相参积累方法 | |
CN102662183B (zh) | Gps信号捕获方法与系统 | |
JP5072694B2 (ja) | 目標検出装置 | |
CN107450055B (zh) | 基于离散线性调频傅立叶变换的高速机动目标检测方法 | |
CN103323829A (zh) | 基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法 | |
CN110518936B (zh) | 高超声速飞行器大动态多普勒快速捕获方法及通信系统 | |
CN108562883A (zh) | 一种多载波雷达系统的最大似然距离估计算法 | |
CN110943758A (zh) | 基于多普勒频率变化率和频率联合搜索的二次捕获方法 | |
CN108415013A (zh) | 基于调频斜率精搜索的低信噪比信号参数提取方法 | |
CN107132513B (zh) | 基于相关距离的距离扩展目标检测方法 | |
CN106788584B (zh) | 一种适应于长码字应用的改进pmf-fft捕获方法 | |
CN113534115B (zh) | 一种主被动联合处理的声纳目标高精度线谱提取方法 | |
CN114545342A (zh) | 利用多通道侦察接收机的雷达脉冲信号参数测量方法 | |
CN113972951B (zh) | 一种低信噪比高动态信号载波多普勒频偏的梳状搜索方法 | |
CN111901266B (zh) | 面向低信噪比和高动态场景的载波多普勒频偏捕获方法 | |
CN112714088A (zh) | 动态多普勒频偏估计方法、装置及存储介质 | |
CN116299166B (zh) | 一种用于chirp信号的低复杂度融合测距方法 | |
CN112968850A (zh) | 微弱bpsk信号载波捕获方法 | |
CN105656511B (zh) | 一种适应于有频偏和低信噪比环境下的差分相关捕获方法 | |
CN109001671B (zh) | 一种跳频信号的目标检测和参数估计方法及装置 | |
CN114280532B (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 |