CN109104390B - 一种高速信号的捕获和跟踪方法及装置 - Google Patents

一种高速信号的捕获和跟踪方法及装置 Download PDF

Info

Publication number
CN109104390B
CN109104390B CN201811061247.9A CN201811061247A CN109104390B CN 109104390 B CN109104390 B CN 109104390B CN 201811061247 A CN201811061247 A CN 201811061247A CN 109104390 B CN109104390 B CN 109104390B
Authority
CN
China
Prior art keywords
speed signal
frequency offset
pilot
signal
effective high
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
Application number
CN201811061247.9A
Other languages
English (en)
Other versions
CN109104390A (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.)
Ruixinfeng Aerospace Technology Beijing Co ltd
Original Assignee
Beijing Rinfon Technology Co ltd
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 Beijing Rinfon Technology Co ltd filed Critical Beijing Rinfon Technology Co ltd
Priority to CN201811061247.9A priority Critical patent/CN109104390B/zh
Publication of CN109104390A publication Critical patent/CN109104390A/zh
Application granted granted Critical
Publication of CN109104390B publication Critical patent/CN109104390B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • H04L27/2659Coarse or integer frequency offset determination and synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2614Peak power aspects
    • H04L27/262Reduction thereof by selection of pilot symbols
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • H04L27/266Fine or fractional frequency offset determination and synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2673Details of algorithms characterised by synchronisation parameters
    • H04L27/2675Pilot or known symbols
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • H04L2027/0024Carrier regulation at the receiver end
    • H04L2027/0026Correction of carrier offset

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

本发明涉及一种高速信号的捕获和跟踪方法及装置,该方法包括以下步骤:对高速信号进行抽取,计算从高速信号中抽取的两个相邻的接收序列的差分相关值,差分相关值最大点为导频序列起始点,以捕获有效高速信号;接收有效高速信号,对有效高速信号进行FFT处理,以获取有效高速信号的粗频偏并校准;对于校准后的有效高速信号,从导频序列起始点的相邻预设范围内抽取接收序列,与预存序列进行互相关计算,互相关峰值最大点为导频序列精确起始点;对有效高速信号进行解扰解扩,通过扫频进行精确频偏估计,并计算相位偏移,以进行精确频偏校正和相偏校正。本发明的方法及装置采用层层递进方式快速跟踪高速移动物体,并且准确解调。

Description

一种高速信号的捕获和跟踪方法及装置
技术领域
本发明涉及载波同步领域,尤其涉及一种高速信号的捕获和跟踪方法及装置。
背景技术
在高速移动的情况下提供良好的网络覆盖一直是移动通信组网的难题。首先,由于高速运动所导致的多普勒频偏会使接收端功率产生变化,移动速度越快,功率越高,产生的频偏也越大,严重影响基站和终端的接收性能。其次,高速移动带来的切换频繁,会使终端接入成功率低、掉话率高,严重影响网络性能。特别地,移动速度剧烈变动的情况下,信号更是无法快速接入网络。
因此,提供一种高速信号的捕获和跟踪方法及装置。
发明内容
鉴于上述问题,提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的高速信号的捕获和跟踪方法及装置,使得能够在移动速度剧烈变动的情况下有效捕获、跟踪和解调高速移动信号。
根据本发明的一个方面,提供一种高速信号的捕获和跟踪方法,包括以下步骤:
接收高速信号,并对高速信号进行抽取,计算从高速信号中抽取的两个相邻的接收序列的差分相关值,差分相关值最大点为导频序列起始点,以捕获有效高速信号,其中,抽取的两个相邻的接收序列为导频序列的两个对称序列;
接收有效高速信号,对有效高速信号进行FFT处理,以获取有效高速信号的粗频偏,并对该有效高速信号进行粗频偏校准;
对于粗频偏校准后的有效高速信号,从导频序列起始点的相邻预设范围内抽取接收序列,与预存序列进行互相关计算,互相关峰值最大点为导频序列精确起始点;
根据导频序列精确起始点对有效高速信号进行解扰解扩,通过扫频进行精确频偏估计,并计算相位偏移,以对有效高速信号进行精确频偏校正和相偏校正,用于下一步的解调。
进一步地,有效高速信号的粗频偏计算公式如下:
Figure BDA0001797154300000021
其中,ΔfINT为粗频偏,Findex为粗频偏索引位置,FFTsize为FFT的大小,fs为采样速率。
进一步地,通过一次求取FFT求平均功率获取粗频偏索引位置,具体包括以下步骤:
第一步,对导频信号去除调制信息后,通过以下公式获取FFT功率图:
fft_pss=FFT(r_pilot.*pss*)
其中,fft_pss为导频信号的频域点,r_pilot.*pss*为导频信号的时域点;
第二步,根据以下公式找出功率最大的位置点:
[maxvalue,maxpos]=max(|fft_pss|)
其中,[maxvalue,maxpos]为功率最大的位置点,fft_pss为导频信号的频域点;
第三步,根据以下条件判断功率最大的位置点是否是粗频偏索引位置:
{maxvalue/mean(|fft_pss|)}>TH
其中,maxvalue为功率最大的位置点的功率值,mean(|fft_pss|)为导频信号的所有频域点的平均功率值,TH为门限;
若条件成立,则功率最大的位置点为粗频偏索引位置。
进一步地,有效高速信号的粗频偏计算公式如下:
Figure BDA0001797154300000031
其中,
Figure BDA0001797154300000032
为粗频偏,M为调制阶数,FFT(r′k)为有效序列r′k的频谱值。
进一步地,FFT(r′k)的计算步骤包括:
计算频偏点的频谱图,并对计算的频谱图进行多次功率平均,得到FFT(r′k)。
进一步地,通过以下公式计算差分相关值:
Figure BDA0001797154300000033
其中,cn为差分相关值,rn+k
Figure BDA0001797154300000034
为两个相邻窗长为Nd的接收序列,k为采样序号。
进一步地,互相关计算公式如下:
Figure BDA0001797154300000041
其中,ck为互相关函数,rn+k为接收序列,sn为预存序列,k为采样序号。
进一步地,精确频偏估计计算公式如下:
Figure BDA0001797154300000042
其中,
Figure BDA0001797154300000043
为精频偏,M为调制阶数,f为采样点功率,Tb为码元周期,k为采样序号。
根据本发明的另一方面,提供一种实现上述高速信号的捕获和跟踪方法的装置,包括:
信号粗捕获模块,用于接收高速信号,并对高速信号进行抽取,计算从高速信号中抽取的两个相邻的接收序列的差分相关值,差分相关值最大点为导频序列起始点,以捕获有效高速信号,其中,抽取的两个相邻的接收序列为导频序列的两个对称序列;
信号粗频偏估计及校正模块,用于接收有效高速信号,对有效高速信号进行FFT处理,以获取有效高速信号的粗频偏,并对该有效高速信号进行粗频偏校准;
信号精捕获模块,用于对于粗频偏校准后的有效高速信号,从导频序列起始点的相邻预设范围内抽取接收序列,与预存序列进行互相关计算,互相关峰值最大点为导频序列精确起始点;
信号细频偏估计及校正模块,用于根据导频序列精确起始点对有效高速信号进行解扰解扩,通过扫频进行精确频偏估计,以对有效高速信号进行精确频偏校正;
信号相偏估计及校正模块,用于根据精确频偏估计计算相位偏移,对有效高速信号进行相偏校正,用于下一步的解调。
进一步地,在信号粗频偏估计及校正模块中,有效高速信号的粗频偏计算公式如下:
Figure BDA0001797154300000051
其中,ΔfINT为粗频偏,Findex为粗频偏索引位置,FFTsize为FFT的大小,fs为采样速率。
本发明与现有技术相比具有以下的优点:
本发明的高速信号的捕获和跟踪方法及装置对于超大频偏单载波信号,能够通过自相关进行粗同步,通过FFT进行粗频偏估计和校正,通过互相关进行精同步,并通过扫频方法进行精频偏估计和校正,还通过精频偏计算并校正相偏,从而以层层递进方式快速捕获并跟踪高速移动物体。
附图说明
以下结合附图和实施例对本发明作进一步说明。
图1是本发明的高速信号的捕获和跟踪方法步骤图;
图2为在SNR=-2,频偏为110KHz时自相关功率图;
图3是在FFTsize=2048,SNR=-2时通过FFT求得的频偏数值图;
图4为本发明获取粗频偏索引位置的一个具体实施方式;
图5为SNR=5时接收的导频信号去除调制信息后FFT显示图;
图6为在SNR=-2时,精确同步的相关峰值功率图;
图7为通过AWGN信道对频偏为11500Hz,扩频为1,信噪比5dB,误比特率为0.0056667的信号解调后的星座图;
图8为在AWGN条件下仿真信号解调效果图;
图9是本发明的高速信号的捕获和跟踪装置框图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语),具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语,应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非被特定定义,否则不会用理想化或过于正式的含义来解释。
图1是本发明的高速信号的捕获和跟踪方法步骤图,如图1所示,本发明提供的高速信号的捕获和跟踪方法,包括以下步骤:
接收高速信号,并对高速信号进行抽取,计算从高速信号中抽取的两个相邻的接收序列的差分相关值,差分相关值最大点为导频序列起始点,以捕获有效高速信号,其中,抽取的两个相邻的接收序列为导频序列的两个对称序列;
接收有效高速信号,对有效高速信号进行FFT处理,以获取有效高速信号的粗频偏,并对该有效高速信号进行粗频偏校准;
对于粗频偏校准后的有效高速信号,从导频序列起始点的相邻预设范围内抽取接收序列,与预存序列进行互相关计算,互相关峰值最大点为导频序列精确起始点;
根据导频序列精确起始点对有效高速信号进行解扰解扩,通过扫频进行精确频偏估计,并计算相位偏移,以对有效高速信号进行精确频偏校正和相偏校正,用于下一步的解调。
本实施例中,高速信号的例子可以如下:
例如移动速度V=3450m/s(10倍音速),加速度a=40g(40倍重力加速度),f0=10GHZ(射频本振),此时的多普勒频移计算如下:
v=3450*1;%%移动速度(m/s)
a=9.8*40;%%加速度(m/s^2)
c=3*1e8;%光速
f0=10*1e9;%%本振
fd=v/c*f0;%%115KHZ
fdd=a/c*f0;%一秒多普勒动态是13K,那么1ms仅仅13HZ,影响很小
fs=40*1e6;%采样功率40M,
ts=1/fs;%%2.5*10^(-8)秒=25ns
IPoint=4;%%M=4内插倍数
freq_pn_Code=fs/IPoint;%伪码速率为10MHz
信道中的功率模型如下:
freq_change=f0*(v+1/2*a*ts)/c;
expf0=exp(1i.*2.*pi.*freq_change./fs.*[0:highN-1]));
该高速信号架构如下:一个时隙内包括一个导频块和多个数据块,组成一个1ms的时隙信号。导频在时间轴上堆成,前后一样。一个时隙信号的多个数据块为10000个QPSK(Quadrature Phase Shift Keyin,正交相移键控)点,一个时隙信号的导频块为512个导频点,导频块包括两个相同的子导频块,每个子导频块为256个导频点。
在本实施例中,接收高速信号时,根据发送信号中频发射特性对等选择接收机处理模式,具体地,若发射机采用重复发送QPSK信号实现内插,接收机则采用矩形滤波器匹配滤波;若发射机采用RRC(Radio Resource Control,无线资源控制)内插滤波,接收机则也采用RRC滤波抽取,采用和发射机对等的中频处理,能够增强接收机处理的灵敏度。
首先,进行对接收的高速信号进行自相关操作,采用单CHIP(码片)一级的自相关处理,采用FIFO(First Input First Output,先入先出队列)处理方式,仅需要一个输入一个输出的功率差,得到总体相关的功率,超过一定门限即认为相关上。具体过程如下:
接收端起初并不知道信号的起始点在哪,利用前导码前后相同的特性进行差分相关。差分相关利用延迟Nd个样值的两个相同训练符号。
差分相关与直接相关法不同,进行差分相关运算的是利用在接收序列rn中的导频序列的对称特性,即时域上形成一个重复的信号,因此将rn
Figure BDA0001797154300000081
进行点乘,并对每个n,将统计窗长Nd,每一个窗内的数据进行求和,统计窗按照逐个采样点移位,相当于在上一个共轭相乘序列上进行滑动,这样,下一个在当前统计窗中只要加上顺延的下一个数据,减去当前窗第一的数据就可以方便的获得下一个统计窗的差分相关值,得到对应的差分相关值cn。该方法差分相关函数c(n)只是与导频序列在接收序列中的位置有关,寻找到cn的最大值,也就找到了最佳定时时刻。
差分相关值如下:
Figure BDA0001797154300000082
其中,cn为差分相关值,rn+k
Figure BDA0001797154300000083
为两个相邻窗长为Nd的接收序列,
Figure BDA0001797154300000084
Figure BDA0001797154300000085
的反转共轭,k为采样序号。判决变量为cn达到最大值的时刻即为最佳定时时刻。
根据上述高速信号架构计算差分相关值,过程如下:
接收信号一时隙(1ms)时间内N=10000采样点(降采样后)加512采样点的数据,n=1…N+512;
接收数据r(n)与他的位移r(n+Nd)接收数据做共轭相乘,获得差分相关序列rx_coru(n)。
rx_coru(n)=r(n)conj(r(n+Nd)
计算统计窗长Nd内差分序列的平均值。
k=0时,
Figure BDA0001797154300000091
且k≠0时,rx_cor_pu(k)=rx_cor_pu(k)+rx_coru(k-1)-rx_coru(k+Nd)
为了克服噪声的影响,可以将多个时隙数据进行平均,这样可以获得更好的检测性能。图2为在SNR=-2,频偏为110KHz时自相关功率图,如图2所示,自相关操作不受频偏的影响,但是最佳功率点前后变动不大,只能估计一定范围内的同步处理,故此适合初始粗同步。
粗同步校准之后,可以启动粗频偏校准。在信号的传输过程中,由于多普勒频移和振荡器的精确度等因素,使得接收信号的载波和本地载波不完全同步,有一定的偏差,并导致相位急剧变化,严重影响了解调器的性能,因此必须在接收机中进行载波同步。输入的中频信号经数字下变频后,变到零频,由于收发载波之间不完全匹配,因此这一零频不是真正意义上的零频,而是收发载波间有一个较小的频差,需要进行载波频偏估计和载波校正,以消除载波偏差的影响。
在实际的通信系统中,接收信号经过下变频变到基带后,总是存在一个小频偏Δf的,载波同步的目的就是要估计出这个Δf值。实际的系统中Δf都是有一定的范围的,如分布在-fmax~+fmax之间,如果在-fmax~+fmax之间对Δf采用一定的算法进行搜索,定义一个似然度,只有当搜索的功率离实际的频偏最近时,似然度值最大,因此便可以得到频偏的估值。
粗同步之后获取时隙的起始位置,由于时隙的切换开关已知,从这个时隙的起始位置,信号进入极坐标,去除调制相位后,开始启动FFT操作的频偏估计,利用功率平滑技术,可以在很低的SNR下精确估计出整数频偏估计,频偏误差颗粒度就是FFT的分辨率。具体过程如下:
采用基于频偏估计与补偿的载波同步算法,该算法首先将频偏估计出来,再对原始信号进行补偿,该操作在解调的过程中始终进行。
经过数字下变频及低通滤波处理后,得到含数据信息和载波频偏的信号,可表示为
Figure BDA0001797154300000101
式中:Δf为残留频偏,Tb为码元周期,n(k)~N(0,2σ2)为加性复高斯白噪声,其同相分量和正交分量的方差均为σ2,θ是初始相位差,Ak是由于噪声影响而得的瞬时幅值,αk是由于噪声影响而附加的相位噪声,φk是被调制的相位,φk∈{0,π}。若不考虑噪声的影响,有
Figure BDA0001797154300000102
M为调制阶数,对于BPSK调制,M=2,对于QPSK调制,M=4。可见调制信息被去掉了。对于本突发系统,初始频偏估计时是利用接收的所有数据(包括前导头和业务数据)进行。
定义r′k
Figure BDA0001797154300000111
n′(k)是一个等效的噪声,可近似认为服从高斯分布。可见r′k是一个被高斯白噪声污染的单载波的形式,此时对此信号进行FFT运算,此时由于调制信息被去掉,剩下的仅剩下频偏信息和噪声。为了进一步增强FFT的抗噪声能力,采用多次功率平滑的FFT叠加算法。
具体地,有效高速信号的粗频偏计算公式如下:
Figure BDA0001797154300000112
其中,
Figure BDA0001797154300000113
为粗频偏,M为调制阶数,FFT(r′k)为有效序列r′k的频谱值。
FFT(r′k)的计算步骤包括:
计算频偏点的频谱图,并对计算的频谱图进行多次功率平均,得到FFT(r′k)。具体如下:
(1)计算频偏点的频谱值|FFT(r′k)|
Figure BDA0001797154300000114
fi的值从-fb/2到fb/2,此时功率的解析带宽(功率分辨率)rbw=fb/FFTsize,例如物理层速率是fb=20.48MHZ,一般FFTsize=2048,则得到rbw=10khz
(2)对计算出来的Fi进行多次功率平均,得到F′i,优选地,N=8次;
F′i=Fi+Fi+1+...+Fi+N
图3是在FFTsize=2048,SNR=-2时通过FFT求得的频偏数值图,如图3所示,通过平滑后在很低的SNR(信噪比)下仍可以准确的估计频偏。
(3)搜索功率最大的功率点,同时由于前面为了去掉功率扩大了M倍,最终结果必须缩小到原先的1/M,这样就找到真正的频偏数值。
Figure BDA0001797154300000121
还有另一种准确寻找粗频偏的方法:发射端发送了导频信息,接收端能够准确找到接收信号导频的起始位置,采用以下方案能够更好更快的估计出粗频偏的位置,具体实施方案如下:
为了节省FFT求平均功率的次数,提供一种FFT求一次平均功率就获取到粗频偏点的方法,可以首先通过对导频信号去除调制信息后,通过一次求取FFT求平均功率即可在很低的SNR下获取粗频偏索引位置,如图4所示,具体包括以下步骤:
第一步,对导频信号去除调制信息后,通过以下公式获取FFT功率图:
fft_pss=FFT(r_pilot.*pss*)
其中,fft_pss为导频信号的频域点,r_pilot.*pss*为导频信号的时域点;
此时去除本地调制信号后的导频信号,功率偏移在FFT后能够非常明显看到,
第二步,根据以下公式找出功率最大的位置点:
[max value,maxpos]=max(|fft_pss|)
其中,[maxvalue,maxpos]为功率最大的位置点,fft_pss为导频信号的频域点;
第三步,根据以下条件判断功率最大的位置点是否是粗频偏索引位置:
{maxvalue/mean(|fft_pss|)}>TH
其中,maxvalue为功率最大的位置点的功率值,mean(|fft_pss|)为导频信号的所有频域点的平均功率值,TH为门限;
若条件成立,则功率最大的位置点为粗频偏索引位置。
实现的具体算法如下:
Figure BDA0001797154300000131
为了避免虚警,相关峰值的绝对功率和相对功率比值都要满足要求,才算准确找到相关峰值点。这一步如果不满足要求,说明峰值功率点没有找对。
图5为SNR=5时接收的导频信号去除调制信息后FFT显示图,如图5所示,此时,相对功率比值PAR达到18,大于检测门限10的要求,说明此时导频寻找正确。
通过索引、FFT的大小以及此时的采样速率就可以得到频偏大小,计算公式如下:
Figure BDA0001797154300000141
其中,ΔfINT为粗频偏,Findex为粗频偏索引位置,FFTsize为FFT的大小,fs为采样速率。
该算法具体实施如下,信号的方向分为正负两个方向,也就是可以计算出频偏的大小有正负。
Figure BDA0001797154300000142
Findex=13,fs=10MHZ,FFTsize=512,则计算法的频偏大小为234.375KHZ
采用此算法的条件是需要准确找到接收信号导频的起始位置。
粗频偏校准之后,可启动精确的相关同步,FFT的功率分辨率就是FFT估计后的残留频偏rbw(解析带宽),rbw*SYN_LEN/fs/2<1/2这个残留频偏不会对本次精确同步造成影响,或者说影响可以忽略。
完成初始频偏测量校准和粗同步之后,可以启动精同步处理,由于小频偏积分对相位的影响减少,不会影响信号的精同步操作,fcoarse=rbw/2=5KHZ
Figure BDA0001797154300000151
Lp=256,rbw=10KHZ,fs=10MHZ,计算得到
Figure BDA0001797154300000152
因此在1次导频互相关周期内相干累加时由于频偏而累积的最大相位为
Figure BDA0001797154300000153
由此而造成的信号能量的损失为
Figure BDA0001797154300000154
精同步处理如下:
互相关计算公式如下:
Figure BDA0001797154300000155
其中,ck为互相关函数,rn+k为接收序列,sn为长度L的预存序列,
Figure BDA0001797154300000156
为sn的反转共轭,k为采样序号,在k=pss_τ的位置时候达到最大值,所以寻找到ck的最大值,也就找到了最佳定时时刻。
根据上述高速信号架构计算互相关值,过程如下:
接收时隙(1ms)时间内10000采样点(降采样后)加512采样点的数据;n=1…N+512
用pss_timeu(n)码与接收数据做相关,每次移动1个采样点的步长,这样得到采样点的相关结果。
计算相关功率,找出最大功率点对应的pss_timeu(n)码,最大功率点位置即为pss_timeu(n)的起始位置。
相关结果:
Figure BDA0001797154300000161
图6为在SNR=-2时,精确同步的相关峰值功率图,如图6所示,精确同步能够精确到具体采样点的位置。
由于相关处理都是m序列,仅仅是1和-1相乘,故此无需乘法运算,只需要加减法运算。并且第一步中粗同步已经确定了大致的导频位置。此时仅仅需要在粗同步位置POS前后偏移M个点来进行相关。POS-M:POS+M来进行相关。此时如果相关峰值功率过低,说明上面粗同步也是一种误判,精同步不仅可以精确定位导频的起始位置,而且还可以反推粗同步的准确性。为了达到小数同步,ADC采集后的信号进行滤波后,不进行抽取。精细同步时虽然是降采样处理,但是移动的点数却是1/M CHIP的速度移动,此时通过精确同步之后,信号精度达到了1/M CHIP的精度。
精确同步之后启动精确频偏估计,首先精确提取出导频信号,这这个导频信号进行解扩,扩频比用户可以自由选择。例如导频码长度512,扩频长度32,则可以解扩出16个数据符号,由于发送的导频解扰之后数值是常数,所以可以避免扩频码边界的选择问题。解扩之后的数据进行扫频法进行精确功率选择,例如进行10次扫频,功率误差缩小到rbw/10。
在精同步和精频偏校准之后,还需要在解调的过程中对频偏进行跟踪,否则由于频偏跟踪误差而累积的相位差就有可能超出相邻两个星座点之间的相位间隔,从而导致解调结果翻转、错误。
精频偏估计采用基于最大似然估计的算法,将数据分成若干连续的段,估计每一段里面的频偏估计值,再对频偏变化率进行修正,从而完成载波同步跟踪。跟踪算法采用并行处理。
M为调制阶数,对于BPSK调制,M=2,对于QPSK调制,M=4。可见调制信息被去掉了。对于本突发系统,初始估计时是利用已知的前导头进行,因此有M=1。
定义r′k
Figure BDA0001797154300000171
n′(k)是一个等效的噪声,可近似认为服从高斯分布。可见r′k是一个被高斯白噪声污染的单载波的形式,由此可以得到r′k功率和初相的最大似然估计表达式为
Figure BDA0001797154300000172
Figure BDA0001797154300000173
其中,
Figure BDA0001797154300000174
为精频偏,
Figure BDA0001797154300000175
为初相偏,M为调制阶数,f为采样点功率,Tb为码元周期,k为采样序号。
综上,本发明的高速信号的捕获和跟踪方法能够准确跟踪超大频偏单载波BPSK(Binary Phase Shift Keying,二进制相移键控)/QPSK(Quadrature Phase Shift Keyin,正交相移键控)/8PSK(8Phase Shift Keying 8,移相键控)信号,分别通过自相关和互相关获取粗同步和精确同步,并分别通过FFT和扫频方法完成粗频偏和精确功率的估计和校正,采用层层递进方式快速跟踪高速移动物体,甚至能够捕获超过20倍音速超高速飞行器的功率偏移,并且准确解调。图7为通过AWGN信道对频偏为11500Hz,扩频为1,信噪比5dB,误比特率为0.0056667的信号解调后的星座图,如图7所示,信号的相位和频偏完全被修正。图8为在AWGN条件下仿真信号解调效果图,如图8所示,即使在低信噪比下,误比特率也很低。
图9是本发明的高速信号的捕获和跟踪装置框图,如图9所示,本发明提供的高速信号的捕获和跟踪装置,包括:
信号粗捕获模块,用于接收高速信号,并对高速信号进行抽取,计算从高速信号中抽取的两个相邻的接收序列的差分相关值,差分相关值最大点为导频序列起始点,以捕获有效高速信号,其中,抽取的两个相邻的接收序列为导频序列的两个对称序列;
信号粗频偏估计及校正模块,用于接收有效高速信号,对有效高速信号进行FFT处理,以获取有效高速信号的粗频偏,并对该有效高速信号进行粗频偏校准;
信号精捕获模块,用于对于粗频偏校准后的有效高速信号,从导频序列起始点的相邻预设范围内抽取接收序列,与预存序列进行互相关计算,互相关峰值最大点为导频序列精确起始点;
信号细频偏估计及校正模块,用于根据导频序列精确起始点对有效高速信号进行解扰解扩,通过扫频进行精确频偏估计,以对有效高速信号进行精确频偏校正;
信号相偏估计及校正模块,用于根据精确频偏估计计算相位偏移,对有效高速信号进行相偏校正,用于下一步的解调。
本发明的高速信号的捕获和跟踪装置对于超大频偏单载波信号,能够通过自相关进行粗同步,通过FFT进行粗频偏估计和校正,通过互相关进行精同步,并通过扫频方法进行精频偏估计和校正,还通过精频偏计算并校正相偏,从而以层层递进方式快速捕获并跟踪高速移动物体。
进一步地,在信号粗捕获模块中,通过以下公式计算差分相关值:
Figure BDA0001797154300000191
其中,cn为差分相关值,rn+k
Figure BDA0001797154300000192
为两个相邻窗长为Nd的接收序列,k为采样序号。
进一步地,在信号粗频偏估计及校正模块中,有效高速信号的粗频偏计算公式如下:
Figure BDA0001797154300000193
其中,
Figure BDA0001797154300000194
为粗频偏,M为调制阶数,FFT(r′k)为有效序列r′k的频谱值。
进一步地,在信号精捕获模块中,互相关计算公式如下:
Figure BDA0001797154300000195
其中,ck为互相关函数,rn+k为接收序列,
Figure BDA0001797154300000196
为预存序列,k为采样序号。
进一步地,在信号细频偏估计及校正模块中,精确频偏估计计算公式如下:
Figure BDA0001797154300000197
其中,
Figure BDA0001797154300000201
为精频偏,M为调制阶数,f为采样点功率,Tb为码元周期,k为采样序号。
对于装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
对于方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明实施例并不受所描述的动作顺序的限制,因为依据本发明实施例,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作并不一定是本发明实施例所必须的。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (6)

1.一种高速信号的捕获和跟踪方法,其特征在于,包括以下步骤:
接收高速信号,并对高速信号进行抽取,计算从高速信号中抽取的两个相邻的接收序列的差分相关值,差分相关值最大点为导频序列起始点,以捕获有效高速信号,其中,抽取的两个相邻的接收序列为导频序列的两个对称序列;
接收有效高速信号,对有效高速信号进行FFT处理,以获取有效高速信号的粗频偏,并对该有效高速信号进行粗频偏校准;
对于粗频偏校准后的有效高速信号,从导频序列起始点的相邻预设范围内抽取接收序列,与预存序列进行互相关计算,互相关峰值最大点为导频序列精确起始点;
根据导频序列精确起始点对有效高速信号进行解扰解扩,通过扫频进行精确频偏估计,并计算相位偏移,以对有效高速信号进行精确频偏校正和相偏校正,用于下一步的解调。
2.根据权利要求1所述的高速信号的捕获和跟踪方法,其特征在于,有效高速信号的粗频偏计算公式如下:
Figure FDA0002786402570000011
其中,ΔfINT为粗频偏,Findex为粗频偏索引位置,FFTsize为FFT的大小,fs为采样速率。
3.根据权利要求2所述的高速信号的捕获和跟踪方法,其特征在于,通过一次求取FFT求平均功率获取粗频偏索引位置,具体包括以下步骤:
第一步,对导频信号去除调制信息后,通过以下公式获取FFT功率图:
fft_pss=FFT(r_pilot.*pss*)
其中,fft_pss为导频信号的频域点,r_pilot.*pss*为导频信号的时域点;
第二步,根据以下公式找出功率最大的位置点:
[max value,max pos]=max(|fft_pss|)
其中,[max value,max pos]为功率最大的位置点,fft_pss为导频信号的频域点;
第三步,根据以下条件判断功率最大的位置点是否是粗频偏索引位置:
{max value/mean(|fft_pss|)}>TH
其中,max value为功率最大的位置点的功率值,mean(|fft_pss|)为导频信号的所有频域点的平均功率值,TH为门限;
若条件成立,则功率最大的位置点为粗频偏索引位置。
4.根据权利要求2所述的高速信号的捕获和跟踪方法,其特征在于,通过以下公式计算差分相关值:
Figure FDA0002786402570000021
其中,cn为差分相关值,rn+k
Figure FDA0002786402570000022
为两个相邻窗长为Nd的接收序列,k为采样序号;L为接收序列长度。
5.一种实现权利要求1所述高速信号的捕获和跟踪方法的装置,其特征在于,包括:
信号粗捕获模块,用于接收高速信号,并对高速信号进行抽取,计算从高速信号中抽取的两个相邻的接收序列的差分相关值,差分相关值最大点为导频序列起始点,以捕获有效高速信号,其中,抽取的两个相邻的接收序列为导频序列的两个对称序列;
信号粗频偏估计及校正模块,用于接收有效高速信号,对有效高速信号进行FFT处理,以获取有效高速信号的粗频偏,并对该有效高速信号进行粗频偏校准;
信号精捕获模块,用于对于粗频偏校准后的有效高速信号,从导频序列起始点的相邻预设范围内抽取接收序列,与预存序列进行互相关计算,互相关峰值最大点为导频序列精确起始点;
信号细频偏估计及校正模块,用于根据导频序列精确起始点对有效高速信号进行解扰解扩,通过扫频进行精确频偏估计,以对有效高速信号进行精确频偏校正;
信号相偏估计及校正模块,用于根据精确频偏估计计算相位偏移,对有效高速信号进行相偏校正,用于下一步的解调。
6.根据权利要求5所述的高速信号的捕获和跟踪装置,其特征在于,在信号粗频偏估计及校正模块中,有效高速信号的粗频偏计算公式如下:
Figure FDA0002786402570000031
其中,ΔfINT为粗频偏,Findex为粗频偏索引位置,FFTsize为FFT的大小,fs为采样速率。
CN201811061247.9A 2018-09-12 2018-09-12 一种高速信号的捕获和跟踪方法及装置 Active CN109104390B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811061247.9A CN109104390B (zh) 2018-09-12 2018-09-12 一种高速信号的捕获和跟踪方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811061247.9A CN109104390B (zh) 2018-09-12 2018-09-12 一种高速信号的捕获和跟踪方法及装置

Publications (2)

Publication Number Publication Date
CN109104390A CN109104390A (zh) 2018-12-28
CN109104390B true CN109104390B (zh) 2021-01-12

Family

ID=64866098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811061247.9A Active CN109104390B (zh) 2018-09-12 2018-09-12 一种高速信号的捕获和跟踪方法及装置

Country Status (1)

Country Link
CN (1) CN109104390B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110113074B (zh) * 2019-05-17 2021-02-26 中国电子科技集团公司第五十四研究所 一种基于fpga的高速伪码并行捕获与跟踪方法
CN110601718B (zh) * 2019-09-05 2022-11-04 珠海市杰理科技股份有限公司 扩频通信捕获方法、装置及系统
CN110868237B (zh) * 2019-11-21 2021-08-24 Oppo广东移动通信有限公司 一种信号同步方法、装置以及计算机存储介质
CN110995632B (zh) * 2019-11-29 2023-03-21 深圳市统先科技股份有限公司 卫星通信带宽复用电路
CN111510409B (zh) * 2020-04-16 2020-12-29 清华大学 利用bpsk数据估计天基机会信号多普勒的方法和系统
CN112039578B (zh) * 2020-09-09 2021-05-04 北京航空航天大学 基于fpga的星载ads-b接收机的信号处理模块
CN112600784B (zh) * 2020-11-23 2022-11-22 中国电子科技集团公司第二十研究所 一种基于二次差分相关的大频偏位同步方法
CN112968719B (zh) * 2021-02-25 2022-08-16 重庆两江卫星移动通信有限公司 一种多分辨率通信信号同步捕获方法、装置、设备及介质
CN113872910B (zh) * 2021-10-11 2024-05-28 南京中科晶上通信技术有限公司 基于fpga的载波同步方法、装置及设备
CN114285709B (zh) * 2021-12-31 2023-04-25 北京中科晶上科技股份有限公司 接收信号的相位跟踪方法、装置和信号处理系统
CN116938657B (zh) * 2023-09-15 2023-12-29 武汉船舶通信研究所(中国船舶集团有限公司第七二二研究所) 一种dsss-oqpsk信号解调方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101394391A (zh) * 2008-11-03 2009-03-25 华北电力大学 基于四维混沌系统的ofdm同步方法
CN102801479A (zh) * 2012-08-14 2012-11-28 中兴通讯股份有限公司 一种扫频频偏估计方法和装置
CN104536016A (zh) * 2014-11-05 2015-04-22 北京大学 一种gnss新体制信号捕获装置及方法
CN104852873A (zh) * 2014-02-18 2015-08-19 中兴通讯股份有限公司 一种频偏估计方法、装置及帧同步子系统
EP3101860A1 (en) * 2015-06-03 2016-12-07 Imagination Technologies Limited Sampling frequency offset calculation
CN105245303B (zh) * 2015-08-28 2017-08-29 北京理工大学 一种高速突发解调同步系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101394391A (zh) * 2008-11-03 2009-03-25 华北电力大学 基于四维混沌系统的ofdm同步方法
CN102801479A (zh) * 2012-08-14 2012-11-28 中兴通讯股份有限公司 一种扫频频偏估计方法和装置
CN104852873A (zh) * 2014-02-18 2015-08-19 中兴通讯股份有限公司 一种频偏估计方法、装置及帧同步子系统
CN104536016A (zh) * 2014-11-05 2015-04-22 北京大学 一种gnss新体制信号捕获装置及方法
EP3101860A1 (en) * 2015-06-03 2016-12-07 Imagination Technologies Limited Sampling frequency offset calculation
CN105245303B (zh) * 2015-08-28 2017-08-29 北京理工大学 一种高速突发解调同步系统

Also Published As

Publication number Publication date
CN109104390A (zh) 2018-12-28

Similar Documents

Publication Publication Date Title
CN109104390B (zh) 一种高速信号的捕获和跟踪方法及装置
EP2016728B1 (en) Delay-doppler channel response demodulation method and apparatus
US7623562B2 (en) Initial synchronization acquiring device and method for parallel processed DS-CDMA UWB system and DS-CDMA system&#39;s receiver using the same
US7684473B2 (en) Receiver for wireless communication network with extended range
US7778363B2 (en) Synchronization timing detecting apparatus, receiving apparatus, and synchronization timing detecting method
EP0892528B1 (en) Carrier recovery for DSSS signals
CN107911329B (zh) 一种信号分析仪ofdm信号解调方法
EP1191704B1 (en) Method and apparatus for acquiring slot timing and frequency offset correction
KR100230029B1 (ko) 코히어런트 수신기에 대한 오프셋 주파수 추정을 위한 방법 및 장치
CN109756968B (zh) 单载波扩频系统的精确同步定时方法及精确同步定时装置
CN106603451B (zh) 一种基于延时自相关的高动态多普勒频偏及频偏变化率估计方法
CN101969321A (zh) 基于fft的直接序列扩频系统的大频偏二次捕获方法
CN105743612A (zh) 实时盲解调大频移短时突发信号的方法
JP2008530951A (ja) 予め符号化された部分応答信号用の復調器および受信器
KR101009304B1 (ko) 속도 응답 시간 추적
CN110278169A (zh) 一种时域频偏估计算法
WO2001013537A1 (en) Doppler spread estimation using channel autocorrelation function hypotheses
US7346098B2 (en) Communication receiver
CN113890591B (zh) 一种低轨星座系统终端载波同步方法及载波同步解调装置
EP1873924A1 (en) Method for determining a pulse position in a signal
GB2524464A (en) Frequency error estimation
US8238408B2 (en) Method for processing a sampled rectified ultra wide band signal
EP1793546B1 (en) Pilot-assisted Doppler frequency estimation
KR101683937B1 (ko) Umts 통신 네트워크용 수신기 내의 주파수 오프셋의 추정을 계산하기 위한 프로세스, umts 통신 시스템용 수신기 및 umts 통신 네트워크용 휴대폰
KR100665100B1 (ko) 디지털 초협대역 단말기용 주파수 옵셋 자동 보상 장치 및그 방법과 그를 이용한 수신 시스템

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
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Room 109-111, 1 / F, 17 / F, Zhongguancun Software Park, 8 Dongbeiwang West Road, Haidian District, Beijing, 100193

Patentee after: Ruixinfeng Aerospace Technology (Beijing) Co.,Ltd.

Country or region after: China

Address before: Room 109-111, 1 / F, 17 / F, Zhongguancun Software Park, 8 Dongbeiwang West Road, Haidian District, Beijing, 100193

Patentee before: BEIJING RINFON TECHNOLOGY Co.,Ltd.

Country or region before: China