CN114838809A - Audio signal measuring method for self-adaptively improving frequency measurement precision - Google Patents

Audio signal measuring method for self-adaptively improving frequency measurement precision Download PDF

Info

Publication number
CN114838809A
CN114838809A CN202210282390.0A CN202210282390A CN114838809A CN 114838809 A CN114838809 A CN 114838809A CN 202210282390 A CN202210282390 A CN 202210282390A CN 114838809 A CN114838809 A CN 114838809A
Authority
CN
China
Prior art keywords
data
frequency
calculating
value
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.)
Pending
Application number
CN202210282390.0A
Other languages
Chinese (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/CN114838809A/en
Publication of CN114838809A publication Critical patent/CN114838809A/en
Pending legal-status Critical Current

Links

Images

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

The invention provides an audio signal measuring method for improving frequency measuring precision in a self-adaptive mode, and aims to provide an audio signal measuring method for improving the calculation resolution of FFT processing in a self-adaptive mode under the condition that the sampling rate and the actual FFT processing calculation data quantity are not changed, so that the signal frequency detecting precision is improved. The invention comprises the following procedures: firstly, setting the sampling rate of an analog-digital converter, setting a deviation proportion k and setting the data extraction number m to be 1, performing windowing processing on n acquired data during initial acquisition, performing FFT (fast Fourier transform) operation processing to acquire corresponding data in a complex form, performing a calibration algorithm according to the acquired parameters and the data to calculate the main frequency f1 of an input signal, further calculating the multiple relation x between the main frequency f1 of the input signal and the current frequency resolution fp, performing self-adaptive adjustment on the data extraction number m according to the comparison result of the multiple relation x and the set deviation proportion k, re-sampling after adjustment to acquire a frequency value f2 with higher precision, and outputting a measurement result after comparison. The invention is applied to the technical field of audio signal measurement.

Description

Audio signal measuring method for self-adaptively improving frequency measurement precision
Technical Field
The invention is applied to the technical field of audio signal measurement, and particularly relates to an audio signal measurement method for adaptively improving frequency measurement precision.
Background
In life, various electronic devices with audio input, such as mobile phones, computers and the like, and various electrical appliances have the function of audio acquisition. The audio acquisition modules need to be verified in terms of their functions and performance before being shipped out of the factory. Common parameters for evaluating the audio acquisition function are: the frequency of the input signal, the amplitude of the input signal, noise floor, harmonic components, etc. The range of sound frequency which can be distinguished by human ears is 20-20000 Hz, so that the frequency of signals needing to be tested by the audio module is also in the range. According to the nyquist sampling theorem, as long as the sampling frequency is greater than or equal to twice the highest frequency of the effective signal, the sampling value can contain all the information of the original signal, and the sampled signal can be restored to the original signal without distortion. During the test, the sampling frequencies of the analog-to-digital converters are 48KHz, 96KHz and 192 KHz. In order to verify whether the frequency of the signal acquired by the audio module is accurate, the input signal needs to be converted from the time domain to the frequency domain, so that various indexes of the signal can be calculated more conveniently. The FFT (fast fourier transform) converts the time domain signal into the frequency domain signal, and the input parameters include: signal under test, sampling frequency.
The resolution calculation formula of the FFT process is: fp = Fs/N, fp is frequency resolution, Fs is sampling rate, and N is the data number of the acquired signals. If the accuracy of measuring the signal frequency is to be improved, the signal sampling rate needs to be reduced, or the number of data of the acquired signal needs to be increased.
For the uncertainty of the frequency of the front-end input signal, it cannot be determined how much the sampling rate is set to be more appropriate, and when the hardware determines, the chip of the analog-to-digital converter also determines, so that the sampling rate that the analog-to-digital converter can select to use is fixed. If the data volume of the acquired signals is simply increased, more resources and cost are consumed.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides the audio signal measuring method which can adaptively improve the calculation resolution of the FFT processing under the condition that the sampling rate and the actual FFT processing calculation data volume are not changed, thereby improving the signal frequency detection precision.
The technical scheme adopted by the invention is as follows: the audio signal measuring method of the present invention comprises the steps of:
step S1, setting a fixed sampling rate Fs and a deviation ratio k for the analog-to-digital converter for sampling;
step S2, setting a data extraction number m =1 during initial sampling, and not processing the data volume sampled by the analog-to-digital converter;
step S3, performing primary signal acquisition, reading N data sampled by the analog-to-digital converter and storing the data in a memory;
step S4, after finishing the storage, executing windowing processing on the N data, then performing FFT operation processing to generate N data in complex form, and calculating the main frequency f1 of the input signal according to the sampling rate Fs of the analog-to-digital converter, the N data in complex form, the data extraction number m and the calibration algorithm of the window function;
step S5, setting the multiple relation between the main frequency f1 of the input signal and the current frequency resolution fp as x, and x = f1/fp, fp = Fs/N, calculating the numerical value of the current multiple relation x according to the main frequency f1 of the input signal, and comparing the value of the multiple relation x with the value of the set deviation proportion k;
step S6, when x is less than k, adjusting a data extraction number m according to a primary frequency f1 of an input signal obtained by primary sampling, the number N of primary sampling data, a set deviation proportion k and a current frequency resolution fp, taking a value of the data extraction number m according to an intermediate variable N, wherein N = (fp ×/(N × f 1), simultaneously executing secondary signal acquisition according to the adjusted data extraction number m, obtaining new N data, executing windowing processing on the new N data, and then performing FFT operation processing to generate new N complex form data;
step S7, calculating the main frequency f2 of the input signal according to the sampling rate Fs of the analog-digital converter, the new N complex form data, the data extraction number m and the calibration algorithm of the window function;
and step S8, calculating the numerical value of the current multiple relation x according to the adjusted main frequency f2 of the input signal, comparing the value of the multiple relation x with the set value of the deviation proportion k, indicating that the frequency resolution reaches the set range when x > = k, and outputting the measurement result, otherwise, returning to the step S6.
According to the scheme, the audio signal acquired for the first time is analyzed, and whether the second processing analysis is needed or not is judged according to the analysis result. When evaluating the extracted data value, comparing the multiple relation x of the main frequency f1 of the input signal and the current frequency resolution fp with a set deviation proportion k, wherein k is a user set value, and when the k value is larger, the frequency resolution fp is higher. Whether the current detection frequency meets the resolution required by a user can be automatically evaluated, and finally, the frequency calculation precision can be improved. And the user is not required to adjust the sampling rate of the analog-to-digital converter and the number of collected data according to the frequency of the input signal. Wherein, the selection of the window function in the steps S4 and S7 is set according to the user' S requirement.
Preferably, the window function in step S4 and step S7 is a nuttall fourth-order third-order window function.
According to the scheme, the window function is selected according to the requirement of a user, the sidelobe attenuation of thd + n can be < -80dB by selecting nuttall four-term third-order window functions, and the measurement accuracy is further ensured.
Further preferably, the calculating frequency of the calibration algorithm of the nuttall quadrinomial third-order window function includes the following steps:
step C1, the data is processed by FFT operation to be N complex form data x + yj, the modulo p of the N complex data is obtained by calculation,
Figure DEST_PATH_IMAGE001
step C2, finding the data with the maximum module value
Figure DEST_PATH_IMAGE002
Data second largest in sum
Figure DEST_PATH_IMAGE003
And the data corresponds to the sequential position in the N data
Figure DEST_PATH_IMAGE004
And
Figure DEST_PATH_IMAGE005
step C3, calculating calibration coefficient according to the parameters
Figure DEST_PATH_IMAGE006
Figure DEST_PATH_IMAGE007
Step C4, calculating calibration coefficient according to the parameters
Figure DEST_PATH_IMAGE009
Figure DEST_PATH_IMAGE010
In which
Figure DEST_PATH_IMAGE011
Figure DEST_PATH_IMAGE012
Figure DEST_PATH_IMAGE013
All four-term third-order calculation of Nuttall
Figure 32238DEST_PATH_IMAGE009
A constant of the coefficient;
step C5, calculating the frequency f,
Figure DEST_PATH_IMAGE014
according to the scheme, the frequency calibration algorithms of the window functions with different FFT processing are different, the calibration algorithm suitable for the frequency f can be selected according to the window functions to calculate the frequency f, and the calibration algorithm used in the scheme is a bimodal spectral line calibration method.
Preferably, in step S6, the value of the data extraction number m should satisfy m > = n and be a divisor of the current frequency resolution fp.
According to the scheme, when the data extraction number m is set, the sampling rate after data extraction is guaranteed to be an integer by limiting that the extraction number must be a common divisor of the sampling rate.
Preferably, the method for performing secondary signal acquisition according to the adjusted data extraction number m in step S6 specifically includes: reading m data when reading signals each time, storing the last read data in the m data read each time into a memory, and obtaining new N signal data after finishing storing the N data.
According to the scheme, the data volume actually read is N x m through the adjusted data extraction number m, the data finally obtained in each signal reading process is obtained, and more accurate sampling data are obtained through stable numerical values obtained after multiple times of reading.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
As shown in fig. 1, in the present embodiment, the present invention includes the following steps:
step S1, setting a fixed sampling rate Fs and a deviation ratio k for the analog-to-digital converter for sampling;
step S2, setting a data extraction number m =1 during initial sampling, and not processing the data volume sampled by the analog-to-digital converter;
step S3, performing primary signal acquisition, reading N data sampled by the analog-to-digital converter and storing the data in a memory;
step S4, after finishing the storage, executing windowing processing on the N data, then performing FFT operation processing to generate N data in complex form, and calculating the main frequency f1 of the input signal according to the sampling rate Fs of the analog-to-digital converter, the N data in complex form, the data extraction number m and the calibration algorithm of the window function;
step S5, setting the multiple relation between the main frequency f1 of the input signal and the current frequency resolution fp as x, and x = f1/fp, fp = Fs/N, calculating the numerical value of the current multiple relation x according to the main frequency f1 of the input signal, and comparing the value of the multiple relation x with the value of the set deviation proportion k;
and step S6, when x is less than k, adjusting the data extraction number m according to the primary frequency f1 of the input signal obtained by primary sampling, the number N of the primary sampling data, the set deviation proportion k and the current frequency resolution fp, taking the value of the data extraction number m according to the intermediate variable N, wherein N = (fp ×/(N × f 1), simultaneously executing secondary signal acquisition according to the adjusted data extraction number m, obtaining new N data, executing windowing processing on the new N data, and then performing FFT operation processing to generate new N complex form data. The method for performing secondary signal acquisition according to the adjusted data extraction number m specifically comprises the following steps: reading m data when reading signals each time, storing the last read data in the m data read each time into a memory, and obtaining new N signal data after finishing storing the N data;
step S7, calculating the main frequency f2 of the input signal according to the sampling rate Fs of the analog-digital converter, the new N complex form data, the data extraction number m and the calibration algorithm of the window function;
and step S8, calculating the numerical value of the current multiple relation x according to the adjusted main frequency f2 of the input signal, comparing the value of the multiple relation x with the set value of the deviation proportion k, indicating that the frequency resolution reaches the set range when x > = k, and outputting the measurement result, otherwise, returning to the step S6.
In this embodiment, the window function in step S4 and step S7 selects nuttall quadrinomial third-order window function, and then the sidelobe attenuation of thd + n can be < -80 dB.
In addition, the window functions with corresponding performance can be freely selected according to the audio signal parameters to be measured in the steps S4 and S7, and an exemplary window function is shown in table 1:
Figure DEST_PATH_IMAGE015
in this embodiment, the calculating frequency of the calibration algorithm of the nuttall fourth-order window function includes the following steps:
step C1, the data is processed by FFT operation to be N complex form data x + yj, and the data is obtained by calculationThe modulo p of the N complex data,
Figure 277275DEST_PATH_IMAGE001
(ii) a Wherein x in the complex form data x + yj represents a real part of the complex form data, and y represents an imaginary part of the complex form data;
step C2, finding the data with the maximum module value
Figure 525853DEST_PATH_IMAGE002
Data second largest in sum
Figure 921063DEST_PATH_IMAGE003
And the data corresponds to the sequential position in the N data
Figure 735435DEST_PATH_IMAGE004
And
Figure 387258DEST_PATH_IMAGE005
step C3, calculating calibration coefficient according to the parameters
Figure 654292DEST_PATH_IMAGE006
Figure 322033DEST_PATH_IMAGE007
Step C4, calculating calibration coefficient according to the parameters
Figure 990912DEST_PATH_IMAGE009
Figure 984276DEST_PATH_IMAGE010
Wherein
Figure 66501DEST_PATH_IMAGE011
Figure 803513DEST_PATH_IMAGE012
Figure 326899DEST_PATH_IMAGE013
All four-term third-order calculation of Nuttall
Figure 350218DEST_PATH_IMAGE009
The constant of the coefficient(s) is (are),
Figure DEST_PATH_IMAGE016
Figure DEST_PATH_IMAGE017
Figure DEST_PATH_IMAGE018
step C5, calculating the frequency f,
Figure 795106DEST_PATH_IMAGE014
in this embodiment, the value of the data extraction number m in step S6 needs to satisfy m > = n and is a divisor of the current frequency resolution fp.
The first embodiment is as follows:
the frequency of an input audio signal is 20Hz, and the sampling number N of data is 8192;
Figure DEST_PATH_IMAGE019
comparison of conventional measurement methods:
Figure DEST_PATH_IMAGE020
example two:
the frequency of an input audio signal is 30Hz, and the sampling number N of data is 8192;
Figure DEST_PATH_IMAGE021
comparison of conventional measurement methods:
Figure DEST_PATH_IMAGE022
example three:
the frequency of an input audio signal is 40Hz, and the sampling number N of data is 8192;
Figure DEST_PATH_IMAGE023
comparison of conventional measurement methods:
Figure DEST_PATH_IMAGE024
and (4) conclusion: the sampling rate of an analog-digital converter is set through application, a deviation proportion k is set, the number m of data extraction is set to be 1, windowing processing is carried out on n acquired data during initial acquisition, then FFT operation processing is carried out to acquire corresponding data in a complex form, a calibration algorithm is executed according to the acquired parameters and the acquired data to calculate the main frequency f1 of an input signal, the multiple relation x of the main frequency f1 of the input signal and the current frequency resolution fp is further calculated, the number m of data extraction is adaptively adjusted according to the comparison result of the multiple relation x and the set deviation proportion k, secondary sampling and sampling data screening are carried out through the adjusted number of data extraction, a measured value is calculated according to a function after parameters are automatically adjusted to acquire a measured value, and the adaptive improvement of the measurement resolution is achieved. The method achieves the self-adaptive improvement of the calculation resolution of the FFT processing under the condition that the sampling rate and the actual FFT processing calculation data volume are not changed, and improves the signal frequency detection precision under the condition that the consumed resources and the cost are the same.

Claims (5)

1. An audio signal measuring method for adaptively improving frequency measurement accuracy is characterized by comprising the following steps:
step S1, setting a fixed sampling rate Fs and a deviation ratio k for the analog-to-digital converter for sampling;
step S2, setting a data extraction number m =1 during initial sampling, and not processing the data volume sampled by the analog-to-digital converter;
step S3, performing primary signal acquisition, reading N data sampled by the analog-to-digital converter and storing the data in a memory;
step S4, after finishing the storage, executing windowing processing on the N data, then performing FFT operation processing to generate N data in complex form, and calculating the main frequency f1 of the input signal according to the sampling rate Fs of the analog-to-digital converter, the N data in complex form, the data extraction number m and the calibration algorithm of the window function;
step S5, setting the multiple relation between the main frequency f1 of the input signal and the current frequency resolution fp as x, and x = f1/fp, fp = Fs/N, calculating the numerical value of the current multiple relation x according to the main frequency f1 of the input signal, and comparing the value of the multiple relation x with the value of the set deviation proportion k;
step S6, when x is less than k, adjusting a data extraction number m according to a primary frequency f1 of an input signal obtained by primary sampling, the number N of primary sampling data, a set deviation proportion k and a current frequency resolution fp, taking a value of the data extraction number m according to an intermediate variable N, wherein N = (fp ×/(N × f 1), simultaneously executing secondary signal acquisition according to the adjusted data extraction number m, obtaining new N data, executing windowing processing on the new N data, and then performing FFT operation processing to generate new N complex form data;
step S7, calculating the main frequency f2 of the input signal according to the sampling rate Fs of the analog-digital converter, the new N complex form data, the data extraction number m and the calibration algorithm of the window function;
and step S8, calculating the numerical value of the current multiple relation x according to the adjusted main frequency f2 of the input signal, comparing the value of the multiple relation x with the set value of the deviation proportion k, indicating that the frequency resolution reaches the set range when x > = k, and outputting the measurement result, otherwise, returning to the step S6.
2. The method as claimed in claim 1, wherein the window function in steps S4 and S7 is nuttall quaternary window function.
3. The audio signal measuring method for adaptively improving the frequency measurement accuracy according to claim 2, wherein the calculating the frequency by the calibration algorithm of the nuttally quadrinomial third-order window function comprises the following steps:
step C1, the data is processed by FFT operation to be N complex form data x + yj, the modulo p of the N complex data is obtained by calculation,
Figure 601036DEST_PATH_IMAGE001
step C2, finding the data with the maximum module value
Figure 443090DEST_PATH_IMAGE002
Data second largest in sum
Figure 103878DEST_PATH_IMAGE003
And the data corresponds to the sequential position in the N data
Figure 528037DEST_PATH_IMAGE004
And
Figure 897970DEST_PATH_IMAGE005
step C3, calculating calibration coefficient according to the parameters
Figure 492899DEST_PATH_IMAGE006
Figure 239270DEST_PATH_IMAGE007
Step C4, calculating calibration coefficient according to the parameters
Figure 767203DEST_PATH_IMAGE008
Figure 305107DEST_PATH_IMAGE009
Wherein
Figure 669224DEST_PATH_IMAGE010
Figure 202973DEST_PATH_IMAGE011
Figure 336145DEST_PATH_IMAGE012
All four-term third-order calculation of Nuttall
Figure 562727DEST_PATH_IMAGE008
A constant of the coefficient;
step C5, calculating the frequency f,
Figure 414140DEST_PATH_IMAGE013
4. the audio signal measuring method for adaptively improving the frequency measurement accuracy according to claim 1, wherein: in step S6, the value of the data extraction number m needs to satisfy m > = n and is a divisor of the current frequency resolution fp.
5. The audio signal measuring method according to claim 1, wherein the method for performing secondary signal acquisition according to the adjusted data extraction number m in step S6 specifically comprises: reading m data when reading signals each time, storing the last read data in the m data read each time into a memory, and obtaining new N signal data after finishing storing the N data.
CN202210282390.0A 2022-03-22 2022-03-22 Audio signal measuring method for self-adaptively improving frequency measurement precision Pending CN114838809A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210282390.0A CN114838809A (en) 2022-03-22 2022-03-22 Audio signal measuring method for self-adaptively improving frequency measurement precision

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210282390.0A CN114838809A (en) 2022-03-22 2022-03-22 Audio signal measuring method for self-adaptively improving frequency measurement precision

Publications (1)

Publication Number Publication Date
CN114838809A true CN114838809A (en) 2022-08-02

Family

ID=82562509

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210282390.0A Pending CN114838809A (en) 2022-03-22 2022-03-22 Audio signal measuring method for self-adaptively improving frequency measurement precision

Country Status (1)

Country Link
CN (1) CN114838809A (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6298363B1 (en) * 1997-11-13 2001-10-02 The Johns Hopkins University Adaptive windowing of FFT data for increased resolution and sidelobe rejection
CN101219064A (en) * 2006-09-11 2008-07-16 株式会社东芝 Ultrasonic diagnostic device and a method of controlling the ultrasonic diagnostic device
CN103207319A (en) * 2013-03-12 2013-07-17 广东电网公司电力科学研究院 Harmonic wave measurement method of electricity signal of digital substation under non-synchronous sampling condition
CN104391282A (en) * 2014-11-24 2015-03-04 天津大学 Method for improving imaging quality by spectrum correction
CN105307095A (en) * 2015-09-15 2016-02-03 中国电子科技集团公司第四十一研究所 Method for high-resolution audio frequency measurement based on FFT (Fast Fourier Transform)
CN105513608A (en) * 2015-07-23 2016-04-20 中国电子科技集团公司第四十一研究所 Audio signal analysis method
CN110377873A (en) * 2019-08-12 2019-10-25 山东超越数控电子股份有限公司 A kind of signal spectrum analysis method and device based on FFT
CN111308199A (en) * 2020-03-12 2020-06-19 国网湖南省电力有限公司 Double spectral line interpolation DFT harmonic wave analysis method, system and medium based on spectrum resolution self-adaption
US20210022027A1 (en) * 2019-07-17 2021-01-21 SiTune Corporation Spectrum monitoring
CN112541157A (en) * 2020-11-30 2021-03-23 西安精密机械研究所 Signal frequency accurate estimation method
CN113611319A (en) * 2021-04-07 2021-11-05 珠海市杰理科技股份有限公司 Wind noise suppression method, device, equipment and system based on voice component

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6298363B1 (en) * 1997-11-13 2001-10-02 The Johns Hopkins University Adaptive windowing of FFT data for increased resolution and sidelobe rejection
CN101219064A (en) * 2006-09-11 2008-07-16 株式会社东芝 Ultrasonic diagnostic device and a method of controlling the ultrasonic diagnostic device
CN103207319A (en) * 2013-03-12 2013-07-17 广东电网公司电力科学研究院 Harmonic wave measurement method of electricity signal of digital substation under non-synchronous sampling condition
CN104391282A (en) * 2014-11-24 2015-03-04 天津大学 Method for improving imaging quality by spectrum correction
CN105513608A (en) * 2015-07-23 2016-04-20 中国电子科技集团公司第四十一研究所 Audio signal analysis method
CN105307095A (en) * 2015-09-15 2016-02-03 中国电子科技集团公司第四十一研究所 Method for high-resolution audio frequency measurement based on FFT (Fast Fourier Transform)
US20210022027A1 (en) * 2019-07-17 2021-01-21 SiTune Corporation Spectrum monitoring
CN110377873A (en) * 2019-08-12 2019-10-25 山东超越数控电子股份有限公司 A kind of signal spectrum analysis method and device based on FFT
CN111308199A (en) * 2020-03-12 2020-06-19 国网湖南省电力有限公司 Double spectral line interpolation DFT harmonic wave analysis method, system and medium based on spectrum resolution self-adaption
CN112541157A (en) * 2020-11-30 2021-03-23 西安精密机械研究所 Signal frequency accurate estimation method
CN113611319A (en) * 2021-04-07 2021-11-05 珠海市杰理科技股份有限公司 Wind noise suppression method, device, equipment and system based on voice component

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
WANG RUNJIE等: "An Adaptive Cutoff Frequency Selection Approach for Fast Fourier Transform Method and Its Application into Short-Term Traffic Flow Forecasting", ISPRS INTERNATIONAL JOURNAL OF GEO-INFORMATION, vol. 9, no. 12, 31 December 2020 (2020-12-31) *
徐健等: "一种基于FFT的高分辨率音频频率测量方法", 电子质量, no. 2, 29 February 2020 (2020-02-29) *
沙长源;蔡萍;: "多抽样率原理在舰船振动信号分析中的应用与研究", 电子设计工程, no. 20, 20 October 2015 (2015-10-20) *
贾朱植;杨理践;祝洪宇;张吉龙;: "基于高分辨率谱估计的早期转子断条故障诊断", 仪器仪表学报, no. 02, 15 February 2017 (2017-02-15) *

Similar Documents

Publication Publication Date Title
CN104685903B (en) The apparatus and method measured for generating audio disturbances
US9363599B2 (en) Systems and methods for protecting a speaker
CN109313909B (en) Method, device, apparatus and system for evaluating consistency of microphone array
RU2596592C2 (en) Spatial audio processor and method of providing spatial parameters based on acoustic input signal
US7602925B2 (en) Audio feedback processing system
US9357300B2 (en) Systems and methods for protecting a speaker
US20070121955A1 (en) Room acoustics correction device
CN109831733A (en) Test method, device, equipment and the storage medium of audio broadcast performance
US20210306782A1 (en) Method and system of audio device performance testing
US20140321655A1 (en) Sensitivity Calibration Method and Audio Device
Robinson et al. Comparative study of the performance of smartphone-based sound level meter apps, with and without the application of a ½” IEC-61094-4 working standard microphone, to IEC-61672 standard metering equipment in the detection of various problematic workplace noise environments
US9143863B2 (en) Apparatus, systems and methods for inaudibly identifying an audio accessory using spectral shaping
CN104869519B (en) A kind of method and system of test microphone background noise
EP3605535A1 (en) Audio device and computer readable program
JP5627440B2 (en) Acoustic apparatus, control method therefor, and program
Zhu et al. Influence of sound source characteristics in determining objective speech intelligibility metrics
CN114838809A (en) Audio signal measuring method for self-adaptively improving frequency measurement precision
CN112530450A (en) Sample-precision delay identification in the frequency domain
US10965265B2 (en) Method and device for adjusting audio signal, and audio system
EP2763433B1 (en) Apparatus, systems and methods for inaudibly identifying an audio accessory using spectral shaping
CN112927705B (en) Frequency response calibration method and related product
TWI740206B (en) Correction system and correction method of signal measurement
JP2018032931A (en) Acoustic signal processing device, program and method
TWI475870B (en) Testing apparatus and method for telephone apparatus
EP3837862A1 (en) Method of operating an ear level audio system and an ear level audio system

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