CN111829650A - Algorithm for testing frequency sweep parameters by using single-frequency signal combination - Google Patents
Algorithm for testing frequency sweep parameters by using single-frequency signal combinationDownload PDFInfo
- Publication number
- CN111829650A CN111829650ACN202010748510.2ACN202010748510ACN111829650ACN 111829650 ACN111829650 ACN 111829650ACN 202010748510 ACN202010748510 ACN 202010748510ACN 111829650 ACN111829650 ACN 111829650A
- Authority
- CN
- China
- Prior art keywords
- frequency
- filter
- value
- filters
- preset data
- 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
Links
- 238000012360testing methodMethods0.000titleclaimsabstractdescription37
- 238000010408sweepingMethods0.000claimsabstractdescription10
- 238000012545processingMethods0.000claimsdescription18
- 230000015572biosynthetic processEffects0.000claimsdescription15
- 238000012937correctionMethods0.000claimsdescription15
- 238000000034methodMethods0.000claimsdescription15
- 238000003786synthesis reactionMethods0.000claimsdescription15
- 238000001228spectrumMethods0.000claimsdescription14
- 238000004364calculation methodMethods0.000claimsdescription7
- 238000000605extractionMethods0.000claimsdescription3
- 230000009286beneficial effectEffects0.000abstractdescription5
- 210000005069earsAnatomy0.000abstractdescription5
- 230000011218segmentationEffects0.000abstractdescription2
- 238000010586diagramMethods0.000description3
- 238000005259measurementMethods0.000description3
- 230000007704transitionEffects0.000description2
- 238000000354decomposition reactionMethods0.000description1
- 230000000694effectsEffects0.000description1
- 238000005516engineering processMethods0.000description1
- 238000011156evaluationMethods0.000description1
- 230000006870functionEffects0.000description1
- 238000013507mappingMethods0.000description1
- 238000000691measurement methodMethods0.000description1
- 238000005065miningMethods0.000description1
- 238000005070samplingMethods0.000description1
- 238000004904shorteningMethods0.000description1
Images
Classifications
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
- G10L25/18—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
- G10L25/51—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
Technical Field
The invention relates to an acoustic testing method, in particular to an algorithm for testing frequency sweep parameters by using single-frequency signal combination.
Background
The measurement of the early frequency characteristics is achieved by a point-by-point mapping method. In the whole measurement process, the amplitude of a signal input to the measured network is kept unchanged, correspondingly output voltages under different frequencies are recorded, and an amplitude-frequency characteristic curve of the network can be drawn on coordinate paper according to obtained data. Obviously, this method is not only cumbersome and time consuming, but also may miss some details on the measured curve due to the insufficiently dense measurement frequency intervals, making the resulting curve less accurate.
In the sweep measurement method, a constant-amplitude sweep signal is added to the input end of a tested circuit, and then an oscilloscope is used for displaying the change of the amplitude of the signal after the signal passes through the tested circuit. Because the frequency of the sweep frequency signal is continuously changed, the amplitude-frequency characteristic of the tested circuit can be directly displayed on the screen of the oscilloscope.
In summary, the existing frequency characteristic test has the following disadvantages: 1. the length of the frequency sweeping signal of the full frequency band is too long; 2. the common single-frequency set is not beneficial to calculation and has longer time; 3. the test does not conform to the hearing habits of the human ear.
Disclosure of Invention
The invention provides an algorithm for testing frequency sweep parameters by using a single-frequency signal combination, which solves the technical problems that the too long time of frequency characteristic testing in the related technology is not beneficial to calculation and the testing does not conform to human auditory sense.
According to one aspect of the present invention, there is provided an algorithm for testing sweep frequency parameters using a single frequency signal combination, comprising the steps of:
s100, generating a sweep frequency signal of a single-frequency signal combination for output;
s200, after the played sweep frequency signal is acquired, dividing the sweep frequency signal by contrasting the generated single-frequency signal combination, and calculating the divided signal in a single-frequency signal mode;
and S300, displaying the calculated acoustic parameters of each single-frequency signal on an xy two-dimensional coordinate graph with the frequency point as an x axis to obtain a result of testing the frequency sweeping parameters.
By adopting the technical scheme, the frequency sweep signal of the single-frequency signal combination is generated for outputting, the same segmentation signal is calculated in a single-frequency signal mode, and the expected frequency sweep signal is set by inputting the preset data, so that the technical problems that the existing frequency characteristic test full-frequency band frequency sweep signal is too long in time, a common single-frequency set is not beneficial to calculation and is very long in time, and the test does not accord with the hearing habit of human ears are solved.
According to another aspect of the present invention, there is provided a system for testing sweep frequency parameters using a single frequency signal combination, comprising:
an input unit for inputting preset data;
the processing unit receives preset data input by the input unit and processes the preset data to obtain a coefficient table;
sending the coefficient table obtained by processing to a variable filter unit and a display unit;
the display unit is used for displaying and processing the coefficient table;
a variable filter unit having n-1 filters, setting parameters of the filters being determined by filter parameters in the coefficient table; the circuit is used for sending the frequency sweeping signal to an output amplifier and outputting the frequency sweeping signal to the circuit;
the back-extraction unit is used for acquiring the signals output by the circuit and sending the signals to the graphic processor;
and the graphic processor is used for processing the signals output by the circuit to obtain acoustic parameters and displaying the acoustic parameters on an xy two-dimensional coordinate graph with the frequency point as an x axis.
Further, the preset data comprises a plurality of groups of frequencies and gains; the user sets the desired swept frequency signal by inputting gain data at a particular frequency point.
Further, the method for processing the coefficient table by the processing unit comprises the following steps:
s1, receiving preset data input by the input unit, and calculating an output waveform amplitude spectrum of a filter;
s2, calculating a synthesis characteristic of the full frequency band of the filter for obtaining the output waveform amplitude spectrum in the step S1, thereby obtaining a synthesis characteristic of n-1 filters, and calculating an error value of the synthesis characteristic of each of the n-1 filters and preset data, so as to obtain an error of the frequency band corresponding to each filter;
s3, determining the types of n-1 filters, including whether three gains of preset data corresponding to three frequency points of a frequency band corresponding to the filters meet preset conditions or not, and determining the filter type as a notch filter if the preset conditions are met; if the preset condition is not met, determining the filter type as a peak filter;
the preset conditions are as follows: the absolute value of the maximum value in all gains in the preset data is smaller than the absolute value of the minimum value or the absolute value of the maximum value and the minimum value are both negative values;
s4, correcting the center frequency of the filter of the notch filter type in n-1 filters, calculating the frequency with the maximum difference value with the preset data in the synthesis characteristic of the filter, judging whether the result of subtracting the center frequency of the filter is a positive value or a negative value, and if the result is the positive value, subtracting a preset correction value from the center frequency of the filter to obtain the corrected center frequency;
if the value is negative, adding a preset correction value to the center frequency of the filter to obtain a corrected center frequency;
s5, correcting the center frequency of the filter of the peak filter type in n-1 filters, calculating the frequency with the maximum difference value with the preset data in the synthesis characteristic of the filter, judging whether the result of subtracting the center frequency of the filter is a positive value or a negative value, and if the result is the positive value, adding a preset correction value to the center frequency of the filter to obtain the corrected center frequency;
if the value is negative, subtracting a preset correction value from the center frequency of the filter to obtain a corrected center frequency;
s6, calculating filter parameters of all filters of n-1 filters, including calculating peak factors of the filters through the correction center frequency of the filters and the bandwidth of the corresponding frequency band, and taking the gain corresponding to the frequency in the preset data with the same frequency of the middle frequency point of the frequency band corresponding to the filters as the filter gain of the filters; taking the corrected center frequency, the peak factor and the filter gain as filter coefficients;
and recording the filter parameters of all the filters of the n-1 filters obtained by calculation to obtain a coefficient table.
Further, the step S1 of receiving the preset data input by the input unit, and calculating the output waveform amplitude spectrum of a filter includes the following steps:
s111, judging whether all gains in preset data meet preset conditions or not; the preset conditions are as follows: the absolute value of the maximum value in all gains in the preset data is smaller than the absolute value of the minimum value or the absolute value of the maximum value and the minimum value are both negative values;
s112, if the requirement is met, the filter type is determined to be a notch filter, the center frequency is determined to be the frequency with the maximum value, and the filter gain is determined to be the maximum value in all gains of the preset data;
s113, if the data do not meet the requirement, determining the filter type to be a peak filter, determining the center frequency to be the frequency of taking the minimum value, and determining the gain of the filter to be the minimum value in all gains of the preset data;
s114, calculating the difference value between the maximum value and the minimum value in all gains in the preset data as a gain difference;
s115, judging whether the gain difference is less than 3 db;
s116, if the gain difference is larger than 3db, taking the difference between the maximum value and the minimum value of all frequencies in the preset data as a temporary bandwidth; calculating a crest factor of the filter using the center frequency and the temporary bandwidth;
s117, if the gain difference is less than 3db, setting the crest factor to a predetermined value without calculating the crest factor;
and S118, calculating an output waveform amplitude spectrum of the filter through the center frequency, the crest factor and the filter gain.
Further, the predetermined value in the step S117 is 0.25 or 0.5.
Furthermore, the n-1 filters correspond to a frequency range, the frequency range is divided into a plurality of frequency bands, each frequency band comprises three frequency points of preset data, each frequency band corresponds to one filter of the n-1 filters, and an error value of a synthetic characteristic of the filter of each frequency band and the preset data is calculated to obtain an error of each frequency band.
Further, the correction amount in steps S4 and S5 is 0.08 times the bandwidth of the error band of the frequency band corresponding to the filter.
The invention has the beneficial effects that:
firstly, shortening the time length of a sweep signal to shorten the test time length and simultaneously ensure that a frequency band can be covered;
secondly, smooth transition from low frequency to high frequency according with the hearing habits of human ears;
thirdly, ensuring that the interference of other frequency bands is avoided; to ensure the test stability.
Drawings
Fig. 1 is a flowchart of an algorithm for testing sweep parameters using a single frequency signal combination according to an embodiment of the present invention;
FIG. 2 is a block diagram of a system for testing frequency sweep parameters using single frequency signal combinations in accordance with an embodiment of the present invention;
FIG. 3 is a flow diagram of a method for a processing unit to process a derived coefficient table according to an embodiment of the invention;
FIG. 4 is a flow chart of the calculation of the output waveform amplitude spectrum of a filter according to an embodiment of the present invention;
fig. 5 is a schematic diagram of an implementation of an algorithm for testing sweep parameters by using a single-frequency signal combination according to an embodiment of the present invention.
Detailed Description
The subject matter described herein will now be discussed with reference to example embodiments. It should be understood that these embodiments are discussed only to enable those skilled in the art to better understand and thereby implement the subject matter described herein, and are not intended to limit the scope, applicability, or examples set forth in the claims. Changes may be made in the function and arrangement of elements discussed without departing from the scope of the disclosure. Various examples may omit, substitute, or add various procedures or components as needed. For example, the described methods may be performed in an order different from that described, and various steps may be added, omitted, or combined. In addition, features described with respect to some examples may also be combined in other examples.
In this embodiment, an algorithm for testing sweep frequency parameters by using a single frequency signal combination is provided, as shown in fig. 1 and 5, a flowchart according to the present invention is shown, as shown in fig. 1 and 5, and the method includes: s100, generating a sweep frequency signal of a single-frequency signal combination for output; generating a sinusoidal single-frequency signal of y ═ Asin (ω x + φ), decomposing the signal of the whole frequency band into a plurality of sinusoidal single-frequency signals of y ═ Asin (ω x + φ) according to the decomposition of an octave, setting the maximum duration 10s and the minimum period 10T according to T ═ 1/f under the condition of setting the sampling rate, and adjusting the period or duration of the generated single-frequency signals to ensure that the duration of frequency sweeping signals and the test stability achieve the effect of relative balance;
s200, after the played sweep frequency signal is acquired, dividing the sweep frequency signal by contrasting the generated single-frequency signal combination, and calculating the divided signal in a single-frequency signal mode; after data are acquired through back mining, dividing the data according to the number of the data by contrasting the generated single-frequency signal combination according to a one-to-one correspondence principle, and then independently calculating each single-frequency signal;
and S300, displaying the calculated acoustic parameters of each single-frequency signal on an xy two-dimensional coordinate graph with the frequency point as an x axis to obtain a result of testing the frequency sweeping parameters.
The algorithm of this embodiment can be implemented by various forms of systems, and as shown in fig. 2, this embodiment provides a system for testing sweep parameters by using a single-frequency signal combination, including:
aninput unit 1 for inputting preset data, optionally but not limited to: a keyboard, a touch screen;
the preset data comprises a plurality of groups of frequencies (Hz) and gains (db); a user sets a desired frequency sweep signal by inputting gain data at a specific frequency point;
the processing unit 2 receives the preset data input by theinput unit 1 and processes the preset data to obtain a coefficient table;
and sends the coefficient table obtained by the processing to thevariable filter unit 4 and thedisplay unit 5;
adisplay unit 5 for displaying the processed coefficient table;
avariable filter unit 4 having n-1 filters, setting parameters of each filter including a center frequency (f0), a filter crest factor (Q), and a filter gain (G), the setting parameters of the filter being determined by filter parameters in a coefficient table;
thevariable filter unit 4 sends the sweep frequency signal to theoutput amplifier 6 and outputs the sweep frequency signal to the circuit;
theextraction unit 7 is used for collecting the signals output by the circuit and sending the signals to thegraphic processor 8;
and thegraphic processor 8 is used for processing the signals output by the circuit to obtain acoustic parameters and displaying the acoustic parameters on an xy two-dimensional coordinate graph with the frequency point as an x axis.
As shown in fig. 3, the method for processing the coefficient table by the processing unit 2 includes the following steps:
s1, receiving preset data input by theinput unit 1, and calculating an output waveform amplitude spectrum of a filter;
as shown in fig. 4, step S1 specifically includes the following steps: s111, judging whether all gains in preset data meet preset conditions or not; the preset conditions are as follows: the absolute value of the maximum value in all gains in the preset data is smaller than the absolute value of the minimum value or the absolute value of the maximum value and the minimum value are both negative values; (the gain of the preset data may be a positive value or a negative value)
S112, if the requirement is met, the filter type is determined to be a notch filter, the center frequency is determined to be the frequency with the maximum value, and the filter gain is determined to be the maximum value in all gains of the preset data;
s113, if the data do not meet the requirement, determining the filter type to be a peak filter, determining the center frequency to be the frequency of taking the minimum value, and determining the gain of the filter to be the minimum value in all gains of the preset data;
s114, calculating the difference value between the maximum value and the minimum value in all gains in the preset data as a gain difference; for example, the maximum value of all gains of the preset data is +9dB, the minimum value is +4dB, the filter type is a peak filter, and the gain difference is 5.
S115, judging whether the gain difference is less than 3 db;
s116, if the gain difference is larger than 3db, taking the difference between the maximum value and the minimum value of all frequencies in the preset data as a temporary bandwidth; calculating a crest factor of the filter using the center frequency and the temporary bandwidth;
s117, if the gain difference is less than 3db, setting the peak factor to a predetermined value without calculating the peak factor (peak factor Q is 0.25 or 0.5);
s118, calculating an output waveform amplitude spectrum of the filter through the center frequency, the peak factor and the filter gain; the output waveform amplitude spectrum of the filter can be converted into a frequency-gain table and stored in a preset table;
correspondingly obtaining the gain of the frequency in the preset data by the output waveform amplitude spectrum of the filter;
s2, calculating the synthesis characteristic of the full frequency band (20Hz to 20000Hz) of the filter which obtains the output waveform amplitude spectrum in the step S1, thereby obtaining the synthesis characteristics of other n-1 filters, and calculating the error value of the synthesis characteristic of each filter in the n-1 filters and preset data, so as to obtain the error of the frequency band corresponding to each filter;
dividing a frequency range to be set into a plurality of frequency bands, wherein each frequency band comprises three frequency points of preset data, each frequency band corresponds to one filter, and calculating the synthetic characteristic of the filter of each frequency band and the error value of the preset data to obtain the error of each frequency band; for example, a band is composed of 70Hz-130Hz-230Hz, the synthetic characteristics are 70Hz-3db, 130Hz-5db, 230Hz-7db, the preset data are 70Hz-4db, 130Hz-6db, 230Hz-8db, and the error of the band is the sum of the differences of 4db-3db, 6db-5db, 8db-7db, namely 3. Of course the error may also be negative.
And the large error of the evaluation only in frequency points is reduced by evaluating in each frequency band.
S3, determining the types of n-1 filters, including whether three gains of preset data corresponding to three frequency points of a frequency band corresponding to the filters meet preset conditions or not, and determining the filter type as a notch filter if the preset conditions are met; if the preset condition is not met, determining the filter type as a peak filter;
the preset conditions are as follows: the absolute value of the maximum value in all gains in the preset data is smaller than the absolute value of the minimum value or the absolute value of the maximum value and the minimum value are both negative values;
s4, correcting the center frequency of the filter of the notch filter type in n-1 filters, calculating the frequency with the maximum difference value with the preset data in the synthesis characteristic of the filter, judging whether the result of subtracting the center frequency of the filter is a positive value or a negative value, and if the result is the positive value, subtracting a preset correction value from the center frequency of the filter to obtain the corrected center frequency;
if the value is negative, adding a preset correction value to the center frequency of the filter to obtain a corrected center frequency;
s5, correcting the center frequency of the filter of the peak filter type in n-1 filters, calculating the frequency with the maximum difference value with the preset data in the synthesis characteristic of the filter, judging whether the result of subtracting the center frequency of the filter is a positive value or a negative value, and if the result is the positive value, adding a preset correction value to the center frequency of the filter to obtain the corrected center frequency;
if the value is negative, subtracting a preset correction value from the center frequency of the filter to obtain a corrected center frequency;
s6, calculating filter parameters of all filters of n-1 filters, including calculating peak factors of the filters through the correction center frequency of the filters and the bandwidth of the corresponding frequency band, and taking the gain corresponding to the frequency in the preset data with the same frequency of the middle frequency point of the frequency band corresponding to the filters as the filter gain of the filters; taking the corrected center frequency, the peak factor and the filter gain as filter coefficients;
recording the filter parameters of all the filters of the n-1 filters obtained by calculation to obtain a coefficient table;
the correction value is the bandwidth of the error frequency band of the frequency band corresponding to the filter, namely 0.08;
according to the method, the expected frequency sweep signal is set by inputting the preset data, the interference of redundant frequency bands is reduced, and the duration of the frequency sweep signal is shortened by adjusting the frequency sweep signal; the low-frequency to high-frequency smooth transition which accords with the auditory habits of human ears can be tested through the setting of the preset data, namely the low-frequency band and the high-frequency band of the sweep frequency signal can be adjusted according with the auditory habits of human ears.
Claims (8)
1. An algorithm for testing swept frequency parameters using a combination of single frequency signals, comprising the steps of:
s100, generating a sweep frequency signal of a single-frequency signal combination for output;
s200, after the played sweep frequency signal is acquired, dividing the sweep frequency signal by contrasting the generated single-frequency signal combination, and calculating the divided signal in a single-frequency signal mode;
and S300, displaying the calculated acoustic parameters of each single-frequency signal on an xy two-dimensional coordinate graph with the frequency point as an x axis to obtain a result of testing the frequency sweeping parameters.
2. A system for testing swept frequency parameters using a combination of single frequency signals, comprising:
an input unit for inputting preset data;
the processing unit receives preset data input by the input unit and processes the preset data to obtain a coefficient table;
sending the coefficient table obtained by processing to a variable filter unit and a display unit;
the display unit is used for displaying and processing the coefficient table;
a variable filter unit having n-1 filters, setting parameters of the filters being determined by filter parameters in the coefficient table; the circuit is used for sending the frequency sweeping signal to an output amplifier and outputting the frequency sweeping signal to the circuit;
the back-extraction unit is used for acquiring the signals output by the circuit and sending the signals to the graphic processor;
and the graphic processor is used for processing the signals output by the circuit to obtain acoustic parameters and displaying the acoustic parameters on an xy two-dimensional coordinate graph with the frequency point as an x axis.
3. A system for testing swept frequency parameters using single frequency signal combinations as claimed in claim 2, wherein the preset data comprises multiple sets of frequencies and gains; the user sets the desired swept frequency signal by inputting gain data at a particular frequency point.
4. A system for testing swept frequency parameters using single frequency signal combinations according to claim 2, wherein the method for processing the resulting coefficient table by the processing unit comprises the steps of:
s1, receiving preset data input by the input unit, and calculating an output waveform amplitude spectrum of a filter;
s2, calculating a synthesis characteristic of the full frequency band of the filter for obtaining the output waveform amplitude spectrum in the step S1, thereby obtaining a synthesis characteristic of n-1 filters, and calculating an error value of the synthesis characteristic of each of the n-1 filters and preset data, so as to obtain an error of the frequency band corresponding to each filter;
s3, determining the types of n-1 filters, including whether three gains of preset data corresponding to three frequency points of a frequency band corresponding to the filters meet preset conditions or not, and determining the filter type as a notch filter if the preset conditions are met; if the preset condition is not met, determining the filter type as a peak filter;
the preset conditions are as follows: the absolute value of the maximum value in all gains in the preset data is smaller than the absolute value of the minimum value or the absolute value of the maximum value and the minimum value are both negative values;
s4, correcting the center frequency of the filter of the notch filter type in n-1 filters, calculating the frequency with the maximum difference value with the preset data in the synthesis characteristic of the filter, judging whether the result of subtracting the center frequency of the filter is a positive value or a negative value, and if the result is the positive value, subtracting a preset correction value from the center frequency of the filter to obtain the corrected center frequency;
if the value is negative, adding a preset correction value to the center frequency of the filter to obtain a corrected center frequency;
s5, correcting the center frequency of the filter of the peak filter type in n-1 filters, calculating the frequency with the maximum difference value with the preset data in the synthesis characteristic of the filter, judging whether the result of subtracting the center frequency of the filter is a positive value or a negative value, and if the result is the positive value, adding a preset correction value to the center frequency of the filter to obtain the corrected center frequency;
if the value is negative, subtracting a preset correction value from the center frequency of the filter to obtain a corrected center frequency;
s6, calculating filter parameters of all filters of n-1 filters, including calculating peak factors of the filters through the correction center frequency of the filters and the bandwidth of the corresponding frequency band, and taking the gain corresponding to the frequency in the preset data with the same frequency of the middle frequency point of the frequency band corresponding to the filters as the filter gain of the filters; taking the corrected center frequency, the peak factor and the filter gain as filter coefficients;
and recording the filter parameters of all the filters of the n-1 filters obtained by calculation to obtain a coefficient table.
5. A system for testing swept frequency parameters by using single frequency signal combination according to claim 4, wherein the step S1 is to receive preset data inputted from the input unit, and the calculating of the output waveform amplitude spectrum of a filter comprises the following steps:
s111, judging whether all gains in preset data meet preset conditions or not; the preset conditions are as follows: the absolute value of the maximum value in all gains in the preset data is smaller than the absolute value of the minimum value or the absolute value of the maximum value and the minimum value are both negative values;
s112, if the requirement is met, the filter type is determined to be a notch filter, the center frequency is determined to be the frequency with the maximum value, and the filter gain is determined to be the maximum value in all gains of the preset data;
s113, if the data do not meet the requirement, determining the filter type to be a peak filter, determining the center frequency to be the frequency of taking the minimum value, and determining the gain of the filter to be the minimum value in all gains of the preset data;
s114, calculating the difference value between the maximum value and the minimum value in all gains in the preset data as a gain difference;
s115, judging whether the gain difference is less than 3 db;
s116, if the gain difference is larger than 3db, taking the difference between the maximum value and the minimum value of all frequencies in the preset data as a temporary bandwidth; calculating a crest factor of the filter using the center frequency and the temporary bandwidth;
s117, if the gain difference is less than 3db, setting the crest factor to a predetermined value without calculating the crest factor;
and S118, calculating an output waveform amplitude spectrum of the filter through the center frequency, the crest factor and the filter gain.
6. A system for testing swept frequency parameters using single frequency signal combination according to claim 1, wherein the predetermined value in step S117 is 0.25 or 0.5.
7. The system for testing frequency sweep parameters by using single frequency signal combination as claimed in claim 4, wherein the n-1 filters correspond to a frequency range, the frequency range is divided into a plurality of frequency bands, each frequency band comprises three frequency points of preset data, each frequency band corresponds to one filter of the n-1 filters, and an error value of the synthesized characteristic of the filter of each frequency band and the preset data is calculated to obtain an error of each frequency band.
8. A system for testing swept frequency parameters using single frequency signal combinations according to claim 4, wherein the amount of correction in steps S4 and S5 is 0.08 times the bandwidth of the error frequency band of the frequency band corresponding to the filter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010748510.2ACN111829650B (en) | 2020-07-30 | 2020-07-30 | Algorithm for testing frequency sweep parameters by using single-frequency signal combination |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010748510.2ACN111829650B (en) | 2020-07-30 | 2020-07-30 | Algorithm for testing frequency sweep parameters by using single-frequency signal combination |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111829650Atrue CN111829650A (en) | 2020-10-27 |
| CN111829650B CN111829650B (en) | 2021-09-03 |
Family
ID=72920366
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010748510.2AActiveCN111829650B (en) | 2020-07-30 | 2020-07-30 | Algorithm for testing frequency sweep parameters by using single-frequency signal combination |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111829650B (en) |
Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5541866A (en)* | 1991-11-28 | 1996-07-30 | Kabushiki Kaisha Kenwood | Device for correcting frequency characteristic of sound field |
| US20020134612A1 (en)* | 2001-03-26 | 2002-09-26 | Khan Tawassul A. | Mapping reservoir characteristics using earth's nonlinearity as a seismic attribute |
| US20090296265A1 (en)* | 2008-05-30 | 2009-12-03 | Masashi Kisaka | Hard-disk drive |
| CN101715200A (en)* | 2009-11-04 | 2010-05-26 | 武汉虹翼信息有限公司 | Method for analyzing situation of mobile network by drive test data of sweep signal generator |
| CN101711969A (en)* | 2009-10-30 | 2010-05-26 | 江苏大学 | Device and method for biological processing with pulsed multi-frequency ultrasonic waves in sweep-frequency mode and use thereof |
| CN103823136A (en)* | 2014-02-24 | 2014-05-28 | 北京航空航天大学 | Multifrequency quick sensitivity test frequency point combination method based on power balance |
| CN104545940A (en)* | 2014-12-31 | 2015-04-29 | 深圳先进技术研究院 | Auditory feedback regulating method and system based on stimulus frequency otoacoustic emission |
| CN104935258A (en)* | 2014-03-18 | 2015-09-23 | 苏州普源精电科技有限公司 | A Sweeping Signal Generator That Can Generate Multiple Frequency Standards |
| US20150289073A1 (en)* | 2014-01-30 | 2015-10-08 | Invensense, Inc. | Electrical testing and feedthrough cancellation for an acoustic sensor |
| CN105842566A (en)* | 2016-05-19 | 2016-08-10 | 中国第汽车股份有限公司 | Automatic digital filter amplitude frequency characteristic curve mapping system and mapping method thereof |
| CN205898894U (en)* | 2016-08-17 | 2017-01-18 | 南京信息工程大学 | A Digital Storage Scanner |
| CN108107350A (en)* | 2018-01-05 | 2018-06-01 | 桂林电子科技大学 | A kind of long-range wide-band amplitude versus frequency characte measuring instrument |
| CN108112056A (en)* | 2016-11-24 | 2018-06-01 | 联芯科技有限公司 | The frequency sweeping method and device of user terminal |
| CN108153242A (en)* | 2017-12-23 | 2018-06-12 | 西安交通大学 | A kind of feed system servo controller and its control method and parameter setting method |
| CN108988979A (en)* | 2018-07-10 | 2018-12-11 | 深圳磊诺科技有限公司 | A kind of full frequency band frequency sweeping method, device and small base station |
| CN109039497A (en)* | 2017-06-09 | 2018-12-18 | 中兴通讯股份有限公司 | Simple signal combination distribution method, data emitting method, method of reseptance and equipment |
| CN109188453A (en)* | 2018-09-11 | 2019-01-11 | 哈尔滨工业大学 | Dynamic frequency scanning interfeerometry ranging system and distance measuring method based on locking phase gamma correction |
| CN109508511A (en)* | 2018-12-24 | 2019-03-22 | 中国航空工业集团公司西安飞机设计研究所 | Frequency sweeping method in frequency response analysis measurement |
| CN109975822A (en)* | 2019-04-16 | 2019-07-05 | 重庆大学 | Merge the optical-fiber laser dynamic distance system of frequency sweep and single-frequency |
| CN110687213A (en)* | 2019-08-26 | 2020-01-14 | 天津大学 | A Selection Method of Optimal Excitation Frequency for Nonlinear Ultrasonic Mixing |
| CN111024819A (en)* | 2019-12-23 | 2020-04-17 | 复旦大学 | A kind of jewelry and jewelry acoustic resonance spectrum imitation method |
| CN111261188A (en)* | 2020-01-20 | 2020-06-09 | 中国电力科学研究院有限公司 | Method and device for determining noise spectrum of high voltage transformer |
- 2020
- 2020-07-30CNCN202010748510.2Apatent/CN111829650B/enactiveActive
Patent Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5541866A (en)* | 1991-11-28 | 1996-07-30 | Kabushiki Kaisha Kenwood | Device for correcting frequency characteristic of sound field |
| US20020134612A1 (en)* | 2001-03-26 | 2002-09-26 | Khan Tawassul A. | Mapping reservoir characteristics using earth's nonlinearity as a seismic attribute |
| US20090296265A1 (en)* | 2008-05-30 | 2009-12-03 | Masashi Kisaka | Hard-disk drive |
| CN101711969A (en)* | 2009-10-30 | 2010-05-26 | 江苏大学 | Device and method for biological processing with pulsed multi-frequency ultrasonic waves in sweep-frequency mode and use thereof |
| CN101715200A (en)* | 2009-11-04 | 2010-05-26 | 武汉虹翼信息有限公司 | Method for analyzing situation of mobile network by drive test data of sweep signal generator |
| US20150289073A1 (en)* | 2014-01-30 | 2015-10-08 | Invensense, Inc. | Electrical testing and feedthrough cancellation for an acoustic sensor |
| CN103823136A (en)* | 2014-02-24 | 2014-05-28 | 北京航空航天大学 | Multifrequency quick sensitivity test frequency point combination method based on power balance |
| CN104935258A (en)* | 2014-03-18 | 2015-09-23 | 苏州普源精电科技有限公司 | A Sweeping Signal Generator That Can Generate Multiple Frequency Standards |
| CN104545940A (en)* | 2014-12-31 | 2015-04-29 | 深圳先进技术研究院 | Auditory feedback regulating method and system based on stimulus frequency otoacoustic emission |
| CN105842566A (en)* | 2016-05-19 | 2016-08-10 | 中国第汽车股份有限公司 | Automatic digital filter amplitude frequency characteristic curve mapping system and mapping method thereof |
| CN205898894U (en)* | 2016-08-17 | 2017-01-18 | 南京信息工程大学 | A Digital Storage Scanner |
| CN108112056A (en)* | 2016-11-24 | 2018-06-01 | 联芯科技有限公司 | The frequency sweeping method and device of user terminal |
| CN109039497A (en)* | 2017-06-09 | 2018-12-18 | 中兴通讯股份有限公司 | Simple signal combination distribution method, data emitting method, method of reseptance and equipment |
| CN108153242A (en)* | 2017-12-23 | 2018-06-12 | 西安交通大学 | A kind of feed system servo controller and its control method and parameter setting method |
| CN108107350A (en)* | 2018-01-05 | 2018-06-01 | 桂林电子科技大学 | A kind of long-range wide-band amplitude versus frequency characte measuring instrument |
| CN108988979A (en)* | 2018-07-10 | 2018-12-11 | 深圳磊诺科技有限公司 | A kind of full frequency band frequency sweeping method, device and small base station |
| CN109188453A (en)* | 2018-09-11 | 2019-01-11 | 哈尔滨工业大学 | Dynamic frequency scanning interfeerometry ranging system and distance measuring method based on locking phase gamma correction |
| CN109508511A (en)* | 2018-12-24 | 2019-03-22 | 中国航空工业集团公司西安飞机设计研究所 | Frequency sweeping method in frequency response analysis measurement |
| CN109975822A (en)* | 2019-04-16 | 2019-07-05 | 重庆大学 | Merge the optical-fiber laser dynamic distance system of frequency sweep and single-frequency |
| CN110687213A (en)* | 2019-08-26 | 2020-01-14 | 天津大学 | A Selection Method of Optimal Excitation Frequency for Nonlinear Ultrasonic Mixing |
| CN111024819A (en)* | 2019-12-23 | 2020-04-17 | 复旦大学 | A kind of jewelry and jewelry acoustic resonance spectrum imitation method |
| CN111261188A (en)* | 2020-01-20 | 2020-06-09 | 中国电力科学研究院有限公司 | Method and device for determining noise spectrum of high voltage transformer |
Non-Patent Citations (5)
| Title |
|---|
| BIN LIU: "《Analysis of sound and vibration interaction on a crack and its use in high-frequency parameter selection for vibro-acoustic modulation testing》", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》* |
| ION CORNEL MITULETU: "《A method for an accurate estimation of natural frequencies using swept-sine acoustic excitation 》", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》* |
| 杨文伟: "《基于压电陶瓷的PVA-ECC混凝土柱裂缝损伤监测》", 《沈阳建筑大学学报(自然科学版)》* |
| 牛鹏宇: "《基于Matlab的低频频率特性测试虚拟仪器设计》", 《测控技术》* |
| 薛媛元: "《低频频率特性测试仪虚拟仪器设计》", 《测控技术》* |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111829650B (en) | 2021-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5450049B2 (en) | Method and apparatus in an acoustic system | |
| CN101241150B (en) | Apparatus, method for processing signal and method for generating signal | |
| US4607216A (en) | Apparatus for measurement by digital spectrum analyzer | |
| US10461711B2 (en) | Method and apparatus for outputting audio signal, method for controlling volume | |
| JP3640641B2 (en) | Method and apparatus for generating a calibration sound field | |
| US9955261B2 (en) | Method and apparatus for adjusting a cross-over frequency of a loudspeaker | |
| CN113116333B (en) | Method and system for detecting brain impedance distribution difference data | |
| CN101374374A (en) | Audio test device and method | |
| CN110557711B (en) | Earphone testing method and earphone | |
| CN111829650B (en) | Algorithm for testing frequency sweep parameters by using single-frequency signal combination | |
| CN107124678B (en) | A Measurement System of Audio Harmonic Distortion | |
| US4421121A (en) | Method and apparatus for obtaining a non-cephalic referential electroencephalogram | |
| US9807519B2 (en) | Method and apparatus for analyzing and visualizing the performance of frequency lowering hearing aids | |
| CN114760567A (en) | Method, device and equipment for calibrating sound quality effect of earphone and storage medium | |
| US20210249028A1 (en) | Method and an audio processing unit for detecting a tone | |
| CN118413797A (en) | Abnormal sound detection method, system and equipment | |
| CN118335113A (en) | Abnormal sound detection method, system, equipment and storage medium | |
| US9693171B2 (en) | Sound field measuring device, method, and program | |
| CN112082641B (en) | A correction processing method for vibration measurement data of vibration sensor | |
| CN108877902A (en) | Ultrasound image brightness adjusting method and regulating system | |
| CN109391875B (en) | Audio debugging method and device, storage medium and processor | |
| JP3272991B2 (en) | How to correct sound field characteristics | |
| EP1547242A2 (en) | Device and method for determining the level of an input signal intended to be applied to a receiving system | |
| WO2021135235A1 (en) | Automatic debugging method and system for loudspeaker | |
| CN118828341B (en) | Visualization of three-dimensional audio data as frequency-angle 3D model and data processing method thereof |
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 |