CN105510706B - 一种高精度欠采样测频方法 - Google Patents
一种高精度欠采样测频方法 Download PDFInfo
- Publication number
- CN105510706B CN105510706B CN201511021008.7A CN201511021008A CN105510706B CN 105510706 B CN105510706 B CN 105510706B CN 201511021008 A CN201511021008 A CN 201511021008A CN 105510706 B CN105510706 B CN 105510706B
- Authority
- CN
- China
- Prior art keywords
- frequency
- channel
- signal
- fin
- initial phase
- 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
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
一种高精度欠采样测频方法,首先对信号进行信道折叠,将输入信号根据按照采样频率均匀折叠在多个信道通路上,随后通过相位控制将各信道的信号按照傅里叶变换的计算长度对信道信号进行截取,通过使用不同采样频率的AD进行采样得到各个通道的采样值,然后对各个通道采集信号进行傅里叶变换,得到傅里叶变换结果后分别通过相位判断得到对应谱线位置,再使用谱线校正方法对谱线结果进行校正,得到各个通道频率结果,最后得到各个信道频率进而得到最后的估计频率。本发明方法通过使用多个低速AD替代了原有的需要高速AD的场合,相比直接使用这种位数的高速AD,提高了采样的分辨率。
Description
技术领域
本发明涉及一种频率测量技术,特别是一种高精度欠采样测频方法。
背景技术
传统的高频采样技术主要依赖于高速AD的技术以及信道化技术对输入信号进行处理,这种处理方式依赖与AD的发展水平,高速AD通常采样位数较低,无法获得高分辨率的信号,不适用于高频采样。
反辐射导引头等被动接收系统中通常需要对输入频率进行有效的估计之后进行相关操作,因为受制于采样定律(被采样信号的频率需要低于采样信号频率的2倍)以及AD的采样速率,所以在外界输入未知频率的高频信号时候,接收系统对于高频信号的估计会有困难,产生估计错误或者估计缓慢等情况。这是因为接收系统首先需要对信号进行准确的范围估计,通过变频完成将输入信号变换到AD可以采集的频率域内的操作,最后通过结合AD的计算结果和变频参数才可以对原信号的频率参数进行估计。
在实际使用过程中,通常无法准确地对输入信号范围进行估计,因此,需要多次对信号进行变频操作来确认输入信号的范围,如果变频操作不恰当还会产生错误的结果。而且因为变频没有目标指向性,在变频操作过程中会耗费大量时间。同时,根据采样定律,若想对较大范围频率信号进行采样,需要使用高速的AD已完成相关的频率估计。首先,高速的AD会带来位数的降低,使采样信号的分辨率有所降低;其次,受电子制造技术影响,采购高速的AD会导致成本会指数上升,是设备的成本大幅增多;随后,当前最高速级别的AD采样速率基本集中在5G以内,而对应的输入信号范围已覆盖几十G;最后,高速AD会带来布板布线的困难,EMI/EMC极难控制,增加了设计及制作难度。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种通过使用多个低速AD电路代替高速AD电路的提高了采样分辨率的高精度欠采样测频方法。
本发明的技术解决方案是:一种高精度欠采样测频方法,包括如下步骤:
(1)选取M个AD电路,并将M个AD电路的采样频率分别记为f1,f2,f3,…,fM;
(2)将频率为fin的待采样信号按照fin_i=fin-n*fset/N的方式划分为N个信道分配信号,第i个信道分配信号的频率为fin_i,其中,n=0,1,2,…,N-1,i=1,2,3,…,N,N为正整数,fset的取值范围为[0,lcm(f1,f2,…,fM)],lcm表示最小公倍数;
(3)对第i个信道分配信号进行相位控制,令第i个信道分配信号初始相位为θi,其中,i为正整数且i的初值为1;
(4)使用M个AD电路对第i个信道分配信号进行采样得到采样数据,并记为Di1、Di2、Di3,…,DiM,对采样数据Di1、Di2、Di3,…,DiM进行傅里叶变换得到傅里叶变换幅度谱、傅里叶变换相位谱,进而得到采样数据Di1,Di2,Di3,…,DiM的频率集合fi1,fi2,fi3,…,fiM,初始相位集合θi1,θi2,θi3,…,θiM,其中,fi1,fi2,fi3,…,fiM和θi1,θi2,θi3,…,θiM均为二维数组;
(5)对频率集合fi1,fi2,fi3,…,fiM和初始相位集合θi1,θi2,θi3,…,θiM分别使用内插法进行校正得到校正频率集合f’i1,f’i2,f’i3,…,f’iM和校正初始相位集合θ’i1,θ’i2,θ’i3,…,θ’iM;
(6)使用信道初始输入相位θi依次筛选校正初始相位集合θ’i1,θ’i2,θ’i3,…,θ’iM,得到校正初始相位集合中每个校正初始相位的二维数组中与θi最接近的1个校正初始相位数值,进而得到校正频率集合中M个校正初始相位数值分别对应的频率得到采样数据Di1,Di2,Di3,…,DiM的频率集合分别为其中,k,m,p,l均为自然数;
(7)在频率集合中的每个频率集中寻找与其他频率集中频率相同或最接近的1个频率,进而得到M个频率,并分别记为进而得到第i个信道频率fin_i为
其中,包括M个频率集;
(8)i=i+1,重复步骤(3)至步骤(7)直至i=N,得到N个信道的信道频率fin_i;
(9)分别获取第i个信道的信道信噪比SNRi,进而得到待采样信号的频率fin为
所述的M≥2。
所述的θi为0°。
所述的M=3。
本发明与现有技术相比的优点在于:
(1)本发明方法通过使用多个低速AD替代了原有的需要高速AD的试验场景,相比直接使用这种位数的高速AD,提高了采样的分辨率;
(2)本发明方法通过使用内插校正算法校正了傅里叶变换的栅栏效应,相比现有的直接采样计算有着更为精确的频率结果;
(3)本发明方法克服了现有技术因高速器件成本随电路频率的升高而快速升高的缺陷,通过使用多个低速AD大大降低了测频成本和硬件设计成本;
(4)本发明方法通过利用信噪比的方式完成频率的估计,与现有的直接高速AD采样技术相比,能够获得与信噪比相关的频率估计。
附图说明
图1为本发明一种高精度欠采样测频方法原理流程图。
具体实施方式
本发明针对现有的测频技术不能确定高频信号频率的不足,通过使用多个特定频率的低速采样AD完成对超出采样定理频率信号的频率确定工作,解决了某些适用场景中需要对高频信号频率进行精确测量的问题,本发明高精度欠采样测频方法首先对信号进行信道折叠,将输入信号根据按照采样频率均匀折叠在多个信道通路上,随后通过相位控制将各信道的信号按照FFT(快速傅里叶变换)的计算长度对信道信号进行截取,通过使用不同采样频率的AD进行采样得到各个通道的采样值,再对各个通道采集信号进行FFT/DFT(离散傅里叶变换),得到傅里叶变换结果后分别通过相位判断得到对应谱线位置,随后使用candan或Jacobsen等谱线校正方法对谱线结果进行校正,得到各个通道频率结果,最后得到各个信道频率进而得到最后的估计频率,下面结合附图对本发明方法进行详细说明。
如图1所示,本发明测频方法包括如下步骤:
步骤1:将待采样信号(S(t),其频率为fin)输入到按照设定频率(fset,fset<lcm(f1,f2,f3),lcm表示最小公倍数)均匀划分为N个信道的采集系统中,每个信道分别进行折叠变频处理(fin_i=fin-n*fset/N,n=0,1..N-1,i=1..N)后得到N个分配信道信号(Si表示,其频率为fin_i,i=1..N),其中,N为正整数,f1、f2、f3分别为本发明选取的3个AD的采样频率,一般不少于2个AD,本发明以3个为例。
步骤2:在第1个信道中对分配信号进行输入相位控制,保证信道分配信号的初始输入相位为确定值(θi,一般保证初始相位为0°)。其他信道均为同样操作,此后步骤3~5,均以第1个信道为例,所述操作在每一信道中均为同样的操作过程。
步骤3:对信道分配信号使用多个不同频率的AD进行采样得到采样数据,并记为D11、D12、D13,对采样数据分别进行傅里叶变换得到对应的幅度谱、相位谱,进而得到采样数据的主频率幅度集合f11、f12、f13,采样数据的主频率初始相位集合θ11、θ12、θ13,其中,f11、f12、f13、θ11、θ12、θ13均为二维数组。
步骤4:对傅里叶变换后的结果f11、f12、f13和θ11、θ12、θ13使用内插法进行校正,计算得到校正结果集f11’、f12’、f13’和θ11’、θ12’、θ13’,然后利用信道初始输入相位θ1对照θ11’,θ12’,θ13’获得真实值的集合坐标(即分别从θ11’或者θ12’或者θ13’的二维数组中挑选与θ1最接近的数值),将坐标结果带入频率集合f11’、f12’、f13’获得对应的D11、D12、D13的频率结果因此可以构建D11、D12、D13的频率谱分别为:k、m、p为自然数。
步骤5:根据中国余数定理对频率谱进行整合,即在频率谱中寻找最小偏差相近值,分别为因频率精度与AD采
样频率正相关,所以用AD采样频率构建权值计算信道频率的输出值,其计算式为
步骤6:将各个信道频率利用各个信道信噪比SNRi作为权值对各信道频率进行融合,其融合计算式为完成欠采样测频。
下面结合实施例对本发明方法进行详细说明,假设某个信道内信道搬移得到的输入信号的频率为31354,输入信号通过相位控制模块对输入相位进行控制。然后在各分支通道中选取AD采样速率为f1=30K、f2=40K和f3=5K。随后将采样的N=10000个数据进行直接FFT计算,可以得到对应的FFT结算频率值,分别为[8644,31356],[1353,28647],[2255,47745];随后使用Candan内插方法对结果进行估计和修正,得到对应的修正结果为[8646.0,31354.0],[1354.0,28646.0],[18646,31354.0],随后通过频谱矫正的方法对频谱值进行选取,得到各采样AD对应的实际结果:31354.0,1354.0,31354.0,随后通过余数定理对频率进行筛选和权值计算,得到单一通道的频率估计结果31354Hz,各通道结果可以通过信噪比构建的权值进行频率融合结算,最后得到对应的频率结果。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (3)
1.一种高精度欠采样测频方法,其特征在于包括如下步骤:
(1)选取M个AD电路,并将M个AD电路的采样频率分别记为f1,f2,f3,…,fM,其中,M≥2;
(2)将频率为fin的待采样信号按照fin_i=fin-n*fset/N的方式划分为N个信道分配信号,第i个信道分配信号的频率为fin_i,其中,n=0,1,2,…,N-1,i=1,2,3,…,N,N为正整数,fset的取值范围为[0,lcm(f1,f2,…,fM)],lcm表示最小公倍数;
(3)对第i个信道分配信号进行相位控制,令第i个信道分配信号初始相位为θi,其中,i为正整数且i的初值为1;
(4)使用M个AD电路对第i个信道分配信号进行采样得到采样数据,并记为Di1、Di2、Di3,…,DiM,对采样数据Di1、Di2、Di3,…,DiM进行傅里叶变换得到傅里叶变换幅度谱、傅里叶变换相位谱,进而得到采样数据Di1,Di2,Di3,…,DiM的频率集合fi1,fi2,fi3,…,fiM,初始相位集合θi1,θi2,θi3,…,θiM,其中,fi1,fi2,fi3,…,fiM和θi1,θi2,θi3,…,θiM均为二维数组;
(5)对频率集合fi1,fi2,fi3,…,fiM和初始相位集合θi1,θi2,θi3,…,θiM分别使用内插法进行校正得到校正频率集合f’i1,f’i2,f’i3,…,f’iM和校正初始相位集合θ’i1,θ’i2,θ’i3,…,θ’iM;
(6)使用信道初始输入相位θi依次筛选校正初始相位集合θ’i1,θ’i2,θ’i3,…,θ’iM,得到校正初始相位集合中每个校正初始相位的二维数组中与θi最接近的1个校正初始相位数值,进而得到校正频率集合中M个校正初始相位数值分别对应的频率得到采样数据Di1,Di2,Di3,…,DiM的频率集合分别为其中,k,m,p,l均为自然数;
(7)在频率集合中的每个频率集中寻找与其他频率集中频率相同或最接近的1个频率,进而得到M个频率,并分别记为进而得到第i个信道频率fin_i为
其中,包括M个频率集;
(8)i=i+1,重复步骤(3)至步骤(7)直至i=N,得到N个信道的信道频率fin_i;
(9)分别获取第i个信道的信道信噪比SNRi,进而得到待采样信号的频率fin为
2.根据权利要求1所述的一种高精度欠采样测频方法,其特征在于:所述的θi为0°。
3.根据权利要求1或2所述的一种高精度欠采样测频方法,其特征在于:所述的M=3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511021008.7A CN105510706B (zh) | 2015-12-30 | 2015-12-30 | 一种高精度欠采样测频方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511021008.7A CN105510706B (zh) | 2015-12-30 | 2015-12-30 | 一种高精度欠采样测频方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105510706A CN105510706A (zh) | 2016-04-20 |
CN105510706B true CN105510706B (zh) | 2018-12-14 |
Family
ID=55718827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511021008.7A Active CN105510706B (zh) | 2015-12-30 | 2015-12-30 | 一种高精度欠采样测频方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510706B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3896465A4 (en) * | 2019-01-22 | 2022-01-12 | Mitsubishi Electric Corporation | FREQUENCY DETECTION CIRCUIT |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107192877B (zh) * | 2017-05-19 | 2019-08-02 | 上海航天控制技术研究所 | 一种基于快速傅里叶变换的相位判读方法 |
CN107561357B (zh) * | 2017-07-05 | 2020-08-14 | 中国电子科技集团公司第三十八研究所 | 一种基于信道化的高精度瞬时测频方法和装置 |
CN108776261B (zh) * | 2018-05-16 | 2020-05-12 | 电子科技大学 | 一种大带宽跳频窄带信号多通道欠采样测频方法 |
CN112129983B (zh) * | 2020-09-25 | 2024-05-10 | 中北大学 | 一种基于等时间间隔等效取样的波形恢复数据处理方法 |
CN113138313B (zh) * | 2021-03-19 | 2022-04-08 | 北京航空航天大学 | 一种基于余数匹配的光学欠采样频率恢复方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630148A (zh) * | 2013-11-01 | 2014-03-12 | 中国科学院物理研究所 | 信号取样平均仪和信号取样平均方法 |
CN103941087A (zh) * | 2014-04-09 | 2014-07-23 | 天津大学 | 欠采样速率下的高频余弦信号的频率测量方法及其装置 |
CN104007316A (zh) * | 2014-05-29 | 2014-08-27 | 天津大学 | 一种欠采样速率下的高精度频率测量方法及其测量仪 |
-
2015
- 2015-12-30 CN CN201511021008.7A patent/CN105510706B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103630148A (zh) * | 2013-11-01 | 2014-03-12 | 中国科学院物理研究所 | 信号取样平均仪和信号取样平均方法 |
CN103941087A (zh) * | 2014-04-09 | 2014-07-23 | 天津大学 | 欠采样速率下的高频余弦信号的频率测量方法及其装置 |
CN104007316A (zh) * | 2014-05-29 | 2014-08-27 | 天津大学 | 一种欠采样速率下的高精度频率测量方法及其测量仪 |
Non-Patent Citations (3)
Title |
---|
信噪比定量估计及其在加权叠加中的应用;刘文辉等;《吐哈油气》;20041231;第9卷(第4期);369-374 * |
基于中国余数定理的欠采样下余弦信号的频率估计;黄翔东等;《物理学报》;20141015;第63卷(第19期);198403-3,198403-5,198403-6 * |
鲁棒的闭式中国余数定理及其在欠采样频率估计中的应用;王文杰等;《信号处理》;20130930;第29卷(第9期);1206-1211 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3896465A4 (en) * | 2019-01-22 | 2022-01-12 | Mitsubishi Electric Corporation | FREQUENCY DETECTION CIRCUIT |
Also Published As
Publication number | Publication date |
---|---|
CN105510706A (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105510706B (zh) | 一种高精度欠采样测频方法 | |
CN103455811B (zh) | 室内无线定位指纹采集方法及系统 | |
CN108414833B (zh) | 一种信号分量频率的精确估计方法 | |
CN103457603B (zh) | 一种基于平均频谱测试adc动态参数的方法 | |
CN109214371B (zh) | 一种基于快速傅氏变换的正弦信号滤波系统 | |
CN105787219B (zh) | 一种利用临近频点采样建立传导干扰耦合通道多元线性回归模型的方法 | |
CN104833937A (zh) | 一种基于mir-rsd高精度余弦窗插值fft算法的谐波测量通道校准方法 | |
CN107144738B (zh) | 一种基于直通线的多端口自动夹具损耗和相位补偿方法 | |
CN114025379B (zh) | 一种宽带多信号检测方法、装置和设备 | |
CN105137373B (zh) | 一种指数信号的去噪方法 | |
CN103543426A (zh) | 一种网络分析仪分波段校准内插补偿方法 | |
CN104378085A (zh) | 一种射频阻抗匹配调试方法 | |
CN105373708B (zh) | 一种基于参数优化的改进广义s变换的时频分析方法 | |
CN109444547A (zh) | 基于二端口网络的rfid芯片阻抗测量方法及装置 | |
CN106645942A (zh) | 一种低成本高精度嵌入式信号采集分析系统和方法 | |
CN106597392B (zh) | 一种sar系统相位误差补偿方法 | |
CN109655775A (zh) | 一种任意波形发生器幅度扫频多尺度校准方法和装置 | |
CN109583575B (zh) | 基于深度学习提高仪器矢量信号分析性能的处理方法 | |
CN101308175A (zh) | 相位谱分析仪 | |
CN106124081A (zh) | 永磁同步电机精确多点实时测温方法及系统 | |
CN111082834B (zh) | 一种基于啁啾信号正交解调的射频时延快速测量装置 | |
CN111611686A (zh) | 一种通信信号时频域的检测方法 | |
CN103092815A (zh) | 校准装置及其对监测装置中的传递函数进行校准的方法 | |
CN107203495B (zh) | 一种s参数圆形插值方法 | |
CN110716167A (zh) | 任意波形发生器幅度扫频自适应分辨率校准方法及系统 |
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 |