CN109412644A - 一种直扩msk信号的多普勒频率估计方法 - Google Patents

一种直扩msk信号的多普勒频率估计方法 Download PDF

Info

Publication number
CN109412644A
CN109412644A CN201811064884.1A CN201811064884A CN109412644A CN 109412644 A CN109412644 A CN 109412644A CN 201811064884 A CN201811064884 A CN 201811064884A CN 109412644 A CN109412644 A CN 109412644A
Authority
CN
China
Prior art keywords
signal
frequency
code
doppler
direct expansion
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.)
Granted
Application number
CN201811064884.1A
Other languages
English (en)
Other versions
CN109412644B (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201811064884.1A priority Critical patent/CN109412644B/zh
Publication of CN109412644A publication Critical patent/CN109412644A/zh
Application granted granted Critical
Publication of CN109412644B publication Critical patent/CN109412644B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/70754—Setting of search window, i.e. range of code offsets to be searched

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种直扩MSK信号的多普勒频率估计方法，将直扩MSK信号构造成一种近似直扩BPSK信号形式，接着在大多普勒低信噪比情况下，采用交叉模糊函数的三角拟合对多普勒频偏值进行优化。具体为：利用中频匹配滤波器将接收到直扩MSK信号转换为近似直扩BPSK信号；对近似直扩BPSK信号利用基于FFT的伪码并行捕获算法进行码相位搜索；通过交叉模糊函数的三角拟合对多普勒频偏值进行优化，得到精确的多普勒频率估计值，并利用频率精估计前后的振幅、频率信息提高捕获概率。本发明能够有效提高多普勒频偏估计精度。

Description

一种直扩MSK信号的多普勒频率估计方法
技术领域
本发明属于直扩信号同步技术领域，具体而言为一种直扩MSK信号的多普 勒频率估计方法。
背景技术
目前，扩频***多采用BPSK、QPSK调制方式，但这两种调制方式无法适 用于存在严重非线性失真、多普勒频移与多径衰落的应用领域。直扩MSK通信 ***同时具有扩频***的低截获性、多用户随机选址能力、抗干扰性能强等优点 和MSK信号的包络恒定、频谱利用率高、能量集中、旁瓣衰减快、带外辐射功 率低、对非线性失真不敏感等优点，在战术数据链、民用航空地空数据链、导弹 制导指令传输、卫星通信等领域得到了广泛应用。所以，直扩MSK信号在直扩 BPSK/QPSK信号无法适用的领域仍然具有很好的应用前景。
G.J.R.Povey等人首先提出基于数字部分匹配滤波器与FFT相结合 (PMF-FFT)的捕获模型，这种方法虽然在一定程度上缓解了多普勒频偏对伪码 捕获性能的影响，并且实现了伪码相位和载波频偏的二维捕获，但是这种方法主 要适应于MPSK信号，并且多普勒频偏的捕获范围较小，高动态、低信噪比环 境下仍然不适用。
发明内容
本发明的目的在于提供一种直扩MSK信号的多普勒频率估计方法，解决高 动态、低信噪比环境下扩频信号的伪码相位和多普勒频偏二维捕获不准确的问题。
实现本发明目的的技术解决方案为：一种直扩MSK信号的多普勒频率估计 方法，包括以下步骤：
步骤1、利用中频匹配滤波器将接收到直扩MSK信号转换为近似直扩BPSK 信号；
步骤2、对近似直扩BPSK信号利用基于FFT的伪码并行捕获算法进行码相 位搜索，在频域实现本地伪码与接收伪码的圆周相关运算，得到所***相位上的 相关值；
步骤3、对频域的非相关积分输出进行检测，对非相干积分超过恒虚警检测 阈值的结果，通过对交叉模糊函数的三角拟合进行多普勒频偏值fd优化，得到精 确的多普勒估计值。
优选地，步骤1中采用中频匹配滤波器h(t)接收直扩MSK信号，滤波器h(t) 的冲激响应为：
其中，Tc＝1/Rc为扩频码周期，Rc为扩频码速率；fc为载波频率。
优选地，步骤1中得到的近似直扩BPSK信号表达式为：
其中，A为信号幅度，为初始相位，Ts为接收信号的采样间隔，d(·)为传 输符号，c(·)为PN码，τ为接收信号与本地伪码之间的码相位误差，fd为多普 勒频偏。
优选地，步骤2中，接收信号与本地PN码相关的结果为：
其中，为本地PN码的副本，为码相位差估计值， 为第k个搜索单元的多普勒频率估计值，K为频率搜索的 单元总数，fsp为频率搜索步进，N是相关过程总的样本数。
优选地，步骤3中第k个频率搜索单元的非相关积分输出的结果为：
其中，TN＝NTs，L为非相干积分长度，Rc(·) 为PN码的归一化自相关函数。
优选地，步骤3中，假设第k个频率搜索区中的非相干积分结果Sk超过 恒虚警检测的阈值，保存该超过阈值的点和相邻频率搜索中相同位置的两个 点，根据幅值大小的顺序，将这三个点定义为Pmax，Pmid和Pmin，其频率和幅 度分别定义为fmax，fmid和fmin和Amax，Amid和Amin，采用三角拟合找到交叉 模糊函数的精确顶点位置。
优选地，精确的多普勒估计值为：
其中，且α∈[-0.5,0]。
本发明与现有技术相比，其显著优点为：1)本发明通过中频匹配滤波器将 直扩MSK信号转换为直扩BPSK信号，降低了接收机同步的难度；2)本发明 利用对交叉模糊函数的三角拟合，对多普勒频率值进行精确的估计；3)本发明 有效提高接收机的捕获概率。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是基于FFT的伪码并行捕获原理框图。
图2是基于FFT的伪码并行捕获方法流程图。
图3是信号正确检测与虚警时的三角拟合方案示意图，(a)是正确检测时的 三角拟合方案示意图，(b)是发生虚警时的三角拟合方案示意图。
图4是一种直扩MSK信号的多普勒频率估计方法的流程图。
具体实施方式
如图4所示，一种直扩MSK信号的多普勒频率估计方法，首先将直扩MSK 信号转换成近似直扩BPSK信号，即将中频直扩MSK信号经冲激响应为h(t)的 中频匹配滤波器进行处理：
其中，Tc＝1/Rc为扩频码码片周期，Rc为扩频码速率；fc为 载波频率。根据冲激响应，可求得匹配滤波器的频率响应为：
用凸优化算法求得匹配滤波器的滤波器系数。将中频直扩MSK信号经此匹 配滤波器，可得输出结果为：
其中，A为信号幅度，为初始相位，Ts为接收信号的采样间隔，d(·)为 传输符号，c(·)为PN码，τ为接收信号与本地伪码之间的码相位误差，fd为多 普勒频偏。
如图1所示，基于FFT的并行捕获算法基本思想是两个离散信号时域的圆 周相关等价于其频域信号的共轭相乘，因此利用FFT和IFFT计算可以在频域 实现本地伪码与接收伪码的圆周相关运算，通过3次FFT计算即可计算得到 所***相位上的相关值。接收信号与本地PN码相关的结果为
其中，为本地PN码的副本，为码相位差估计值， 为第k个搜索单元的多普勒频率估计值，K为频率搜索的 单元总数，fsp为频率搜索步进，N是相关过程总的样本数。
对频域的非相关积分输出进行检测，第k个频率搜索单元的非相关积分输出 的结果为：
其中，TN＝NTs，L为非相干积分长度，Rc(·) 为PN码的归一化自相关函数。
假设第k个频率搜索区中的非相干积分结果Sk'超过恒虚警(constant falsealarmrate,CFAR)检测的阈值，保存该超过阈值的点和相邻频率搜索中相同位置 的两个点。根据幅值大小的顺序，将这三个点定义为Pmax，Pmid和Pmin，其频率 和幅度分别定义为fmax，fmid和fmin和Amax，Amid和Amin，采用三角拟合来找到 交叉模糊函数(cross-ambiguityfunction，CAF)的精确顶点位置。拟合情况如图 2所示。首先，评估Pmax和Pmin的线性方程。其次，评估在Pmid上拟合的线性方 程，该方程的线性系数与第一步中导出的系数符号相反。最后，根据等腰三角形 的特征，获得CAF的顶点作为拟合在两条线上的等腰三角形的顶点，其被定义 为点Pv，频率fv和幅度Av表示为：
如图3a所示，如果信号被正确检测且fmin＞fmax＞fmid(fmid＞fmax＞fmin的分析情况类似)，那么三个点的幅度值可表示为：
Amax＝A|sin c(Δfd,k'TN)|Rc(Δτ)
Amid＝A|sin c(Δfd,k'-1TN)|Rc(Δτ) (8)
Amin＝A|sin c(Δfd,k'+1TN)|Rc(Δτ)
由式(8)，fv与Av表示为
Av＝Amax+(Amid-Amin)/2 (10)
比较Pmax、Pv的频率与幅值，可得
Av≥Amax (12)
频率精估计后的剩余频差为：
其中，为频率精估计前的剩余频差且为 了保证Pmax、Pmid和Pmin均在CAF的主瓣内且三者幅值大小可明显区分，设 置频率搜索步进fsp＝2/(3TN)。精确估计后频率估计值为：
其中，且α∈[-0.5,0]。因此，Δf的范围为
(8×10-5)·fsp≤Δf≤0.5·fsp (15)
由式(15)可知，三角拟合后多普勒频率估计精度大大提高。
相反，如果没有正确检测信号，即发生虚警时，如图3b所示，三点的 频率与幅值关系将改变。式(11)(12)不再适用但可用来排除虚警。
在频率搜索中，邻近正确频率搜索单元的虚警概率高于其他搜索单元。
因此，本发明中三角拟合的最大虚警概率为：
其中，Pd是非相干积分的检测概率，为非相 干积分的虚警概率，Vt是检测门限。由于Pd＜1，Pfa＜1，可得因此，在相同的虚警概率要求下，可以将CFAR检测门限设置得较低，提高 了检测概率。根据这一特点，该方法不仅可以提高多普勒频率的估计精度， 而且可以提高捕获概率。
频率精估计前的剩余频差的范围为[0,fsp/2]，频率精估计后剩余频差范围缩 小为[(8×10-5)·fsp,0.5·fsp]，可见三角拟合后多普勒频率估计精度大大提高。

Claims (7)

1.一种直扩MSK信号的多普勒频率估计方法，其特征在于，包括以下步骤：
步骤1、利用中频匹配滤波器将接收到直扩MSK信号转换为近似直扩BPSK信号；
步骤2、对近似直扩BPSK信号利用基于FFT的伪码并行捕获算法进行码相位搜索，在频域实现本地伪码与接收伪码的圆周相关运算，得到所***相位上的相关值；
步骤3、对频域的非相关积分输出进行检测，对非相干积分超过恒虚警检测阈值的结果，通过对交叉模糊函数的三角拟合进行多普勒频偏值fd优化，得到精确的多普勒估计值。
2.根据权利要求1所述的一种直扩MSK信号的多普勒频率估计方法，其特征在于，步骤1中采用中频匹配滤波器h(t)接收直扩MSK信号，滤波器h(t)的冲激响应为：
其中，Tc＝1/Rc为扩频码周期，Rc为扩频码速率；fc为载波频率。
3.根据权利要求1所述的一种直扩MSK信号的多普勒频率估计方法，其特征在于，步骤1中得到的近似直扩BPSK信号表达式为：
其中，A为信号幅度，为初始相位，Ts为接收信号的采样间隔，d(·)为传输符号，c(·)为PN码，τ为接收信号与本地伪码之间的码相位误差，fd为多普勒频偏。
4.根据权利要求1所述的一种直扩MSK信号的多普勒频率估计方法，其特征在于，步骤2中，接收信号与本地PN码相关的结果为：
其中，为本地PN码的副本，为码相位差估计值，为第k个搜索单元的多普勒频率估计值，K为频率搜索的单元总数，fsp为频率搜索步进，N是相关过程总的样本数。
5.根据权利要求1所述的一种直扩MSK信号的多普勒频率估计方法，其特征在于，步骤3中第k个频率搜索单元的非相关积分输出的结果为：
其中，TN＝NTs，L为非相干积分长度，Rc(·)为PN码的归一化自相关函数。
6.根据权利要求1所述的一种直扩MSK信号的多普勒频率估计方法，其特征在于，步骤3中，假设第k个频率搜索区中的非相干积分结果Sk超过恒虚警检测的阈值，保存该超过阈值的点和相邻频率搜索中相同位置的两个点，根据幅值大小的顺序，将这三个点定义为Pmax，Pmid和Pmin，其频率和幅度分别定义为fmax，fmid和fmin和Amax，Amid和Amin，采用三角拟合找到交叉模糊函数的精确顶点位置。
7.根据权利要求1所述的一种直扩MSK信号的多普勒频率估计方法，其特征在于，精确的多普勒估计值为：
其中，且α∈[-0.5,0]。
CN201811064884.1A 2018-09-13 2018-09-13 一种直扩msk信号的多普勒频率估计方法 Active CN109412644B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811064884.1A CN109412644B (zh) 2018-09-13 2018-09-13 一种直扩msk信号的多普勒频率估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811064884.1A CN109412644B (zh) 2018-09-13 2018-09-13 一种直扩msk信号的多普勒频率估计方法

Publications (2)

Publication Number Publication Date
CN109412644A true CN109412644A (zh) 2019-03-01
CN109412644B CN109412644B (zh) 2021-02-12

Family

ID=65464745

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811064884.1A Active CN109412644B (zh) 2018-09-13 2018-09-13 一种直扩msk信号的多普勒频率估计方法

Country Status (1)

Country Link
CN (1) CN109412644B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110113074A (zh) * 2019-05-17 2019-08-09 中国电子科技集团公司第五十四研究所 一种基于fpga的高速伪码并行捕获与跟踪方法
CN110912844A (zh) * 2019-11-28 2020-03-24 上海瀚讯信息技术股份有限公司 一种基于大数据分析的信道估计优化方法
CN110958035A (zh) * 2019-11-12 2020-04-03 天津津航计算技术研究所 一种基于三角包络拟合的频率二次捕获装置及方法
CN113640842A (zh) * 2021-08-25 2021-11-12 北京理工大学 一种基于智能多普勒搜索的直扩信号捕获方法
CN113890563A (zh) * 2021-09-24 2022-01-04 天津津航计算技术研究所 一种直扩信号时频二维内插精估方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6879647B1 (en) * 2000-09-29 2005-04-12 Northrop Grumman Corporation Radio receiver AM-MSK processing techniques
CN105141340A (zh) * 2015-07-24 2015-12-09 南京理工大学 直扩msk信号的全数字接收方法
CN105790788A (zh) * 2016-04-28 2016-07-20 南京理工大学 一种直扩msk信号的伪码-多普勒联合捕获方法
CN107493117A (zh) * 2016-06-12 2017-12-19 南京理工大学 一种高动态下直扩msk信号的二维联合捕获方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6879647B1 (en) * 2000-09-29 2005-04-12 Northrop Grumman Corporation Radio receiver AM-MSK processing techniques
CN105141340A (zh) * 2015-07-24 2015-12-09 南京理工大学 直扩msk信号的全数字接收方法
CN105790788A (zh) * 2016-04-28 2016-07-20 南京理工大学 一种直扩msk信号的伪码-多普勒联合捕获方法
CN107493117A (zh) * 2016-06-12 2017-12-19 南京理工大学 一种高动态下直扩msk信号的二维联合捕获方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
朱雯等: "一种改进的快速频偏捕获算法", 《信息化研究》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110113074A (zh) * 2019-05-17 2019-08-09 中国电子科技集团公司第五十四研究所 一种基于fpga的高速伪码并行捕获与跟踪方法
CN110958035A (zh) * 2019-11-12 2020-04-03 天津津航计算技术研究所 一种基于三角包络拟合的频率二次捕获装置及方法
CN110912844A (zh) * 2019-11-28 2020-03-24 上海瀚讯信息技术股份有限公司 一种基于大数据分析的信道估计优化方法
CN110912844B (zh) * 2019-11-28 2022-03-22 上海瀚讯信息技术股份有限公司 一种基于大数据分析的信道估计优化方法
CN113640842A (zh) * 2021-08-25 2021-11-12 北京理工大学 一种基于智能多普勒搜索的直扩信号捕获方法
CN113640842B (zh) * 2021-08-25 2023-12-05 北京理工大学 一种基于智能多普勒搜索的直扩信号捕获方法
CN113890563A (zh) * 2021-09-24 2022-01-04 天津津航计算技术研究所 一种直扩信号时频二维内插精估方法

Also Published As

Publication number Publication date
CN109412644B (zh) 2021-02-12

Similar Documents

Publication Publication Date Title
CN109412644A (zh) 一种直扩msk信号的多普勒频率估计方法
CN107493117B (zh) 一种高动态下直扩msk信号的二维联合捕获方法
CN101729089B (zh) 通信***的发射机、接收机及其同步方法
CN105790788B (zh) 一种直扩msk信号的伪码-多普勒联合捕获方法
CN102338878B (zh) 一种多径抑制boc码跟踪方法以及码跟踪环
CN106879068B (zh) 一种强多径环境下信号的到达时间估计方法
CN101777933A (zh) 机群链路的加密跳码扩频信号生成与捕获***
CN105743612B (zh) 实时盲解调大频移短时突发信号的方法
CN103905085B (zh) 一种猝发混合扩频水声隐蔽通信方法
CN101567727A (zh) 一种差分循环移位扩频水声通信方法
EP2015463A1 (en) System and method for obtaining frequency and time synchronization in a wideband communication system
CN109088838A (zh) 一种高动态下直扩dpsk信号的伪码—多普勒快速捕获方法
CN105141340A (zh) 直扩msk信号的全数字接收方法
CN104570017B (zh) 一种针对北斗导航卫星***D1导航电文中的Neumann‑Hoffman码规避与剥离方法
KR100809020B1 (ko) 이동 통신 시스템에서의 단말기의 초기 동기 획득 장치 및그 방법
CN105553506B (zh) 一种长码扩频信号的快速捕获方法及装置
CN102170314A (zh) 一种双曲调频扩频水声通信方法
CN104199067A (zh) 多径环境下全球卫星导航***接收机无模糊处理方法
CN107370705A (zh) 高动态弱连续相位调制信号捕获中fft的优化方法
CN110247867A (zh) 水声多普勒估计方法及装置、水声通信方法及***
CN103439718B (zh) 一种高阶boc调制信号的无模糊跟踪单元
CN103873421B (zh) 一种多径信道下的符号同步方法
CN102841359B (zh) 一种抗电文翻转的直扩信号载波伪码二维捕获方法
CN102243309A (zh) Gnss互相关干扰抑制方法和装置
CN102798871A (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