CN115296966A - 一种适用于高动态场景的多普勒码偏估计方法 - Google Patents
一种适用于高动态场景的多普勒码偏估计方法 Download PDFInfo
- Publication number
- CN115296966A CN115296966A CN202210916481.5A CN202210916481A CN115296966A CN 115296966 A CN115296966 A CN 115296966A CN 202210916481 A CN202210916481 A CN 202210916481A CN 115296966 A CN115296966 A CN 115296966A
- Authority
- CN
- China
- Prior art keywords
- doppler
- error
- data
- value
- change
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000004891 communication Methods 0.000 abstract description 7
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 238000001914 filtration Methods 0.000 abstract description 2
- 230000010287 polarization Effects 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009123 feedback regulation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- 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
- H04J—MULTIPLEX COMMUNICATION
- H04J3/00—Time-division multiplex systems
- H04J3/02—Details
- H04J3/06—Synchronising arrangements
- H04J3/0602—Systems characterised by the synchronising information used
- H04J3/0617—Systems characterised by the synchronising information used the synchronising signal being characterised by the frequency or phase
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0008—Modulated-carrier systems arrangements for allowing a transmitter or receiver to use more than one type of modulation
-
- 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)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种适用于高动态场景的多普勒码偏估计方法,涉及通信系统中的高动态同步技术领域。本发明包括低通滤波、采样率变换、定时误差估计、差值计算、码偏计算等步骤,可实现高动态场景下的多普勒码偏估计功能,为时钟频率补偿提供依据,解决了接收机在高动态时的同步难题。本发明适用于通信中任意速率的调相、调幅调相信号,且不需要导频等辅助信息,可对随机信号进行多普勒码偏盲估计,可适应任意多普勒大小,特别适用于通信系统的高动态场景。
Description
技术领域
本发明涉及通信领域,特别是指一种适用于高动态场景的多普勒码偏估计方法。
背景技术
随着我国航空航天技术的快速发展,低轨卫星、高速飞行器等高机动平台的应用越来越多,因此高动态场景下的通信测控技术也越来越受重视。由于高机动平台的轨道存在一定变化,运动方向和速度也一直在改变,要想快速完成时钟同步存在很大困难,因此需要预先对码偏进行估计和补偿,为后续时钟同步打下基础。目前对于高动态的分析,大多是只对多普勒频移进行了分析和计算,很少对码偏进行分析和估计。而且,现有的多普勒估计方法很多依赖于导频等辅助信息,不能在没有辅助信息的情况下进行盲估计。因此对多普勒码偏的盲估计的研究很有必要。
发明内容
有鉴于此,本发明提供一种适用于高动态场景的多普勒码偏估计方法。该方法既能高精度恢复数据,又可以恢复同步时钟,还不需要反馈调整。
本发明的目的是这样实现的:
一种适用于高动态场景的多普勒码偏估计方法,包括以下步骤:
(1)接收低通滤波后的零中频信号r(t);
(2)将r(t)转换为4倍过采样数据r(k);
(3)采用平方法对r(k)进行定时误差估计,估计长度为L个符号,每4*L个数据作为一段,每段数据输出一个定时误差t_error(i);其中,L为整数;
(4)计算相邻两段数据的定时误差的差值t1:
t1(i)=t_error(i+1)-t_error(i)
存储t1,存储长度为M;其中,M为整数,i为数据段的序号;
(5)对t1进行累加求和,再除以间隔的符号个数,得到每个符号的定时变化,即为归一化码偏:
其中,N为累加的段数,N取整数值。
进一步地,步骤(3)中,平方法定时误差估计的方式为:
其中,t_error(i)表示第i段数据的定时误差,tan-1(x)表示对x取反正切值。
进一步地,步骤(5)中N的确定方式为:
1)寻找非零的t1,记录其位置;
2)计算相邻非零值的位置差,即间隔几段数据,记录为Δk;
3)取Δk的最大值或者最小值或者均值作为N值;具体来说:
多普勒变化较小的场景,N取Δk的最大值;多普勒变化较大的场景,N取Δk的最小值;多普勒变化几乎为零的场景,N取Δk的均值;
N确定后,相邻两个非零t1所对应的数据,将其定时误差的变化按线性变化考虑,用误差值的变化除以间隔符号数,完成平均码偏的计算。
本发明相比背景技术具有如下优点:
1、本发明可工作在不同的高动态场景,可适应任意多普勒大小,不需要导频等辅助信息,可对信号进行多普勒码偏盲估计。
2、本发明能够适应多种调制方式,且能够工作于低信噪比的情况。
3、本发明采用纯数字方式实现,可嵌入其它任何数字模块中,便于工程应用。
附图说明
图1是本发明方法的流程图。
图2是步骤(5)码偏计算的流程图。
具体实施方式
下面对本发明做进一步的详细说明。
如图1所示,一种适用于高动态场景的多普勒码偏估计方法,包括以下步骤:
(1)接收低通滤波后的零中频信号r(t);
(2)将r(t)转换为4倍过采样数据r(k);
(3)采用平方法对r(k)进行定时误差估计,估计长度为L个符号。每4*L个数据作为一段,每段数据输出一个定时误差t_error(i);其中,L为整数。
平方法定时误差估计的方式为:
其中,t_error(i)表示第i段数据的定时误差,tan-1(x)表示对x取反正切值。
(4)计算相邻两段数据的定时误差的差值t1,t1(i)=t_error(i+1)-t_error(i),并存储t1,存储长度为M;其中,M为整数;
(5)计算差值的均值。如图2所示,通过对t1的累加求和,再除以间隔的符号个数,得到每个符号的定时变化,即为归一化码偏;N为累加长度,即间隔的符号个数,N取整数值。
步骤(5)中N的确定方式为:
1)寻找非零的t1,记录其位置;
2)然后计算相邻非零值的位置差,即间隔几段数据,记录为Δk;
3)取Δk的最大值或者最小值或者均值作为N值。
多普勒变化较小的场景,N适合取Δk的最大值;多普勒变化较大的场景,N适合取Δk的最小值;多普勒变化几乎为零的场景,N适合取Δk的均值。
总之,本发明包括低通滤波、采样率变换、定时误差估计、差值计算、码偏计算等步骤,可实现高动态场景下的多普勒码偏估计功能,为时钟频率补偿提供依据,解决了接收机在高动态时的同步难题。本发明适用于通信中任意速率的调相、调幅调相信号,且不需要导频等辅助信息,可对随机信号进行多普勒码偏盲估计,可适应任意多普勒大小,特别适用于通信系统的高动态场景。
以上所述仅为本发明的一种具体实施方式,但本发明的保护范围并不局限于此。任何熟悉本技术领域的技术人员根据本发明的技术方案及其发明构思作出的等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (3)
1.一种适用于高动态场景的多普勒码偏估计方法,其特征在于,包括以下步骤:
(1)接收低通滤波后的零中频信号r(t);
(2)将r(t)转换为4倍过采样数据r(k);
(3)采用平方法对r(k)进行定时误差估计,估计长度为L个符号,每4*L个数据作为一段,每段数据输出一个定时误差t_error(i);其中,L为整数;
(4)计算相邻两段数据的定时误差的差值t1:
t1(i)=t_error(i+1)-t_error(i)
存储t1,存储长度为M;其中,M为整数,i为数据段的序号;
(5)对t1进行累加求和,再除以间隔的符号个数,得到每个符号的定时变化,即为归一化码偏:
其中,N为累加的段数,N取整数值。
3.根据权利要求1所述的一种适用于高动态场景的多普勒码偏估计方法,其特征在于,步骤(5)中N的确定方式为:
1)寻找非零的t1,记录其位置;
2)计算相邻非零值的位置差,即间隔几段数据,记录为Δk;
3)取Δk的最大值或者最小值或者均值作为N值;具体来说:
多普勒变化较小的场景,N取Δk的最大值;多普勒变化较大的场景,N取Δk的最小值;多普勒变化几乎为零的场景,N取Δk的均值;
N确定后,相邻两个非零t1所对应的数据,将其定时误差的变化按线性变化考虑,用误差值的变化除以间隔符号数,完成平均码偏的计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210916481.5A CN115296966B (zh) | 2022-08-01 | 2022-08-01 | 一种适用于高动态场景的多普勒码偏估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210916481.5A CN115296966B (zh) | 2022-08-01 | 2022-08-01 | 一种适用于高动态场景的多普勒码偏估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115296966A true CN115296966A (zh) | 2022-11-04 |
CN115296966B CN115296966B (zh) | 2024-06-25 |
Family
ID=83826953
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210916481.5A Active CN115296966B (zh) | 2022-08-01 | 2022-08-01 | 一种适用于高动态场景的多普勒码偏估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115296966B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1284220A (zh) * | 1997-11-20 | 2001-02-14 | 鹰眼技术公司 | 低功率的基于卫星的地面定位系统 |
US6366599B1 (en) * | 1998-03-16 | 2002-04-02 | Trimble Navigation Limited | Fast acquisition of spread-spectrum signals by dynamically varying spacing of search bins |
US20060170589A1 (en) * | 2005-01-28 | 2006-08-03 | Samsung Electronics Co., Ltd. | Apparatus and method for maintaining time synchronization in AGPS receiver |
CN101242197A (zh) * | 2007-02-06 | 2008-08-13 | 上海微小卫星工程中心 | 一种基于数字插值的多普勒模拟实现方法 |
KR20110125581A (ko) * | 2010-05-13 | 2011-11-21 | 한국전자통신연구원 | 항법 신호 처리 장치 및 방법 |
CN102759739A (zh) * | 2011-04-26 | 2012-10-31 | 中国科学院微电子研究所 | 快速联合捕获装置及其方法 |
CN105099498A (zh) * | 2014-05-20 | 2015-11-25 | 中国科学院国家天文台 | 扩频信号捕获系统和方法 |
CN111555775A (zh) * | 2020-04-22 | 2020-08-18 | 中国电子科技集团公司第五十四研究所 | 一种适用于跳频抗干扰系统高动态环境下的同步方法 |
-
2022
- 2022-08-01 CN CN202210916481.5A patent/CN115296966B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1284220A (zh) * | 1997-11-20 | 2001-02-14 | 鹰眼技术公司 | 低功率的基于卫星的地面定位系统 |
US6366599B1 (en) * | 1998-03-16 | 2002-04-02 | Trimble Navigation Limited | Fast acquisition of spread-spectrum signals by dynamically varying spacing of search bins |
US20060170589A1 (en) * | 2005-01-28 | 2006-08-03 | Samsung Electronics Co., Ltd. | Apparatus and method for maintaining time synchronization in AGPS receiver |
CN101242197A (zh) * | 2007-02-06 | 2008-08-13 | 上海微小卫星工程中心 | 一种基于数字插值的多普勒模拟实现方法 |
KR20110125581A (ko) * | 2010-05-13 | 2011-11-21 | 한국전자통신연구원 | 항법 신호 처리 장치 및 방법 |
CN102759739A (zh) * | 2011-04-26 | 2012-10-31 | 中国科学院微电子研究所 | 快速联合捕获装置及其方法 |
CN105099498A (zh) * | 2014-05-20 | 2015-11-25 | 中国科学院国家天文台 | 扩频信号捕获系统和方法 |
CN111555775A (zh) * | 2020-04-22 | 2020-08-18 | 中国电子科技集团公司第五十四研究所 | 一种适用于跳频抗干扰系统高动态环境下的同步方法 |
Non-Patent Citations (2)
Title |
---|
于海亮;胡小平;唐康华;: "惯性信息辅助高动态接收机捕获问题的研究", 航天控制, no. 05, 15 October 2008 (2008-10-15) * |
李冰蟾;毛波;: "卫星信号模拟器中多普勒频移的研究与实现", 科技致富向导, no. 36, 25 December 2014 (2014-12-25) * |
Also Published As
Publication number | Publication date |
---|---|
CN115296966B (zh) | 2024-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109617570B (zh) | 一种无数据辅助的宽带跳频直扩信号全数字同步方法 | |
US9774441B2 (en) | Clock and data recovery using receiver clock spread spectrum modulation and offset compensation | |
CN104215981B (zh) | 一种接收机高动态环境下自适应跟踪方法 | |
CN105684373B (zh) | 具有频率偏移校正和匹配滤波器组解码的cpfsk接收机 | |
CN101841509B (zh) | 采样时钟偏移估计与补偿的方法及接收装置 | |
CN111624402B (zh) | 微弱pm信号载波频率精确估计方法 | |
CN106603451B (zh) | 一种基于延时自相关的高动态多普勒频偏及频偏变化率估计方法 | |
CN111884984B (zh) | 快变载波多普勒频移捕获系统 | |
CN106302296A (zh) | 高动态窄带信号频率跟踪方法 | |
CN115296966A (zh) | 一种适用于高动态场景的多普勒码偏估计方法 | |
CN113900127A (zh) | 基于强跟踪平方根容积卡尔曼滤波的卫星信号跟踪环路 | |
CN110932720B (zh) | 卫星双向时间比对信号的鉴频方法及系统 | |
US8422593B2 (en) | Parallel automatic frequency offset estimation method and apparatus | |
CN114244674B (zh) | 一种超宽带基带接收机频偏估计方法及装置 | |
KR101629493B1 (ko) | 다중 파일럿 신호간의 진폭 및 위상 차 추출 장치 및 방법 | |
CN113965436A (zh) | 一种多普勒频偏及多普勒频偏加速度联合估计方法 | |
CN116319211B (zh) | 一种用于qam信号的多阶卡尔曼载波跟踪方法、跟踪环路及信号接收机 | |
US8571150B2 (en) | Frequency offset compensation apparatus | |
JPH10257117A (ja) | 受信デジタル信号の搬送周波数と受信局部発振器の周波数との間の周波数差を推定する装置、および対応する受信機 | |
CN111245500A (zh) | 一种低轨卫星高动态突发信号快速参数估计方法 | |
CN118487643A (zh) | 一种基于数据辅助的高动态自适应载波跟踪装置 | |
CN116505978B (zh) | 一种蓝牙信号的处理方法、装置、电子设备及存储介质 | |
US7336750B1 (en) | Optimal one-shot phase and frequency estimation for timing acquisition | |
CN115189862B (zh) | 一种高精度同步数据同步时钟恢复方法 | |
CN111190206A (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 |