CN110068727A - 一种基于Candan-Rife综合内插的单频信号频率估计方法 - Google Patents
一种基于Candan-Rife综合内插的单频信号频率估计方法 Download PDFInfo
- Publication number
- CN110068727A CN110068727A CN201910279608.5A CN201910279608A CN110068727A CN 110068727 A CN110068727 A CN 110068727A CN 201910279608 A CN201910279608 A CN 201910279608A CN 110068727 A CN110068727 A CN 110068727A
- Authority
- CN
- China
- Prior art keywords
- frequency
- discrete
- interpolation
- data sequence
- candan
- 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 74
- 238000001228 spectrum Methods 0.000 claims abstract description 54
- 238000004364 calculation method Methods 0.000 claims abstract description 3
- 238000005070 sampling Methods 0.000 claims description 32
- 238000001514 detection method Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims 3
- 238000003786 synthesis reaction Methods 0.000 claims 3
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 7
- 238000012937 correction Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mathematical Physics (AREA)
- Complex Calculations (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于Candan‑Rife综合内插的单频信号频率估计方法,包括:获取待处理的单频信号采样数据序列,计算数据序列的离散傅里叶变换和幅度谱;搜索数据序列的幅度谱最大值所对应的离散频率索引;利用数据序列的离散傅里叶变换,计算Candan插值相对频率偏差;利用数据序列的幅度谱,计算Rife插值相对频率偏差;对离散频率索引左、右相邻两离散频率索引所对应幅度谱的比较结果,结合Rife和Candan插值相对频率偏差,估计得到综合相对频率偏差;根据综合相对频率偏差和离散频率索引,利用插值公式估计出单频信号的频率。本发明在不同频段采用不同的估计方法,可在不增加计算量的前提下,提高信号频率估计的精度和稳健性,适合对单频信号频率实时、高精度估计。
Description
技术领域
本发明涉及一种基于Candan-Rife综合内插的单频信号频率估计方法,属于信号处理技术领域。
背景技术
由于单频信号生成简单,且具有较高的速度分辨率,在无线通讯、雷达和声纳等领域得到了广泛的应用,对被噪声污染的单频信号的频率进行精确估计一直是信号处理研究的热点之一。
长期以来,国内外学者对单频信号的频率估计进行了深入、系统的研究,提出了很多估计算法,归纳起来主要包括:(1)基于参数模型的现代谱估计;(2)最大似然估计;(3)基于时域相位估计和(4)基于离散傅里叶变换(Discrete Fourier Transform,DFT)频谱校正的方法。其中基于现代谱估计的方法具有频率分辨率高且不依赖于数据长度等优点,但该类方法对谱模型误差敏感,且运算量大;最大似然法(Maximum Likelihood Estimation,MLE)方差性能最好,接近克拉美-罗下界(Cramer-Rao Low Bound,CRLB),但因MLE是一种非线性最小二乘拟合问题,需要多维搜索才能得到频率估计结果,其运算量太大,难于实现实时处理;基于时域相位的方法如线性回归频率估计器(Linear Prediction,LP)和通用加权线性预测估计器(Generalized Weighted Linear Predictor,GWLP),在信噪比足够大的情况下可以达到CRLB,运算量小但对信噪比要求高,存在频率估计范围小且估计性能不一致等问题。
基于DFT频谱校正的方法,物理意义明确,可以利用快速傅里叶变换(FFT)快速实现、实时性好,而且具有较高的信噪比增益和对算法参数不敏感等优点,是一个综合性能最佳的方法,因此得到了广泛的应用和研究。基于DFT频谱校正的方法目前主要有能量重心法、FFT+FT谱连续细化分析傅立叶变换法和插值法有三种校正方法。其中插值法由于具有实现简单,运算量小等优点,得到的研究最为广泛。插值法主要包括两点插值法和三点插值法,两点插值法的典型代表是Rife插值法,该算法利用单频信号幅度谱主瓣内的两个谱线的对相对频率偏差进行计算,从而估计出信号频率;三点内插法的典型代表是Candan插值法,该算法利用单频信号离散傅里叶变换主瓣附近的三根谱线计算相对频率偏差,从而估计出信号频率。分析表明:当单频信号相对频率偏差较大时,Rife插值法估计精度较高,而当相对频率偏差较小时,Rife插值法估计精度较低;与之相反,当相对频率偏差较小时,Candan插值法估计精度高,而相对频率偏差较大时,Candan插值法估计精度低,也就是说当前算法的估计精度对相对频率偏差有较高的依赖性,稳健性差。
发明内容
本发明所要解决的技术问题在于克服现有技术的不足,提供一种基于Candan-Rife综合内插的单频信号频率估计方法,该方法可在不增加运算量的条件下,频率估计精度较常规Candan和Rife直接插值法对相对频率偏差的依赖性大大降低,稳健性高,可以满足雷达和声纳信号处理的精度要求和实时性需求,工程实用性强。
本发明具体采用以下技术方案解决上述技术问题:
一种基于Candan-Rife综合内插的单频信号频率估计方法,包括以下步骤:
步骤(1)、获取待处理的单频信号采样数据序列x(n),n=0,1,…,N-1,其中N为检测到的单频信号脉宽长度所对应的采样点个数,且N取值为2的整数次幂及N≥4;
步骤(2)、根据单频信号采样数据序列x(n)计算得到数据序列的离散傅里叶变换X(l)和幅度谱P(k);
步骤(3)、搜索数据序列的幅度谱P(k)最大值所对应的离散频率索引kp;
步骤(4)、利用数据序列的离散傅里叶变换X(l)计算Candan插值相对频率偏差
步骤(5)、利用数据序列的幅度谱P(k)计算Rife插值相对频率偏差
步骤(6)通过对离散频率索引kp左、右相邻两离散频率索引所对应幅度谱的比较结果,结合相对频率偏差和估计得到综合相对频率偏差
步骤(7)、根据综合相对频率偏差和离散频率索引kp,利用插值公式估计出单频信号的频率
进一步地,作为本发明的一种优选技术方案,所述步骤(1)获取待处理的单频信号采样数据序列x(n),具体为:
从传感器接收N个采样点的实时采集数据作为待处理的单频信号采样数据序列x(n);或从存储器中提取从检测到信号时刻起始的N个采样点的数据作为待处理的单频信号采样数据序列x(n)。
进一步地,作为本发明的一种优选技术方案,所述步骤(2)计算得到数据序列的离散傅里叶变换X(l)和幅度谱P(k),具体为:
步骤(21)、计算单频信号采样数据序列x(n)的离散傅里叶变换X(l):
其中,l为X(l)的离散频率索引,j表示虚数单位,即
步骤(22)、依据离散傅里叶变换X(l)计算单频信号采样数据序列x(n)的脉冲信号幅度谱P(k):
且k=0,1,2…,N/2
其中,k为P(k)的离散频率索引,X(k)为x(n)的离散傅里叶变换,||代表取模值运算。
进一步地,作为本发明的一种优选技术方案,所述步骤(3)搜索数据序列的幅度谱P(k)最大值所对应的离散频率索引kp,采用公式:
其中,表示在1≤k≤N/2-1范围内搜索幅度谱P(k)的最大值所对应的离散频率索引。
进一步地,作为本发明的一种优选技术方案,所述步骤(4)计算Candan插值相对频率偏差采用公式:
其中,Re{}代表取实部运算;X(kp)代表离散频率索引kp对应的离散傅里叶变换;X(kp-1)代表离散频率索引kp-1对应的离散傅里叶变换;X(kp+1)代表离散频率索引kp+1对应的离散傅里叶变换。
进一步地,作为本发明的一种优选技术方案,所述步骤(5)计算Rife插值相对频率偏差采用公式:
其中,P(kp)代表离散频率索引kp对应的幅度谱;P(kp-1)代表离散频率索引kp-1对应的幅度谱;P(kp+1)代表离散频率索引kp+1对应的幅度谱。
进一步地,作为本发明的一种优选技术方案,所述步骤(6)估计得到综合相对频率偏差具体为:
步骤(61)、提取离散频率索引kp左、右相邻两离散频率索引kp-1和kp+1所对应的幅度谱结果Al和Ar:
Al=P(kp-1) Ar=P(kp+1)
步骤(62)、分别提取Al和Ar中的最大值和最小值Amax和Amin:
Amax=max(Al,AR) Amin=min(Al,AR)
其中,max()代表取大值运算,min()代表取小值运算;
步骤(63)、根据相对频率偏差和的最大值和最小值Amax及Amin,估计得到综合相对频率偏差
其中,η为设定的特征门限,η为大于0且小于1的数。
进一步地,作为本发明的一种优选技术方案,所述步骤(7)中利用插值公式估计出单频信号的频率具体为:
其中,Δf为长度为N的离散傅里叶变换的频率分辨率,Δf=fs/N,fs为采样频率
本发明采用上述技术方案,能产生如下技术效果:
本发明的基于Candan-Rife综合内插的单频信号频率估计方法,与现有的方法相比,能产生以下有益效果:
(1)本发明的估计方法通过对相对频率偏差分布范围进行预估,在不同的频段采用不同的估计方法:由于当相对频率偏差δ较小时,Candan插值估计精度优于Rife插值,而当δ较小时,Candan插值估计精度劣于Rife插值,因此仅在δ较小时采用Candan插值法,比对所有δ均利用Candan插值法对δ的稳健性高;
(2)本发明的估计方法通过对相对频率偏差分布范围进行预估,在不同的频段采用不同的估计方法:当相对频率偏差δ较大时,信号DFT主瓣内两根谱线的幅值均较大,抗噪声能力强,且主瓣内次大谱线的幅值远大于旁瓣谱线的幅值,利用Rife插值法出现插值方向错误的概率很小,因此仅在δ较大时采用Rife插值法,比对所有δ均利用Rife插值法对δ的稳健性高;
(3)本发明的估计方法较原Candan插值和Rife插值,可在不增加计算量的前提下,提高频率估计的精度和稳健性,适用于对实时性和频率估计精度要求高的场合。
附图说明
图1为本发明单频信号频率估计方法的流程示意图。
图2为本发明中实施例1的仿真单频信号幅度谱图。
图3为本发明中实施例2的仿真单频信号幅度谱图。
具体实施方式
下面结合说明书附图对本发明的实施方式进行描述。
如图1所示,本发明设计了一种基于Candan-Rife综合内插的单频信号频率估计方法,该方法具体包括如下步骤:
步骤(1)、获取待处理的单频信号采样数据序列x(n),n=0,1,…,N-1;包括:从传感器接收N个采样点的实时采集数据作为待处理的单频信号采样数据序列x(n),n=0,1,…,N-1;或从存储器中提取从检测到信号时刻起始的N个采样点的数据作为待处理的单频信号采样数据序列x(n),n=0,1,…,N-1,所述的N为检测到的单频信号脉宽长度所对应的采样点个数,取值为2的整数次幂,且N≥4。
步骤(2)、对所述单频信号采样数据序列x(n)做离散傅里叶变换,计算得到数据序列的离散傅里叶变换X(l)和幅度谱P(k),具体包括如下两步:
步骤(21)、计算单频信号采样数据序列x(n)的离散傅里叶变换X(l):
其中,l为离散傅里叶变换X(l)的离散频率索引,j表示虚数单位,即该式通过快速傅里叶变换实现;
步骤(22)、依据离散傅里叶变换X(l)计算x(n)的幅度谱P(k):
且k=0,1,2…,N/2 (2)
其中k为P(k)的离散频率索引,X(k)为x(n)的离散傅里叶变换,||代表取模值运算。
在步骤(22)步中,数据序列x(n)的离散傅里叶变换X(l)即公式(1),是通过快速傅里叶变换实现的,利用快速傅里叶变换可以降低算法的运算量,提高算法的计算效率;公式(2)中k=0,1,2…,N/2是因为实数据序列的离散傅里叶变换关于中心对称,因此k可以只取前N/2+1个点。
步骤(3)、搜索数据序列x(n)的幅度谱P(k)的最大值所对应的离散频率索引kp,采用如下公式:
其中,表示在1≤k≤N/2-1范围内搜索幅度谱P(k)的最大值所对应的离散频率索引;
步骤(4)、利用数据序列的离散傅里叶变换X(l)计算Candan插值相对频率偏差
其中,Re{}代表取实部运算;X(kp)代表离散频率索引kp对应的离散傅里叶变换;X(kp-1)代表离散频率索引kp-1对应的离散傅里叶变换;X(kp+1)代表离散频率索引kp+1对应的离散傅里叶变换。
在步骤(4)中,计算Candan插值相对频率偏差本发明的方法将其作为相对频率偏差的预估值之一,该插值算法在相对频率偏差δ较小时,估计精度较高。
步骤(5)、利用数据序列的幅度谱P(k)计算Rife插值相对频率偏差
其中,P(kp)代表离散频率索引kp对应的幅度谱;P(kp-1)代表离散频率索引kp-1对应的幅度谱;P(kp+1)代表离散频率索引kp+1对应的幅度谱。
在步骤(5)中,计算Rife插值相对频率偏差本发明的方法将其作为相对频率偏差的预估值之一,该插值算法在相对频率偏差δ较大时,估计精度较高。
步骤(6)、通过对离散频率索引kp左、右相邻两离散频率索引所对应幅度谱的比较结果,结合相对频率偏差和估计得到综合相对频率偏差包括如下步骤:
步骤(61)、提取最大值离散频率索引kp左、右相邻两离散频率索引kp-1和kp+1所对应的幅度谱结果Al和Ar:
Al=P(kp-1) Ar=P(kp+1) (6)
步骤(62)、分别提取Al和Ar中的最大值和最小值Amax和Amin:
Amax=max(Al,AR) (7)
Amin=min(Al,AR) (8)
其中,max()代表取大值运算,min()代表取小值运算。
步骤(63)、根据相对频率偏差和的最大值和最小值Amax及Amin,估计出综合相对频率偏差
在步骤(63)中,η为设定的特征门限;特征门限η可以取大于0且小于1之间的任一值,但是通过蒙特卡洛仿真实验得到的优选值为η=0.75;
步骤(7)、根据综合相对频率偏差和离散频率索引kp,利用插值公式估计出单频信号的频率
其中,Δf为长度为N的离散傅里叶变换的频率分辨率,Δf=fs/N,fs为采样频率。
本发明的实施例中,仿真接收单频信号模型为:
其中A为信号幅度,为初始相位,τ为脉冲宽度,f0为信号中心频率,即待估计频率参数。w(t)为均值为0,方差为σ2高斯白噪声,方差σ2的大小由信噪比SNR决定:SNR=10log(A2/2σ2)。
以采样频率fs对上述单频信号进行离散采样可得到单频信号采样数据序列:
其中Nτ=int(fsτ),int()代表四舍五入运算。
为了验证本发明方法能够对相对频率偏差分布范围进行预估,在不同的频段采用不同的估计方法,提高频率估计的精度和稳健性,现列举两个实施例进行验证说明。
实施例1、
本实施例中,仿真信号参数分别设置为:信号幅度A=1,初始相位脉宽τ=0.256s,信号中心频率f0=750Hz,采样频率fs=4000Hz,观测数据序列点数N=512,对应的频率分辨率Δf=fs/N=4000/512=7.8125,相对频率偏差为δ=0,最大峰值频率索引为96,信噪比SNR=0dB,设定特征门限η=0.75。
依据步骤(2),计算所述数据序列x(n)的离散傅里叶变换X(l)和功率谱P(k),P(k)如图2所示。
依据步骤(3),搜索幅度谱P(k)的最大值所对应的离散频率索引kp=96。
依据步骤(4),利用数据序列的离散傅里叶变换X(l)计算Candan插值相对频率偏差
依据步骤(5),利用数据序列的幅度谱P(k)计算Rife插值相对频率偏差
因为P(kp-1)=0.1327,P(kp+1)=0.1005,P(kp-1)>P(kp+1),因此有:
依据步骤(6),通过对相对频率偏差和的分布范围进行预估,进而估计综合相对频率偏差包括如下步骤:
依据步骤(61)提取最大值离散频率索引左、右相邻两离散频率索引kp-1和kp+1所对应的幅度谱结果Al和Ar:
Al=P(kp-1)=0.1327 Ar=P(kp+1)=0.1005
依据步骤(62)分别提取Al和Ar中的最大值和最小值Amax和Amin:
Amax=max(Al,AR)=Al=0.1327
Amin=min(Al,AR)=AR=0.1005
因此有Amin/Amax=0.7575>η,依据步骤(63)得到综合相对频率偏差
依据步骤(7),利用插值公式估计出单频信号的频率
因此,估计频率的相对频率偏差为:
实施例2、
本实施例中,仿真信号参数分别设置为:信号幅度A=1,初始相位脉宽τ=0.512s,信号中心频率f0=350Hz,采样频率fs=4000Hz,观测数据序列点数N=1024,对应的频率分辨率Δf=fs/N=4000/1024=3.9063,相对频率偏差为δ=-0.4,最大峰值频率索引为90,信噪比SNR=0dB,设定特征门限η=0.75。
依据步骤(2),计算所述数据序列x(n)的离散傅里叶变换X(l)和功率谱P(k),P(k)如图3所示。
依据步骤(3),搜索幅度谱P(k)的最大值所对应的离散频率索引kp=90。
依据步骤(4),利用数据序列的离散傅里叶变换X(l)计算Candan插值相对频率偏差
依据步骤(5),利用数据序列的幅度谱P(k)计算Rife插值相对频率偏差
因为P(kp-1)=0.7156,P(kp+1)=0.257,P(kp-1)>P(kp+1),因此有:
依据步骤(6),通过对相对频率偏差分布范围进行预估,进而估计综合相对频率偏差包括如下步骤:
依据步骤(61)提取最大值离散频率索引左、右相邻两离散频率索引kp-1和kp+1所对应的幅度谱结果Al和Ar:
Al=P(kp-1)=0.7156 Ar=P(kp+1)=0.2570
依据步骤(62)分别提取Al和Ar中的最大值和最小值Amax和Amin:
Amax=max(Al,AR)=Al=0.7156
Amin=min(Al,AR)=AR=0.2570
因此有Amin/Amax=0.3591<η,依据步骤(63)得到综合相对频率偏差
依据步骤(7),利用插值估计出单频信号的频率
因此,估计频率的相对频率偏差为:
综上,本发明方法通过对相对频率偏差的分布范围进行预估,在不同频段采用不同的估计方法,仅在δ较大时采用Rife插值法,仅在δ较小时采用Candan插值法,比对所有δ均利用Rife插值法对δ的稳健性高;在不增加计算量的前提下,提高信号频率估计的精度,适合对单频信号频率进行实时、高精度估计,用于对实时性和频率估计精度要求高的场合。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (8)
1.一种基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,包括以下步骤:
步骤(1)、获取待处理的单频信号采样数据序列x(n),n=0,1,…,N-1,其中N为检测到的单频信号脉宽长度所对应的采样点个数,且N取值为2的整数次幂及N≥4;
步骤(2)、根据单频信号采样数据序列x(n)计算得到数据序列的离散傅里叶变换X(l)和幅度谱P(k);
步骤(3)、搜索数据序列的幅度谱P(k)最大值所对应的离散频率索引kp;
步骤(4)、利用数据序列的离散傅里叶变换X(l)计算Candan插值相对频率偏差
步骤(5)、利用数据序列的幅度谱P(k)计算Rife插值相对频率偏差
步骤(6)通过对离散频率索引kp左、右相邻两离散频率索引所对应幅度谱的比较结果,结合相对频率偏差和估计得到综合相对频率偏差
步骤(7)、根据综合相对频率偏差和离散频率索引kp,利用插值公式估计出单频信号的频率
2.根据权利要求1所述基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,所述步骤(1)获取待处理的单频信号采样数据序列x(n),具体为:
从传感器接收N个采样点的实时采集数据作为待处理的单频信号采样数据序列x(n);或从存储器中提取从检测到信号时刻起始的N个采样点的数据作为待处理的单频信号采样数据序列x(n)。
3.根据权利要求1所述基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,所述步骤(2)计算得到数据序列的离散傅里叶变换X(l)和幅度谱P(k),具体为:
步骤(21)、计算单频信号采样数据序列x(n)的离散傅里叶变换X(l):
其中,l为X(l)的离散频率索引,j表示虚数单位,即
步骤(22)、依据离散傅里叶变换X(l)计算单频信号采样数据序列x(n)的脉冲信号幅度谱P(k):
且k=0,1,2…,N/2
其中,k为P(k)的离散频率索引,X(k)为x(n)的离散傅里叶变换;| |代表取模值运算。
4.根据权利要求1所述基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,所述步骤(3)搜索数据序列的幅度谱P(k)最大值所对应的离散频率索引kp,采用公式:
其中,表示在1≤k≤N/2-1范围内搜索幅度谱P(k)的最大值所对应的离散频率索引。
5.根据权利要求1所述基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,所述步骤(4)计算Candan插值相对频率偏差采用公式:
其中,Re{}代表取实部运算;X(kp)代表离散频率索引kp对应的离散傅里叶变换;X(kp-1)代表离散频率索引kp-1对应的离散傅里叶变换;X(kp+1)代表离散频率索引kp+1对应的离散傅里叶变换。
6.根据权利要求1所述基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,所述步骤(5)计算Rife插值相对频率偏差采用公式:
其中,P(kp)代表离散频率索引kp对应的幅度谱;P(kp-1)代表离散频率索引kp-1对应的幅度谱;P(kp+1)代表离散频率索引kp+1对应的幅度谱。
7.根据权利要求1所述基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,所述步骤(6)估计得到综合相对频率偏差具体为:
步骤(61)、提取幅度谱P(k)最大值所对应的离散频率索引左、右相邻两离散频率索引kp-1和kp+1所对应的幅度谱结果Al和Ar:
Al=P(kp-1)Ar=P(kp+1)
步骤(62)、分别提取Al和Ar中的最大值和最小值Amax和Amin:
Amax=max(Al,AR)
Amin=min(Al,AR)
其中,max()代表取大值运算,min()代表取小值运算;
步骤(63)、根据相对频率偏差和的最大值和最小值Amax及Amin,估计得到综合相对频率偏差
其中,η为设定的特征门限,η为大于0且小于1的数。
8.根据权利要求1所述基于Candan-Rife综合内插的单频信号频率估计方法,其特征在于,所述步骤(7)中利用插值公式估计出单频信号的频率具体为:
其中,Δf为长度为N的离散傅里叶变换的频率分辨率,Δf=fs/N,fs为采样频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910279608.5A CN110068727B (zh) | 2019-04-09 | 2019-04-09 | 一种基于Candan-Rife综合内插的单频信号频率估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910279608.5A CN110068727B (zh) | 2019-04-09 | 2019-04-09 | 一种基于Candan-Rife综合内插的单频信号频率估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110068727A true CN110068727A (zh) | 2019-07-30 |
CN110068727B CN110068727B (zh) | 2021-03-30 |
Family
ID=67367201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910279608.5A Active CN110068727B (zh) | 2019-04-09 | 2019-04-09 | 一种基于Candan-Rife综合内插的单频信号频率估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110068727B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112541157A (zh) * | 2020-11-30 | 2021-03-23 | 西安精密机械研究所 | 一种信号频率精确估计方法 |
CN113009432A (zh) * | 2020-02-28 | 2021-06-22 | 加特兰微电子科技(上海)有限公司 | 提高测量精度、提高目标探测精度的方法、装置及设备 |
CN114090950A (zh) * | 2021-11-25 | 2022-02-25 | 成都飞机工业(集团)有限责任公司 | 平稳振动信号的宽频噪声等相位压缩方法及装置 |
CN114252699A (zh) * | 2021-12-24 | 2022-03-29 | 中国电子科技集团公司第十四研究所 | 一种离散傅里叶变换比幅测频方法 |
CN114280366A (zh) * | 2021-12-21 | 2022-04-05 | 中国航天科工集团八五一一研究所 | 基于改进型频率插值算法的正弦信号频率估计方法 |
CN115856424A (zh) * | 2023-03-01 | 2023-03-28 | 西安瀚博电子科技有限公司 | 基于峰邻比值的信号频率与幅度自适应提取方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680948A (zh) * | 2012-05-15 | 2012-09-19 | 东南大学 | 一种线性调频信号调频率和起始频率估计方法 |
CN103197141A (zh) * | 2013-03-21 | 2013-07-10 | 广东电网公司电力科学研究院 | 一种测量电力系统信号频率及谐波参数的方法 |
CN103941089A (zh) * | 2014-04-15 | 2014-07-23 | 南京邮电大学 | 基于dft的正弦信号频率估计方法 |
CN105738696A (zh) * | 2016-04-18 | 2016-07-06 | 天津大学 | 全相位时移相位差频率估计方法及装置 |
CN106443178A (zh) * | 2016-09-08 | 2017-02-22 | 东南大学 | 一种基于IQuinn‑Rife综合的正弦信号频率估计方法 |
CN106546949A (zh) * | 2016-10-28 | 2017-03-29 | 东南大学 | 一种基于频率预估计的双阵元正弦信号来波方向估计方法 |
CN107204840A (zh) * | 2017-07-31 | 2017-09-26 | 电子科技大学 | 基于dft和迭代校正的正弦信号频率估计方法 |
CN108469602A (zh) * | 2018-03-29 | 2018-08-31 | 东南大学 | 一种基于谱特征提取的脉冲信号类型自动判别方法 |
CN110007148A (zh) * | 2019-03-28 | 2019-07-12 | 东南大学 | 一种基于离散频谱相位和幅值综合内插的单频信号频率估计方法 |
-
2019
- 2019-04-09 CN CN201910279608.5A patent/CN110068727B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680948A (zh) * | 2012-05-15 | 2012-09-19 | 东南大学 | 一种线性调频信号调频率和起始频率估计方法 |
CN103197141A (zh) * | 2013-03-21 | 2013-07-10 | 广东电网公司电力科学研究院 | 一种测量电力系统信号频率及谐波参数的方法 |
CN103941089A (zh) * | 2014-04-15 | 2014-07-23 | 南京邮电大学 | 基于dft的正弦信号频率估计方法 |
CN105738696A (zh) * | 2016-04-18 | 2016-07-06 | 天津大学 | 全相位时移相位差频率估计方法及装置 |
CN106443178A (zh) * | 2016-09-08 | 2017-02-22 | 东南大学 | 一种基于IQuinn‑Rife综合的正弦信号频率估计方法 |
CN106546949A (zh) * | 2016-10-28 | 2017-03-29 | 东南大学 | 一种基于频率预估计的双阵元正弦信号来波方向估计方法 |
CN107204840A (zh) * | 2017-07-31 | 2017-09-26 | 电子科技大学 | 基于dft和迭代校正的正弦信号频率估计方法 |
CN108469602A (zh) * | 2018-03-29 | 2018-08-31 | 东南大学 | 一种基于谱特征提取的脉冲信号类型自动判别方法 |
CN110007148A (zh) * | 2019-03-28 | 2019-07-12 | 东南大学 | 一种基于离散频谱相位和幅值综合内插的单频信号频率估计方法 |
Non-Patent Citations (4)
Title |
---|
ÇAĞATAY CANDAN: "Fine resolution frequency estimation from three DFT", 《SIGNAL PROCESS》 * |
SHUXIA HUANG等: "Parameter Estimation of Delay-Doppler Underwater Acoustic Multi-Path Channel Based on Iterative Fractional Fourier Transform", 《IEEE ACCESS》 * |
唐鹏飞等: "基于FFT的正弦信号频率估计综合算法", 《船舶电子对抗》 * |
闵柏成等: "一种高精度频率估计的抗噪声调幅干扰算法", 《船舶电子对抗》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113009432A (zh) * | 2020-02-28 | 2021-06-22 | 加特兰微电子科技(上海)有限公司 | 提高测量精度、提高目标探测精度的方法、装置及设备 |
CN113009432B (zh) * | 2020-02-28 | 2024-05-31 | 加特兰微电子科技(上海)有限公司 | 提高测量精度、提高目标探测精度的方法、装置及设备 |
CN112541157A (zh) * | 2020-11-30 | 2021-03-23 | 西安精密机械研究所 | 一种信号频率精确估计方法 |
CN112541157B (zh) * | 2020-11-30 | 2024-03-22 | 西安精密机械研究所 | 一种信号频率精确估计方法 |
CN114090950A (zh) * | 2021-11-25 | 2022-02-25 | 成都飞机工业(集团)有限责任公司 | 平稳振动信号的宽频噪声等相位压缩方法及装置 |
CN114280366A (zh) * | 2021-12-21 | 2022-04-05 | 中国航天科工集团八五一一研究所 | 基于改进型频率插值算法的正弦信号频率估计方法 |
CN114280366B (zh) * | 2021-12-21 | 2023-10-31 | 中国航天科工集团八五一一研究所 | 基于改进型频率插值算法的正弦信号频率估计方法 |
CN114252699A (zh) * | 2021-12-24 | 2022-03-29 | 中国电子科技集团公司第十四研究所 | 一种离散傅里叶变换比幅测频方法 |
CN114252699B (zh) * | 2021-12-24 | 2023-12-15 | 中国电子科技集团公司第十四研究所 | 一种离散傅里叶变换比幅测频方法 |
CN115856424A (zh) * | 2023-03-01 | 2023-03-28 | 西安瀚博电子科技有限公司 | 基于峰邻比值的信号频率与幅度自适应提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110068727B (zh) | 2021-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110068727B (zh) | 一种基于Candan-Rife综合内插的单频信号频率估计方法 | |
CN110007148B (zh) | 一种基于离散频谱相位和幅值综合内插的单频信号频率估计方法 | |
CN107064629B (zh) | 一种基于频率相对偏差预估的分段综合单频信号频率估计方法 | |
CN103941089B (zh) | 基于dft的正弦信号频率估计方法 | |
CN104237871B (zh) | 一种基于相位补偿的时延差估计方法 | |
CN107102255B (zh) | 单一adc采集通道动态特性测试方法 | |
CN103746722B (zh) | 一种跳频信号跳周期和起跳时间估计方法 | |
CN110333389A (zh) | 基于插值dft的正弦信号频率估计方法 | |
CN108509377B (zh) | 一种基于沿特征提取的脉冲信号到达时间与脉宽估计方法 | |
CN104076200B (zh) | 一种双子段相位差频率估计方法及其装置 | |
CN110196407B (zh) | 一种基于频率预估的单矢量水听器信号来波方向估计方法 | |
CN103926599A (zh) | 基于emd迭代阈值滤波的gnss多径效应抑制方法 | |
CN109471095A (zh) | 基于快速迭代插值的fmcw雷达距离估计算法 | |
CN109379310B (zh) | 一种基于Rife-Quinn综合的MPSK信号载频估计方法 | |
CN106546949A (zh) | 一种基于频率预估计的双阵元正弦信号来波方向估计方法 | |
CN114280366A (zh) | 基于改进型频率插值算法的正弦信号频率估计方法 | |
CN102495285B (zh) | 对称窗函数功率重心估计电力谐波参数的方法 | |
CN112035790A (zh) | 井间定位信号频率估计方法 | |
CN102087313B (zh) | 一种卫星搜救信号的频率估计方法 | |
CN109117698B (zh) | 一种基于最小均方误差准则的噪声背景估计方法 | |
CN115826004B (zh) | 一种基于二维角度及时差联合的三星协同直接定位方法 | |
CN105738698A (zh) | 一种基于中心频移的谐波参数估计算法 | |
CN113468474B (zh) | 基于根Mini-Norm的电网频率估计方法 | |
CN113238184A (zh) | 一种基于非圆信号的二维doa估计方法 | |
Cai | Fast frequency measurement algorithm based on zero crossing method |
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 |