CN106059973B - 频偏估计方法和系统 - Google Patents
频偏估计方法和系统 Download PDFInfo
- Publication number
- CN106059973B CN106059973B CN201610600662.1A CN201610600662A CN106059973B CN 106059973 B CN106059973 B CN 106059973B CN 201610600662 A CN201610600662 A CN 201610600662A CN 106059973 B CN106059973 B CN 106059973B
- Authority
- CN
- China
- Prior art keywords
- target point
- baseband signal
- signal
- fft
- point
- 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
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
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2626—Arrangements specific to the transmitter only
- H04L27/2627—Modulators
- H04L27/2628—Inverse Fourier transform modulators, e.g. inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2657—Carrier synchronisation
- H04L27/266—Fine or fractional frequency offset determination and synchronisation
-
- 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)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明涉及一种频偏估计方法和系统,根据基带信号获取基带去调制信号,对基带去调制信号进行快速傅里叶变换,获得离散频谱,根据离散频谱获取频谱幅值最大的第一目标点以及和第一目标点左右相邻的第二目标点和第三目标点,对这三个点进行分段三次样条插值计算,获得第四目标点,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值。本发明使用分段三次样条插值拟合算法,使得对离散频谱的拟合更加贴近真实曲线的形态,因而获得了更真实的逼近效果,对离散频谱的更加逼近的拟合使得插值出的点位置更加精细与准确,而插值出的点位置精度就决定了频偏估计的精度,本发明的频偏估计精度比一般二次拟合插值算法的频偏估计精度高。
Description
技术领域
本发明涉及数字通信技术领域,特别是涉及一种频偏估计方法和系统。
背景技术
在目前的数字通信系统中,收发端载波频偏的存在会使通信系统的性能急剧恶化。数据辅助的频偏估计算法由于其精度高、估计范围广,因而其受到了更多的关注。数据辅助的频偏估计算法主要分为两大类:基于时域分析的估计方法和基于频域分析的估计方法。在频域算法中,基于FFT(Fast Fourier Transformation,快速傅里叶变换)对频谱进行精确插值的频偏估计算法是比较常见的。此种算法频偏估计范围广,对信噪比不敏感,常常受到众多通信系统的青睐。但是由于FFT自身的栅栏效应,FFT的估计精度和有用数据个数以及FFT变换点数有关,因而对于估计精度要求比较高的系统常常难以满足要求。二次抛物线频域内插频偏估计方法虽然使得估计精度有所提高,但此种方法因为只用到二次拟合,所以估计精度提高程度有限,仍然很难满足很多系统的精度要求。
发明内容
基于此,有必要针对现有的二次抛物线频域内插频偏估计方法的估计精度较低的问题,提供一种频偏估计方法和系统。
一种频偏估计方法,包括以下步骤:
获取数字通信系统信道中的基带信号,根据基带信号获取用于去调制的信号,根据用于去调制的信号对基带信号进行去调制操作,获得基带去调制信号;
对基带去调制信号进行快速傅里叶变换,获得离散频谱,根据离散频谱获取频谱幅值最大的第一目标点以及与第一目标点左右相邻的第二目标点和第三目标点;
对第一目标点、第二目标点和第三目标点进行分段三次样条插值计算,获得第四目标点;
获取快速傅里叶变换点数和基带信号的符号周期,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值。
一种频偏估计系统,包括以下单元:
预处理单元,用于获取数字通信系统信道中的基带信号,根据基带信号获取用于去调制的信号,根据用于去调制的信号对基带信号进行去调制操作,获得基带去调制信号;
频谱处理单元,用于对基带去调制信号进行快速傅里叶变换,获得离散频谱,根据离散频谱获取频谱幅值最大的第一目标点以及和第一目标点左右相邻的第二目标点和第三目标点;
插值处理单元,用于对第一目标点、第二目标点和第三目标点进行分段三次样条插值计算,获得第四目标点;
频偏估计单元,用于获取快速傅里叶变换点数和基带信号的符号周期,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值。
根据上述本发明的频偏估计方法和系统,其是根据基带信号获取基带去调制信号,对基带去调制信号进行快速傅里叶变换,获得离散频谱,根据离散频谱获取频谱幅值最大的第一目标点以及和第一目标点左右相邻的第二目标点和第三目标点,对这三个点进行分段三次样条插值计算,获得第四目标点,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值。本发明使用分段三次样条插值拟合算法,使得对离散频谱的拟合更加贴近真实曲线的形态,因而获得了更加真实的逼近效果,对离散频谱的更加逼近的拟合使得插值出的点位置更加精细与准确,而插值出的点位置精度就决定了频偏估计的精度,因而本发明的频偏估计精度比一般二次拟合插值算法的频偏估计精度高。
附图说明
图1为其中一个实施例的频偏估计方法的流程示意图;
图2为其中一个实施例的第一目标点、第二目标点、第三目标点和第四目标点的位置示意图;
图3为其中一个实施例的第一目标点、第二目标点、第三目标点和第四目标点的位置示意图;
图4为其中一个实施例的快速傅里叶变换公式对应的关系曲线图;
图5为其中一个实施例的频偏估计系统的结构示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步的详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不限定本发明的保护范围。
参见图1所示,为本发明一个实施例中频偏估计方法的流程示意图。该实施例中的频偏估计方法,包括以下步骤:
步骤S101:获取数字通信系统信道中的基带信号,根据基带信号获取用于去调制的信号,根据用于去调制的信号对基带信号进行去调制操作,获得基带去调制信号;
本步骤主要是去掉接收的基带信号中的调制信号,便于简化后续的处理过程;
步骤S102:对基带去调制信号进行快速傅里叶变换,获得离散频谱,根据离散频谱获取频谱幅值最大的第一目标点以及和第一目标点左右相邻的第二目标点和第三目标点;
在本步骤中,对基带去调制信号进行快速傅里叶变换,可以得到与基带去调制信号对应的离散频谱,离散频谱中的纵坐标表示频谱幅值,离散频谱中的横坐标表示频谱位置,根据离散频谱中的频谱位置以及幅值可以获得第一目标点、第二目标点和第三目标点;
步骤S103:对第一目标点、第二目标点和第三目标点进行分段三次样条插值计算,获得第四目标点;
在本步骤中,进行分段三次样条插值计算涉及到第一目标点、第二目标点和第三目标点各自的频谱幅值和频谱位置,计算之后可以获得第四目标点,其中也涉及第四目标点的频谱幅值和频谱位置,频谱幅值对应于频谱的纵坐标,频谱位置对应于频谱的横坐标;
步骤S104:获取快速傅里叶变换点数和基带信号的符号周期,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值。
在本实施例中,根据基带信号获取基带去调制信号,对基带去调制信号进行快速傅里叶变换,获得离散频谱,根据离散频谱获取频谱幅值最大的第一目标点以及和第一目标点左右相邻的第二目标点和第三目标点,对这三个点进行分段三次样条插值计算,获得第四目标点,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值。本发明使用分段三次样条插值拟合算法,使得对离散频谱的拟合更加贴近真实曲线的形态,因而获得了更加真实的逼近效果,对离散频谱的更加逼近的拟合使得插值出的点位置更加精细与准确,而插值出的点位置精度就决定了频偏估计的精度,因而本发明的频偏估计精度比一般二次拟合插值算法的频偏估计精度高。根据精度高的频偏估计值对数字通信系统进行补偿,可以提高通信系统的性能稳定性。
在其中一个实施例中,快速傅里叶变换点数大于基带信号中符号的个数,在对基带去调制信号进行快速傅里叶变换之前,在基带去调制信号的尾部补上预设个数的零,预设个数为快速傅里叶变换点数与基带信号中符号的个数的差值。
在本实施例中,快速傅里叶变换点数一般为2的乘方数,而基带信号中符号的个数不一定恰好等于2的乘方数,快速傅里叶变换点数理论上是要求大于等于待变换的符号的个数的,快速傅里叶变换点数大于基带信号中符号的个数,可以保证能对基带信号中所有符号进行快速傅里叶变换,在基带去调制信号的尾部补上预设个数的零,预设个数为快速傅里叶变换点数与基带信号中符号的个数的差值,可以使信号的符号个数与快速傅里叶变换点数一致,补零就是快速傅里叶变换本身的需要,便于统一进行快速傅里叶变换。
在其中一个实施例中,根据基带信号获取用于去调制的信号的步骤包括以下步骤:
获取基带信号中的调制信号,对调制信号进行调制映射和复共轭操作,获得用于去调制的信号。
在本实施例中,接收的基带信号中包括调制信号,该调制信号主要是用于同步的,同步包含数据的捕获,载波频率和相位估计以及码元位同步,这里主要涉及的是频率,调制信号一般表示为复数形式,因此需要对其做复共轭操作,把复数形式转化为实数形式,获得用于去调制的信号,用于去调制的信号为后续的去调制操作提供依据。
在其中一个实施例中,根据用于去调制的信号对基带信号进行去调制操作的步骤包括以下步骤:
将用于去调制的信号与基带信号进行逐点相乘。
在本实施例中,基带信号和用于去调制的信号都是一个信号序列,由于基带信号的每一个信号中都包括了调制信号,而用于去调制的信号又是调制信号的复共轭形式,因此,将用于去调制的信号与基带信号进行逐点相乘后,就可以去掉基带信号中的调制信号,完成去调制操作。
在其中一个实施例中,对第一目标点、第二目标点和第三目标点进行分段三次样条插值计算,获得第四目标点的步骤包括以下步骤:
判断第二目标点的频谱幅值是否大于第三目标点的频谱幅值,若是,如图2所示,则根据进行分段三次样条插值计算;
式中,A(km)、A(km-1)和A(km+1)分别表示第一目标点的频谱幅值、第二目标点的频谱幅值和第三目标点的频谱幅值,km、km-1、km+1和kmax分别表示第一目标点的频谱位置、第二目标点的频谱位置、第三目标点的频谱位置和第四目标点的频谱位置。
在本实施例中,是在第二目标点的频谱幅值大于第三目标点的频谱幅值的情况下使用三次样条插值拟合算法,使得对离散频谱的拟合更加贴近真实曲线的形态,获得更加真实的逼近效果,而且三次样条插值拟合算法的公式仅涉及加减乘除和开平方运算,计算简单,因此,相比于其他插值算法,三次样条插值拟合算法的计算复杂度并未增加。
在其中一个实施例中,对第一目标点、第二目标点和第三目标点进行分段三次样条插值计算,获得第四目标点的步骤包括以下步骤:
判断第二目标点的频谱幅值是否大于第三目标点的频谱幅值,若否,如图3所示,则根据进行分段三次样条插值计算;
式中,A(km)、A(km-1)和A(km+1)分别表示第一目标点的频谱幅值、第二目标点的频谱幅值和第三目标点的频谱幅值,km、km-1、km+1和kmax分别表示第一目标点的频谱位置、第二目标点的频谱位置、第三目标点的频谱位置和第四目标点的频谱位置。
在本实施例中,是在第二目标点的频谱幅值不大于第三目标点的频谱幅值的情况下使用三次样条插值拟合算法,使得对离散频谱的拟合更加贴近真实曲线的形态,获得更加真实的逼近效果,而且三次样条插值拟合算法的公式仅涉及加减乘除和开平方运算,计算简单,因此,相比于其他插值算法,三次样条插值拟合算法的计算复杂度并未增加。
在上述两个实施例中,插值公式是通过分段三次样条曲线拟合计算得出来的。令x0=km-1,x1=km,x2=km+1,对应的令y0=A(km-1),y1=A(km),y2=A(km+1),则令样条曲线S(x)是一个分段定义的光滑曲线,其分段函数曲线我们用三次多项式进行拟合,其两个分段函数可写为:
Si(x)=ai+bi(x-xi)+ci(x-xi)2+di(x-xi)3,i=0,1
其中a0,a1,b0,b1,c0,c1,d0,d1代表8个未知系数。
通过已知条件可以求解出8个未知数,求解的结果为:
得到了8个未知系数后,便得到了两段区间的光滑曲线函数表达式,根据三次函数的特性,可以求解出两个区间段内的最大值点,也即
求解S'0(x)=0表达式得到:
同时要求:
求解S1'(x)=0表达式得到:
同时要求:
在其中一个实施例中,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值的步骤包括以下步骤:
根据获取频偏估计值,式中,表示频偏估计值,kmax表示第四目标点的频谱位置,Nfft表示快速傅里叶变换点数,T表示基带信号的符号周期。
在本实施例中,根据公式可以获得频偏估计值,由于快速傅里叶变换点数和基带信号的符号周期一般都不会有偏差,因此频偏估计的精度与第四目标点的频谱位置,也就是第四目标点的位置精度有关;对离散频谱的更加逼近的拟合使得插值出的点位置更加精细与准确,而插值出的点位置精度就决定了频偏估计的精度,因而频偏估计精度比一般插值算法的频偏估计精度高。
在其中一个实施例中,接收基带信号的物理传输通道为加性高斯白噪声信道,基带信号的表达式为:
rk=sk×ej(2πΔfkT+φ)+nk,k=1,2,...,N
式中,rk表示基带信号,sk为调制信号的复数形式,Δf为基带信号的频率偏差,T表示基带信号的符号周期,N表示基带信号中符号的个数,Φ表示相位偏差,nk为加性复高斯白噪声信号;
基带去调制信号的快速傅里叶变换的公式为:
式中,A(f)表示基带去调制信号进行快速傅里叶变换后的频谱幅值,Nfft表示快速傅里叶变换点数,sk *是sk的复共轭形式,f表示频率。
在本实施例中,基带信号中还包括相位偏差和加性复高斯白噪声,根据上述公式进行计算时会存在偏差,得到的频谱中的最大幅值点会有所偏差,因此才需要使用进行分段三次样条插值计算,使得对离散频谱的拟合更加贴近真实曲线的形态,获得了更加真实的逼近效果。
优选的,在进行快速傅里叶变换时,Nfft大于N,在基带去调制信号的快速傅里叶变换的公式中,当k的取值大于N时,rk为零,sk*为零;
基带去调制信号的快速傅里叶变换的公式所对应的曲线如图4所示,图中, 表示真正的频偏值,Δf为估计的频偏值(即估计的基带信号的频率偏差)。
在离散频谱中寻找幅值最大的点的公式为:
A(km)=max{A(k)},k=1,2,...,Nfft
由上式可知,在离散频谱中各点的频谱位置与快速傅里叶变换点相对应。
本发明的基本思路是搜索去调制信号的离散频谱图中的峰值,理论上去调制信号的离散频谱图的峰值位置为零频位置,当信号带有频偏时,峰值位置便为频偏值的位置,因而峰值位置即对应于载波频偏。FFT作为频域频偏估计算法,其不但频偏估计范围大,可以达到±50%的符号速率,而且还能够在很低的信噪比下工作。此外由于使用了三次样条插值你拟合算法,实现复杂度低,在实际系统中可以得到广泛应用。
根据上述频偏估计方法,本发明还提供一种频偏估计系统,以下就本发明的频偏估计系统的实施例进行详细说明。
参见图5所示,为本发明一个实施例中频偏估计系统的结构示意图。该实施例中的频偏估计系统包括以下单元:
预处理单元210,用于获取数字通信系统信道中的基带信号,根据基带信号获取用于去调制的信号,根据用于去调制的信号对基带信号进行去调制操作,获得基带去调制信号;
频谱处理单元220,用于对基带去调制信号进行快速傅里叶变换,获得离散频谱,根据离散频谱获取频谱幅值最大的第一目标点以及和第一目标点左右相邻的第二目标点和第三目标点;
插值处理单元230,用于对第一目标点、第二目标点和第三目标点进行分段三次样条插值计算,获得第四目标点;
频偏估计单元240,用于获取快速傅里叶变换点数和基带信号的符号周期,根据第四目标点、快速傅里叶变换点数和基带信号的符号周期获取频偏估计值。
在其中一个实施例中,快速傅里叶变换点数大于基带信号中符号的个数,频谱处理单元240还用于在对基带去调制信号进行快速傅里叶变换之前,在基带去调制信号的尾部补上预设个数的零,预设个数为快速傅里叶变换点数与基带信号中符号的个数的差值。
在其中一个实施例中,预处理单元210获取基带信号中的调制信号,对调制信号进行调制映射和复共轭操作,获得用于去调制的信号。
在其中一个实施例中,预处理单元210将用于去调制的信号与基带信号进行逐点相乘,获得基带去调制信号。
在其中一个实施例中,插值处理单元230判断第二目标点的频谱幅值是否大于第三目标点的频谱幅值,若是,则根据
进行分段三次样条插值计算;
式中,A(km)、A(km-1)和A(km+1)分别表示第一目标点的频谱幅值、第二目标点的频谱幅值和第三目标点的频谱幅值,km、km-1、km+1和kmax分别表示第一目标点的频谱位置、第二目标点的频谱位置、第三目标点的频谱位置和第四目标点的频谱位置。
在其中一个实施例中,插值处理单元230判断第二目标点的频谱幅值是否大于第三目标点的频谱幅值,若否,则根据
进行分段三次样条插值计算;
式中,A(km)、A(km-1)和A(km+1)分别表示第一目标点的频谱幅值、第二目标点的频谱幅值和第三目标点的频谱幅值,km、km-1、km+1和kmax分别表示第一目标点的频谱位置、第二目标点的频谱位置、第三目标点的频谱位置和第四目标点的频谱位置。
在其中一个实施例中,频偏估计单元240根据获取频偏估计值,式中,表示频偏估计值,kmax表示第四目标点的频谱位置,Nfft表示快速傅里叶变换点数,T表示基带信号的符号周期。
在其中一个实施例中,接收基带信号的物理传输通道为加性高斯白噪声信道,基带信号的表达式为:
rk=sk×ej(2πΔfkT+φ)+nk,k=1,2,...,N
式中,rk表示基带信号,sk为调制信号的复数形式,Δf为基带信号的频率偏差,T表示基带信号的符号周期,N表示基带信号中符号的个数,Φ表示相位偏差,nk表示加性复高斯白噪声信号;
基带去调制信号的快速傅里叶变换的公式为:
式中,A(f)表示基带去调制信号进行快速傅里叶变换后的频谱幅值,Nfft表示快速傅里叶变换点数,sk *是sk的复共轭形式,f表示频率。
本发明的频偏估计系统与本发明的频偏估计方法一一对应,在上述频偏估计方法的实施例阐述的技术特征及其有益效果均适用于频偏估计系统的实施例中。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种频偏估计方法,其特征在于,包括以下步骤:
获取数字通信系统信道中的基带信号,根据所述基带信号获取用于去调制的信号,根据所述用于去调制的信号对所述基带信号进行去调制操作,获得基带去调制信号;
对所述基带去调制信号进行快速傅里叶变换,获得离散频谱,根据所述离散频谱获取频谱幅值最大的第一目标点以及与所述第一目标点左右相邻的第二目标点和第三目标点;
对所述第一目标点、所述第二目标点和所述第三目标点进行分段三次样条插值计算,获得第四目标点;
获取快速傅里叶变换点数和所述基带信号的符号周期,根据所述第四目标点、所述快速傅里叶变换点数和所述基带信号的符号周期获取频偏估计值;
所述对所述第一目标点、所述第二目标点和所述第三目标点进行分段三次样条插值计算,获得第四目标点的步骤包括以下步骤:
判断所述第二目标点的频谱幅值是否大于所述第三目标点的频谱幅值,若是,则根据进行分段三次样条插值计算;若否,则根据进行分段三次样条插值计算;
式中,A(km)、A(km-1)和A(km+1)分别表示所述第一目标点的频谱幅值、所述第二目标点的频谱幅值和所述第三目标点的频谱幅值,km、km-1、km+1和kmax分别表示所述第一目标点的频谱位置、所述第二目标点的频谱位置、所述第三目标点的频谱位置和所述第四目标点的频谱位置。
2.根据权利要求1所述的频偏估计方法,其特征在于,所述快速傅里叶变换点数大于所述基带信号中符号的个数,在对所述基带去调制信号进行快速傅里叶变换之前,在所述基带去调制信号的尾部补上预设个数的零,所述预设个数为所述快速傅里叶变换点数与所述基带信号中符号的个数的差值。
3.根据权利要求1所述的频偏估计方法,其特征在于,所述根据所述基带信号获取用于去调制的信号的步骤包括以下步骤:
获取所述基带信号中的调制信号,对所述调制信号进行调制映射和复共轭操作,获得所述用于去调制的信号。
4.根据权利要求1所述的频偏估计方法,其特征在于,所述根据所述用于去调制的信号对所述基带信号进行去调制操作的步骤包括以下步骤:
将所述用于去调制的信号与所述基带信号进行逐点相乘。
5.根据权利要求1所述的频偏估计方法,其特征在于,所述根据所述第四目标点、快速傅里叶变换点数和所述基带信号的符号周期获取频偏估计值的步骤包括以下步骤:
根据获取频偏估计值,式中,表示所述频偏估计值,kmax表示所述第四目标点的频谱位置,Nfft表示快速傅里叶变换点数,T表示所述基带信号的符号周期。
6.根据权利要求3所述的频偏估计方法,其特征在于,所述数字通信系统信道为加性高斯白噪声信道,所述基带信号的表达式为:
rk=sk×ej(2πΔfkT+φ)+nk,k=1,2,...,N
式中,rk表示基带信号,sk为所述调制信号的复数形式,Δf为基带信号的频率偏差,T表示所述基带信号的符号周期,N表示所述基带信号中符号的个数,Φ表示相位偏差,nk表示加性高斯白噪声信号;
所述基带去调制信号的快速傅里叶变换的公式为:
式中,A(f)表示所述基带去调制信号进行快速傅里叶变换后的频谱幅值,Nfft表示快速傅里叶变换点数,sk *是sk的复共轭形式,f表示频率。
7.一种频偏估计系统,其特征在于,包括以下单元:
预处理单元,用于获取数字通信系统信道中的基带信号,根据所述基带信号获取用于去调制的信号,根据所述用于去调制的信号对所述基带信号进行去调制操作,获得基带去调制信号;
频谱处理单元,用于对所述基带去调制信号进行快速傅里叶变换,获得离散频谱,根据所述离散频谱获取频谱幅值最大的第一目标点以及与所述第一目标点左右相邻的第二目标点和第三目标点;
插值处理单元,用于对所述第一目标点、所述第二目标点和所述第三目标点进行分段三次样条插值计算,获得第四目标点;
频偏估计单元,用于获取快速傅里叶变换点数和所述基带信号的符号周期,根据所述第四目标点、所述快速傅里叶变换点数和所述基带信号的符号周期获取频偏估计值;
插值处理单元判断所述第二目标点的频谱幅值是否大于所述第三目标点的频谱幅值,若是,则根据进行分段三次样条插值计算;若否,则根据进行分段三次样条插值计算;
式中,A(km)、A(km-1)和A(km+1)分别表示所述第一目标点的频谱幅值、所述第二目标点的频谱幅值和所述第三目标点的频谱幅值,km、km-1、km+1和kmax分别表示所述第一目标点的频谱位置、所述第二目标点的频谱位置、所述第三目标点的频谱位置和所述第四目标点的频谱位置。
8.根据权利要求7所述的频偏估计系统,其特征在于,所述快速傅里叶变换点数大于所述基带信号中符号的个数,所述频谱处理单元还用于在对所述基带去调制信号进行快速傅里叶变换之前,在所述基带去调制信号的尾部补上预设个数的零,所述预设个数为所述快速傅里叶变换点数与所述基带信号中符号的个数的差值。
9.根据权利要求7所述的频偏估计系统,其特征在于,预处理单元获取所述基带信号中的调制信号,对所述调制信号进行调制映射和复共轭操作,获得所述用于去调制的信号。
10.根据权利要求7所述的频偏估计系统,其特征在于,预处理单元将所述用于去调制的信号与所述基带信号进行逐点相乘,获得基带去调制信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610600662.1A CN106059973B (zh) | 2016-07-26 | 2016-07-26 | 频偏估计方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610600662.1A CN106059973B (zh) | 2016-07-26 | 2016-07-26 | 频偏估计方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106059973A CN106059973A (zh) | 2016-10-26 |
CN106059973B true CN106059973B (zh) | 2019-10-18 |
Family
ID=57418661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610600662.1A Active CN106059973B (zh) | 2016-07-26 | 2016-07-26 | 频偏估计方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106059973B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107395287B (zh) * | 2016-05-16 | 2019-12-13 | 深圳市中兴微电子技术有限公司 | 一种频偏估计方法和装置 |
CN108881087A (zh) * | 2018-09-13 | 2018-11-23 | 哈尔滨工业大学(深圳) | 一种用于k–k接收的频偏估计方法 |
CN109067680B (zh) * | 2018-09-19 | 2021-09-07 | 深圳市鼎阳科技股份有限公司 | 一种基带信号的载波频偏估计方法及其装置 |
CN111122402B (zh) * | 2019-12-10 | 2022-06-07 | 北京蛙鸣华清环保科技有限公司 | 基于离散监测点数据的污染物路况图生成方法和系统 |
CN112444788B (zh) * | 2020-11-09 | 2023-11-24 | 厦门大学 | 补零信号的频率估计方法 |
CN115913851A (zh) * | 2022-12-06 | 2023-04-04 | 南宁师范大学 | 一种基于三次样条插值的载波相位估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101702701A (zh) * | 2009-11-20 | 2010-05-05 | 西安电子科技大学 | 极低信噪比下频率偏移的估计与补偿方法 |
CN102195905A (zh) * | 2011-03-31 | 2011-09-21 | 重庆重邮信科通信技术有限公司 | 一种频偏估计方法及装置 |
CN103023831A (zh) * | 2012-12-19 | 2013-04-03 | 中国船舶重工集团公司第七二二研究所 | 一种适用于突发波形的载波频偏估计方法 |
-
2016
- 2016-07-26 CN CN201610600662.1A patent/CN106059973B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101702701A (zh) * | 2009-11-20 | 2010-05-05 | 西安电子科技大学 | 极低信噪比下频率偏移的估计与补偿方法 |
CN102195905A (zh) * | 2011-03-31 | 2011-09-21 | 重庆重邮信科通信技术有限公司 | 一种频偏估计方法及装置 |
CN103023831A (zh) * | 2012-12-19 | 2013-04-03 | 中国船舶重工集团公司第七二二研究所 | 一种适用于突发波形的载波频偏估计方法 |
Non-Patent Citations (2)
Title |
---|
基于傅里叶系数插值的短波信道迭代频偏估计;游行远;《华南理工大学学报(自然科学版)》;20150131;第2页左栏 * |
突发数据传输中的同步技术研究;张媛;《中国优秀硕士论文全文数据库》;20120731;正文第4.4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN106059973A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106059973B (zh) | 频偏估计方法和系统 | |
CN103856431B (zh) | 基于快速傅里叶变换的mpsk调制下开环载波同步算法 | |
WO2014132599A1 (ja) | 受信装置、位相誤差推定方法、及び位相誤差補正方法 | |
CN106546817B (zh) | 一种具有反馈功能的频率估计和能量估计方法 | |
JP3271504B2 (ja) | 周波数推定回路およびそれを用いたafc回路 | |
CN108710027B (zh) | 通道间相位差、幅度差的高精度测量方法 | |
CN104363194A (zh) | 基于波形变换的psk调制识别方法 | |
CN109067680B (zh) | 一种基带信号的载波频偏估计方法及其装置 | |
CN107623647B (zh) | 一种基于离散导频辅助的载波同步方法 | |
CN111131098B (zh) | 基于水声ofdm信号自相关输出峰插值的多普勒因子估计方法 | |
CN103152294B (zh) | 基于信号消除进行噪声估计的方法及系统 | |
CN101242390B (zh) | 基于已知序列相干自相关的载波频偏估计算法及其实现装置 | |
WO2021143644A1 (zh) | 一种载波相位跟踪方法及装置 | |
US20140064418A1 (en) | Dc removal from multicarrier signals | |
CN109660481B (zh) | 一种正交幅度调制解调方法及其星座图优化方法 | |
CN102315835B (zh) | 一种成形滤波器滚降系数估计方法 | |
CN101969416A (zh) | 一种信道系数估计方法及装置 | |
CN111539323A (zh) | 一种循环前缀线性调频信号的频率估计方法与装置 | |
CN105406921A (zh) | 一种mfsk水声通信方法及系统 | |
CN114866393A (zh) | 基于前导序列的信令检测方法及装置 | |
Li et al. | An improved DFT frequency offset estimation algorithm with high accuracy | |
CN110247864A (zh) | 一种sc-fde系统的无线时变信道估计方法 | |
CN101197802A (zh) | 正交频分多路信号载波的大频率偏差估计与纠正方法 | |
CN103248591A (zh) | 基于频谱重心的粗频偏估计方法 | |
CN113824666A (zh) | 一种卫星移动ofdm通信系统信道估计方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |