CN103234624A - 一种高精度的频率估计系统 - Google Patents

一种高精度的频率估计系统 Download PDF

Info

Publication number
CN103234624A
CN103234624A CN2012102949498A CN201210294949A CN103234624A CN 103234624 A CN103234624 A CN 103234624A CN 2012102949498 A CN2012102949498 A CN 2012102949498A CN 201210294949 A CN201210294949 A CN 201210294949A CN 103234624 A CN103234624 A CN 103234624A
Authority
CN
China
Prior art keywords
frequency
signal
input
frequency estimation
module
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
Application number
CN2012102949498A
Other languages
English (en)
Other versions
CN103234624B (zh
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.)
Sundak Semiconductor Technology Shanghai Co ltd
Original Assignee
Suzhou University
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 Suzhou University filed Critical Suzhou University
Priority to CN201210294949.8A priority Critical patent/CN103234624B/zh
Publication of CN103234624A publication Critical patent/CN103234624A/zh
Application granted granted Critical
Publication of CN103234624B publication Critical patent/CN103234624B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提出了一种高精度的频率估计系统，由信号采集模块，信号处理模块，信号显示和存储模块，配置模块四个部分组成，其中，信号处理模块读取信号采集模块传送的数据，根据配置模块配置的信息按公式进行输入信号的频率估计，首先将输入数字信号利用参考频率进行
Figure DEST_PATH_IMAGE002
和
Figure DEST_PATH_IMAGE004
两个正交分量处理，并对正交分量
Figure DEST_PATH_IMAGE002A
和进行低通滤波，然后求其一周期相位差。多次测量求其相位差的均值，进而高精度的估计输入单音频率，并将估计的结果输出到显示和存储模块。本发明的高精度的频率估计系统，此数字方式频率估计方法，只需较小的时间复杂度和空间复杂度，理论上可以达到任意精度。

Description

一种高精度的频率估计系统
技术领域
本发明涉及一种用于估计高精度的频率信号的系统。
 
背景技术
在实际应用中，常常需要对一些信号的频率进行估计，用于单频率信号的检测。例如在ITU-T制订的V.34数据传输协议的调制解调器应用时，需要对单音信号进行检测估计，并根据单音信号特征进行信道线路探测。而在实际生活中，高精度的频率估计方法可以应用于钢琴等按键乐器的调音等应用中，对按键音采集后进行频率估计，辅助调音师进行准确调音。
 
发明内容
为了解决上述问题，本发明提供一种高精度的频率估计系统。
为实现上述发明目的，本发明所采用的技术方案是：
本发明是一种高精度的频率估计系统，所述频率估计系统由信号采集模块，信号处理模块，信号显示和存储模块，配置模块四个部分组成。
信号采集模块：主要实现单音信号的采集，支持模拟信号输入采集和数字信号输入采集两种模式，将采集的数据传送给信号处理模块进行频率估计；
信号处理模块：读取信号采集模块传送的数据，根据配置模块配置的信息进行输入信号的频率估计，首先将输入数字信号利用参考频率进行                                                
Figure 696447DEST_PATH_IMAGE002
和
Figure 153973DEST_PATH_IMAGE004
两个正交分量处理，并对正交分量和
Figure 843711DEST_PATH_IMAGE004
进行低通滤波，然后求其一周期相位差；多次测量求其相位差的均值，进而高精度的估计输入单音频率，并将估计的结果输出到显示和存储模块；
信号显示和存储模块：用于存储和显示信号处理模块输出的频率估计结果；
配置模块：主要包括参考频率设置和特定功能配置的功能。
优选的，所述采集数据的输入，输入接口包括模拟输入接口和数字输入接口两种模式，模拟输入接口接收来自线路上或麦克风的单音信号，经AD转换后转化为数字信号；数字接口可直接接收数字线路上的数字信号，两种模式的输入通过多路选择开关MUX进行功能选择，通过IO接口输送给处理器。
优选的，所述频率估计系统的处理步骤为，处理器在系统开始工作时首先读取存储器中存储的配置信息，然后进入相应的信号处理模块对输入的数字信号进行高精度频率估计，处理器可以通过显示接口输出频率估计结果，可以通过存储接口输送到存储器。
优选的，所述配置模块中，采用参考频率设置时，系统根据设置的参考频率进行高精度频率估计；采用特定功能配置时，系统将主要对特定功能所要求的频率进行估计，其方法是对于输入样本将该组所有频率依次用做参考频率
Figure 2012102949498100002DEST_PATH_IMAGE007
，计算输入单音与参考频率之差，选输入单音与参考频率最小的
Figure DEST_PATH_IMAGE011
估计输入单音频率。
相对于现有技术，本发明的高精度的频率估计系统，具有以下有益的技术效果：
1.本发明的高精度频率估计方法，理论上可以达到任意精度。
2.采用数字方式频率估计方法，需要较小的时间复杂度和空间复杂度。
3.与幅度无关的单音信号频率估计方法，幅度变化不会影响频率估计的精度。
4.支持两种频率估计模式，包括单频率估计方式和特定功能频率估计模式（多频率同时进行估计模式）；
5.支持模拟输入和数字输入两种信号输入。
附图说明
图1为本发明的一种高精度的频率估计系统的系统方案示意图。
图2为本发明的一种高精度的频率估计系统的硬件原理框图。
 
具体实施方式
以下结合附图详细描述本发明的高精度的频率估计系统的体系结构，但不构成对本发明的限制。
1.高精度的频率估计方法理论基础
假设输入的单音信号频率为
Figure DEST_PATH_IMAGE013
，相位为
Figure DEST_PATH_IMAGE015
，幅值为
Figure DEST_PATH_IMAGE017
，则输入单音信号
Figure DEST_PATH_IMAGE019
可以表示为公式1。对乘以互为正交的2个谐波信号得到
Figure 247721DEST_PATH_IMAGE020
和
Figure 508938DEST_PATH_IMAGE004
两个正交分量（见公式2和公式3），其中
Figure DEST_PATH_IMAGE021
为进行频率估计的参考频率，一般选取被测单音信号的理想频率。对
Figure 53183DEST_PATH_IMAGE002
和
Figure 859596DEST_PATH_IMAGE004
进行积化和差转换，得到公式4、5、6，其中
Figure 676243DEST_PATH_IMAGE009
为输入单音与参考频率之差。对
Figure 226304DEST_PATH_IMAGE020
和正交分量进行低通滤波处理，滤除高频分量后，其输出信号
Figure DEST_PATH_IMAGE023
和
Figure DEST_PATH_IMAGE025
如公式7、8所示。取
Figure DEST_PATH_IMAGE027
，其中
Figure DEST_PATH_IMAGE029
是采样周期，则低通滤波器输出的正交分量信号表示为
Figure DEST_PATH_IMAGE031
和
Figure DEST_PATH_IMAGE033
，即第n个采样时刻的低通滤波器输出的正交分量，表示为公式9、10。由和
Figure 872026DEST_PATH_IMAGE033
表达式可知，第n个采样时刻低通滤波器输出信号的相位
Figure DEST_PATH_IMAGE035
可以表示为公式11，由公式9和10推出
Figure 22515DEST_PATH_IMAGE036
可以表示为公式12。根据公式11、12，可以计算相差一个周期T时间的滤波器相位差
Figure 479036DEST_PATH_IMAGE038
，表示为公式13、14。在实际应用中采用多次测量求平均值（数学期望）减少误差的方法，利用公式14和公式15计算出
Figure 876519DEST_PATH_IMAGE040
。联立公式13和15，可以估计出输入频率和参考频率
Figure 559621DEST_PATH_IMAGE007
之间的频率差
Figure 873578DEST_PATH_IMAGE009
（见公式16），进而通过公式17估计输入单音信号频率。
由公式1～17可知，单音信号频率
Figure DEST_PATH_IMAGE041
的估计精度取决于
Figure 533546DEST_PATH_IMAGE009
的精度，即取决于
Figure 229101DEST_PATH_IMAGE040
，而取决于低通滤波器输出信号的相位
Figure 399499DEST_PATH_IMAGE035
的精度。由公式12可知，低通滤波器输出信号的相位
Figure 899751DEST_PATH_IMAGE036
的精度取决于反正切函数的精度，理论上精度可以达到任意希望精度。
在实际应用过程中，采用数字信号进行单音信号频率估计，与FFT相比在空间复杂度和时间复杂度上具有以下优势。
1)精度
采用FFT方法，其频率精度为
Figure DEST_PATH_IMAGE043
（为采样频率，
Figure DEST_PATH_IMAGE047
为采样点数），当
Figure 5853DEST_PATH_IMAGE048
取8000Hz时，频率分辨率要达到约1Hz，需要
Figure DEST_PATH_IMAGE049
＝8192，这将大大增加系统的时间复杂度和空间复杂度，降低实时性，而且1Hz的精度并不高。
选用公式1～17的方法，理论上可以达到任意高的精度，且对时间复杂度和空间复杂度影响不大。
2)空间复杂度
采用FFT进行频率估计，需要存储一定数量的输入信号样本。由于其估计精度为
Figure 353789DEST_PATH_IMAGE050
，存储的输入信号样本越多其估计精度越高。
而采用公式1～17所述方法，仅需要存储若干个间隔为T时间的样本，即可以实现频率估计，空间复杂度远远小于FFT方法。
3)时间复杂度
采用点的FFT方法，理论上需要个复数加运算和
Figure 554461DEST_PATH_IMAGE054
复数乘运算，为了提高频率分辨率，增大N，将同时增加时间复杂度。
采用公式1～17所述方法，如果选择N阶低通滤波器，其时间复杂度基本在100个乘加运算以内，远远小于FFT方法。
Figure 420218DEST_PATH_IMAGE056
                                   （1）
Figure 766885DEST_PATH_IMAGE058
                                  （2）
Figure 992462DEST_PATH_IMAGE060
                                  （3）
Figure 895827DEST_PATH_IMAGE062
                        （4）
Figure 139726DEST_PATH_IMAGE064
                        （5）
Figure 142448DEST_PATH_IMAGE066
                                         （6）
Figure 42271DEST_PATH_IMAGE068
                                     （7）
Figure 670699DEST_PATH_IMAGE070
                                    （8）
Figure 844804DEST_PATH_IMAGE072
                                （9）
Figure 18427DEST_PATH_IMAGE074
                                  （10）
Figure 467863DEST_PATH_IMAGE076
                                         （11）
Figure 650714DEST_PATH_IMAGE078
                                       （12）
                   （13）
                    （14）
                                      （15）
Figure DEST_PATH_IMAGE086
                                        （16）
Figure DEST_PATH_IMAGE088
                                        （17）
2.一种高精度的频率估计方法系统设计方案
为实现高精度的频率估计，采用如图1所示的系统设计方案。整个系统主要由信号采集模块，信号处理模块，信号显示和存储模块，配置模块四个部分组成。
1)信号采集模块
主要实现单音信号的采集，支持模拟信号输入采集和数字信号输入采集两种模式，将采集的数据传送给信号处理模块进行频率估计。
2)信号处理模块
信号处理模块读取信号采集模块传送的数据，根据配置模块配置的信息按公式（1～17）的方式进行输入信号的频率估计，首先将输入数字信号利用参考频率进行
Figure DEST_PATH_IMAGE089
和
Figure 29709DEST_PATH_IMAGE004
两个正交分量（见公式2和公式3）处理，并对正交分量
Figure 451594DEST_PATH_IMAGE020
和
Figure 687404DEST_PATH_IMAGE004
进行低通滤波，然后求其一周期相位差。多次测量求其相位差的均值，进而高精度的估计输入单音频率，实现高精度的频率估计，并将估计的结果输出到显示和存储模块。
3)信号显示和存储模块
信号显示和存储模块用于存储和显示信号处理模块输出的频率估计结果，根据应用领域不同，该显示和存储模块可以使用不同的形式。
4)配置模块
配置模块主要是对系统功能进行配置工作，主要包括参考频率设置和特定功能配置。采用参考频率设置时，系统根据设置的参考频率进行高精度频率估计。采用特定功能配置等功能时，系统将主要对特定功能所要求的频率进行估计，其方法是对于输入样本将该组所有频率依次用做参考频率
Figure 511134DEST_PATH_IMAGE021
，计算各自的
Figure 541407DEST_PATH_IMAGE009
，选最小
Figure 447659DEST_PATH_IMAGE011
估计输入单音频率。用户可以根据自己的需求进行特定功能配置。
3.一种高精度的频率估计系统硬件设计方案
根据高精度频率估计系统设计方案，采用图2所示的硬件设计方案。输入接口主要包括模拟输入接口和数字输入接口两种模式，其中模拟输入接口接收来自线路上或麦克风的单音信号，经AD转化为数字信号，数字接口可以直接接收数字线路上的数字信号，两种输入模式可以通过多路选择开关MUX进行功能选择，通过IO接口输送给处理器。处理器在系统开始工作时首先读取存储器中存储的配置信息，然后进入信号处理模块对输入的数字信号进行高精度频率估计，处理器可以通过显示接口输出频率估计结果，可以通过存储接口存储估计结果到存储器。用户可以通过键盘和触摸屏进行高精度频率估计系统配置，配置信息存储于系统的存储器中。
4.一种高精度频率估计系统的软件设计方案
系统启动后，软件根据用户输入决定系统进入相应的处理模块。当用户选择进入配置模式时，用户可以选择“单频率估计模式”和“特定功能配置模式”两种模式。当用户选择“单频率估计模式”时，用户输入参考频率或者从特定功能模式下存储的频率中选择参考频率（例如在钢琴键音中选择一个调音目标音，例如中央“dou”，用户可以不知道该目标频率的确切值，系统自动提取对应的频率作为参考频率），并存储参考频率到系统中。当用户选择“特定功能配置模式”时，用户输入一组参考频率和对应标识信息（如钢琴按键音与其频率）并存储于系统配置信息中，用户可以根据实际应用配置多组特定功能，每次开始频率估计时仅有一种特定功能起作用。当开始频率估计时，如果是“单频率估计模式”，将以用户输入的参考频率进行高精度频率估计。如果是“特定功能频率估计模式”，将以该特定功能所配置一组频率分别作为参考频率对输入信号进行频率估计，根据配置的阈值信息输出频率估计信息。其方法是对于输入单音将该组所有频率依次用做参考频率
Figure 221580DEST_PATH_IMAGE007
，计算各自的
Figure 899817DEST_PATH_IMAGE009
，选最小
Figure 100991DEST_PATH_IMAGE011
估计输入单音频率。用户可以根据需要终止频率估计。
上述实施例只是为了说明本发明的技术构思及特点，其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所作出的等效的变化或修饰，都应涵盖在本发明的保护范围内。

Claims (4)

1.一种高精度的频率估计系统，其特征在于，所述频率估计系统由信号采集模块，信号处理模块，信号显示和存储模块，配置模块四个部分组成,
信号采集模块：主要实现单音信号的采集，支持模拟信号输入采集和数字信号输入采集两种模式，将采集的数据传送给信号处理模块进行频率估计；
信号处理模块：读取信号采集模块传送的数据，根据配置模块配置的信息进行输入信号的频率估计，首先将输入数字信号利用参考频率进行                                               和两个正交分量处理，并对正交分量
Figure DEST_PATH_IMAGE002A
和进行低通滤波，然后求其一周期相位差；多次测量求其相位差的均值，进而高精度的估计输入单音频率，并将估计的结果输出到显示和存储模块；
信号显示和存储模块：用于存储和显示信号处理模块输出的频率估计结果；
配置模块：主要包括参考频率设置和特定功能配置的功能。
2.根据权利要求1所述的高精度的频率估计系统，其特征在于，所述采集数据的输入接口包括模拟输入接口和数字输入接口两种模式，模拟输入接口接收来自线路上或麦克风的单音信号，经AD转换后转化为数字信号；数字接口可直接接收数字线路上的数字信号，两种模式的输入通过多路选择开关MUX进行功能选择，通过IO接口输送给处理器。
3.根据权利要求1所述的高精度的频率估计系统，其特征在于，所述频率估计系统的处理步骤为，处理器在系统开始工作时首先读取存储器中存储的配置信息，然后进入信号处理模块对输入的数字信号进行频率估计，处理器可以通过显示接口输出频率估计结果，可以通过存储接口存储频率估计结果到存储器。
4.根据权利要求1所述的高精度的频率估计系统，其特征在于，所述配置模块中，采用参考频率设置时，系统根据设置的参考频率进行频率估计；采用特定功能配置时，系统将主要对特定功能所要求的频率进行估计，其方法是对于输入样本将该组所有频率依次用做参考频率
Figure 2012102949498100001DEST_PATH_IMAGE006
，计算输入单音与各个参考频率之差
Figure DEST_PATH_IMAGE008
，选输入单音与参考频率最小的
Figure DEST_PATH_IMAGE010
估计输入单音频率。
 
CN201210294949.8A 2012-08-20 2012-08-20 一种高精度的频率估计系统 Active CN103234624B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210294949.8A CN103234624B (zh) 2012-08-20 2012-08-20 一种高精度的频率估计系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210294949.8A CN103234624B (zh) 2012-08-20 2012-08-20 一种高精度的频率估计系统

Publications (2)

Publication Number Publication Date
CN103234624A true CN103234624A (zh) 2013-08-07
CN103234624B CN103234624B (zh) 2015-01-07

Family

ID=48882675

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210294949.8A Active CN103234624B (zh) 2012-08-20 2012-08-20 一种高精度的频率估计系统

Country Status (1)

Country Link
CN (1) CN103234624B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107864071A (zh) * 2017-11-02 2018-03-30 江苏物联网研究发展中心 一种面向主动安全的数据动态采集方法、装置及系统
CN112151065A (zh) * 2019-06-28 2020-12-29 力同科技股份有限公司 单音信号频率检测方法、装置、设备及计算机存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101057150A (zh) * 2004-09-13 2007-10-17 菲尔特罗尼克公开有限公司 用于测量接收的信号的频率的方法和设备
CN101374134A (zh) * 2008-07-03 2009-02-25 天津大学 一种适用于时域同步正交频分复用接收机的精确频偏估计方法
US20090131006A1 (en) * 2007-11-20 2009-05-21 Mediatek Inc. Apparatus, integrated circuit, and method of compensating iq phase mismatch
CN101963659A (zh) * 2010-10-25 2011-02-02 哈尔滨工程大学 利用三维正交干涉仪实现宽频段测向的方法
US7912157B2 (en) * 2008-04-04 2011-03-22 Newport Media, Inc. Symmetric pilot processing for robust timing offset and frequency offset estimation in ISDB-T and ISDB-TSB receivers

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101057150A (zh) * 2004-09-13 2007-10-17 菲尔特罗尼克公开有限公司 用于测量接收的信号的频率的方法和设备
US20090131006A1 (en) * 2007-11-20 2009-05-21 Mediatek Inc. Apparatus, integrated circuit, and method of compensating iq phase mismatch
US7912157B2 (en) * 2008-04-04 2011-03-22 Newport Media, Inc. Symmetric pilot processing for robust timing offset and frequency offset estimation in ISDB-T and ISDB-TSB receivers
CN101374134A (zh) * 2008-07-03 2009-02-25 天津大学 一种适用于时域同步正交频分复用接收机的精确频偏估计方法
CN101963659A (zh) * 2010-10-25 2011-02-02 哈尔滨工程大学 利用三维正交干涉仪实现宽频段测向的方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107864071A (zh) * 2017-11-02 2018-03-30 江苏物联网研究发展中心 一种面向主动安全的数据动态采集方法、装置及系统
CN112151065A (zh) * 2019-06-28 2020-12-29 力同科技股份有限公司 单音信号频率检测方法、装置、设备及计算机存储介质
CN112151065B (zh) * 2019-06-28 2024-03-15 力同科技股份有限公司 单音信号频率检测方法、装置、设备及计算机存储介质

Also Published As

Publication number Publication date
CN103234624B (zh) 2015-01-07

Similar Documents

Publication Publication Date Title
CN1904634B (zh) 用来加速全球定位系统信号的关联性处理的方法与装置
CN103592843B (zh) 一种时间戳电路和实现方法
CN109979469B (zh) 信号处理方法、设备及存储介质
CN105527516B (zh) 一种通道校准方法及电子设备
CN106506101B (zh) 一种接收机幅频响应自动化校准方法及装置
CN203287435U (zh) 一种基于stm32f107vct6的微电网谐波与间谐波检测装置
CN104545870A (zh) 用于心率检测的床垫和心率检测方法
CN110045379A (zh) 一种距离测量方法、相关设备及系统
CN103728494A (zh) 一种谐振频率测量及跟踪系统和方法
CN106291105B (zh) 一种基于数字零中频的扫频仪
CN102435842A (zh) 双模式频标比对与频率稳定度测量装置
CN107084662A (zh) 基于fpga的正交信号实时处理方法
CN100562762C (zh) 用于定位全球定位系统相关峰值信号的方法和系统
CN110247866A (zh) 一种基于doa估计的dmwc频谱感知相位校准方法
CN109005002A (zh) 传感器数据处理装置、传感器系统以及用于在所述传感器系统的范畴内确定换算参数的方法
CN103234624B (zh) 一种高精度的频率估计系统
US20230379200A1 (en) Methods and apparatus for demodulating digital signals
US20230099986A1 (en) Pattern generator and built-in-self test device including the same
CN103234623A (zh) 一种高精度的频率估计方法
CN112242844A (zh) 信号处理方法及其系统
US10866079B2 (en) Position sensing device
RU2313101C1 (ru) Анализатор плотности вероятности фазы сигнала
CN108414001A (zh) 非均匀采样正弦波形失真度的确定方法
CN108549618B (zh) 动模量、阻尼比计算方法及装置
CN116988782B (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
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210715

Address after: 199 Ren'ai Road, industrial park, Qingpu District, Shanghai, 201799

Patentee after: Sundec semiconductor technology (Shanghai) Co.,Ltd.

Address before: 215000 199 Ren Yan Road, Suzhou Industrial Park, Jiangsu

Patentee before: SOOCHOW University

TR01 Transfer of patent right
CP01 Change in the name or title of a patent holder

Address after: 199 Ren'ai Road, industrial park, Qingpu District, Shanghai, 201799

Patentee after: Sundak Semiconductor Technology (Shanghai) Co.,Ltd.

Address before: 199 Ren'ai Road, industrial park, Qingpu District, Shanghai, 201799

Patentee before: Sundec semiconductor technology (Shanghai) Co.,Ltd.

CP01 Change in the name or title of a patent holder