CN103837740A - 一种高精度数字瞬时测频方法及装置 - Google Patents

一种高精度数字瞬时测频方法及装置 Download PDF

Info

Publication number
CN103837740A
CN103837740A CN201310728431.5A CN201310728431A CN103837740A CN 103837740 A CN103837740 A CN 103837740A CN 201310728431 A CN201310728431 A CN 201310728431A CN 103837740 A CN103837740 A CN 103837740A
Authority
CN
China
Prior art keywords
frequency
signal
conversion
digital
carried out
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.)
Pending
Application number
CN201310728431.5A
Other languages
English (en)
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.)
Beijing Aerospace Measurement and Control Technology Co Ltd
Original Assignee
Beijing Aerospace Measurement and Control 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 Aerospace Measurement and Control Technology Co Ltd filed Critical Beijing Aerospace Measurement and Control Technology Co Ltd
Priority to CN201310728431.5A priority Critical patent/CN103837740A/zh
Publication of CN103837740A publication Critical patent/CN103837740A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

本发明公开了一种高精度数字瞬时测频方法及装置。其中,该方法包括:采用射频接收前端进行射频信号变频、滤波处理,将X波段信号变换为中频信号;将中频信号进行AD变换,通过数字信号处理器件将变换后的中频信号进行下变频处理,得到中频数字信号;将中频数字信号变换为基带数字信号;采用基于FFT和信号平方的双信号估算法,通过FPGA对基带数字信号进行计算处理,得到测频结果。解决了相关技术中X波段脉冲信号的测频精度不高的问题,本发明采用基于FFT和信号平方的双信号频率估计算法,结合FFT测频算法相对快速和一维平方信号谱峰搜索后插值估计频率间隔的精确性,利用FPGA硬件实现相关算法,保证了测频的快速性和准确性。

Description

一种高精度数字瞬时测频方法及装置
技术领域
本发明涉及G01R2300频率测量测试领域,特别是涉及一种高精度数字瞬时测频方法及装置。 
背景技术
在20世纪五六十年代,瞬时测频技术产生,主要分为模拟瞬时测频和数字瞬时测频两种。其中,模拟瞬时测频设备简单易于实现,但测频精度不高,且稳定性差。数字瞬时测频具有稳定、精确度高、实时性好等特点。 
数字瞬时测频即用数字信号处理的方法对信号的载频进行测量,要求数字测频必须可以对大带宽的信号进行高速、高精度的测量。数字测频算法种类繁多,根据算法所利用的信号不同信息分量,数字测频算法可分为:利用幅度信息的测频算法、利用相位信息的测频算法、利用频率信息的测频算法、利用功率谱信息的测频算法4类。 
利用频率信息进行测频,包括在时域测频和频域测频两类算法。在时域上对频率进行推算的典型算法为频率推算法,而在频域上对频率进行推算的典型算法为傅里叶变换法。 
傅里叶变换法是将一段信号进行傅里叶变换,将时域信息转化为频域信息,在频域上的峰值点及对应该信号的频率,其测频精度依赖于采样点数的多少。 
对于X波段的载频快速测量而言,由于X波段脉冲信号是高频信号,其调制后的信号最小脉宽为0.2us,在信号工作比小于0.000004时,要求在0.4us时间内完成相对精确的频率测量,基于上述各个因素,最终的测频精度很难满 足≤±200KHz的要求,即测频精度不高。 
针对相关技术中X波段脉冲信号的测频精度不高的问题,目前尚未提出有效的解决方案。 
发明内容
针对相关技术中X波段脉冲信号的测频精度不高的问题,本发明提供了一种高精度数字瞬时测频方法及装置,用以解决上述技术问题。 
根据本发明的一个方面,本发明提供了一种高精度数字瞬时测频方法,其中,该方法包括:采用射频接收前端进行射频信号变频、滤波处理,将X波段信号变换为中频信号;将所述中频信号进行模数AD变换,通过数字信号处理器件将变换后的中频信号进行下变频处理,得到中频数字信号;将所述中频数字信号变换为基带数字信号;采用基于快速傅氏变换(Fast Fourier Transformation,简称为FFT)和信号平方的双信号估算法,通过现场可编程门阵列(Field-Programmable Gate Array,简称为FPGA)对所述基带数字信号进行计算处理,得到测频结果。 
优选地,通过所述FPGA对所述基带数字信号进行计算处理,包括以下步骤:所述FPGA根据同步触发对所述基带数字信号进行采集,并采用快进慢出的方式读取所述信号的数据;对读取的所述数据进行有限长单位冲激响应滤波器(Finite Impulse Response,简称为FIR)带通滤波处理,并进行缓存;将滤波后的数据进行快速傅里叶变换,得到相应的频谱信息,去除所述频谱信息中的镜频部分,保留正频率部分,提取最大值为频谱值,以及其位置;根据采样率,频谱值,以及其位置进行映射运算得到相应的频率值,确定频率的位置;采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果。 
优选地,采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果,包括以下步骤:根据信号中两频率分量的幅度比,通过一维平方信号谱峰搜索,得到所述信号的和频与差频分量的估计值;并利用插值技术得到频 率结果。 
优选地,在得到所述频率结果后,所述方法还包括:按照约定的通讯协议进行频率转换,将转换后的频率结果输出。 
优选地,所述方法还包括:在测频过程中进行状态管理和故障上报。 
根据本发明的另一方面,本发明还提供了一种高精度数字瞬时测频装置,其中,该装置包括:中频变换模块,用于采用射频接收前端进行射频信号变频、滤波处理,将X波段信号变换为中频信号;模数变换模块,用于将所述中频信号进行模数AD变换,通过数字信号处理器件将变换后的中频信号进行下变频处理,得到中频数字信号;测频模块,用于将所述中频数字信号变换为基带数字信号;采用基于快速傅氏变换FFT和信号平方的双信号估算法,通过现场可编程门阵列FPGA对所述基带数字信号进行计算处理,得到测频结果。 
优选地,所述测频模块包括:数据采集单元,用于根据同步触发对所述基带数字信号进行采集,并采用快进慢出的方式读取所述信号的数据;滤波缓存单元,用于对读取的所述数据进行FIR带通滤波处理,并进行缓存;时频变换单元,用于将滤波后的数据进行快速傅里叶变换,得到相应的频谱信息,去除所述频谱信息中的镜频部分,保留正频率部分,提取最大值为频谱值,以及其位置;频率计算单元,用于根据采样率,频谱值,以及其位置进行映射运算得到相应的频率值,确定频率的位置;双信号估算单元,用于采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果。 
优选地,所述双信号估算单元包括:双信号估算子单元,用于根据信号中两频率分量的幅度比,通过一维平方信号谱峰搜索,得到所述信号的和频与差频分量的估计值;并利用插值技术得到频率结果。 
优选地,所述装置还包括:结果输出模块,用于按照约定的通讯协议进行频率转换,将转换后的频率结果输出。 
优选地,所述装置还包括:管理模块,用于在测频过程中进行状态管理和故障上报。 
本发明采用射频接收前端进行射频信号的变频、滤波,将X频段信号变换到中频;利用宽带ADC完成变换后经过专用的数字信号处理器件,进行数字下变频处理降低数据流速率;再将中频数字信号变换成基带数字信号,最后送入FPGA进行计算处理,完成频率的快速测量。解决了相关技术中X波段脉冲信号的测频精度不高的问题,本发明采用基于FFT和信号平方的双信号频率估计算法,结合FFT测频算法相对快速和一维平方信号谱峰搜索后插值估计频率间隔的精确性,利用FPGA硬件实现相关算法,保证了测频的快速性和准确性。 
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。 
附图说明
图1是根据本发明实施例的高精度数字瞬时测频方法的流程图; 
图2是根据本发明实施例的高精度数字瞬时测频装置的结构框图; 
图3是根据本发明实施例的数字测频组件功能框图; 
图4是根据本发明实施例的数字测频组件的射频接收前端原理框图; 
图5是根据本发明实施例的数字测频组件的硬件电路组成框图。 
具体实施方式
为了解决相关技术中X波段脉冲信号的测频精度不高的问题,本发明提供了一种高精度数字瞬时测频方法及装置,以下结合附图以及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不限定本发明。 
本实施例提供了一种高精度数字瞬时测频方法,图1是根据本发明实施例的高精度数字瞬时测频方法的流程图,如图1所示,该方法包括以下步骤(步骤S102-步骤S106): 
步骤S102,采用射频接收前端进行射频信号变频、滤波处理,将X波段信号变换为中频信号。 
步骤S104,将中频信号进行AD变换,通过数字信号处理器件将变换后的中频信号进行下变频处理,得到中频数字信号。 
步骤S106,将中频数字信号变换为基带数字信号;采用基于快速傅氏变换FFT和信号平方的双信号估算法,通过FPGA对基带数字信号进行计算处理,得到测频结果。 
本实施例采用射频接收前端进行射频信号的变频、滤波,将X频段信号变换到中频;利用宽带ADC完成变换后经过专用的数字信号处理器件,进行数字下变频处理降低数据流速率,再将中频数字信号变换成基带数字信号,最后送入FPGA进行计算处理,完成频率的快速测量。解决了相关技术中X波段脉冲信号的测频精度不高的问题,本发明采用基于FFT和信号平方的双信号频率估计算法,结合FFT测频算法相对快速和一维平方信号谱峰搜索后插值估计频率间隔的精确性,利用FPGA硬件实现相关算法,保证了测频的快速性和准确性。 
在上述步骤S106中,通过FPGA对基带数字信号进行计算处理,可以通过以下优选实施方式实现:FPGA根据同步触发对基带数字信号进行采集,并采用快进慢出的方式读取信号的数据;对读取的数据进行FIR带通滤波处理,并进行缓存;将滤波后的数据进行快速傅里叶变换,得到相应的频谱信息,去除频谱信息中的镜频部分,保留正频率部分,提取最大值为频谱值,以及其位置;根据采样率,频谱值,以及其位置进行映射运算得到相应的频率值,确定频率的位置;采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果。 
对于采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果,可以通过以下优选实施方式实现:根据信号中两频率分量的幅度比,通过一维平方信号谱峰搜索,得到信号的和频与差频分量的估计值;并利用插值技 术得到频率结果。 
在得到频率结果后,该方法还包括:按照约定的通讯协议进行频率转换,将转换后的频率结果输出。 
该方法还包括:在测频过程中进行状态管理和故障上报。 
下面对上述步骤S106中,采用基于快速傅氏变换FFT和信号平方的双信号估算法,通过现场可编程门阵列FPGA对基带数字信号进行计算处理,得到测频结果的具体过程,进行介绍。该过程可通过以下几个阶段实现。 
(1)采集控制 
根据同步触发对信号进行采集,连续采集120us,最多包含24个脉冲,数据传输采用双沿DDR方式,降低数据传输时钟频率。FPGA端采用IDDR内置双沿读取模块读取数据,保证数据读取正确率。对接收到的双沿数据进行排序恢复,恢复为一列连续数据。 
(2)数据缓存 
对输入数据进行缓存,数据率高达800MB/s,缓存控制需采取扩位降速,并通过严格的布局布线规则,保证FPGA数据缓存正确。 
(3)滤波 
数据读取时采用快进慢出的方式进行,高速写入,读取时根据计算需要,慢速读出数据进行FIR带通滤波,提高信号的信杂比,减少干扰,提高计算准确度。 
(4)时频变换 
将滤波后的数据进行快速傅里叶变换,得到相应的频谱信息,去除镜频部分,保留正频率部分,提取最大值及其位置。 
(5)频率计算 
根据采样率,频谱值,及其位置进行映射运算得到相应的频率值。 
FFT是把信号映射到正交的负指数空间,仿佛FFT把信号映射到正交的方波空间。比较其系数的大小就可以估计其所在的频率空间。 
设信号为s(n)=cos(Ωn+φ)+e(n) n=1,2.....N 
其中Ω为数字频率,N为样本空间,e(n)是均值为0,方差为σ2的高斯白噪声。 
则取其符号为: 
x(n)=sign[s(n)]=sign[cos(Ωn+φ)+e(n)] 
Φ Nc kn = sin g { cos [ 2 Π ( k + 1 ) ( n - 1 ) / N }
φ Ns kn = sign { sin [ 2 π ( k + 1 ) ( n - 1 ) / N }
Xc ( k ) = Σ n = 1 N x ( n ) φ Nc kn
Xs ( k ) = Σ n = 1 N x ( n ) φ Ns kn
Pxcs(k)=Xc(k)2+Xs(k)2
找出Pxcs(k)的最大值,就可确定频率的位置。 
(6)信号平方与FFT结合的双信号频率估算 
本发明基于FFT测频技术,发明了一种信号平方与FFT结合的双信号频率估算算法。根据信号中两频率分量的幅度比,只需一次一维平方信号谱峰搜索,就可以得到该信号的和频与差频分量的估计值,并利用插值技术提高测频精度。该算法能够精确地估计频率间隔小的双信号频率,且容易扩展到复信号,FPGA硬件实现简单。 
无噪声下的双实正弦信号观测序列为: 
x(n)=a1sin(ω1n)+a2sin(ω2n)  n=0,1…..,N-1 
a1和a2为信号频率分量的确定性幅度,角频率ω1=2πf1/fs,ω2=2πf2/fs,fs为采样频率,N为FFT长度。 
双实正弦信号的平方序列s(n)定义为: 
s(n)=x2(n)=0.5(a1 2+a2 2)+0.5a1 2cos(2ω1n)+0.5a2 2cos(2ω2n)+a1a2[cos(ω-n)-cos(ω+n)] 
其中ω-=ω12,ω+=ω12; 
假设s(n)的DFT为S(k) 
S ( k ) = Σ n = 0 N - 1 s ( n ) exp ( - j 2 πnk / N ) , k = 0,1 . . . . . , N - 1
注意到s(n)中的0.5(a1 2+a2 2)为直流分量的幅度,它只影响|S(0)|,不影响其他的谱线值。因此,我们可以假设谱峰对应的位置k≥1,谱峰搜索和插值都不会受到直流分量的影响。 
理论上,在频谱|S(k)|中,有两个相等最强频率分量对应ω-和ω+,两个弱频率分量对应2ω1和2ω2,现假设两强频率分量的频率位置为k-和k+,分别对应于具有最高谱峰值的ω-和ω+,即是 
| S ( k - ) | = | S ( k + ) | = Max 1 ≤ k ≤ N - 1 { | S ( k ) | }
可以通过一维谱峰搜索得到最高谱峰位置,计算位置就可估计出ω-和ω+,即 
ω ^ - = 2 π k - / N , ω ^ + = 2 π k + / N
其中
Figure DEST_PATH_GDA0000486889840000084
Figure DEST_PATH_GDA0000486889840000085
表示ω-和ω+的估计,则 
f ^ 1 = ( ω ^ + + ω ^ - ) fs 4 π f ^ 2 = ( ω ^ + - ω ^ - ) fs 4 π
其中
Figure DEST_PATH_GDA0000486889840000087
分别表示对f1、f2的估计。 
由以上推导可知,当信号的两个频率间隔很小时,该算法仍能准确区分,即该算法可以准确区分频率间隔小的双信号。 
(7)格式转换输出 
得到频率结果后,按照约定的通讯协议进行转换,将结果输出。 
(8)状态管理 
对AD,FPGA,电平转换等主要器件进行状态管理和故障上报,通过灯来指示电源及FPGA状态,AD,通讯及其它状态由FPGA进行检测汇总并上报。 
对应于上述实施例介绍的高精度数字瞬时测频方法,本实施例提供了一种高精度数字瞬时测频装置,用以实现上述实施例。图2是根据本发明实施例的高精度数字瞬时测频装置的结构框图,如图2所示,该装置包括:中频变换模块10、模数变换模块20和测频模块30。下面对该结构进行详细介绍。 
中频变换模块10,用于采用射频接收前端进行射频信号变频、滤波处理,将X波段信号变换为中频信号; 
模数变换模块20,连接至中频变换模块10,用于将中频信号进行AD变换,通过数字信号处理器件将变换后的中频信号进行下变频处理,得到中频数字信号; 
测频模块30,连接至模数变换模块20,用于将中频数字信号变换为基带数字信号;采用基于FFT和信号平方的双信号估算法,通过FPGA对基带数字信号进行计算处理,得到测频结果。 
本实施例采用射频接收前端进行射频信号的变频、滤波,将X频段信号变换到中频;利用宽带ADC完成变换后经过专用的数字信号处理器件,进行数字下变频处理降低数据流速率,再将中频数字信号变换成基带数字信号,最后送入FPGA进行计算处理,完成频率的快速测量。解决了相关技术中X波段脉冲信号的测频精度不高的问题,本发明采用基于FFT和信号平方的双信号频率估计算法,结合FFT测频算法相对快速和一维平方信号谱峰搜索后插值估计频率间隔的精确性,利用FPGA硬件实现相关算法,保证了测频的快速性和准确性。 
优选地,上述测频模块30包括:数据采集单元,用于根据同步触发对基带数字信号进行采集,并采用快进慢出的方式读取信号的数据;滤波缓存单元,用于对读取的数据进行FIR带通滤波处理,并进行缓存;时频变换单元,用于将滤波后的数据进行快速傅里叶变换,得到相应的频谱信息,去除频谱信息中 的镜频部分,保留正频率部分,提取最大值为频谱值,以及其位置;频率计算单元,用于根据采样率,频谱值,以及其位置进行映射运算得到相应的频率值,确定频率的位置;双信号估算单元,用于采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果。 
优选地,上述双信号估算单元包括:双信号估算子单元,用于根据信号中两频率分量的幅度比,通过一维平方信号谱峰搜索,得到信号的和频与差频分量的估计值;并利用插值技术得到频率结果。 
优选地,上述装置还可以包括:结果输出模块,用于按照约定的通讯协议进行频率转换,将转换后的频率结果输出。管理模块,用于在测频过程中进行状态管理和故障上报。 
图3是根据本发明实施例的数字测频组件功能框图,如图3所示,本实施例中的数字测频组件以功能划分包括三部分:用于射频信号变换、位于模数变换之前的射频接收前端、宽带模数转换器(ADC)以及位于ADC之后的高速数字信号处理单元(FPGA)。 
图4是根据本发明实施例的数字测频组件的射频接收前端原理框图,如图4所示,射频接收前端主要由限幅器、电子开关、放大器、混频器、滤波器组成,完成放大、变频、滤波等功能,即将接收到的射频信号变换到中频,但是中频频率的选择和后面的采样是息息相关的。对射频模拟信号的采样数字化可有多种方法,本实施例使用的是基于Nyquist采样定理的低通采样及带通采样。 
图5是根据本发明实施例的数字测频组件的硬件电路组成框图,如图5所示,射频接收前端输出的中频信号经ADC进行模数转换后,输入到高速数字信号处理单元(FPGA)进行计算。 
在频域上对频率进行推算使用的是傅里叶变换法,频率推算法又称为曲线拟合法,是在假设信号为一个单一载频的情况下导出的一种算法。傅里叶变换法是将一段信号进行傅里叶变换,将时域信息转换为频域信息,在频域上的峰值点即对应该信号的频率,其测频精度依赖于采样点数的多少。 
假设被测量的信号具有如下形式: 
Figure DEST_PATH_GDA0000486889840000111
如果采样周期为T,则信号变为如下序列: 
则可推出:si+1-si-1=2sicos(ωT) 
采样率fs=1/T,则信号频率与采样频率之间的关系为: 
ω = arccos ( s i + 1 + s i - 1 2 s i ) × f s
在已知采样率情况下,只需要3个采样点就可以得出信号的频率。 
从以上的描述中可知,在本发明的技术方案中,高精度数字瞬时测频组件,将X波段信号变换到中频,经ADC变换后送入专用数字信号处理器件,采用基于FFT和信号平方的双信号估算法,快速、准确进行测量。具体方法如下:(1)采用射频接收前端进行射频信号变频、滤波,将X波段信号变换到中频;(2)将中频信号进行AD变换,将变换后的结果送入数字信号处理器件,进行下变频处理;(3)将中频数字信号变换成基带数字信号,送入FPGA进行计算处理,完成频率测量。进一步地,使用频带压缩法将X波段信号变换到中频。采用下变频技术降低数据流速率,将中频信号变换为基带信号。采用FFT和信号平方的双信号频率估计算法,保证测频的快速性和准确性。 
尽管为示例目的,已经公开了本发明的优选实施例,本领域的技术人员将意识到各种改进、增加和取代也是可能的,因此,本发明的范围应当不限于上述实施例。 

Claims (10)

1.一种高精度数字瞬时测频方法,其特征在于,所述方法包括:
采用射频接收前端进行射频信号变频、滤波处理,将X波段信号变换为中频信号;
将所述中频信号进行模数AD变换,通过数字信号处理器件将变换后的中频信号进行下变频处理,得到中频数字信号;
将所述中频数字信号变换为基带数字信号;采用基于快速傅氏变换FFT和信号平方的双信号估算法,通过现场可编程门阵列FPGA对所述基带数字信号进行计算处理,得到测频结果。
2.如权利要求1所述的方法,其特征在于,通过所述FPGA对所述基带数字信号进行计算处理,包括以下步骤:
所述FPGA根据同步触发对所述基带数字信号进行采集,并采用快进慢出的方式读取所述信号的数据;
对读取的所述数据进行有限长单位冲激响应滤波器FIR带通滤波处理,并进行缓存;
将滤波后的数据进行快速傅里叶变换,得到相应的频谱信息,去除所述频谱信息中的镜频部分,保留正频率部分,提取最大值为频谱值,以及其位置;
根据采样率,频谱值,以及其位置进行映射运算得到相应的频率值,确定频率的位置;
采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果。
3.如权利要求2所述的方法,其特征在于,采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果,包括以下步骤:
根据信号中两频率分量的幅度比,通过一维平方信号谱峰搜索,得到所述信号的和频与差频分量的估计值;并利用插值技术得到频率结果。
4.如权利要求3所述的方法,其特征在于,在得到所述频率结果后,所述方法还包括:
按照约定的通讯协议进行频率转换,将转换后的频率结果输出。
5.如权利要求1至4中任一项所述的方法,其特征在于,所述方法还包括:
在测频过程中进行状态管理和故障上报。
6.一种高精度数字瞬时测频装置,其特征在于,所述装置包括:
中频变换模块,用于采用射频接收前端进行射频信号变频、滤波处理,将X波段信号变换为中频信号;
模数变换模块,用于将所述中频信号进行模数AD变换,通过数字信号处理器件将变换后的中频信号进行下变频处理,得到中频数字信号;
测频模块,用于将所述中频数字信号变换为基带数字信号;采用基于快速傅氏变换FFT和信号平方的双信号估算法,通过现场可编程门阵列FPGA对所述基带数字信号进行计算处理,得到测频结果。
7.如权利要求6所述的装置,其特征在于,所述测频模块包括:
数据采集单元,用于根据同步触发对所述基带数字信号进行采集,并采用快进慢出的方式读取所述信号的数据;
滤波缓存单元,用于对读取的所述数据进行有限长单位冲激响应滤波器FIR带通滤波处理,并进行缓存;
时频变换单元,用于将滤波后的数据进行快速傅里叶变换,得到相应的频谱信息,去除所述频谱信息中的镜频部分,保留正频率部分,提取最大值为频谱值,以及其位置;
频率计算单元,用于根据采样率,频谱值,以及其位置进行映射运算得到相应的频率值,确定频率的位置;
双信号估算单元,用于采用基于快速傅氏变换FFT和信号平方的双信号估算法,得到频率结果。
8.如权利要求7所述的装置,其特征在于,所述双信号估算单元包括:
双信号估算子单元,用于根据信号中两频率分量的幅度比,通过一维平方信号谱峰搜索,得到所述信号的和频与差频分量的估计值;并利用插值技术得到频率结果。
9.如权利要求8所述的装置,其特征在于,所述装置还包括:
结果输出模块,用于按照约定的通讯协议进行频率转换,将转换后的频率结果输出。
10.如权利要求6至9中任一项所述的装置,其特征在于,所述装置还包括:
管理模块,用于在测频过程中进行状态管理和故障上报。
CN201310728431.5A 2013-12-25 2013-12-25 一种高精度数字瞬时测频方法及装置 Pending CN103837740A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310728431.5A CN103837740A (zh) 2013-12-25 2013-12-25 一种高精度数字瞬时测频方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310728431.5A CN103837740A (zh) 2013-12-25 2013-12-25 一种高精度数字瞬时测频方法及装置

Publications (1)

Publication Number Publication Date
CN103837740A true CN103837740A (zh) 2014-06-04

Family

ID=50801445

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310728431.5A Pending CN103837740A (zh) 2013-12-25 2013-12-25 一种高精度数字瞬时测频方法及装置

Country Status (1)

Country Link
CN (1) CN103837740A (zh)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104104635A (zh) * 2014-07-15 2014-10-15 中国电子科技集团公司第四十一研究所 一种比相法瞬时测频编码方法
CN105486920A (zh) * 2015-11-27 2016-04-13 北京振兴计量测试研究所 一种基于数字信道化的窄脉冲频率测量方法
CN105891599A (zh) * 2016-03-31 2016-08-24 福州大学 一种基于改进的dodf-wspd的电网频率跟踪方法
CN106027156A (zh) * 2016-04-29 2016-10-12 成都卓力致远科技有限公司 一种基于光学模数转换的微波信号测频方法及装置
CN106405221A (zh) * 2016-08-17 2017-02-15 中国电子科技集团公司第四十研究所 一种快速高精度测量信号频率功率参数的装置和方法
CN106841776A (zh) * 2017-03-28 2017-06-13 北京润科通用技术有限公司 一种测频机及测量信号频率的方法
CN107976572A (zh) * 2017-12-18 2018-05-01 中国科学院深圳先进技术研究院 一种信道化高精度实时频谱分析方法及其***
CN108667484A (zh) * 2018-03-26 2018-10-16 西南电子技术研究所(中国电子科技集团公司第十研究所) 非相干扩频数字收发信机瞬时测频与解调方法
WO2018188228A1 (zh) * 2017-04-13 2018-10-18 中国电子科技集团公司第二十四研究所 高精度频率测量***及方法
CN108777599A (zh) * 2018-08-15 2018-11-09 上海频语电子技术有限公司 一种微波宽带瞬时测频接收***
CN109194598A (zh) * 2018-11-16 2019-01-11 上海工程技术大学 一种通用psk调制解调***
CN110794208A (zh) * 2019-11-14 2020-02-14 中电科仪器仪表有限公司 一种基于数字相位推算的宽带快速调制域测频方法
CN112946580A (zh) * 2021-01-14 2021-06-11 无锡国芯微电子***有限公司 一种多处理器协同辐射源频率参数估计装置及方法
CN113075450A (zh) * 2021-02-22 2021-07-06 中国电子科技集团公司第二十九研究所 基于采样率解析宽带频率压缩信号射频和中频频率的方法
CN113406386A (zh) * 2021-06-23 2021-09-17 中国电子科技集团公司第二十九研究所 一种基于数字下变频的信号频率精确估计方法
CN113884761A (zh) * 2021-09-13 2022-01-04 中国电子科技集团公司第二十九研究所 一种提高测频精度的方法、装置、设备及介质
CN113965217A (zh) * 2021-02-01 2022-01-21 湖南国科雷电子科技有限公司 一种双通道、单通道s/c/x三波段宽带单比特数字测频接收机
CN114034927A (zh) * 2021-11-01 2022-02-11 南京国电南自电网自动化有限公司 一种跟频插值采样的信号测量方法及***
CN114184837A (zh) * 2021-12-09 2022-03-15 电子科技大学 一种基于Cordic算法的瞬时测频方法
CN116879627A (zh) * 2023-09-04 2023-10-13 中国电子科技集团公司第二十九研究所 一种纳秒级非相参窄脉冲序列的测频***

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1804642A (zh) * 2006-01-23 2006-07-19 天津市德力电子仪器有限公司 中频全数字化频谱仪
CN201464557U (zh) * 2009-03-16 2010-05-12 邵啸 一种中频频谱监测装置
US20120038369A1 (en) * 2010-08-13 2012-02-16 Tektronix, Inc. Time-domain measurements in a test and measurement instrument
JP2012042462A (ja) * 2010-08-13 2012-03-01 Tektronix Inc 試験測定機器及びトリガ方法
CN103067104A (zh) * 2012-12-27 2013-04-24 上海创远仪器技术股份有限公司 基于数字本振对射频信号高速扫频频谱测量的***及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1804642A (zh) * 2006-01-23 2006-07-19 天津市德力电子仪器有限公司 中频全数字化频谱仪
CN201464557U (zh) * 2009-03-16 2010-05-12 邵啸 一种中频频谱监测装置
US20120038369A1 (en) * 2010-08-13 2012-02-16 Tektronix, Inc. Time-domain measurements in a test and measurement instrument
JP2012042462A (ja) * 2010-08-13 2012-03-01 Tektronix Inc 試験測定機器及びトリガ方法
CN103067104A (zh) * 2012-12-27 2013-04-24 上海创远仪器技术股份有限公司 基于数字本振对射频信号高速扫频频谱测量的***及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张昌菊: "双信号快速测频技术及FPGA实现", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104104635B (zh) * 2014-07-15 2017-12-01 中国电子科技集团公司第四十一研究所 一种比相法瞬时测频编码方法
CN104104635A (zh) * 2014-07-15 2014-10-15 中国电子科技集团公司第四十一研究所 一种比相法瞬时测频编码方法
CN105486920A (zh) * 2015-11-27 2016-04-13 北京振兴计量测试研究所 一种基于数字信道化的窄脉冲频率测量方法
CN105891599A (zh) * 2016-03-31 2016-08-24 福州大学 一种基于改进的dodf-wspd的电网频率跟踪方法
CN106027156A (zh) * 2016-04-29 2016-10-12 成都卓力致远科技有限公司 一种基于光学模数转换的微波信号测频方法及装置
CN106027156B (zh) * 2016-04-29 2019-02-05 成都卓力致远科技有限公司 一种基于光学模数转换的微波信号测频方法及装置
CN106405221B (zh) * 2016-08-17 2019-01-18 中国电子科技集团公司第四十一研究所 一种快速高精度测量信号频率功率参数的装置和方法
CN106405221A (zh) * 2016-08-17 2017-02-15 中国电子科技集团公司第四十研究所 一种快速高精度测量信号频率功率参数的装置和方法
CN106841776A (zh) * 2017-03-28 2017-06-13 北京润科通用技术有限公司 一种测频机及测量信号频率的方法
CN106841776B (zh) * 2017-03-28 2019-08-06 北京润科通用技术有限公司 一种测频机及测量信号频率的方法
WO2018188228A1 (zh) * 2017-04-13 2018-10-18 中国电子科技集团公司第二十四研究所 高精度频率测量***及方法
CN107976572A (zh) * 2017-12-18 2018-05-01 中国科学院深圳先进技术研究院 一种信道化高精度实时频谱分析方法及其***
CN108667484A (zh) * 2018-03-26 2018-10-16 西南电子技术研究所(中国电子科技集团公司第十研究所) 非相干扩频数字收发信机瞬时测频与解调方法
CN108667484B (zh) * 2018-03-26 2020-03-31 西南电子技术研究所(中国电子科技集团公司第十研究所) 非相干扩频数字收发信机瞬时测频与解调方法
CN108777599A (zh) * 2018-08-15 2018-11-09 上海频语电子技术有限公司 一种微波宽带瞬时测频接收***
CN109194598B (zh) * 2018-11-16 2020-12-11 上海工程技术大学 一种通用psk调制解调***
CN109194598A (zh) * 2018-11-16 2019-01-11 上海工程技术大学 一种通用psk调制解调***
CN110794208A (zh) * 2019-11-14 2020-02-14 中电科仪器仪表有限公司 一种基于数字相位推算的宽带快速调制域测频方法
CN112946580A (zh) * 2021-01-14 2021-06-11 无锡国芯微电子***有限公司 一种多处理器协同辐射源频率参数估计装置及方法
CN113965217B (zh) * 2021-02-01 2022-12-30 湖南国科雷电子科技有限公司 一种双通道、单通道s/c/x三波段宽带单比特数字测频接收机
CN113965217A (zh) * 2021-02-01 2022-01-21 湖南国科雷电子科技有限公司 一种双通道、单通道s/c/x三波段宽带单比特数字测频接收机
CN113075450A (zh) * 2021-02-22 2021-07-06 中国电子科技集团公司第二十九研究所 基于采样率解析宽带频率压缩信号射频和中频频率的方法
CN113075450B (zh) * 2021-02-22 2023-04-25 中国电子科技集团公司第二十九研究所 基于采样率解析宽带频率压缩信号射频和中频频率的方法
CN113406386A (zh) * 2021-06-23 2021-09-17 中国电子科技集团公司第二十九研究所 一种基于数字下变频的信号频率精确估计方法
CN113884761A (zh) * 2021-09-13 2022-01-04 中国电子科技集团公司第二十九研究所 一种提高测频精度的方法、装置、设备及介质
CN113884761B (zh) * 2021-09-13 2023-04-25 中国电子科技集团公司第二十九研究所 一种提高测频精度的方法、装置、设备及介质
CN114034927A (zh) * 2021-11-01 2022-02-11 南京国电南自电网自动化有限公司 一种跟频插值采样的信号测量方法及***
CN114184837A (zh) * 2021-12-09 2022-03-15 电子科技大学 一种基于Cordic算法的瞬时测频方法
CN114184837B (zh) * 2021-12-09 2022-10-18 电子科技大学 一种基于Cordic算法的瞬时测频方法
CN116879627A (zh) * 2023-09-04 2023-10-13 中国电子科技集团公司第二十九研究所 一种纳秒级非相参窄脉冲序列的测频***
CN116879627B (zh) * 2023-09-04 2023-11-21 中国电子科技集团公司第二十九研究所 一种纳秒级非相参窄脉冲序列的测频***

Similar Documents

Publication Publication Date Title
CN103837740A (zh) 一种高精度数字瞬时测频方法及装置
CN101523224B (zh) 实时示波器上的实时谱触发***
CN104122444B (zh) 全数字中频频谱分析仪及频谱分析方法
CN105512431A (zh) 一种基于相位噪声数学模型的相位噪声测量方法
CN106546817B (zh) 一种具有反馈功能的频率估计和能量估计方法
CN107239611B (zh) 一种矢量信号分析装置及方法
CN104035012B (zh) 一种局部放电特高频信号检测的调理电路
CN111884984B (zh) 快变载波多普勒频移捕获***
CN109407501B (zh) 一种基于相关信号处理的时间间隔测量方法
CN102868402A (zh) 一种测试模数转换器主要性能指标的测试方法
CN106027179A (zh) 一种基于综合互素分析的宽带频谱感知方法及其装置
Song et al. Real-time multi-gigahertz sub-Nyquist spectrum sensing system for mmwave
CN102495280B (zh) 一种抗噪音宽带频率测量方法及锁相频率计
CN114531329B (zh) 一种多路msk信号载波频率估计方法、***及应用
CN113242114B (zh) 一种猝发式窄脉冲信号的检测***及方法
CN108900445B (zh) 一种信号符号率估计的方法及装置
CN103684464A (zh) 一种相关型微波辐射计中频信号欠采样处理方法
CN110261830A (zh) 一种针对雷达回波模拟器的性能校准器
EP3055704B1 (en) Analog to information converter
CN103902133A (zh) 电磁触控接收装置和方法
CN103905360A (zh) 加入“判极”操作的非合作bpsk信号解码方法
CN106777505A (zh) 基于频偏识别的欠采样信号的鲁棒的频率估计方法及装置
Qiu et al. NYFR output pulse radar signal TOA analysis using extended Fourier transform and its TOA estimation
CN116360235A (zh) 一种基于SerDes的TDC实现装置
CN105743830B (zh) 一种mfsk信号的解调方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140604