CN111829650B - 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 combination Download PDFInfo
- Publication number
- CN111829650B CN111829650B CN202010748510.2A CN202010748510A CN111829650B CN 111829650 B CN111829650 B CN 111829650B CN 202010748510 A CN202010748510 A CN 202010748510A CN 111829650 B CN111829650 B CN 111829650B
- 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.)
- Active
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 35
- 238000010408 sweeping Methods 0.000 claims abstract description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 20
- 238000003786 synthesis reaction Methods 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 5
- 210000005069 ears Anatomy 0.000 abstract description 5
- 230000011218 segmentation Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000007704 transition Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
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
The invention discloses an algorithm for testing frequency sweep parameters by using a single-frequency signal combination, which relates to the technical field of acoustic testing and comprises the following steps: generating a frequency sweep signal of the single frequency signal combination for output; after the played sweep frequency signal is collected, the sweep frequency signal is divided by contrasting the generated single-frequency signal combination, and the divided signal is calculated in a single-frequency signal mode; then 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 invention generates the sweep frequency signal of the single-frequency signal combination for outputting, calculates the same segmentation signal in a single-frequency signal mode, and sets the expected sweep frequency signal by inputting the preset data, thereby solving the technical problems that the existing frequency characteristic test full-frequency band sweep frequency signal has overlong duration, a common single-frequency set is not beneficial to calculating and has stronger duration at the same time, and the test does not accord with the auditory habits of human ears.
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:
an input 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 the input unit 1 and processes the preset data to obtain a coefficient table;
and sends the coefficient table obtained by the processing to the variable filter unit 4 and the display unit 5;
a display unit 5 for displaying the processed coefficient table;
a variable 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;
the variable filter unit 4 sends the sweep frequency signal to the output amplifier 6 and outputs the sweep frequency signal to the circuit;
the extraction unit 7 is used for collecting the signals output by the circuit and sending the signals to the graphic processor 8;
and the graphic 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 the input 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 (6)
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;
inputting preset data, and processing the preset data to obtain a coefficient table;
the method for processing the preset data to obtain the coefficient table 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;
recording the filter parameters of all the filters of the n-1 filters obtained by calculation to obtain a coefficient table;
sending the processed coefficient table to a variable filter unit with n-1 filters, and sending the frequency sweep signal to an output amplifier and outputting the frequency sweep signal to a circuit by the variable filter unit;
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;
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;
the preset data comprises a plurality of groups of frequencies and gains; a user sets a desired frequency sweep signal by inputting gain data at a specific frequency point;
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.
3. A system for testing swept frequency parameters by using single frequency signal combination according to claim 2, wherein the step S1 is implemented by receiving preset data input from the input unit, and the calculating an 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.
4. A system for testing swept frequency parameters using single frequency signal combination according to claim 3, wherein the predetermined value in step S117 is 0.25 or 0.5.
5. A system for testing swept frequency parameters using single frequency signal combination according to claim 2, 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.
6. A system for testing swept frequency parameters using single frequency signal combinations as claimed in claim 2, 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.2A CN111829650B (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.2A CN111829650B (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 |
---|---|
CN111829650A CN111829650A (en) | 2020-10-27 |
CN111829650B true CN111829650B (en) | 2021-09-03 |
Family
ID=72920366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010748510.2A Active CN111829650B (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 (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101715200A (en) * | 2009-11-04 | 2010-05-26 | 武汉虹翼信息有限公司 | Method for analyzing situation of mobile network by drive test data of sweep signal generator |
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 | 苏州普源精电科技有限公司 | Sweep signal generator capable of generating a plurality of frequency markers |
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 | 南京信息工程大学 | Digital storage sweep generator |
CN108112056A (en) * | 2016-11-24 | 2018-06-01 | 联芯科技有限公司 | The frequency sweeping method and device of user terminal |
CN108107350A (en) * | 2018-01-05 | 2018-06-01 | 桂林电子科技大学 | A kind of long-range wide-band amplitude versus frequency characte measuring instrument |
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 |
CN109508511A (en) * | 2018-12-24 | 2019-03-22 | 中国航空工业集团公司西安飞机设计研究所 | Frequency sweeping method in frequency response analysis measurement |
CN110687213A (en) * | 2019-08-26 | 2020-01-14 | 天津大学 | Method for selecting optimal excitation frequency of nonlinear ultrasonic frequency mixing |
CN111024819A (en) * | 2019-12-23 | 2020-04-17 | 复旦大学 | Jewelry and jewelry acoustic resonance spectrum imitation method |
CN111261188A (en) * | 2020-01-20 | 2020-06-09 | 中国电力科学研究院有限公司 | High-voltage transformer noise spectrum determination method and device |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1993011647A1 (en) * | 1991-11-28 | 1993-06-10 | Kabushiki Kaisha Kenwood | Device for correcting frequency characteristic of sound field |
US6631783B2 (en) * | 2001-03-26 | 2003-10-14 | Nonlinear Seismic Imaging, Inc. | Mapping reservoir characteristics using earth's nonlinearity as a seismic attribute |
JP2009289373A (en) * | 2008-05-30 | 2009-12-10 | Hitachi Global Storage Technologies Netherlands Bv | Hard-disk device |
CN101711969B (en) * | 2009-10-30 | 2012-08-29 | 江苏大学 | Device and method for biological processing with pulsed multi-frequency ultrasonic waves in sweep-frequency mode and use thereof |
US9661433B2 (en) * | 2014-01-30 | 2017-05-23 | Invensense, Inc. | Electrical testing and feedthrough cancellation for an acoustic sensor |
CN109188453A (en) * | 2018-09-11 | 2019-01-11 | 哈尔滨工业大学 | Dynamic frequency scanning interfeerometry ranging system and distance measuring method based on locking phase gamma correction |
CN109975822B (en) * | 2019-04-16 | 2023-03-28 | 重庆大学 | Fiber laser dynamic ranging system integrating frequency sweep and single frequency |
-
2020
- 2020-07-30 CN CN202010748510.2A patent/CN111829650B/en active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101715200A (en) * | 2009-11-04 | 2010-05-26 | 武汉虹翼信息有限公司 | Method for analyzing situation of mobile network by drive test data of sweep signal generator |
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 | 苏州普源精电科技有限公司 | Sweep signal generator capable of generating a plurality of frequency markers |
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 | 南京信息工程大学 | Digital storage sweep generator |
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 |
CN109508511A (en) * | 2018-12-24 | 2019-03-22 | 中国航空工业集团公司西安飞机设计研究所 | Frequency sweeping method in frequency response analysis measurement |
CN110687213A (en) * | 2019-08-26 | 2020-01-14 | 天津大学 | Method for selecting optimal excitation frequency of nonlinear ultrasonic frequency mixing |
CN111024819A (en) * | 2019-12-23 | 2020-04-17 | 复旦大学 | Jewelry and jewelry acoustic resonance spectrum imitation method |
CN111261188A (en) * | 2020-01-20 | 2020-06-09 | 中国电力科学研究院有限公司 | High-voltage transformer noise spectrum determination method and device |
Non-Patent Citations (5)
Title |
---|
《A method for an accurate estimation of natural frequencies using swept-sine acoustic excitation 》;Ion Cornel Mituletu;《Mechanical Systems and Signal Processing》;20180724;全文 * |
《Analysis of sound and vibration interaction on a crack and its use in high-frequency parameter selection for vibro-acoustic modulation testing》;Bin Liu;《Mechanical Systems and Signal Processing》;20200329;全文 * |
《低频频率特性测试仪虚拟仪器设计》;薛媛元;《测控技术》;20090318;第31页-第34页 * |
《基于Matlab的低频频率特性测试虚拟仪器设计》;牛鹏宇;《测控技术》;20111118;全文 * |
《基于压电陶瓷的PVA-ECC混凝土柱裂缝损伤监测》;杨文伟;《沈阳建筑大学学报(自然科学版)》;20200715;第628页-第636页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111829650A (en) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5450049B2 (en) | Method and apparatus in an acoustic system | |
US5677987A (en) | Feedback detector and suppressor | |
Hudde | Measurement of the eardrum impedance of human ears | |
US4607216A (en) | Apparatus for measurement by digital spectrum analyzer | |
US20020154785A1 (en) | Adjusting a loudspeaker to its acoustic environment: the ABC system | |
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 | |
EP2814031B1 (en) | Method and apparatus for outputting audio signal, and method of adjusting volume of the audio signal | |
KR20120080593A (en) | An auditory test and compensation method | |
CN113116333B (en) | Method and system for detecting brain impedance distribution difference data | |
CN109684908B (en) | Signal filtering method for thromboelastography | |
JP2017090888A (en) | Method for modeling characteristic of instrument | |
CN111829650B (en) | Algorithm for testing frequency sweep parameters by using single-frequency signal combination | |
CN110557711B (en) | Earphone testing method and earphone | |
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 | |
CN105594231B (en) | Sound field measurement apparatus and sound field measuring method | |
CN107124678B (en) | Audio harmonic distortion measuring system | |
CN109391875B (en) | Audio debugging method and device, storage medium and processor | |
CN114598974A (en) | Bone conduction earphone equalization method based on distortion product otoacoustic emission | |
JP3272991B2 (en) | How to correct sound field characteristics | |
CN115580806B (en) | Headset noise reduction method based on automatic weight calculation of filter and noise reduction headset | |
CN115424636B (en) | Voice quality measurement method and device based on voice transmission index | |
US20210249028A1 (en) | Method and an audio processing unit for detecting a tone | |
JPH08294195A (en) | Signal processing method |
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 |