CN107064628A - 高精度频率测量系统及方法 - Google Patents
高精度频率测量系统及方法 Download PDFInfo
- Publication number
- CN107064628A CN107064628A CN201710240279.4A CN201710240279A CN107064628A CN 107064628 A CN107064628 A CN 107064628A CN 201710240279 A CN201710240279 A CN 201710240279A CN 107064628 A CN107064628 A CN 107064628A
- Authority
- CN
- China
- Prior art keywords
- frequency
- module
- signal
- measured value
- frequency measurement
- 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
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
- G01R23/04—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage adapted for measuring in circuits having distributed constants
-
- 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
- G01R23/163—Spectrum analysis; Fourier analysis adapted for measuring in circuits having distributed constants
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种高精度频率测量系统及方法,该系统包括:模数转换模块,用于将模拟中频信号转换为数字中频信号;混频模块,用于产生两路正交的本地载波将数字中频信号变频到数字基带信号;抽取滤波模块,用于对数字基带信号进行低通滤波与抽取,以降低数据速率;傅立叶变换模块,用于将短数据进行离散傅立叶变换得到频域信号;测频模块,基于频域信号中幅度最大值及相邻的两个计算值,采用三点插值测频得到第一频率测量值;扫描模块,以第一频率测量值为中心,在扫描范围内采用步进扫描,按照傅立叶变换方式逐点计算幅度最大值,得到扫描的第二频率测量值;选择器,用于选择第一、二频率测量值其中一个测量的结果。提高了频率测量精度。
Description
技术领域
本发明涉及信号处理技术领域,特别是涉及一种基于数字信号处理的短数据高精度频率测量系统及方法。
背景技术
参数估计是信号与信息处理学科的重要组成部分,也是近年来非常活跃、发展迅速的一个研究领域。在时间序列中,信号频率是一个重要的信号参数,测量淹没在噪声中的正弦波频率是现代信号处理中最有实用价值的技术之一,是测试所有谱估计性能的基础,也是信号处理技术的基础。该技术在雷达、电子对抗、声纳等领域有着广泛的应用,频率测量技术的进步必然推动以上应用领域的发展,随着现代通信及信息处理技术的迅猛发展,对频率测量技术的研究必然会提出更高的要求。
然而,现有的经典法(即,频率测量方法)都是以傅立叶变换(DFT)为基础的,快速傅立叶变换(FFT)的发现,促使经典法得到了更广泛的应用。但是经典法的局限性在于其测量分辨率与数据长度成正比,要提高分辨率必须增加数据长度。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种高精度频率测量系统及方法,用于解决现有技术中测量短数据精度不高的问题。
为实现上述目的及其他相关目的,本发明提供一种高精度频率测量系统,包括:
模数转换模块,接收模拟中频信号,用于将所述模拟中频信号转换为数字中频信号;
混频模块,其输入端连接所述模数转换模块的输出端,用于产生两路正交的本地载波将数字中频信号变频到数字基带信号;
抽取滤波模块,其输入端连接所述混频模块的输出端,用于对所述数字基带信号进行低通滤波与抽取处理,以降低数据速率;
傅立叶变换模块,其输入端连接于所述抽取滤波模块的输出端,用于将所述短数据进行离散傅立叶变换得到频域信号;
测频模块,其输入端连接所述傅立叶变换模块的输出端,基于傅立叶变换的频域信号中幅度最大值及相邻的两个计算值,采用三点插值测频得到第一频率测量值;
扫描模块,其输入端连接所述测频模块的输出端,以所述第一频率测量值为中心,在扫描范围内采用步进扫描,按照傅立叶变换方式逐点计算幅度最大值,得到扫描的第二频率测量值;
选择器,其输入端分别连接所述测频模块、及扫描模块的输出端,用于选择第一频率测量值与第二频率测量值中任意一个作为频率测量的结果。
本发明的另一目的在于提供一种高精度频率测量方法,包括:
将模拟中频信号进行模数转换生成数字中频信号;
利用混频模块产生两路正交的本地载波,将数字中频信号变频到数字基带信号;
对所述数字基带信号进行低通滤波与抽取处理,以降低数据速率;
将所述短数据进行离散傅立叶变换得到其对应的频域信号;
基于傅立叶变换中频域信号的幅度最大值及相邻的两个计算值,采用三点插值测频得到第一频率测量值;
以所述第一频率测量值为中心,在扫描范围内采用小频率步进,按照傅立叶变换方式进行计算,根据幅度最大值得到扫描的第二频率测量值;
选择第一频率测量值与第二频率测量值中任意一个作为频率测量的结果。
如上所述,本发明的高精度频率测量系统及方法,具有以下有益效果:
采用三点插值和精细扫描的方式突破了数据长度对频率测量精度的限制,即使通过短数据也能获取得很高的频率测量精度;同时,通过设置旁路控制和二选一电路,不仅增强了整个系统的灵活性,还避免了资源浪费。
附图说明
图1显示为本发明提供的一种高精度频率测量系统结构框图;
图2显示为图1中高精度频率测量系统中的抽取滤波模块的结构框图;
图3显示为图1中高精度频率测量系统中的频率信号的DFT幅值样本;
图4显示为本发明提供的一种高精度频率测量方法流程图。
元件标号说明:
1 模数转换模块
2 混频模块
3 抽取滤波模块
4 傅立叶变换模块
5 测频模块
6 扫描模块
7 选择器
8 参数配置模块
9 时钟模块
S1~S7 步骤1至步骤7
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
请参阅图1,本发明提供一种高精度频率测量系统,包括:
模数转换模块1,接收模拟中频信号,用于将所述模拟中频信号转换为数字中频信号;
混频模块2,其输入端连接所述模数转换模块1的输出端,用于产生两路正交的本地载波将数字中频信号变频到数字基带信号;
根据系统需求可采用过采样或带通采样等采样策略,当采用过采样策略时,载波频率与模拟中频频率相同,当采用带通采样策略时,载波频率需与带通采样后的信号中心频率一致。
另外,所述混频模块包括频率源与乘法器,所述频率源采用直接频率合成器21(DDS)实现,产生两路正交的本地载波;通过两路乘法器分别对应连接第一混频电路22、第二混频电路23,将数字中频信号下变频到数字基带信号。
抽取滤波模块3,其输入端连接所述混频模块2的输出端,用于对所述数字基带信号进行低通滤波与抽取处理,以降低数据速率;
其中,第一抽取滤波模块31与第二抽取滤波模块32对数字基带信号进行低通滤波和抽取,一方面滤除高频噪声;另一方面降低数据率;抽取倍数的选择需保证抽取后的信号频谱不发生混叠。
傅立叶变换模块4,其输入端连接于所述抽取滤波模块3的输出端,用于将所述短数据进行离散傅立叶变换得到频域信号;
其中,采用离散傅里叶变换电路(DFT/FFT),DFT变换电路采用快速傅里叶变换(FFT)算法,对短时间内经过混频、抽取和滤波后的短数据(速率降低的数据)进行离散傅里叶变换。
测频模块5,其输入端连接所述傅立叶变换模块4的输出端,基于傅立叶变换的频域信号中幅度最大值及相邻的两个计算值,采用三点插值测频得到第一频率测量值;
其中,所述测频模块包括幅值排序电路与三点插值电路,所述幅值排序电路,用于根据离散傅立叶变换所得的幅值大小进行排序得到幅度最大值;所述三点插值电路,用于根据幅度最大值以及其相邻两个计算值为依据,利用三点插值算法,进行多种三点插值算法运算得到第一频率测量值。
扫描模块6,其输入端连接所述测频模块5的输出端,以所述第一频率测量值为中心,在扫描范围内采用步进扫描,按照傅立叶变换方式逐点计算幅度最大值,得到扫描的第二频率测量值;
其中,所述扫描模块为精细扫描电路,用于根据三点插值测频得到第一频率测量值为中心,在扫描范围内采用小频率步进,按照傅立叶变换思路进行计算幅值最大点,得到对应的为第二频率测量值。
选择器7,其输入端分别连接所述测频模块5、扫描模块6的输出端,用于选择第一频率测量值与第二频率测量值中任意一个作为频率测量的结果。
所述选择器为二选一电路,其中,在基于三点插值测频模块和扫描模块的精细测量结果中选择一个作为系统的测量结果,该功能增强了系统的灵活性,可根据实际需求进行合适的选择以尽可能的降低系统复杂度。
在本实施例中,通过设置测频模块与扫描模块,相对于以傅立叶变换(DFT)为基础的扫频测量方法,突破了数据长度对频率测量精度的限制,即使通过短数据也能获取得很高的频率测量精度。
具体地,所述高精度频率测量系统还包括:参数配置模块7,其用于根据外部输入的配置信息对所述傅立叶变换模块中数据长度、抽取滤波模块的抽取倍数、滤波系数与旁路选择电路进行参数配置;
时钟模块8,其用于根据外部输入的配置信息产生各个模块所需的时钟信号。
在本实施例中,所述参数配置模块与时钟模块8均是外部输入的配置信息进行配置,提高了整个系统的灵活度。
请参阅图2,为图1中高精度频率测量系统中的抽取滤波模块的结构框图,包括:
所述抽取滤波电路包括积分梳妆滤波器31CIC、半带滤波器32HB(Half-BandFilter)、FIR(Finite Impulse Response)滤波器33、可变抽取器34与多个旁路选择电路A,其中,所述积分梳妆滤波器31、半带滤波器32、FIR滤波器33与可变抽取器34的首位依次相连,且所述积分梳妆滤波器31、半带滤波器32、FIR滤波器33与可变抽取器34均对应并联一个旁路选择电路A。
具体地,在本实施例中,所示CIC滤波器31和HB滤波器32可以对高数据率的信号进行快速抽取,使数据率快速降下来,由于CIC滤波器31的系数都是1,因此,只有加减运算,硬件实现非常简单,但其过渡带和阻带衰减特性不是很好。HB滤波器32的抽取因子固定为2,其滤波器系数近一半为零,可以节省一半的乘法器,非常适合采样率降一半的应用需求。FIR滤波器33的主要用途是对信道进行整形滤波,可变抽取电路34则可以进一步降低数据率。设置旁路选择电路使得系统具有更高的灵活性,满足多种应用需求。
请参阅图3,为图1中高精度频率测量系统中的频率信号的DFT幅值样本,包括:
一个单频点信号的DFT变换的幅值样本,其中,Fk表示DFT运算结果中幅值最大的点所对应的频率,Fk+1和Fk-1是其相邻的两个计算频率,Fpeak表示信号的真实频率,在理想情况下Fpeak位于Fk+1和Fk-1之间,在经典频率估计算法中,通常根据傅里叶变换的结果将Fk作为频率估计的结果,其频率估计精度最高只能达到DFT的物理分辨率,受到数据长度的影响。本发明中利用三点插值算法得出一个小数修正项δ,用以表示Fpeak与Fk之间的距离,最终得到更为精确的信号频率估计值Fpeak,可以完成三点插值的算法包括Jacobsen,Quinn和Macleod等人提出算法:
其中,式(1)至式(4)中Xk、Xk+1和Xk-1分别表示Fk、Fk+1和Fk-1对应的DFT计算结果,Re表示取实部,P和Q表示可变常数,用于调节不同窗函数带来的影响。三点插值的频率估计结果通过下式得到:
Fpeak=Fk+δfs/N (5)
其中,式(5)中fs为采样频率,N为参与DFT运算的点数,小数修正项δ,Fk表示DFT运算结果中幅值最大的点所对应的频率,Fpeak表示信号的真实频率。
在本实施例中,通过三点插值算法提高了频率测量系统的精度,打破了经典法的局限性,即测量分辨率与数据长度成正比,要提高分辨率必须增加数据长度的手段。
请参阅图4,为本发明提供的一种高精度频率测量方法流程图,包括:
步骤S1,将模拟中频信号进行模数转换生成数字中频信号;
其中,具体可采用模数转换器对输入的模拟中频信号进行转换,如ADC(模数转换电路)。
步骤S2,利用混频模块产生两路正交的本地载波,将数字中频信号变频到数字基带信号;
其中,具体可采用混频模块对数字中频信号进行处理,得到数字基带信号,在此不一一赘述。
步骤S3,对所述数字基带信号进行低通滤波与抽取处理,以降低数据速率;
其中,具体可采用抽取滤波模块,对所述数字基带信号进行滤波和处理,得到短数据。
步骤S4,将所述短数据进行离散傅立叶变换得到其对应的频域信号;
其中,具体可采用傅立叶变换模块进行离散傅立叶变换处理。
步骤S5,基于傅立叶变换中频域信号的幅度最大值及相邻的两个计算值,采用三点插值测频得到第一频率测量值;
其中,具体可采用测频模块进行处理。
步骤S6,以所述第一频率测量值为中心,在扫描范围内采用小频率步进,按照傅立叶变换方式进行计算,根据幅度最大值得到扫描的第二频率测量值;
其中,具体可采用扫描模块进行处理。
步骤S7,选择第一频率测量值与第二频率测量值中任意一个作为频率测量的结果。
其中,具体可采用二选一电路,选择其中一个作为频率测量的结果。
所述步骤6的具体方式,详述如下:
确定扫描范围和扫描步进m,在扫描范围内以第一频率测量值为中心;按照以下公式进行扫描计算:
式(6)中,x[n]为数据序列,X[m]中幅值最大的点所对应的频率即为测量结果,m为小数,其取值代表所扫描的频率点。
综上所述,本发明采用三点插值和精细扫描的方式突破了数据长度对频率测量精度的限制,即使通过短数据也能获取得很高的频率测量精度;同时,通过设置旁路控制和二选一电路,不仅增强了整个系统的灵活性,还避免了资源浪费。所以,本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (10)
1.一种高精度频率测量系统,其特征在于,包括:
模数转换模块,接收模拟中频信号,用于将所述模拟中频信号转换为数字中频信号;
混频模块,其输入端连接所述模数转换模块的输出端,用于产生两路正交的本地载波以将所述数字中频信号变频到数字基带信号;
抽取滤波模块,其输入端连接所述混频模块的输出端,用于对所述数字基带信号进行低通滤波与抽取处理以降低数据速率;
傅立叶变换模块,其输入端连接于所述抽取滤波模块的输出端,用于将所述短数据进行离散傅立叶变换得到频域信号;
测频模块,其输入端连接所述傅立叶变换模块的输出端,基于傅立叶变换的频域信号中幅度最大值及相邻的两个计算值,采用三点插值测频得到第一频率测量值;
扫描模块,其输入端连接所述测频模块的输出端,以所述第一频率测量值为中心,在扫描范围内采用步进扫描,按照傅立叶变换方式逐点计算幅度最大值,得到扫描的第二频率测量值;
选择器,其输入端分别连接所述测频模块及扫描模块的输出端,用于选择第一频率测量值与第二频率测量值中任意一个作为频率测量的结果。
2.根据权利要求1所述的高精度频率测量系统,其特征在于,所述混频模块包括频率源与乘法器,所述频率源采用直接频率合成器实现,产生两路正交的本地载波;通过两路乘法器分别对应连接第一混频电路、第二混频电路,将数字中频信号下变频到数字基带信号。
3.根据权利要求1所述的高精度频率测量系统,其特征在于,所述抽取滤波电路包括积分梳妆滤波器、半带滤波器、FIR滤波器、可变抽取器与多个旁路选择电路,其中,所述积分梳妆滤波器、半带滤波器、FIR滤波器与可变抽取器的首位依次相连,且所述积分梳妆滤波器、半带滤波器、FIR滤波器与可变抽取器均对应并联一个旁路选择电路。
4.根据权利要求1所述的高精度频率测量系统,其特征在于,所述测频模块包括幅值排序电路与三点插值电路,所述幅值排序电路,用于根据离散傅立叶变换所得的幅值大小进行排序得到幅度最大值;所述三点插值电路,用于根据幅度最大值以及其相邻两个计算值为依据,利用三点插值算法计算得到第一频率测量值。
5.根据权利要求1所述的高精度频率测量系统,其特征在于,所述扫描模块为精细扫描电路,用于根据三点插值测频得到第一频率测量值为中心,在扫描范围内采用小频率步进,按照傅立叶变换方式进行计算,得到对应的幅值最大点为第二频率测量值。
6.根据权利要求1所述的高精度频率测量系统,其特征在于,所述选择器为二选一电路。
7.根据权利要求1所述的高精度频率测量系统,其特征在于,还包括:参数配置模块,其用于根据外部输入的配置信息对所述傅立叶变换模块中数据长度、抽取滤波模块的抽取倍数、滤波系数与旁路选择电路进行参数配置。
8.根据权利要求1所述的高精度频率测量系统,其特征在于,还包括:时钟模块,其用于根据外部输入的配置信息产生各个模块所需的时钟信号。
9.一种高精度频率测量方法,其特征在于,包括:
将模拟中频信号进行模数转换生成数字中频信号;
利用混频模块产生两路正交的本地载波,将数字中频信号变频到数字基带信号;
对所述数字基带信号进行低通滤波与抽取处理,以降低数据速率;
将所述短数据进行离散傅立叶变换得到其对应的频域信号;
基于傅立叶变换中频域信号的幅度最大值及相邻的两个计算值,采用三点插值测频得到第一频率测量值;
以所述第一频率测量值为中心,在扫描范围内采用小频率步进,按照傅立叶变换方式进行计算,根据幅度最大值得到扫描的第二频率测量值;
选择第一频率测量值与第二频率测量值中任意一个作为频率测量的结果。
10.根据权利要求9所述的高精度频率测量方法,其特征在于,所述在扫描范围内采用小频率步进,采用傅立叶变换方式进行计算,根据幅度最大值得到扫描的第二频率测量值的步骤,包括:
确定扫描范围和扫描步进m,在扫描范围内以第一频率测量值为中心;按照以下公式进行扫描计算:
<mrow>
<mi>X</mi>
<mo>&lsqb;</mo>
<mi>m</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>x</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>N</mi>
</mfrac>
</msup>
</mrow>
式中,x[n]为数据序列,X[m]中幅值最大的点所对应的频率即为测量结果,m为小数,其取值代表所扫描的频率点。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710240279.4A CN107064628B (zh) | 2017-04-13 | 2017-04-13 | 高精度频率测量系统及方法 |
US16/493,184 US20200011911A1 (en) | 2017-04-13 | 2017-07-21 | High-precision frequency measuring system and method |
PCT/CN2017/093839 WO2018188228A1 (zh) | 2017-04-13 | 2017-07-21 | 高精度频率测量系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710240279.4A CN107064628B (zh) | 2017-04-13 | 2017-04-13 | 高精度频率测量系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064628A true CN107064628A (zh) | 2017-08-18 |
CN107064628B CN107064628B (zh) | 2019-08-16 |
Family
ID=59601185
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710240279.4A Active CN107064628B (zh) | 2017-04-13 | 2017-04-13 | 高精度频率测量系统及方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20200011911A1 (zh) |
CN (1) | CN107064628B (zh) |
WO (1) | WO2018188228A1 (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490255A (zh) * | 2018-04-16 | 2018-09-04 | 中国电子科技集团公司第二十九研究所 | 一种基于两级fft的数字宽带高精度测频方法及系统 |
WO2018188228A1 (zh) * | 2017-04-13 | 2018-10-18 | 中国电子科技集团公司第二十四研究所 | 高精度频率测量系统及方法 |
CN109142867A (zh) * | 2018-10-25 | 2019-01-04 | 闽南理工学院 | 基于改进滑窗离散傅里叶变换的谐波检测方法、设备 |
CN109782076A (zh) * | 2018-12-28 | 2019-05-21 | 北京航天测控技术有限公司 | 一种用于微波通讯装备的中频信号频率测试方法 |
CN109917405A (zh) * | 2019-03-04 | 2019-06-21 | 中国电子科技集团公司第十一研究所 | 一种激光测距方法及系统 |
CN110007122A (zh) * | 2018-01-05 | 2019-07-12 | 罗德施瓦兹两合股份有限公司 | 用于自动设置示波器的信号分析电路和方法 |
CN110068730A (zh) * | 2019-05-07 | 2019-07-30 | 中国科学院电子学研究所 | 双频段频谱数据采集方法及装置 |
CN110596457A (zh) * | 2019-09-26 | 2019-12-20 | 中国科学院微电子研究所 | 数字扫频系统及方法 |
CN110609172A (zh) * | 2019-08-15 | 2019-12-24 | 上海微波技术研究所(中国电子科技集团公司第五十研究所) | 脉冲调制信号载波频率测量系统及方法 |
CN111220846A (zh) * | 2020-03-10 | 2020-06-02 | 星汉时空科技(北京)有限公司 | 一种高速采样全数字化频率稳定度的测试设备及方法 |
CN111562438A (zh) * | 2020-04-28 | 2020-08-21 | 重庆金美通信有限责任公司 | 一种基于fft和相位差的正弦信号频率估计方法和装置 |
CN114838809A (zh) * | 2022-03-22 | 2022-08-02 | 珠海市运泰利自动化设备有限公司 | 一种自适应提高频率测量精度的音频信号测量方法 |
CN115856424A (zh) * | 2023-03-01 | 2023-03-28 | 西安瀚博电子科技有限公司 | 基于峰邻比值的信号频率与幅度自适应提取方法 |
CN117268288A (zh) * | 2023-07-26 | 2023-12-22 | 北京大学长三角光电科学研究院 | 光学衍射层析成像激光扫描方法、装置和电子设备 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022091329A1 (ja) * | 2020-10-30 | 2022-05-05 | 三菱電機株式会社 | 周波数検出器 |
CN113075451A (zh) * | 2021-02-23 | 2021-07-06 | 北京鸿普惠信息技术有限公司 | 一次调频中正序分量补偿角度偏移提高频率精度的方法 |
CN113406386B (zh) * | 2021-06-23 | 2023-04-25 | 中国电子科技集团公司第二十九研究所 | 一种基于数字下变频的信号频率精确估计方法 |
CN114244384B (zh) * | 2021-12-20 | 2023-06-16 | 浙江嘉科电子有限公司 | 一种用于零中频接收机的连续波峰值捕获电路及方法 |
CN116112104B (zh) * | 2023-01-30 | 2024-01-12 | 成都九洲迪飞科技有限责任公司 | 一种基于fpga的频谱监测方法及前端接收机 |
CN118191414B (zh) * | 2024-05-15 | 2024-07-23 | 北京航空航天大学 | 一种基于机床自激励的数控机床频率特性测试方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07239351A (ja) * | 1994-02-25 | 1995-09-12 | Nippon Signal Co Ltd:The | 周波数分析装置 |
CN101179549A (zh) * | 2007-12-14 | 2008-05-14 | 清华大学 | 采用三点加权插值算法的通信信号载频估计方法 |
CN101261292A (zh) * | 2008-04-14 | 2008-09-10 | 湖南大学 | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 |
CN105137175A (zh) * | 2015-08-06 | 2015-12-09 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4896102A (en) * | 1988-06-13 | 1990-01-23 | Scientific-Atlanta, Inc. | Spectrum analyzer |
JPH0552883A (ja) * | 1991-08-23 | 1993-03-02 | Nec Corp | 周波数検出器 |
CN101603985B (zh) * | 2009-07-15 | 2011-11-16 | 北京航空航天大学 | 高准确度正弦信号测量方法 |
CN103837740A (zh) * | 2013-12-25 | 2014-06-04 | 北京航天测控技术有限公司 | 一种高精度数字瞬时测频方法及装置 |
CN106443122B (zh) * | 2016-08-18 | 2020-05-15 | 中国电子科技集团公司第四十一研究所 | 一种宽频带大动态信号高精度测量装置及方法 |
CN106405235B (zh) * | 2016-10-17 | 2019-09-10 | 深圳市鼎阳科技有限公司 | 一种频谱分析仪及其数据处理方法 |
CN106443180B (zh) * | 2016-11-09 | 2023-03-07 | 优利德科技(中国)股份有限公司 | 测量闸门自适应的高精度频率测量方法及测量装置 |
CN107064628B (zh) * | 2017-04-13 | 2019-08-16 | 中国电子科技集团公司第二十四研究所 | 高精度频率测量系统及方法 |
-
2017
- 2017-04-13 CN CN201710240279.4A patent/CN107064628B/zh active Active
- 2017-07-21 US US16/493,184 patent/US20200011911A1/en not_active Abandoned
- 2017-07-21 WO PCT/CN2017/093839 patent/WO2018188228A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07239351A (ja) * | 1994-02-25 | 1995-09-12 | Nippon Signal Co Ltd:The | 周波数分析装置 |
CN101179549A (zh) * | 2007-12-14 | 2008-05-14 | 清华大学 | 采用三点加权插值算法的通信信号载频估计方法 |
CN101261292A (zh) * | 2008-04-14 | 2008-09-10 | 湖南大学 | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 |
CN105137175A (zh) * | 2015-08-06 | 2015-12-09 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018188228A1 (zh) * | 2017-04-13 | 2018-10-18 | 中国电子科技集团公司第二十四研究所 | 高精度频率测量系统及方法 |
CN110007122A (zh) * | 2018-01-05 | 2019-07-12 | 罗德施瓦兹两合股份有限公司 | 用于自动设置示波器的信号分析电路和方法 |
CN108490255A (zh) * | 2018-04-16 | 2018-09-04 | 中国电子科技集团公司第二十九研究所 | 一种基于两级fft的数字宽带高精度测频方法及系统 |
CN109142867A (zh) * | 2018-10-25 | 2019-01-04 | 闽南理工学院 | 基于改进滑窗离散傅里叶变换的谐波检测方法、设备 |
CN109782076A (zh) * | 2018-12-28 | 2019-05-21 | 北京航天测控技术有限公司 | 一种用于微波通讯装备的中频信号频率测试方法 |
CN109917405A (zh) * | 2019-03-04 | 2019-06-21 | 中国电子科技集团公司第十一研究所 | 一种激光测距方法及系统 |
CN110068730A (zh) * | 2019-05-07 | 2019-07-30 | 中国科学院电子学研究所 | 双频段频谱数据采集方法及装置 |
CN110068730B (zh) * | 2019-05-07 | 2020-08-07 | 中国科学院电子学研究所 | 双频段频谱数据采集方法及装置 |
CN110609172A (zh) * | 2019-08-15 | 2019-12-24 | 上海微波技术研究所(中国电子科技集团公司第五十研究所) | 脉冲调制信号载波频率测量系统及方法 |
CN110596457A (zh) * | 2019-09-26 | 2019-12-20 | 中国科学院微电子研究所 | 数字扫频系统及方法 |
CN110596457B (zh) * | 2019-09-26 | 2022-04-29 | 中国科学院微电子研究所 | 数字扫频系统及方法 |
CN111220846A (zh) * | 2020-03-10 | 2020-06-02 | 星汉时空科技(北京)有限公司 | 一种高速采样全数字化频率稳定度的测试设备及方法 |
CN111220846B (zh) * | 2020-03-10 | 2022-04-19 | 星汉时空科技(北京)有限公司 | 一种高速采样全数字化频率稳定度的测试设备及方法 |
CN111562438A (zh) * | 2020-04-28 | 2020-08-21 | 重庆金美通信有限责任公司 | 一种基于fft和相位差的正弦信号频率估计方法和装置 |
CN111562438B (zh) * | 2020-04-28 | 2022-04-08 | 重庆金美通信有限责任公司 | 一种基于fft和相位差的正弦信号频率估计方法和装置 |
CN114838809A (zh) * | 2022-03-22 | 2022-08-02 | 珠海市运泰利自动化设备有限公司 | 一种自适应提高频率测量精度的音频信号测量方法 |
CN115856424A (zh) * | 2023-03-01 | 2023-03-28 | 西安瀚博电子科技有限公司 | 基于峰邻比值的信号频率与幅度自适应提取方法 |
CN117268288A (zh) * | 2023-07-26 | 2023-12-22 | 北京大学长三角光电科学研究院 | 光学衍射层析成像激光扫描方法、装置和电子设备 |
CN117268288B (zh) * | 2023-07-26 | 2024-03-08 | 北京大学长三角光电科学研究院 | 光学衍射层析成像激光扫描方法、装置和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
WO2018188228A1 (zh) | 2018-10-18 |
CN107064628B (zh) | 2019-08-16 |
US20200011911A1 (en) | 2020-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107064628A (zh) | 高精度频率测量系统及方法 | |
CN101825660B (zh) | 欠采样下的正弦信号频率的高效测量方法及实施装置 | |
CN103973324B (zh) | 一种宽带数字接收机及其实时频谱处理方法 | |
CN103457603B (zh) | 一种基于平均频谱测试adc动态参数的方法 | |
CN103941087B (zh) | 欠采样速率下的高频余弦信号的频率测量方法及其装置 | |
CN105510706B (zh) | 一种高精度欠采样测频方法 | |
CN110020452B (zh) | 电磁兼容性测试仿真分析方法及其系统 | |
CN102043091B (zh) | 数字化高精度相位检测器 | |
CN104375111A (zh) | 对密集频谱进行快速高精度细化校正的方法 | |
CN111800091A (zh) | 一种数字下变频的实现方法及存储介质 | |
CN104007316A (zh) | 一种欠采样速率下的高精度频率测量方法及其测量仪 | |
CN108650048A (zh) | 一种高精度数字阵列多通道延时补偿方法 | |
US7564386B2 (en) | Pre-processing data samples from parallelized data converters | |
CN104483619A (zh) | 一种基于虚拟仪器的频率特性测试系统 | |
CN101063695B (zh) | 无功功率计算电路和方法 | |
CN106972833A (zh) | 一种可任意倍数重采样的数字上变频器 | |
CN103941280B (zh) | 基于冲激响应不变法的数字核脉冲高斯成形方法 | |
CN107659290B (zh) | 带宽扩展滤波器及其设计方法 | |
CN108919202A (zh) | 一种基于认知机制的非均匀动态滤波器组的构建方法 | |
CN105548711B (zh) | 一种多频信息滤波递推解调方法 | |
CN104991119A (zh) | 一种消除伪峰、谱泄漏效应的互素谱分析方法及其装置 | |
CN102072987A (zh) | 短区间正弦信号的相位估计法及其实验装置 | |
CN108008188A (zh) | 一种通道间相位差、幅度差的测量方法 | |
CN114114231A (zh) | 一种基于fft-czt的雷达测距方法 | |
CN107026301A (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 |