JPH1194642A - Diagnostic device for contribution of source of sound or vibration - Google Patents

Abstract

PROBLEM TO BE SOLVED: To evaluate a quantitative contribution rate in real time in the state where a residual remains. SOLUTION: A vibration sensor is mounted on each of in-vehicle principal panels serving as sources of sound or vibration, and microphones are installed in positions near the ears of a front-seat driver which serve as points of evaluation. An adaptive filter computing part 6 performs digital filtering process while changing a filter factor so that the difference between the sum of the digital filtering process values and a microphone output is minimized. A multiplier circuit part 8 computes the square of an input signal, and a contribution rate computing circuit part 9 adds together the outputs of the multiplier circuit part 8 in a predetermined time (i) (time (k)-(i) to time (k)) to compute a mean square, and computes the ratio Yaab /(Yk +Ek V) of the mean square to the sum of the total Yk V of the digital filtering process values and a residual Ek V as a contribution rate for display on a contribution rate display part 5.

Description

DETAILED DESCRIPTION OF THE INVENTION

[0001]

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a sound source / vibration source contribution diagnosis apparatus, and more particularly, to a transmission system in which one or more sound sources or vibration sources cause sound or vibration at an evaluation point. For example, in the fields of general mechanical devices, automobiles, building structures, OA equipment, railways, etc., visual and quantitative diagnosis of sound sources or vibration sources that greatly contribute to the evaluation points, and efficient improvement The present invention relates to a sound source / vibration source contribution diagnosis device that can be used for performing a diagnosis.

[0002]

2. Description of the Related Art Conventionally, as a method of analyzing the contribution of a sound source or a vibration source to an evaluation point, a correlation component is obtained from a transfer function between each sound source or a vibration source and the evaluation point in a frequency domain, and each correlation component is determined by the correlation component. A method for determining the degree of contribution of a sound source or a vibration source to an evaluation point has been proposed.

Further, recently, in the time domain, a method of diagnosing the contribution of a sound source or a vibration source using an adaptive digital filter (Japanese Patent Laid-Open No. 5-26722), and a method of calculating and displaying the contribution rate (Japanese Unexamined Patent Publication No. 7-243
No. 90) has also been proposed.

[0004]

However, in the above-described conventional frequency domain contribution analysis method, only amplitude information is used, and phase information is not taken into account. Therefore, it is difficult to analyze if there is a strong correlation between inputs. become. Further, the analysis is time-consuming, results are not immediately obtained, and long-term steady data is required.

In the conventional time domain contribution analysis method (Japanese Unexamined Patent Publication No. Hei 5-26722), it is possible to perform analysis in an actual operation state, that is, in real time, but as a result, only a time waveform is obtained. Since the degree of contribution is determined based on the magnitude of the amplitude of the time waveform, only a qualitative and approximate degree of contribution can be obtained in real time.

In order to solve this problem, a method of representing the instantaneous quantitative degree of contribution (contribution ratio) (JP-A-7-24)
390) is proposed, but in this method, a deviation between the sum of outputs from the respective digital filters and a target value,
That is, when the residual becomes minimum (0), the contribution rate is calculated based on the sum of the outputs from the digital filters. For example, when the sound source / vibration source is not specified, When detecting the position where it is located (replace the sensor with the position considered to be the sound source / vibration source, measure and calculate each time, and determine the position of the smallest of the calculated residuals of the sound source / vibration source ) And the number of sensors is smaller than the number of sound sources / vibration sources. There is no disclosure of a calculation method in a case where the residual is left with a size that cannot be ignored, and it is insufficient for practical use immediately.

[0007] The present invention has been made in view of such a point,
The aim is to provide a more practical sound source / vibration source contribution diagnosis device by enabling real-time and quantitative contribution rate evaluation even when residuals remain. And

[0008]

In order to achieve the above object, the present invention is directed to a sensor for detecting sound or vibration at a predetermined portion as an analog value, and an analog-to-digital converter for converting the analog value to a digital value. Converting means; a filter for performing digital filtering processing using the digital value and the filter coefficient; changing means for changing the filter coefficient so that a deviation between the filter output and a target value is minimized; In the state where the value has decreased to a value, the ratio of the average value of the filter output within a predetermined time to the sum of the filter output and the deviation, or the ratio of the root mean square value of the filter output within the predetermined time to the sum And a display means for displaying the contribution rate. A.

According to a second aspect of the present invention, there are provided a plurality of sensors for detecting sound or vibration at a plurality of predetermined parts as analog values, a plurality of analog-digital converting means for converting the analog values into digital values, respectively. A plurality of filters each performing digital filtering processing using the digital value and the filter coefficient, and changing means for changing a filter coefficient in at least one of the filters so that a deviation between a sum of the filter outputs and a target value is minimized. And in a state where the deviation is reduced to a predetermined value,
The ratio of the average value of the filter outputs within the predetermined time to the sum of the sum of the filter outputs and the deviation within a predetermined time period, or the ratio of the root mean square value of the filter outputs within the predetermined time period to the sum It is configured to include a contribution rate calculating means for calculating a contribution rate to an evaluation point and a display means for displaying the contribution rate.

Further, the invention of claim 3 is based on claims 1 and 2
The filter according to the present invention performs digital filtering processing using a plurality of digital values generated by delaying the digital value converted by the analog-digital conversion means and a plurality of filter coefficients, and changing means. However, at least one of the plurality of filter coefficients is changed so that the deviation is minimized.

According to the first and second aspects of the present invention, the sound at the sound source or the vibration at the vibration source is detected as an analog value, and the analog value is converted into a digital value. Digital filtering processing using a filter coefficient that is changed so that the deviation between the filtering processing result and the target value is minimized, and in a state where the deviation is reduced to a predetermined value, that is, when the deviation is close to the minimum value. In this state, the contribution ratio of the sound source or the vibration source to the evaluation point is calculated and displayed based on the digital filtering processing result. This makes it possible to diagnose the contribution of the sound source or the vibration source to the evaluation point when the digital filtering processing result (or the sum of the digital filtering processing results) converges to a predetermined value near the minimum value. The filtering process may be performed using a plurality of digital values generated by delaying the digital values and a plurality of changed filter coefficients of which at least one is changed so that the deviation is minimized. In the invention of the second aspect, the contribution ratio of the sound source or the vibration source to the evaluation point is calculated and displayed in a state where the residual remains, so that the degree of contribution can be determined at a glance even in the state where the residual remains. be able to.

More specifically, in the first aspect of the present invention,
The analog value of sound or vibration at a predetermined portion detected by the sensor is converted into a digital value by analog-digital conversion means for converting the analog value into a digital value. The filter performs digital filtering processing using the filter coefficient and the digital value converted by the analog-digital conversion means. The changing means changes the filter coefficient so that the deviation between the filter output and the target value is minimized, and the contribution rate calculating means evaluates using the filter output within a predetermined time in a state where the deviation is reduced to a predetermined value. Calculate the contribution to the point. This contribution is displayed on the display means.

[0013] As the contribution ratio, when the deviation is reduced to a predetermined value, that is, when the deviation remains to a non-negligible level, the average of the filter output within a predetermined time with respect to the sum of the filter output and the deviation is obtained. Value ratios can be used. In addition, the ratio of the root mean square value of the filter output within a predetermined time to this sum can be used as the contribution rate to the evaluation point.

As described above, according to the first aspect of the present invention, the contribution ratio of one sound source or vibration source including the residual is calculated and displayed. The contribution rate of the source to the evaluation point can be determined at a glance.

According to a second aspect of the present invention, a plurality of sensors detect sound or vibration at a plurality of predetermined portions as analog values, and convert the analog values into digital values by analog-digital conversion means. Each of the plurality of filters performs digital filtering processing using the digital value converted by the analog-digital conversion means and the filter coefficient, and the changing means performs at least the filtering so that the deviation between the sum of the filter outputs and the target value is minimized. The filter coefficient of one filter is changed, and the contribution ratio calculating means calculates the contribution ratio to the evaluation point with the deviation reduced to a predetermined value, as in the first aspect of the invention. This contribution is displayed on the display means.

As the contribution ratio, the ratio of the average value of the filter output within a predetermined time to the sum of the total sum and the deviation of the filter output, or the ratio of the mean value of the filter output within a predetermined time to this sum is used. Can be.

As described above, according to the second aspect of the present invention, the contribution ratios of a plurality of sound sources or vibration sources are calculated and displayed.
Even in the state where the residual remains, the contribution rate of a plurality of sound sources or vibration sources to the evaluation points can be determined at a glance.

Each of the filter according to the first aspect of the present invention and the filter according to the second aspect of the present invention includes a plurality of digital values generated by delaying the digital value converted by the analog-digital conversion means and a plurality of filters. A digital filtering process can be performed using the coefficients. In this case, in the invention of claim 1, the changing means changes at least one of the plurality of filter coefficients so as to minimize the deviation.
The changing means may change at least one of the plurality of filter coefficients in at least one filter.

[0019]

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of the present invention will be described below in detail with reference to the drawings. In the present embodiment, the present invention is applied to an analysis of the contribution of vehicle interior noise.

In general, when considering measures to reduce the noise in the cabin of a driver during driving, the driver or occupant of the driver or the occupant of the sound using the main panels A to G in the cabin shown in FIG. It is very important to determine the contribution to the sound. This is because effective measures can be implemented if measures are taken with vibration damping / sound absorbing materials or the like in the order of the parts (panels) having the largest contribution.

As shown in FIG. 1, the sound source / vibration source contribution diagnosing apparatus according to the present embodiment includes a plurality of vibration sensors 2a each composed of an acceleration sensor for detecting sound of a sound source or vibration of a vibration source by an analog value. To 2 g, and a microphone 3 for detecting sound at the evaluation point. These vibration sensors 2a to 2g are attached to each of the main panels A to G in the vehicle interior, which are sound sources or vibration sources, and the microphone 3 is installed at a position near the ear of the front driver who is the evaluation point.

Each of the vibration sensors 2a to 2g and the microphone 3 are connected to a contribution analyzer 4. The contribution analysis device 4 includes an adaptive filter operation unit 6, a personal computer for calculating the contribution ratio, and the like, and includes a time interval for calculating the contribution ratio, a convergence condition of the adaptive digital filter, a band pass described later. Filter frequency,
And a CRT or the like, to which a change means 7A such as a keyboard for setting analysis conditions such as delay time and the like, and display conditions such as the type, scale, and coordinates of a graph, which will be described later, are connected. And a contribution rate display section 5.

As shown in FIG. 2, the adaptive filter operation unit 6 includes low-pass filters (LPF) 20 to 27, which are anti-aliasing filters connected to the microphone 3 and each of the vibration sensors 2a to 2g,
A / D converters 30 to 37, which are connected to each of 0 to 27 and convert an analog value to a digital value at a predetermined cycle,
A delay circuit (DELA) connected to each of the D converters 30 to 37 and delaying an input signal by a predetermined time and outputting the delayed signal
Y) Adaptive digital filter (ADF) 4 connected to each of 50 to 57 and delay circuits 51 to 57 and performing digital filtering processing at the same cycle as the A / D conversion cycle.
1 to 47, connected to the output terminals of the adaptive digital filters 41 to 47, respectively.
Adder circuit 60 for calculating the sum of 47 outputs, delay circuit 50
Connected band-pass filter (BPF) 40 with and,, and an adder circuit 62 which calculates and outputs a deviation between the adding circuit 60 an output of bandpass filter 40 outputs a residual signal E _k. The output terminal of the adder circuit 60 is connected to each of the adaptive digital filters 41 to 47 and connected to the contribution ratio calculation circuit unit 7, and the output terminals of the adaptive digital filters 41 to 47 are also connected to the contribution ratio calculation circuit unit 7. It is connected to the.

The bandpass filter 40 can be constituted by an analog filter such as Chebyshev or Butterworth, or a digital filter.

Each of the above adaptive digital filters is
Although it is constituted by hardware such as a DPS (Digital Signal Processor), as shown schematically in FIG. 3, a digital value (vector quantity) X generated from an input digital value x _k and a filter coefficient W (Vector quantity)
And a block 64 composed of an FIR type, IIR type or lattice type digital filter for performing a digital filtering process by calculating a product sum with a LMS (least mean square) algorithm which is a general filter coefficient updating algorithm , A filter 66 that updates the filter coefficients so that the residual signal is minimized using the Newton method or the steepest descent method. The adaptive digital filter is preferably configured by hardware as described above in terms of speeding up the operation, but may be configured by software.

Block of the above adaptive digital filter
Numeral 64 denotes an input directory as schematically shown in FIG.
Digital value x_kDigital value x generated by delaying
_k-1 , X _k-2 , ... x_kLL delay circuits that output
D1, D2,..., DL
Updated filter coefficient w_0k, W_1k, ... w_LkAnd Digi
Tal value x_k-1 , X_k-2 , ... x_kLL to calculate the product of
+1 multiplication circuits S0, S1,..., SL
And a sum circuit SUM for calculating the sum of the products
I have.

The matrix X _k is _represented by [x _k , x _k−1 , x _k−2 ,.
, X _kL ], and the matrix W _k is _represented by [w _0k , w _1k , w _2k ,.
, W _Lk ], each of the adaptive digital filters outputs a signal represented by the following equation. ^{_{X k W k T = x k}} w 0k + x k-1 w 1k + x k-2 w 2k + ··· + x kL w Lk ··· (1) where, T represents an transpose, following X _k W _k Let ^T be simply X _k
Expressed as W _k .

The adaptive filter operation unit 6 can be represented schematically as shown in FIG. 4, and adds the digital filtering processing value according to the following equation (2),
, The sum Y _k (vector amount) at k time (current time) is output. In FIG. 4, x _ak , x _bk , x _{ck ,.}
x _gk is the output (scalar amount) of the delay circuits 51 to 57.

Y _k = X _ak W _ak + X _bk W _bk + X _ck W _ck +... + X _gk W _gk (2) where X _ak , X _bk , X _{ck ,.} type digital to correspond to the filter 41 to 47 is the aforementioned matrix X _k bearing the suffix a to g, the adaptive digital filter 4
The digital values x _ak input at time k to 1 to ₄₇ ,
x _bk, x _ck, shows the digital value generated from _{··· x gk, W ak, W} bk, W ck, a · · · W _gk subscript a~g corresponds to adaptive digital filters 41 to 47 The above-mentioned matrix W _k attached to each of the adaptive digital filters 41 to 41
47 shows the filter coefficient at k time 47.

The band at time k is output from the adding circuit 62.
Pass filter 40 output d_kAnd sum Y _kThe residual that is the difference from
Signal E_k= D_k-Y_kIs output. And blocks
At 66, the residual signal E is obtained using the LMS algorithm or the like.
_kIs gradually reduced and converges to the minimum value.
The data coefficient is updated to be the filter coefficient at the time k + 1. This
, The residual signal E_kMay converge to the minimum value of 0
However, it usually converges to a predetermined value near the minimum value. Also,
At least one filter coefficient must be updated according to the conditions.
However, normally, all the filter coefficients are updated.

The contribution ratio calculating circuit 7, the output of the band-pass filter 40, the adaptive digital filter 41 to 47 each output, and a plurality of contributions connected to the output of the adder circuit 62 for outputting a residual signal E _k It is constituted by rate calculating means. As shown in FIG. 5, each of the contribution rate calculation means adds the output of the multiplication circuit section 8 for calculating the square of the input signal and the output of the multiplication circuit section 8 within a predetermined time i (time k-i to time k). And a contribution ratio calculation circuit unit 9 for calculating a contribution ratio using the root mean square and the like, and performing coordinate conversion for performing display in accordance with the type, coordinates and scale of the specified graph. It is composed of a coordinate conversion circuit section 10 for performing the operation.

Assuming that Ya _k = X _ak W _ak , the mean square value Ya _ab of the panel A within the predetermined time i is equal to the predetermined time i
Is expressed as shown in the following equation (3) using the digital filtering processing value in the above. ^{_{Ya ab = (Ya k 2 +}} Ya k-1 2 + Ya k-2 2 + ··· + Ya ki 2) / (i + 1) ··· (3) _{Similarly, Yb k = X bk W bk} , Yc k = X _ck W _ck・ ・ ・ Y
If d _k = X _gk W _gk , the mean square values Yb _ab , Yc _ab ,..., Yg _ab for panels B to G are expressed as follows, and the mean square value for the residual signal E _{ab is} also expressed as follows. ^{_{Yb ab = (Yb k 2 +}} Yb k-1 2 + Yb k-2 2 + ··· + Yb ki 2) / (i + 1) Yc ab = (Yc k 2 + Yc k-1 2 + Yc k-2 2 + ··· ^{_{+ Yc ki 2) / (i}} + 1) ··· Yg ab = (Yg k 2 + Yg k-1 2 + Yg k-2 2 + ··· + Yg ki 2) / (i + 1) E ab = (E k 2 + E k- ^{_{1 2 + E k-2 2}} + ··· + E ki 2) / (i + 1) ··· (4) In addition, k time output from the adding circuit 60 (the current time)
The sum Y _k and for the residual signal E _k, 2 mean square value Yb _ab for panel _{B~G, Yc ab, ···, Yg} ab, E
_If the ratio of _ab is the contribution ratio, the contribution ratio is _Yaab / ( _Yk +
E _k ), Yb _ab / (Y _k + E _k ),... Yg _ab / (Y
_k + _Ek ) and _Eab / ( _Yk + _Ek ).

Next, the operation of the present embodiment will be described. Detection signals, which are analog signals detected by the vibration sensors 2a to 2g, are input to low-pass filters 21 to 27, respectively, and high-frequency components are removed. Output signals of the low-pass filters 21 to 27 are input to A / D converters 31 to 37, respectively, A / D-converted at predetermined cycles, and output as digital signals. These digital signals are input to delay circuits 51 to 57 and delayed, and then input to adaptive digital filters 41 to 47. As described above, digital filtering is performed at the same cycle as the A / D conversion cycle. The output signals Ya _{k to} Yg _k are input to the arithmetic circuit 60, respectively.

A detection signal (target signal) from the microphone 3 is input to a low-pass filter 20, which is an anti-aliasing filter, and the output signal is input to an A / D converter 30, and the digital signal is output at a predetermined period. Is converted and output. This digital signal is supplied to the delay circuit 50
After that, the signal is input to a band-pass filter 40 for limiting an analysis frequency band, and the output signal d _k is input to an adding circuit 62.

Here, by using the above-mentioned delay circuits 51 to 57, it is possible to correct the delay caused by the signal passing through the band pass filter 40 for limiting the analysis frequency band. That is, it is possible to secure the synchronization on the time axis between the signal d _k and the signals Ya _{k to} Yg _k, and it is also possible to search for a sound source or a vibration source. In order to secure the synchronism of the signals on the time axis, the delay circuits 51 to 5
7, but usually, the delay circuits 51 to 51 are used.
In many cases, contribution analysis can be performed without using 57. The delay circuit 50 is used for performance check using a simulation or a model signal, operation check of a circuit, and the like, and usually need not be used.

The output from each block of each of the adaptive digital filters 41 to 47 is input to a contribution rate operation circuit 7, and the contribution rate operation circuit 7 outputs the output signals Ya _{k to} Yg _{k of} each adaptive digital filter. , The sum Y _k , and the residual signal E _k to calculate the contribution ratio as described above.

The contribution ratio is displayed on the contribution ratio display section 5 in real time as a graph as shown in FIGS. 6, 8 and 9 and a bar graph as shown in FIG. In addition,
In each figure, ch1 to ch7 indicate the values of panels A to G, err indicates the magnitude of the residual signal, and desire
Indicates a target signal.

The filter coefficient is changed by the block 66 so that the residual signal _Ek gradually decreases and converges to a predetermined value near the minimum value. Is calculated and displayed including the residual.

By changing the average time interval by changing the predetermined time i in the contribution rate calculator 7, the display time interval of the bar graph display can be adjusted, and the average time interval is shortened for a signal that changes rapidly. By specifying a longer average time interval for a specified and stationary signal, the result can be easily checked. The vertical axis is 50%, 100%
By making the full scale such as% variable, the result can be more easily confirmed.

In the above description, the contribution ratio is calculated using the mean square value. However, the contribution ratio may be calculated as follows using the average value of Ya _{k to} Yg _k . _{(Ya k + Ya k-1} + Ya k-2 + ··· + Ya ki) / {(i + 1) · (Y k + E k)} (Yb k + Yb k-1 + Yb k-2 + ··· + Yb ki) / {(i + 1) · (Y k + E k)} ··· (Yg k + Yg k-1 + Yg k-2 + ··· + Yg ki) / {(i + 1) · (Y k + E k)} Figure Reference numeral 11 denotes a result of a contribution analysis at a constant speed (50 km / h) by a conventional time domain contribution analysis method and apparatus. From this figure, it is possible to judge a panel that has a qualitatively large contribution, but it is not possible to judge the degree of a quantitative contribution. The same applies to the result of the contribution analysis at the time of accelerating traveling. Although the chronological change of the contribution can be determined qualitatively, the degree of the quantitative contribution cannot be determined.

FIG. 7 shows an output result at the time of constant speed running according to the above embodiment. From the bar graph displayed here, the predetermined engine speed (562.70R)
PM, this value can be changed), and the contribution rate of each panel vibration to the vehicle interior noise can be quantitatively and clearly understood.

FIGS. 6 and 8 show output results obtained by the contribution analyzing apparatus of the above embodiment. In the graph of FIG. 6, the horizontal axis represents time, and the vertical axis represents the contribution rate. In the graph of FIG. 8, the horizontal axis represents the engine speed and the vertical axis represents the contribution rate. From these, the time (or rotation speed) dependency of the contribution rate of each panel can be quantitatively determined. From this graph, a bar graph display similar to that shown in FIG. 7 at an arbitrary time (or rotation speed) is also possible. FIG. 9 is a graph when the scale is changed.

In the above description, a case where a large number of signals are input to the adaptive filter operation unit has been described as an example. However, the present invention can also be applied to a single input. In FIG.
A block diagram of a sound source / vibration source contribution diagnosis apparatus when a single signal is input using only the vibration sensor 2a as a sensor is shown.

In this case, the arithmetic circuit 60 becomes unnecessary, and at least one of the plurality of filter coefficients w _0k , w _1k , w _2k ,..., W _Lk provided in the adaptive digital filter is used. Update it.

The contribution ratio in this case is represented by the mean square value Y
When a _ab is used, it is expressed as Ya _ab / (Y _ak + E _k ), and when an average value is used, it is expressed as follows.

(Ya_k+ Ya_k-1 + Ya_k-2 + ... +
Ya_ki) / ｛(I + 1) · (Ya _k+ E_kOthers are the same as above, and the description is omitted.

In this embodiment, especially when the input signals from the vibration sensors 2a to 2g contain many noise components, the same characteristics as the bandpass filter 40 for limiting the frequency band to be analyzed are used. The bandpass filter of
A stage preceding each adaptive digital filter, that is, between the adaptive digital filters 41-47 and the delay circuits 51-57, or between the A / D converters 31-37 and the delay circuits 51-51.
In some cases, a more accurate result can be obtained by inserting the line between the line 57 and the line 57.

When the frequency band to be analyzed changes with time in synchronization with the rotation speed, analyzing the contribution using a tracking filter synchronized with the rotation speed is also an effective analysis method.

The sound of this embodiment is usually detected by a microphone and a signal amplified by a sound level meter amplifier or the like is used. Vibration is usually detected by a displacement sensor, a speed sensor, an acceleration sensor or the like, and amplified by the amplifier. Use the obtained signal. However, it is also possible to substitute an electric signal which is completely correlated with sound or vibration.

In the present embodiment, the frequency characteristic of the contribution ratio can be obtained by changing the filtering output signal to the frequency domain by FFT (Fast Fourier Transform) or the like. In addition, the filtering output signal is averaged within a predetermined time, and by using the obtained frequency spectrum or power spectrum and dividing by the frequency spectrum or power spectrum of the sum of a plurality of filter outputs, the contribution of each input in the frequency domain is obtained. Rate is required.

Further, by changing each digital filter coefficient to a frequency domain by FFT or the like, a transfer characteristic between each input signal (sound or vibration) and the evaluation point signal (sound or vibration) is obtained. Effective information can be obtained for soundproofing and anti-vibration measures or noise reduction design.

In the above description, an example in which sound or vibration is detected as an analog value and then converted into a digital value has been described.
A sensor that directly detects digital values may be used.

As described above, this embodiment diagnoses the contribution of sound or vibration to the evaluation point, including the phase information, by using the adaptive digital filter. Since the data that changes quickly, etc., can be easily handled and the contribution ratio is determined in real time and quantitatively, the utility is extremely high.
Further, the degree of contribution can be determined in real time and quantitatively in comparison with the conventional masking method, frequency domain contribution analysis method, and time domain contribution analysis method, so that the method and apparatus are very effective.

Therefore, soundproofing and antivibration measures can be performed efficiently and accurately, and the contribution rate is determined quantitatively, so that the effect of the measures can be quantitatively predicted to some extent.

Further, in the active control of sound and vibration, the present invention can be applied to a case where a sound source and a vibration source having high causality are selected and a case where the effect of the active control is quantitatively predicted.

In addition to the above-described embodiments, the present invention provides a noise and vibration problem for general equipment, a sound and vibration problem for cars, a source and vibration source identification for railway noise and vibration, a noise and vibration problem for OA equipment, a building structure. It can be applied to various problems such as noise and vibration of objects.

[0057]

As described above, according to the first, second and third aspects of the present invention, in a complicated acoustic system or vibration system,
In the actual operation state, the effect that the contribution rate of the sound source or the vibration source to the evaluation point can be determined at a glance in a real-time and quantitative manner while the residual remains is obtained.

[Brief description of the drawings]

FIG. 1 is a block diagram illustrating an apparatus for analyzing contribution of in-vehicle noise according to an embodiment of the present invention.

FIG. 2 is a block diagram illustrating details of an adaptive filter operation unit in FIG. 1;

FIG. 3 is a block diagram schematically showing the adaptive digital filter of FIG. 1;

FIG. 4 is a schematic diagram of an adaptive filter operation unit in FIG. 1;

FIG. 5 is a block diagram of a contribution rate calculation circuit unit of FIG. 1;

FIG. 6 is a graph of an output result with the horizontal axis representing time according to the present embodiment;

FIG. 7 is a bar graph showing an output result according to the embodiment.

FIG. 8 is a graph of an output result with the horizontal axis representing the engine speed according to the present embodiment.

FIG. 9 is a graph of an output result with the horizontal axis representing the engine speed according to the present embodiment;

FIG. 10 is a schematic diagram showing details of a block 64 in FIG. 3;

FIG. 11 is a graph showing a conventional output result.

FIG. 12 is a block diagram illustrating a device for analyzing the contribution of in-vehicle noise in the case of a single input.

[Explanation of symbols]

 2a to 2g Vibration sensor 3 Microphone 4 Contribution analysis display device 5 Contribution ratio display unit 6 Adaptive filter operation unit 7 Contribution ratio operation circuit unit 20 to 27 Low pass filter 30 to 37 A / D converter 41 to 47 Adaptive digital filter 40 band Pass filter 50-57 Delay circuit

Claims (3)

[Claims] 1. A sensor for detecting sound or vibration at a predetermined portion as an analog value, an analog-to-digital converter for converting the analog value to a digital value, and a digital filtering process using the digital value and a filter coefficient. A filter, changing means for changing the filter coefficient so that a deviation between the filter output and a target value is minimized, and a filter for the sum of the filter output and the deviation in a state where the deviation is reduced to a predetermined value. Contribution ratio calculating means for calculating, as a contribution ratio to an evaluation point, a ratio of an average value of the filter outputs within a predetermined time period or a ratio of a root mean square value of the filter outputs to the sum within the predetermined time period; A display means for displaying a sound source, and a contribution diagnosis apparatus for a sound source / vibration source, comprising: 2. A plurality of sensors for detecting sound or vibration at a plurality of predetermined parts as analog values, a plurality of analog-digital converting means for converting the analog values into digital values, respectively, the digital values and filter coefficients A plurality of filters for performing digital filtering processing using: a changing means for changing a filter coefficient of at least one of the filters so that a deviation between a sum of the filter outputs and a target value is minimized; In a state in which the sum of the filter output and the deviation within a predetermined period is the ratio of the average value of the filter output within the predetermined time to the sum of the deviation, or the ratio of the filter output within the predetermined time to the sum Contribution ratio calculating means for calculating a ratio of a root mean square as a contribution ratio to an evaluation point; Contribution diagnosis apparatus of sound and vibration sources, including a Shimesuru display means. 3. The filter performs a digital filtering process using a plurality of digital values generated by delaying the digital value converted by the analog-to-digital converter and a plurality of filter coefficients. 3. The apparatus according to claim 1, wherein at least one of the plurality of filter coefficients is changed such that the deviation is minimized.