CN106949959B - A kind of method of impulse response signal sound intersity measurement - Google Patents
A kind of method of impulse response signal sound intersity measurement Download PDFInfo
- Publication number
- CN106949959B CN106949959B CN201710160108.0A CN201710160108A CN106949959B CN 106949959 B CN106949959 B CN 106949959B CN 201710160108 A CN201710160108 A CN 201710160108A CN 106949959 B CN106949959 B CN 106949959B
- Authority
- CN
- China
- Prior art keywords
- response signal
- impulse response
- signal
- inner product
- value
- 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
- 230000004044 response Effects 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005259 measurement Methods 0.000 title claims abstract description 16
- 239000000523 sample Substances 0.000 claims abstract description 21
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 238000001228 spectrum Methods 0.000 claims description 15
- 239000000700 radioactive tracer Substances 0.000 claims description 4
- 230000001052 transient effect Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 2
- 230000006698 induction Effects 0.000 claims description 2
- 210000001367 artery Anatomy 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
Abstract
The present invention provides a kind of method of impulse response signal sound intersity measurement, include the following steps: that, with two microphones probe A, B acquisition pulse response signal acoustic pressure data, the acoustic pressure data importing computer that two probes are measured is denoted as P' respectivelyA(t)、P'B(t)；Determine impulse response signal pAi(t) natural frequency valueWith impulse response signal pBi(t) natural frequency valueAccording to impulse response signal pAi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseImpulse response signal pBi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseParameter Reconstruction Impulse response signal, is denoted as respectivelyBy Reconstruction Impulse response signalSubstitute into instantaneous METHOD FOR SOUND POWER CALCULATION formula.The present invention can decompose impulse response signal without using Fourier transformation, but each impulse response signal is extracted from original signal, directly calculate its sound intensity；The sound intensity value for identifying impulse response signal can be clearly measured, the characteristic for accurately judging impulse response sound wave is facilitated.
Description
Technical field
The present invention relates in particular to a kind of methods of impulse response signal sound intersity measurement.
Background technique
The intensity of sound intensity characterization acoustic wave energy.Sound intersity measurement method is widely used in acoustic power level and calculates, vibrates identifing source
Equal fields.Engineer judges each vibration component actual contribution by the corresponding sound intensity value of Frequency point each on observation sound intensity frequency spectrum
Acoustic energy, and then take targetedly noise reduction measure.Noise Intensity Testing Technique is always problem concerned by people.
Currently, most common sound intersity measurement method is two microphones Mutual spectrum, Fourier transformation is that Mutual spectrum surveys the sound intensity
Core procedure, the time-domain sound pressure signal p that Fourier transformation measures two microphones A, channel BA(t)、pB(t) it is transformed into frequency domain
Sound pressure signal PA(ω)、PB(ω), substitution frequency domain METHOD FOR SOUND POWER CALCULATION formula finds out the sound intensity after calculating its cross-spectrum.For stable state just, cosine
Noise waves, using its sound intensity of Mutual spectrum survey calculation, the frequency spectrum of acquisition can effectively and reasonably reflect the feature of sound wave.
But in practical projects, sound wave existence form is numerous, and pulse shock is a kind of common exiting form, and impact produces
Raw impulse response.When calculating the sound intensity of impulse response wave using Mutual spectrum, impulse response signal is decomposed into more by Fourier transformation
A different frequency just, cosine signal, cause the sound intensity value on frequency spectrum to be dispersed on multiple and different Frequency points, in some instances it may even be possible to
Occurs the negative sound intensity on component frequency point, survey calculation result reflects that the impulse response noise sound intensity is not easy enough, and physical features are unknown
It is aobvious.
The present invention proposes that a kind of new sound intersity measurement calculates method.This method is pressed not using Fourier transformation as core
According to impulse response signal, itself inherent feature carrys out constructed fuction, and the items for calculating impulse response signal using inner product approach are special
Parameter is levied, according to the characteristic parameter Reconstruction Impulse response signal of acquisition, reconstruction signal is finally directly substituted into time domain sound intensity integral
Formula calculates, and solves the sound intensity value of impulse response sound wave.The method of the present invention effectively measuring can identify impulse response sound
The feature of wave, explicit physical meaning.
Summary of the invention
It in response to the deficiencies in the existing technology, can the present invention provides a kind of method of impulse response signal sound intersity measurement
To extract each impulse response signal from original signal, its sound intensity is directly calculated, clearly measurement identifies in this way
The sound intensity value of impulse response signal facilitates the characteristic for accurately judging impulse response sound wave.
The present invention achieves the above technical objects by the following technical means.
A kind of method of impulse response signal sound intersity measurement, includes the following steps:
S01: using two microphones probe A, B acquisition pulse response signal acoustic pressure data, and the acoustic pressure data measured is imported and is counted
Calculation machine, is denoted as P' respectivelyA(t)、PB'(t)；
S02: tracer signal P'A(t) and P'B(t) time span is T, to signal P'A(t) and P'B(t) direct computation of DFT is done
Leaf transformation determines that the impulse response signal quantity for including in channel is m by the peak value number of frequency spectrum, wherein i-th of arteries and veins of probe A
It rushes response signal and is set as pAi(t), i-th of impulse response signal of probe B is set as pBi(t), pA(t) and pB(t) it is denoted as respectively:
Wherein 0 < t < T；
S03: impulse response signal p is determinedAi(t) natural frequency valueWith impulse response signal pBi(t) intrinsic frequency
Value
Construct signal:
As initial value i=1, LA0(t)=pA(t)、LB0(t)=pB(t)；
Construct signal LAi-1(t) Fourier transformation is done, determines maximum amplitude and respective frequencies point, determines building signal LAi-1
(t) frequency estimationDetermine frequency search sectionIn frequency search sectionIt is interior to obtain the H subdivision thereof points divided equallyWherein, w=0,
1,…,H；Constructed fuctionBy constructed fuctionWith building signal LAi-1(t) inner product operation is done, inner product module is found
Extreme point obtains response signal p according to the frequency at extreme pointAi(t) natural frequency value
Construct signal LBi-1(t) Fourier transformation is done, maximum amplitude and respective frequencies point is determined, determines constructor LBi-1
(t) frequency estimationDetermine frequency search sectionIn frequency search sectionIt is interior to obtain the H subdivision thereof points divided equallyWherein, w=0,
1,…,H；Constructed fuctionBy constructed fuctionWith building signal LBi-1(t) inner product operation is done, inner product module is found
Extreme point obtains response signal p according to the frequency at extreme pointBi(t) natural frequency value
S04: according to impulse response signal pAi(t) natural frequency valueConstruct the function moved on a timeline:
According to impulse response signal pBi(t) natural frequency valueConstruct the function moved on a timeline:
Wherein: Δ T is the time span of function, and τ is variable and τ ∈ (0, T- Δ T)；
Construct signal LAi-1(t) with constructed fuction cAi(t) inner product operation, constructor LB arei-1(t) with constructed fuction cBi
(t) inner product operation is done, determines inner product modulus maximum, according to determining inner product modulus maximum, so that it is determined that impulse response signal pAi
(t) initial timeSignal pBi(t) initial time
S05: interception building signal LAi-1(t) in the periodOn sound pressure signal, be denoted asInterception building
Signal LBi-1(t) in the periodOn sound pressure signal, be denoted as
Constructed fuction:
Wherein, nAiAnd nBiFor variable；
Constructed fuction normalized is obtained:
yAi(t)=uAi(t)/||uAi(t)||
yBi(t)=uBi(t)/||uBi(t)||
By function yAi(t) and intercept signalInner product operation is carried out, by function yBi(t) and intercept signalInto
Row inner product operation, obtains XAi、XBiExpression formula:
Inner product module is found using thinning method step by step | | XAi| | and | | XBi| | maximum value, according to XAi、XBiExpression formula meter
Calculate attenuation coefficientWith
S06: according to inner product value Xmax Ai、Xmax BiWith constructed fuction uAi(t)、uBi(t) numerical value is calculated:
Wherein:
Xmax Ai、Xmax BiRespectively | | XAi||、||XBi| | obtain corresponding numerical value when maximum value；
According to C is calculatedAiAnd CBiPhase and amplitude:
Wherein:
Arg(CAi) indicate to seek CAiArgument, Arg (CBi) indicate to seek CBiArgument；
|CAi| plural C is sought in expressionAiAbsolute value, | CBi| plural C is sought in expressionBiAbsolute value；
S07: according to impulse response signal pAi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseImpulse response signal pBi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseParameter Reconstruction Impulse response signal, is denoted as respectivelyWherein:
S08: by Reconstruction Impulse response signalSubstitute into instantaneous METHOD FOR SOUND POWER CALCULATION formula:
Wherein, ρ: atmospheric density；
D: the spacing of probe A and B；
The time average for taking the signal transient sound intensity obtains the sound intensity of this impulse response signal:
S09: if when i < m, i=i+1, to LAi-1(t) and LBi-1(t) S03-S08 is executed, calculates remaining m-1
The sound intensity of impulse response signal；
As i=m, circulation terminates.
Further, in the S03 step
Beneficial effects of the present invention:
1. the method for impulse response signal sound intersity measurement of the present invention can decompose pulse without using Fourier transformation
Response signal, but each impulse response signal is extracted from original signal, directly calculate its sound intensity.
2. the method for impulse response signal sound intersity measurement of the present invention, can clearly measure and identify impulse response
The sound intensity value of signal facilitates the characteristic for accurately judging impulse response sound wave.
3. the method for impulse response signal sound intersity measurement of the present invention, can avoid the occurrence of the negative sound intensity.
Detailed description of the invention
Fig. 1 is the flow chart of the method for impulse response signal sound intersity measurement of the present invention.
Fig. 2 is the impulse response signal of embodiment of the present invention.
Fig. 3 is the sound intensity level spectrogram of the present invention obtained using Mutual spectrum measurement embodiment.
Specific embodiment
Present invention will be further explained with reference to the attached drawings and specific examples, but protection scope of the present invention is simultaneously
It is without being limited thereto.
As shown in Figure 1, a kind of method of impulse response signal sound intersity measurement, includes the following steps:
S01: using two microphones probe A, B acquisition pulse response signal acoustic pressure data, and the acoustic pressure data measured is imported and is counted
Calculation machine, is denoted as P' respectivelyA(t)、PB'(t)；
S02: tracer signal pA(t) and pB(t) time span is T, to signal pA(t) and pB(t) discrete fourier change is done
It changes, determines that the impulse response signal quantity for including in channel is m by the peak value number of frequency spectrum, wherein i-th of pulse of probe A is rung
Induction signal is set as pAi(t), i-th of impulse response signal of probe B is set as pBi(t), pA(t) and pB(t) it is denoted as respectively:
Wherein 0 < t < T；
S03: impulse response signal p is determinedAi(t) natural frequency valueWith impulse response signal pBi(t) intrinsic frequency
Value
Construct signal:
As initial value i=1, LA0(t)=pA(t)、LB0(t)=pB(t)；
Construct signal LAi-1(t) Fourier transformation is done, determines maximum amplitude and respective frequencies point, determines building signal LAi-1
(t) frequency estimationDetermine frequency search sectionIn frequency search sectionIt is interior to obtain the H subdivision thereof points divided equallyWherein, w=0,
1,…,H；Constructed fuctionBy constructed fuctionWith building signal LAi-1(t) inner product operation is done, inner product module is found
Extreme point obtains response signal p according to the frequency at extreme pointAi(t) natural frequency value
Construct signal LBi-1(t) Fourier transformation is done, maximum amplitude and respective frequencies point is determined, determines constructor LBi-1
(t) frequency estimationDetermine frequency search sectionIn frequency search sectionIt is interior to obtain the H subdivision thereof points divided equallyWherein, w=0,
1,…,H；Constructed fuctionBy constructed fuctionWith building signal LBi-1(t) inner product operation is done, inner product module is found
Extreme point obtains response signal p according to the frequency at extreme pointBi(t) natural frequency value
S04: according to impulse response signal pAi(t) natural frequency valueConstruct the function moved on a timeline:
According to impulse response signal pBi(t) natural frequency valueConstruct the function moved on a timeline:
Wherein: Δ T is the time span of function, and τ is variable and τ ∈ (0, T- Δ T)；
Construct signal LAi-1(t) with constructed fuction cAi(t) inner product operation, constructor LB arei-1(t) with constructed fuction cBi
(t) inner product operation is done, determines inner product modulus maximum, according to determining inner product modulus maximum, so that it is determined that impulse response signal pAi
(t) initial timeSignal pBi(t) initial time
S05: interception building signal LAi-1(t) in the periodOn sound pressure signal, be denoted asInterception building
Signal LBi-1(t) in the periodOn sound pressure signal, be denoted as
Constructed fuction:
Wherein, nAiAnd nBiFor variable；
Constructed fuction normalized is obtained:
yAi(t)=uAi(t)/||uAi(t)||
yBi(t)=uBi(t)/||uBi(t)||
By function yAi(t) and intercept signalInner product operation is carried out, by function yBi(t) and intercept signalInto
Row inner product operation, obtains XAi、XBiExpression formula:
Inner product module is found using thinning method step by step | | XAi| | and | | XBi| | maximum value, according to XAi、XBiExpression formula meter
Calculate attenuation coefficientWith
S06: according to inner product value Xmax Ai、Xmax BiWith constructed fuction uAi(t)、uBi(t) numerical value is calculated:
Wherein:
Xmax Ai、Xmax BiRespectively | | XAi||、||XBi| | obtain corresponding numerical value when maximum value；
According to C is calculatedAiAnd CBiPhase and amplitude:
Wherein:
Arg(CAi) indicate to seek CAiArgument, Arg (CBi) indicate to seek CBiArgument；
|CAi| plural C is sought in expressionAiAbsolute value, | CBi| plural C is sought in expressionBiAbsolute value；
S07: according to impulse response signal pAi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseImpulse response signal pBi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseParameter Reconstruction Impulse response signal, is denoted as respectivelyWherein:
S08: by Reconstruction Impulse response signalSubstitute into instantaneous METHOD FOR SOUND POWER CALCULATION formula:
Wherein, ρ: atmospheric density；
D: the spacing of probe A and B；
The time average for taking the signal transient sound intensity obtains the sound intensity of this impulse response signal:
S09: if when i < m, i=i+1；To LAi-1(t) and LBi-1(t) S03-S08 is executed, calculates remaining m-1
The sound intensity of impulse response signal；
As i=m, circulation terminates.
In the S03 stepWhen, calculated natural frequency value is more accurate.
Embodiment: as shown in Fig. 2,
Assuming that there are one by 3 impulse responses into the plane sound wave being grouped as:
Wherein, Ai、ni、τi、fi、The amplitude of respectively i-th (i=1,2,3) a impulse response signal, rises attenuation coefficient
Begin time, intrinsic frequency, initial phase.Design parameter value is shown in Table 1.Response signal 2,3 occurs simultaneously, and response signal 1 is below
Occur on time.
| Response signal | Intrinsic frequency (Hz) | Initial time (s) | Amplitude (pa) | Attenuation coefficient | Phase |
| Signal 1 | 450.21 | 1.3 | 20 | 4.2 | 0.2 |
| Signal 2 | 250.82 | 0.3 | 20 | 4.6 | 0.2 |
| Signal 3 | 100.15 | 0.3 | 20 | 5.3 | 0.5 |
Table 1: the characteristic parameter of each impulse response acoustic signals
S01: using two microphones probe A, B acquisition pulse response signal acoustic pressure data, and the acoustic pressure data measured is imported and is counted
Calculation machine, probe A are denoted as P' in rear end in front end, probe B respectivelyA(t)、P'B(t)；
Signal pB(t) than signal pA(t) a tiny time section is laggedWherein d is A, B
The spacing (unit: m) of probe, c are sound propagation velocity (unit: m/s) in air.
S02: tracer signal P'A(t) and P'B(t) time span is that T is 3s, does discrete Fourier transform to signal, is led to
The peak value number for crossing frequency spectrum determines that wherein i-th of impulse response signal of probe A is set as comprising 3 impulse response signals in channel
pAi(t), i-th of impulse response signal of probe B is set as pBi(t), pA(t) and pB(t) it is denoted as respectively:
Wherein 0 < t < 3；
S03: impulse response signal p is determinedAi(t) natural frequency valueWith impulse response signal pBi(t) intrinsic frequency
Value
Construct signal:
As initial value i=1, LA0(t)=pA(t)、LB0(t)=pB(t)；
Construct signal LAi-1(t) Fourier transformation is done, determines maximum amplitude and respective frequencies point, determines building signal LAi-1
(t) frequency estimationDetermine frequency search sectionWherein,It is searched in frequency
Rope sectionIt is interior to obtain the H subdivision thereof points divided equally, as H=100,ThenWherein, w=0,1 ..., 100；Constructed fuctionBy constructed fuctionWith building
Signal LAi-1(t) inner product operation is done, inner product module extreme point is found, according to the frequency at extreme point, obtains response signal pAi(t)
Natural frequency value
Construct signal LBi-1(t) Fourier transformation is done, maximum amplitude and respective frequencies point is determined, determines constructor LBi-1
(t) frequency estimationDetermine frequency search sectionWherein,It is searched in frequency
Rope sectionIt is interior to obtain 100 subdivision thereof points divided equallyIts
In, w=0,1 ..., 100；Constructed fuctionBy constructed fuctionWith building signal LBi-1(t) inner product operation is done,
Inner product module extreme point is found, according to the frequency at extreme point, obtains response signal pBi(t) natural frequency value
S04: according to impulse response signal pAi(t) natural frequency valueConstruct the function moved on a timeline:
According to impulse response signal pBi(t) natural frequency valueConstruct the function moved on a timeline:
Wherein: Δ T is the time span of function, value 0.02s.τ is variable and τ ∈ (0, T- Δ T)；
Construct signal LAi-1(t) with constructed fuction cAi(t) inner product operation, constructor LB arei-1(t) with constructed fuction cBi
(t) inner product operation is done, determines inner product modulus maximum, according to determining inner product modulus maximum, so that it is determined that impulse response signal pAi
(t) initial timeSignal pBi(t) initial time
S05: interception building signal LAi-1(t) in the periodOn sound pressure signal, be denoted asInterception building
Signal LBi-1(t) in the periodOn sound pressure signal, be denoted as
Constructed fuction:
Wherein, nAiAnd nBiFor variable；
Constructed fuction normalized is obtained:
yAi(t)=uAi(t)/||uAi(t)||
yBi(t)=uBi(t)/||uBi(t)||
By function yAi(t) and intercept signalInner product operation is carried out, by function yBi(t) and intercept signalInto
Row inner product operation, obtains XAi、XBiExpression formula:
Inner product module is found using thinning method step by step | | XAi| | and | | XBi| | maximum value, according to XAi、XBiExpression formula meter
Calculate attenuation coefficientWith
S06: according to inner product value Xmax Ai、Xmax BiWith constructed fuction uAi(t)、uBi(t) numerical value is calculated:
Wherein:
Xmax Ai、Xmax BiRespectively | | XAi||、||XBi| | obtain corresponding numerical value when maximum value；
According to C is calculatedAiAnd CBiPhase and amplitude:
Wherein:
Arg(CAi) indicate to seek CAiArgument, Arg (CBi) indicate to seek CBiArgument；
|CAi| plural C is sought in expressionAiAbsolute value, | CBi| plural C is sought in expressionBiAbsolute value；
S07: according to impulse response signal pAi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseImpulse response signal pBi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitial phaseParameter Reconstruction Impulse response signal, recognition result are shown in Table 2, are denoted as respectivelyWherein:
Reconstruction signal is as shown in Figure 2；
Table 2: each impulse response signal characteristic parameter that the method for the present invention is extracted
S08: by Reconstruction Impulse response signalSubstitute into instantaneous METHOD FOR SOUND POWER CALCULATION formula:
Wherein, ρ=1.293: atmospheric density (unit: Kg/m3)；D=0.01: the spacing (unit: m) of probe A and B；
The time average for taking the signal transient sound intensity obtains the sound intensity of this impulse response signal:
S09: if when i < m, i=i+1, to LAi-1(t) and LBi-1(t) step 3- step 8 is executed, remaining m-1 is calculated
The sound intensity of a impulse response signal, calculated result are shown in Table 3；
| Response signal | The theoretical sound intensity (W/m2) | Calculated result (W/m of the present invention2) | Relative error |
| Signal 1 | 0.01813109 | 0.01806788 | 0.3486% |
| Signal 2 | 0.01655058 | 0.01663622 | 0.5174% |
| Signal 3 | 0.01427526 | 0.01433991 | 0.4528% |
Table 3: the method for the present invention METHOD FOR SOUND POWER CALCULATION result
As i=m, circulation terminates.
Comparison calculates the sound intensity of impulse response signal using two microphones Mutual spectrum, and Mutual spectrum has used many description arteries and veins
The sound intensity of response signal is rushed, can not be indicated in the table.It is indicated using third-octave method.As shown in Figure 3, it can be seen that sound
It is dispersed in by force on different Frequency points, or even the negative sound intensity occurs.And the calculated sound intensity value result of this method is accurate, and it is each
Impulse response signal corresponds to each other, explicit physical meaning.
The embodiment is a preferred embodiment of the present invention, but present invention is not limited to the embodiments described above, not
In the case where substantive content of the invention, any conspicuous improvement that those skilled in the art can make, replacement
Or modification all belongs to the scope of protection of the present invention.
Claims (2)
1. a kind of method of impulse response signal sound intersity measurement, which comprises the steps of:
S01: using two microphones probe A, B acquisition pulse response signal acoustic pressure data, and the acoustic pressure data measured is imported computer,
It is denoted as P' respectivelyA(t)、P'B(t)；
S02: tracer signal P'A(t) and P'B(t) time span is T, to signal P'A(t) and P'B(t) discrete fourier change is done
It changes, determines that the impulse response signal quantity for including in channel is m by the peak value number of frequency spectrum, wherein i-th of pulse of probe A is rung
Induction signal is set as pAi(t), i-th of impulse response signal of probe B is set as pBi(t), pA(t) and pB(t) it is denoted as respectively:
Wherein 0 < t < T；
S03: impulse response signal p is determinedAi(t) natural frequency valueWith impulse response signal pBi(t) natural frequency valueConstruct signal:
As initial value i=1, LA0(t)=pA(t)、LB0(t)=pB(t)；
Construct signal LAi-1(t) Fourier transformation is done, determines maximum amplitude and respective frequencies point, determines building signal LAi-1(t)
Frequency estimationDetermine frequency search sectionIn frequency search sectionIt is interior to obtain the H subdivision thereof points divided equallyWherein, w=0,
1,…,H；Constructed fuctionBy constructed fuctionWith building signal LAi-1(t) inner product operation is done, inner product module is found
Extreme point obtains response signal p according to the frequency at extreme pointAi(t) natural frequency value
Construct signal LBi-1(t) Fourier transformation is done, maximum amplitude and respective frequencies point is determined, determines constructor LBi-1(t)
Frequency estimationDetermine frequency search sectionIn frequency search sectionIt is interior to obtain the H subdivision thereof points divided equallyWherein, w=0,
1,…,H；Constructed fuctionBy constructed fuctionWith building signal LBi-1(t) inner product operation is done, inner product module is found
Extreme point obtains response signal p according to the frequency at extreme pointBi(t) natural frequency value
S04: according to impulse response signal pAi(t) natural frequency valueConstruct the function moved on a timeline:
According to impulse response signal pBi(t) natural frequency valueConstruct the function moved on a timeline:
Wherein: Δ T is the time span of function, and τ is variable and τ ∈ (0, T- Δ T)；
Construct signal LAi-1(t) with constructed fuction cAi(t) inner product operation, constructor LB arei-1(t) with constructed fuction cBi(t)
Inner product operation is done, determines inner product modulus maximum, according to determining inner product modulus maximum, so that it is determined that impulse response signal pAi(t)
Initial timeSignal pBi(t) initial time
S05: interception building signal LAi-1(t) in the periodOn sound pressure signal, be denoted asInterception building signal
LBi-1(t) in the periodOn sound pressure signal, be denoted as
Constructed fuction:
Wherein, nAiAnd nBiFor variable；
Constructed fuction normalized is obtained:
yAi(t)=uAi(t)/||uAi(t)||
yBi(t)=uBi(t)/||uBi(t)||
By function yAi(t) and intercept signalInner product operation is carried out, by function yBi(t) and intercept signalCarry out inner product
Operation obtains XAi、XBiExpression formula:
Inner product module is found using thinning method step by step | | XAi| | and | | XBi| | maximum value, according to XAi、XBiExpression formula calculate
Attenuation coefficientWith
S06: according to inner product value XmaxAi、XmaxBiWith constructed fuction uAi(t)、uBi(t) numerical value is calculated:
Wherein:
XmaxAi、XmaxBiRespectively | | XAi||、||XBi| | obtain corresponding numerical value when maximum value；
According to C is calculatedAiAnd CBiPhase and amplitude:
Wherein:
Arg(CAi) indicate to seek CAiArgument, Arg (CBi) indicate to seek CBiArgument；
|CAi| plural C is sought in expressionAiAbsolute value, | CBi| plural C is sought in expressionBiAbsolute value；
S07: according to impulse response signal pAi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeJust
Beginning phaseImpulse response signal pBi(t) initial timeIntrinsic frequencyAttenuation coefficientAmplitudeInitially
PhaseParameter Reconstruction Impulse response signal, is denoted as respectivelyWherein:
S08: by Reconstruction Impulse response signalSubstitute into instantaneous METHOD FOR SOUND POWER CALCULATION formula:
Wherein, ρ: atmospheric density；
D: the spacing of probe A and B；
The time average for taking the signal transient sound intensity obtains the sound intensity of this impulse response signal:
S09: if when i < m, i=i+1, to LAi-1(t) and LBi-1(t) S03-S08 is executed, remaining m-1 pulse is calculated
The sound intensity of response signal；
As i=m, circulation terminates.
2. the method for impulse response signal sound intersity measurement according to claim 1, which is characterized in that in the S03 step
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710160108.0A CN106949959B (en) | 2017-03-17 | 2017-03-17 | A kind of method of impulse response signal sound intersity measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710160108.0A CN106949959B (en) | 2017-03-17 | 2017-03-17 | A kind of method of impulse response signal sound intersity measurement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106949959A CN106949959A (en) | 2017-07-14 |
| CN106949959B true CN106949959B (en) | 2019-10-01 |
Family
ID=59472160
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710160108.0A Active CN106949959B (en) | 2017-03-17 | 2017-03-17 | A kind of method of impulse response signal sound intersity measurement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106949959B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1250559A (en) * | 1997-03-12 | 2000-04-12 | 诺基亚电信公司 | Estimating CIR using correlation with line removal |
| CN102854006A (en) * | 2011-06-22 | 2013-01-02 | 霍尼韦尔国际公司 | Severity analysis apparatus and method for shafts of rotating machinery |
| CN103424183A (en) * | 2013-08-06 | 2013-12-04 | 常州大学 | Method for eliminating abnormal interference on detection for mechanical vibration signals |
| CN103605880A (en) * | 2013-10-25 | 2014-02-26 | 江苏大学 | Closely spaced mode damping ratio precisely-diagnosing method |
| CN103778291A (en) * | 2014-01-16 | 2014-05-07 | 江苏大学 | Variable-period integrated signal damping ratio diagnosis method |
| CN105424167A (en) * | 2015-11-09 | 2016-03-23 | 哈尔滨工程大学 | Underwater self-contained sound intensity instrument and sound intensity measurement method |
| CN106052852A (en) * | 2016-06-01 | 2016-10-26 | 中国电子科技集团公司第三研究所 | Pulse sound signal detection method and device |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5672317A (en) * | 1979-11-16 | 1981-06-16 | Nissan Motor Co Ltd | Detector for knocking oscillation |
| JPH09101199A (en) * | 1995-10-06 | 1997-04-15 | Mitsubishi Heavy Ind Ltd | Rotary vane vibration measurement device |
-
2017
- 2017-03-17 CN CN201710160108.0A patent/CN106949959B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1250559A (en) * | 1997-03-12 | 2000-04-12 | 诺基亚电信公司 | Estimating CIR using correlation with line removal |
| CN102854006A (en) * | 2011-06-22 | 2013-01-02 | 霍尼韦尔国际公司 | Severity analysis apparatus and method for shafts of rotating machinery |
| CN103424183A (en) * | 2013-08-06 | 2013-12-04 | 常州大学 | Method for eliminating abnormal interference on detection for mechanical vibration signals |
| CN103605880A (en) * | 2013-10-25 | 2014-02-26 | 江苏大学 | Closely spaced mode damping ratio precisely-diagnosing method |
| CN103778291A (en) * | 2014-01-16 | 2014-05-07 | 江苏大学 | Variable-period integrated signal damping ratio diagnosis method |
| CN105424167A (en) * | 2015-11-09 | 2016-03-23 | 哈尔滨工程大学 | Underwater self-contained sound intensity instrument and sound intensity measurement method |
| CN106052852A (en) * | 2016-06-01 | 2016-10-26 | 中国电子科技集团公司第三研究所 | Pulse sound signal detection method and device |
Non-Patent Citations (2)
| Title |
|---|
| "汽车振动测试分析系统的研究";陈伟;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20081130;C035-119 * |
| "脉冲激振模态分析法在300MW汽轮机松装叶片静频测量中的应用";王柏仁 等;《湖北电力技术》;19940331(第3期);第2-6页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106949959A (en) | 2017-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9909863B2 (en) | Determination of a room dimension estimate | |
| TWI525334B (en) | Sonic-assisted localization of wireless devices | |
| CN101600144B (en) | Method and system for obtaining a plurality of parameters of electro-acoustic product by adopting continuous logarithmic swept-frequency signal | |
| Sinha et al. | The impulse response of a high-speed jet forced with localized arc filament plasma actuators | |
| CN106643982B (en) | A kind of level measuring method based on resonance frequency of sound wave | |
| CN106255027B (en) | A kind of the sound quality Small Enclosure appraisal procedure and system of non-linear audio system | |
| CN101852638B (en) | Liquid level measurement method based on resonance frequency of sound wave on fixed frequency range | |
| CN101674520A (en) | Method and device for analyzing and testing parameter of vibration type pickup | |
| CN104270208B (en) | Method and device for detecting standing-wave ratio of RRU | |
| CN101426169A (en) | Time-domain tracking filter fast detecting acoustic response parameter of sounding body and system | |
| CN107064846A (en) | The sensitivity detection method and device of live testing apparatus for local discharge | |
| CN101874744B (en) | Ultrasonic guided wave parameter measurement method used for long bone analysis | |
| CN113472390B (en) | Frequency hopping signal parameter estimation method based on deep learning | |
| CN111811617A (en) | Liquid level prediction method based on short-time Fourier transform and convolutional neural network | |
| CN105675126B (en) | A kind of new method for stablizing sound field acoustic pressure for detecting multifrequency multi-source complexity | |
| CN104792450A (en) | Method for indirectly measuring exciting force of equipment acting on mounting base | |
| CN203100824U (en) | Level measurement device based on standing sound wave method and audio card | |
| CN105674065A (en) | Variational mode decomposition-based method for locating leakage point of pipeline by acoustic emission | |
| CN106949959B (en) | A kind of method of impulse response signal sound intersity measurement | |
| Willmarth | Space-time correlations and spectra of wall pressure in a turbulent boundary layer | |
| JP2004061361A (en) | Piping breakage investigating apparatus | |
| CN104568137B (en) | The decision method of specific noise source contribution rate in noise testing | |
| CN107741263B (en) | Liquid level estimation method based on frequency domain information fusion | |
| CN109696660B (en) | Method for accurately measuring amplitude sensitivity and phase of microphone for detecting free sound field | |
| CN201435826Y (en) | Vibrating sound pick-up parameter analyzing and testing device |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201230 Address after: No. 101, Nanxu Avenue, Zhenjiang City, Jiangsu Province, 212000 Patentee after: JIANGSU HUIZHI INTELLECTUAL PROPERTY SERVICES Co.,Ltd. Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301 Patentee before: JIANGSU University Effective date of registration: 20201230 Address after: 224200 No.8 Beihai Road, Dongtai City, Yancheng City, Jiangsu Province Patentee after: Dongtai science and Technology Service Center Address before: No. 101, Nanxu Avenue, Zhenjiang City, Jiangsu Province, 212000 Patentee before: JIANGSU HUIZHI INTELLECTUAL PROPERTY SERVICES Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210517 Address after: 224200 No. 48 Dongjin Avenue, Chengdong New District, Dongtai City, Yancheng City, Jiangsu Province Patentee after: Dongtai Chengdong science and Technology Pioneer Park Management Co.,Ltd. Address before: 224200 No.8 Beihai Road, Dongtai City, Yancheng City, Jiangsu Province Patentee before: Dongtai science and Technology Service Center |
|
| TR01 | Transfer of patent right | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170714 Assignee: Dongtai tepusong Machinery Equipment Co.,Ltd. Assignor: Dongtai Chengdong science and Technology Pioneer Park Management Co.,Ltd. Contract record no.: X2023980043158 Denomination of invention: A Method for Measuring Sound Intensity of Pulse Response Signals Granted publication date: 20191001 License type: Common License Record date: 20231012 |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170714 Assignee: Dongtai Donggao Electronic Information Technology Co.,Ltd. Assignor: Dongtai Chengdong science and Technology Pioneer Park Management Co.,Ltd. Contract record no.: X2023980045160 Denomination of invention: A Method for Measuring Sound Intensity of Pulse Response Signals Granted publication date: 20191001 License type: Common License Record date: 20231102 |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170714 Assignee: Dongtai Gaoxin Mechanical Equipment Co.,Ltd. Assignor: Dongtai Chengdong science and Technology Pioneer Park Management Co.,Ltd. Contract record no.: X2023980046304 Denomination of invention: A Method for Measuring Sound Intensity of Pulse Response Signals Granted publication date: 20191001 License type: Common License Record date: 20231110 |
|
| EE01 | Entry into force of recordation of patent licensing contract |