[go: up one dir, main page]

CN114838809B - 一种自适应提高频率测量精度的音频信号测量方法 - Google Patents

一种自适应提高频率测量精度的音频信号测量方法 Download PDF

Info

Publication number
CN114838809B
CN114838809B CN202210282390.0A CN202210282390A CN114838809B CN 114838809 B CN114838809 B CN 114838809B CN 202210282390 A CN202210282390 A CN 202210282390A CN 114838809 B CN114838809 B CN 114838809B
Authority
CN
China
Prior art keywords
data
frequency
value
signal
extraction number
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
Application number
CN202210282390.0A
Other languages
English (en)
Other versions
CN114838809A (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.)
Intelligent Automation Equipment Zhuhai Co Ltd
Original Assignee
Intelligent Automation Equipment Zhuhai 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 Intelligent Automation Equipment Zhuhai Co Ltd filed Critical Intelligent Automation Equipment Zhuhai Co Ltd
Priority to CN202210282390.0A priority Critical patent/CN114838809B/zh
Publication of CN114838809A publication Critical patent/CN114838809A/zh
Application granted granted Critical
Publication of CN114838809B publication Critical patent/CN114838809B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明为一种自适应提高频率测量精度的音频信号测量方法,旨在提供一种在采样率和实际FFT处理计算数据量不变的情况下自适应提高FFT处理的计算分辨率,从而提高了信号频率检测精度的音频信号测量方法。本发明包括以下流程:首先,设定模数转换器的采样率、设定一偏差比例k以及设定数据抽取数m为1,初次采集时对采集的n个数据执行加窗处理后做FFT运算处理获取相应的复数形式的数据,根据获取的参数和数据执行校准算法计算计算输入信号的主频率f1,进一步计算输入信号的主频率f1与当前的频率分辨率fp的倍数关系x,根据倍数关系x与设定的偏差比例k比对结果对数据抽取数m进行自适应调整,调整后重新采样获取精度更高的频率数值f2,完成比对后输出测量结果。本发明应用于音频信号测量的技术领域。

Description

一种自适应提高频率测量精度的音频信号测量方法
技术领域
本发明应用于音频信号测量的技术领域,特别涉及一种自适应提高频率测量精度的音频信号测量方法。
背景技术
生活中我们常见的手机,电脑等各种带有音频输入的电子设备以及各种电器都有音频采集的功能。这些音频采集模块在出厂之前是需要验证其功能与性能的。评估音频采集功能的常见参数有:输入信号的频率,输入信号的幅值,底噪,谐波分量等等。对于人耳能够分辨的声音频率范围是20~20000Hz,因此我们的音频模块需要测试的信号频率也在此范围。根据奈奎斯特采样定理,只要采样频率大于或等于有效信号最高频率的两倍,采样值就可以包含原始信号的所有信息,被采样的信号就可以不失真地还原成原始信号。在测试过程中,我们常用模数转换器的采样频率有48KHz,96KHz,192KHz。为了验证音频模块采集的信号频率是否准确,我们需要将输入信号从时域转换到频域,这样能够更方便的计算信号的各种指标。FFT处理(快速傅氏变换)就是将时域信号转换成频域信号,其输入参数有:被测试信号,采样频率。
FFT处理的分辨率计算公式为:fp=Fs/N,公式中fp为频率分辨率,Fs为采样率,N为采集信号的数据个数。如果要提高测量信号频率的精度,需要降低信号采样率,或者增加采集信号的数据个数。
对于前端输入信号频率的不确定性,并不能确定采样率设置多少更合适,而且当硬件确定后,模数转换器的芯片也确定了,那么该模数转换器能够选择使用的采样率就固定了。如果只是单纯的增加采集信号的数据量,需要消耗更多的资源和成本。
发明内容
本发明所要解决的技术问题是克服现有技术的不足,提供了一种在采样率和实际FFT处理计算数据量不变的情况下自适应提高FFT处理的计算分辨率,从而提高了信号频率检测精度的音频信号测量方法。
本发明所采用的技术方案是:本发明所述音频信号测量方法包括以下步骤:
步骤S1、为进行采样的模数转换器设定一固定的采样率Fs以及设定偏差比例为k;
步骤S2、初次采样时设置数据抽取数m=1,对模数转换器采样回来的数据量不做任何处理;
步骤S3、执行初次信号采集,读取模数转换器采样回来的N个数据并储存在存储器中;
步骤S4、完成存储后对N个数据执行加窗处理后做FFT运算处理,产生N个的复数形式的数据,并根据模数转换器的采样率Fs、N个复数形式数据、数据抽取数m以及窗函数的校准算法计算输入信号的主频率f1;
步骤S5、设输入信号的主频率f1与当前的频率分辨率fp的倍数关系为x,且x=f1/fp,fp=Fs/N,根据输入信号的主频率f1计算当前倍数关系x的数值,将倍数关系x的值与设定的偏差比例k的值进行比对;
步骤S6、当x<k时,根据初次采样获得输入信号的主频率f1、初次采样数据个数N、设定的偏差比例k以及当前的频率分辨率fp进行数据抽取数m的调整,数据抽取数m的值根据中间变量n进行取值,n=(fp*k)/(N*f1),同时依据调整后的数据抽取数m执行二次信号采集,并获得新的N个数据,对新的N个数据执行加窗处理后做FFT运算处理,产生新的N个复数形式的数据;
步骤S7、根据模数转换器的采样率Fs、新的N个复数形式数据、数据抽取数m以及窗函数的校准算法计算输入信号的主频率f2;
步骤S8、根据调整后输入信号的主频率f2计算当前倍数关系x的数值,将倍数关系x的值与设定的偏差比例k的值进行比对,当x>=k的值时表示频率分辨率达到设定范围并输出测量结果,反之则返回步骤S6。
由上述方案可见,通过将第一次采集的音频信号进行分析,根据分析结果判断是否需要做第二次处理分析。评估抽取数据值时,采用输入信号的主频率f1和当前的频率分辨率fp的倍数关系x与设定的偏差比例k进行对比,其中k为用户设定值,当k值越大, 频率分辨率fp越高。可以自动评估当前检测频率是否满足用户所要求的分辨率,最终可以提高频率计算的精度。而且不需要用户去根据输入信号的频率去调整模数转换器的采样率和采集数据的个数。其中,步骤S4和步骤S7中的窗函数的选择根据用户的需求设定。
一个优选方案是,步骤S4和步骤S7中窗函数为nuttall 四项三阶窗函数。
由上述方案可见,窗函数的选择根据用户的需求设定,通过选择nuttall 四项三阶窗函数目的在于使thd+n的旁瓣衰减能够<-80dB,进而保证测量精度。
进一步的优选方案是,所述nuttall四项三阶窗函数的校准算法计算频率包括以下步骤:
步骤C1、数据经过FFT运算处理后为N个复数形式数据 x+yj,计算获得N个复数数据的模p,
步骤C2、找到模值最大的数据和模值第二大的数据以及数据对应在N个数据中的顺序位置
步骤C3、根据上述参数计算校准系数
步骤C4、根据上述参数计算校准系数,其中均为Nuttall四项三阶的计算系数的常数;
步骤C5、计算频率f,
由上述方案可见,FFT处理不同的窗函数频率校准算法有所不同,可根据窗函数选择适合自己的校准算法进行频率f的计算,本方案中使用的校准算法是双峰谱线校准法。
一个优选方案是,步骤S6中数据抽取数m的取值需满足m>=n且为当前的频率分辨率fp的约数。
由上述方案可见,在设置数据抽取数m时,通过限定抽取数必须是采样率的公约数,进而实现保证数据抽取后的采样率是整数。
一个优选方案是,步骤S6中依据调整后的数据抽取数m执行二次信号采集的方法具体为:每次读取信号时读取m个数据,取当次读取m个数据中最后读出的数据保存至存储器中,完成N次数据保存后获得新的N个信号数据。
由上述方案可见,通过调整后的数据抽取数m使得实际读取的数据量为N*m,取每次读取信号时最后获得的数据,通过取多次读取后的稳定数值,使得获得更精准的采样数据。
附图说明
图1是本发明的流程图。
具体实施方式
如图1所示,在本实施例中,本发明包括以下步骤:
步骤S1、为进行采样的模数转换器设定一固定的采样率Fs以及设定偏差比例为k;
步骤S2、初次采样时设置数据抽取数m=1,对模数转换器采样回来的数据量不做任何处理;
步骤S3、执行初次信号采集,读取模数转换器采样回来的N个数据并储存在存储器中;
步骤S4、完成存储后对N个数据执行加窗处理后做FFT运算处理,产生N个的复数形式的数据,并根据模数转换器的采样率Fs、N个复数形式数据、数据抽取数m以及窗函数的校准算法计算输入信号的主频率f1;
步骤S5、设输入信号的主频率f1与当前的频率分辨率fp的倍数关系为x,且x=f1/fp,fp=Fs/N,根据输入信号的主频率f1计算当前倍数关系x的数值,将倍数关系x的值与设定的偏差比例k的值进行比对;
步骤S6、当x<k时,根据初次采样获得输入信号的主频率f1、初次采样数据个数N、设定的偏差比例k以及当前的频率分辨率fp进行数据抽取数m的调整,数据抽取数m的值根据中间变量n进行取值,n=(fp*k)/(N*f1),同时依据调整后的数据抽取数m执行二次信号采集,并获得新的N个数据,对新的N个数据执行加窗处理后做FFT运算处理,产生新的N个复数形式的数据。其中,依据调整后的数据抽取数m执行二次信号采集的方法具体为:每次读取信号时读取m个数据,取当次读取m个数据中最后读出的数据保存至存储器中,完成N次数据保存后获得新的N个信号数据;
步骤S7、根据模数转换器的采样率Fs、新的N个复数形式数据、数据抽取数m以及窗函数的校准算法计算输入信号的主频率f2;
步骤S8、根据调整后输入信号的主频率f2计算当前倍数关系x的数值,将倍数关系x的值与设定的偏差比例k的值进行比对,当x>=k的值时表示频率分辨率达到设定范围并输出测量结果,反之则返回步骤S6。
在本实施例中,步骤S4和步骤S7中窗函数选择nuttall 四项三阶窗函数,进而实现thd+n的旁瓣衰减能够<-80dB。
另外,步骤S4和步骤S7根据待测量的音频信号参数可自由选取相应性能的窗函数,示例窗函数如表1所示:
在本实施例中,所述nuttall四项三阶窗函数的校准算法计算频率包括以下步骤:
步骤C1、数据经过FFT运算处理后为N个复数形式数据 x+yj,计算获得N个复数数据的模p,;其中,复数形式数据x+yj中的x代表复数形式数据的实部,y代表复数形式数据的虚部;
步骤C2、找到模值最大的数据和模值第二大的数据以及数据对应在N个数据中的顺序位置
步骤C3、根据上述参数计算校准系数
步骤C4、根据上述参数计算校准系数,其中均为Nuttall四项三阶的计算系数的常数,
步骤C5、计算频率f,
在本实施例中,步骤S6中数据抽取数m的取值需满足m>=n且为当前的频率分辨率fp的约数。
实施例一:
输入的音频信号频率为20Hz,数据的采样数量N为8192;
常规测量方法对比:
实施例二:
输入的音频信号频率为30Hz,数据的采样数量N为8192;
常规测量方法对比:
实施例三:
输入的音频信号频率为40Hz,数据的采样数量N为8192;
常规测量方法对比:
结论:通过申请设定模数转换器的采样率、设定一偏差比例k以及设定数据抽取数m为1,在初次采集时对采集的n个数据执行加窗处理后做FFT运算处理获取相应的复数形式的数据,根据获取的参数和数据执行校准算法计算计算输入信号的主频率f1,进一步计算输入信号的主频率f1与当前的频率分辨率fp的倍数关系x,根据倍数关系x与设定的偏差比例k比对结果对数据抽取数m进行自适应调整,进而通过调整后的数据抽取数进行二次采样以及采样数据的筛选,再根据自动调整参数后的函数进行计算获取测量值,实现自适应提高测量分辨率。达到在采样率和实际FFT处理计算数据量不变的情况下自适应提高FFT处理的计算分辨率,保证消耗资源和成本相同的情况下提高信号频率检测精度。

Claims (4)

1.一种自适应提高频率测量精度的音频信号测量方法,其特征在于,它包括以下步骤:步骤S1、为进行采样的模数转换器设定一固定的采样率Fs以及设定偏差比例为k;
步骤S2、初次采样时设置数据抽取数m=1,对模数转换器采样回来的数据量不做任何处理;
步骤S3、执行初次信号采集,读取模数转换器采样回来的N个数据并储存在存储器中;
步骤S4、完成存储后对N个数据执行加窗处理后做FFT运算处理,产生N个的复数形式的数据,并根据模数转换器的采样率Fs、N个复数形式数据、数据抽取数m以及窗函数的校准算法计算输入信号的主频率f1;
步骤S5、设输入信号的主频率f1与当前的频率分辨率fp的倍数关系为x,且x=f1/fp,fp=Fs/N,根据输入信号的主频率f1计算当前倍数关系x的数值,将倍数关系x的值与设定的偏差比例k的值进行比对;
步骤S6、当x<k时,根据初次采样获得输入信号的主频率f1、初次采样数据个数N、设定的偏差比例k以及当前的频率分辨率fp进行数据抽取数m的调整,数据抽取数m的值根据中间变量n进行取值,n=(fp*k)/(N*f1),同时依据调整后的数据抽取数m执行二次信号采集,并获得新的N个数据,对新的N个数据执行加窗处理后做FFT运算处理,产生新的N个复数形式的数据;其中,数据抽取数m的取值需满足m>=n且为当前的频率分辨率fp的约数;
步骤S7、根据模数转换器的采样率Fs、新的N个复数形式数据、数据抽取数m以及窗函数的校准算法计算输入信号的主频率f2;
步骤S8、根据调整后输入信号的主频率f2计算当前倍数关系x的数值,将倍数关系x的值与设定的偏差比例k的值进行比对,当x>=k的值时表示频率分辨率达到设定范围并输出测量结果,反之则返回步骤S6。
2.根据权利要求1所述的一种自适应提高频率测量精度的音频信号测量方法,其特征在于,步骤S4和步骤S7中窗函数为nuttall 四项三阶窗函数。
3.根据权利要求2所述的一种自适应提高频率测量精度的音频信号测量方法,其特征在于,所述nuttall四项三阶窗函数的校准算法计算频率包括以下步骤:
步骤C1.数据经过FFT运算处理后为N个复数形式数据 x+yj,计算获得N个复数数据的模p,
步骤C2.找到模值最大的数据和模值第二大的数据以及数据对应在N个数据中的顺序位置
步骤C3.根据上述参数计算校准系数
步骤C4.根据上述参数计算校准系数,其中均为Nuttall四项三阶的计算系数的常数;
步骤C5.计算频率f,
4.根据权利要求1所述的一种自适应提高频率测量精度的音频信号测量方法,其特征在于,步骤S6中依据调整后的数据抽取数m执行二次信号采集的方法具体为:每次读取信号时读取m个数据,取当次读取m个数据中最后读出的数据保存至存储器中,完成N次数据保存后获得新的N个信号数据。
CN202210282390.0A 2022-03-22 2022-03-22 一种自适应提高频率测量精度的音频信号测量方法 Active CN114838809B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210282390.0A CN114838809B (zh) 2022-03-22 2022-03-22 一种自适应提高频率测量精度的音频信号测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210282390.0A CN114838809B (zh) 2022-03-22 2022-03-22 一种自适应提高频率测量精度的音频信号测量方法

Publications (2)

Publication Number Publication Date
CN114838809A CN114838809A (zh) 2022-08-02
CN114838809B true CN114838809B (zh) 2024-10-18

Family

ID=82562509

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210282390.0A Active CN114838809B (zh) 2022-03-22 2022-03-22 一种自适应提高频率测量精度的音频信号测量方法

Country Status (1)

Country Link
CN (1) CN114838809B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110377873A (zh) * 2019-08-12 2019-10-25 山东超越数控电子股份有限公司 一种基于fft的信号谱线分析方法与装置
CN111308199A (zh) * 2020-03-12 2020-06-19 国网湖南省电力有限公司 基于频谱分辨率自适应的双谱线插值dft谐波分析方法、系统及介质

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1062963C (zh) * 1990-04-12 2001-03-07 多尔拜实验特许公司 用于产生高质量声音信号的解码器和编码器
CA2126810C (en) * 1993-06-29 2004-10-19 Gen Ichimura Audio signal transmitting apparatus and the method thereof
US6298363B1 (en) * 1997-11-13 2001-10-02 The Johns Hopkins University Adaptive windowing of FFT data for increased resolution and sidelobe rejection
JP5300188B2 (ja) * 2006-09-11 2013-09-25 株式会社東芝 超音波診断装置及び超音波診断装置の制御プログラム
CN103207319A (zh) * 2013-03-12 2013-07-17 广东电网公司电力科学研究院 数字化变电站电力信号非同步采样条件下的谐波测量方法
CN104391282A (zh) * 2014-11-24 2015-03-04 天津大学 通过频谱校正来改善成像质量的方法
CN105513608B (zh) * 2015-07-23 2018-12-28 中国电子科技集团公司第四十一研究所 一种音频信号分析方法
CN105307095B (zh) * 2015-09-15 2019-09-10 中国电子科技集团公司第四十一研究所 一种基于fft的高分辨率音频频率测量方法
CN107064628B (zh) * 2017-04-13 2019-08-16 中国电子科技集团公司第二十四研究所 高精度频率测量系统及方法
US11290903B2 (en) * 2019-07-17 2022-03-29 SiTune Corporation Spectrum monitoring
CN112541157B (zh) * 2020-11-30 2024-03-22 西安精密机械研究所 一种信号频率精确估计方法
CN113611319B (zh) * 2021-04-07 2023-09-12 珠海市杰理科技股份有限公司 基于语音成分实现的风噪抑制方法、装置、设备及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110377873A (zh) * 2019-08-12 2019-10-25 山东超越数控电子股份有限公司 一种基于fft的信号谱线分析方法与装置
CN111308199A (zh) * 2020-03-12 2020-06-19 国网湖南省电力有限公司 基于频谱分辨率自适应的双谱线插值dft谐波分析方法、系统及介质

Also Published As

Publication number Publication date
CN114838809A (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
US9363599B2 (en) Systems and methods for protecting a speaker
RU2596592C2 (ru) Пространственный аудио процессор и способ обеспечения пространственных параметров на основе акустического входного сигнала
US9357300B2 (en) Systems and methods for protecting a speaker
CN101674520B (zh) 振动式拾音器参数分析测试方法和装置
CN109831733A (zh) 音频播放性能的测试方法、装置、设备和存储介质
US7203324B2 (en) Audio feedback processing system
CN117171516B (zh) 一种x射线测厚仪数据优化校正方法
US20200082835A1 (en) Methods and apparatus to fingerprint an audio signal via normalization
CN102176667A (zh) 一种基于频率跟踪的工频数字陷波装置与方法
CN114199365A (zh) 一种振动信号处理方法
CN104869519A (zh) 一种测试麦克风本底噪声的方法和系统
JP5627440B2 (ja) 音響装置及びその制御方法、プログラム
CN114838809B (zh) 一种自适应提高频率测量精度的音频信号测量方法
CN106092492B (zh) 一种滤波去噪方法
CN111641422A (zh) 高动态数字接收机自适应抗干扰检测门限确定系统及方法
US11528558B2 (en) Method and an audio processing unit for detecting a tone
WO2018201374A1 (en) Method and device for adjusting audio signal, and audio system
Hall et al. Use of adaptive filtering for improved performance in digital stethoscopes
JPH1198016A (ja) Ad変換器の評価装置
TWI475870B (zh) 電話裝置的測試裝置以及測試方法
US12135233B2 (en) Acoustic noise detection method and system using vibration sensor to detect acoustic noise
JP2012095254A (ja) 音量調整装置、音量調整方法及び音量調整プログラム並びに音響機器
JPH11103252A (ja) Da変換器の評価装置
JP2024057456A (ja) 震度推定装置、震度推定プログラム及び震度推定方法
CN118817064A (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