JPH0915035A - Adaptive signal processing method and equipment therefor, detection method for noise or vibration contribution, transfer characteristics-estimation method for noise or vibration, suppression method to noise or vibration, and simulation method for noise or vibration environments - Google Patents

Adaptive signal processing method and equipment therefor, detection method for noise or vibration contribution, transfer characteristics-estimation method for noise or vibration, suppression method to noise or vibration, and simulation method for noise or vibration environments

Info

Publication number
JPH0915035A
JPH0915035A JP7160259A JP16025995A JPH0915035A JP H0915035 A JPH0915035 A JP H0915035A JP 7160259 A JP7160259 A JP 7160259A JP 16025995 A JP16025995 A JP 16025995A JP H0915035 A JPH0915035 A JP H0915035A
Authority
JP
Japan
Prior art keywords
noise
vibration
signal
adaptive
stage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP7160259A
Other languages
Japanese (ja)
Other versions
JP3393738B2 (en
Inventor
Masanori Watabe
眞徳 渡部
Yasushi Takano
靖 高野
Shinichi Shimoide
新一 下出
Katsuo Oki
克夫 大木
Gichu Ota
義注 太田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP16025995A priority Critical patent/JP3393738B2/en
Publication of JPH0915035A publication Critical patent/JPH0915035A/en
Application granted granted Critical
Publication of JP3393738B2 publication Critical patent/JP3393738B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Abstract

PURPOSE: To provide an adaptive signal processing method which can precisely analyze the impulse response of an acoustic system or a vibration system, and an equipment using the method. CONSTITUTION: The respective signals of noise or vibration which are detected at a plurality of points are inputted in first stage adaptive filters 431-433. The coefficients of the adaptive filters 431-433 are updated so that the difference e1 between the sum of the filter output signals and the noise signal detected at an estimation point where noise or vibration is estimated becomes minimum. Detection signals at a plurality of points are inputted in second stage adaptive filters 434-436, which are so updated that the difference e2 between the sum of the output signals and the difference e1 becomes minimum. The respective adaptive filter inputs are delayed through delay circuits 421-426. Impulse responses of various kinds of acoustic fields can be expressed by installing two stages of adaptive filters and combining the adaptive filters with delay time.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、1個または複数の信号
から評価点における信号を合成するように適応ディジタ
ルフィルタの特性を定め、これにより音響系または振動
系のインパルス応答を解析するようにした適応信号処理
方法及びその装置と、その方法を用いた騒音または振動
寄与率の検出方法、騒音または振動の伝達特性推定方
法、騒音または振動の抑圧方法、及び騒音または振動環
境の模擬方法に関するものである。
BACKGROUND OF THE INVENTION The present invention defines the characteristics of an adaptive digital filter so as to synthesize a signal at an evaluation point from one or a plurality of signals, and thereby analyzes the impulse response of an acoustic system or a vibration system. Adaptive signal processing method and device, and method for detecting noise or vibration contribution rate using the method, noise or vibration transfer characteristic estimation method, noise or vibration suppression method, and noise or vibration environment simulation method Is.

【0002】[0002]

【従来の技術】従来より、音波の空間伝達特性の推定、
機器等からの騒音振動の測定やその低減のために、適応
信号処理を用いた方法が数多く提案されている。そのう
ち、音源または振動源の影響の度合を調べる寄与診断方
法等の従来例として、特開平5−26722号が挙げら
れる。これは 機器等の周辺の複数点で検出される騒
音、もしくは振動の各信号をそれぞれ適応フィルタに入
力し、その適応フィルタの出力信号の和と、騒音もしく
は振動を評価する評価点で検出される騒音もしくは振動
の信号との残差が最小になるように、前記適応フィルタ
の係数を更新し、一定値に収束した時の各適応フィルタ
の出力信号を用いて、評価点における騒音もしくは振動
に対して音源もしくは振動源が寄与している度合を解析
するものである。
2. Description of the Related Art Conventionally, estimation of spatial transfer characteristics of sound waves,
A number of methods using adaptive signal processing have been proposed for measuring and reducing noise and vibration from devices and the like. Among them, Japanese Patent Laid-Open No. 5-26722 is known as a conventional example of a contribution diagnosis method for investigating the degree of influence of a sound source or a vibration source. This is detected by inputting each noise or vibration signal detected at multiple points around the device, etc. to the adaptive filter, and the sum of the output signal of the adaptive filter and the evaluation point for evaluating noise or vibration. The coefficient of the adaptive filter is updated so that the residual difference with the noise or vibration signal is minimized, and the output signal of each adaptive filter when it converges to a constant value is used for noise or vibration at the evaluation point. The degree of contribution of the sound source or the vibration source is analyzed.

【0003】また、音波の空間伝達特性の推定では、信
号発生源からの信号を適応フィルタに入力し、前記適応
フィルタの出力信号と、前記信号発生源からの信号を音
波として空間に放射し、設置した測定点で検出される騒
音情報の信号との差の残差信号が最小になるように、適
応フィルタの係数を更新し、そのフィルタ係数を収束さ
せ、収束した適応フィルタ係数を用いて測定点までの空
間の伝達特性を推定していた。
Further, in estimating the spatial transfer characteristic of a sound wave, a signal from a signal generation source is input to an adaptive filter, and an output signal of the adaptive filter and a signal from the signal generation source are radiated into space as a sound wave, The coefficient of the adaptive filter is updated so that the residual signal of the difference from the noise information signal detected at the installed measurement point is minimized, the filter coefficient is converged, and measurement is performed using the converged adaptive filter coefficient. The transfer characteristics of the space to the point were estimated.

【0004】[0004]

【発明が解決しようとする課題】上記した従来の適応信
号処理方法では、高残響場等の残響時間の長い環境での
場合、インパルス応答が精度よく推定できないため、適
応フィルタの出力信号の和と、騒音もしくは振動の評価
点で検出される騒音もしくは振動の信号との残差が適応
し切れずに、多く残ってしまうという問題がある。さら
に、直接波と反射波の時間間隔が長い場合には、反射波
まで表現しようとすると、適応フィルタのタップ数を多
くとる必要があった。
In the conventional adaptive signal processing method described above, the impulse response cannot be accurately estimated in an environment with a long reverberation time such as a high reverberation field. However, there is a problem that the residual difference from the noise or vibration signal detected at the noise or vibration evaluation point is not fully applied and a large amount remains. Furthermore, when the time interval between the direct wave and the reflected wave is long, it is necessary to increase the number of taps of the adaptive filter to represent the reflected wave.

【0005】本発明の目的は、高残響場等での騒音・振
動の測定や音源・振動源の寄与診断等を正確に行え、構
成の簡単な適応信号処理方法及びその装置と、その方法
を用いた騒音または振動寄与率の検出方法、騒音または
振動の伝達特性推定方法、騒音または振動の抑圧方法、
及び騒音または振動環境の模擬方法を提供することにあ
る。
It is an object of the present invention to provide an adaptive signal processing method and apparatus and a method therefor, which can accurately measure noise and vibration in a highly reverberant field and diagnose contribution of sound source and vibration source. Noise or vibration contribution detection method used, noise or vibration transfer characteristic estimation method, noise or vibration suppression method,
And to provide a method of simulating a noise or vibration environment.

【0006】[0006]

【課題を解決するための手段】本発明は、第1の1個又
は複数個の時系列信号を1個又は複数個の1段目適応フ
ィルタの各々に入力し、該適応フィルタの出力信号の和
と前記第1の時系列信号とは別の第2の時系列信号との
差を表す第1の残差信号が最小になるように、前記1段
目適応フィルタの係数を更新するとともに、前記第1の
時系列信号を1個又は複数個の2段目適応フィルタの各
々に入力し、該適応フィルタの出力信号の和と前記第1
の残差信号との差を表す第2の残差信号が最小になるよ
うに前記2段目適応フィルタを更新することを特徴とす
る適応信号処理方法を開示する。
According to the present invention, a first one or a plurality of time series signals is input to one or a plurality of first stage adaptive filters, and the output signals of the adaptive filters are output. The coefficient of the first-stage adaptive filter is updated so that the first residual signal representing the difference between the sum and the second time-series signal different from the first time-series signal is minimized, The first time series signal is input to each of one or a plurality of second stage adaptive filters, and the sum of the output signals of the adaptive filters and the first
Disclosed is an adaptive signal processing method, characterized in that the second-stage adaptive filter is updated so that the second residual signal representing the difference from the residual signal is updated.

【0007】また本発明は、前記1段目及び2段目の適
応フィルタに入力される前記第1の時系列信号を、その
遅延時間が可変な遅延手段により遅延できるようにした
ことを特徴とする適応信号処理方法を開示する。
Further, the present invention is characterized in that the first time-series signal input to the first-stage and second-stage adaptive filters can be delayed by delay means having a variable delay time. A method for adaptive signal processing is disclosed.

【0008】また本発明は、上記した適応信号処理方法
を用いて、前記第1の時系列信号を、騒音又は振動発生
源周辺に設けた複数の測定点に於て検出した複数の騒音
又は振動の検出信号とし、且つ前記第2の時系列信号を
前記騒音又は振動発生源からの騒音又は振動を評価する
ために設けた評価点に於て検出した騒音又は振動の検出
信号としたときの適応信号処理を行い、該処理終了時の
前記適応フィルタの各出力から、当該出力に対応する前
記測定点に於る検出信号の、前記評価点に於る検出信号
に対する寄与率を求めることを特徴とする騒音又は振動
の寄与率の検出方法を開示する。
According to the present invention, a plurality of noises or vibrations obtained by detecting the first time-series signal at a plurality of measurement points provided around a noise or vibration source using the above adaptive signal processing method. And the second time-series signal is a noise or vibration detection signal detected at an evaluation point provided for evaluating noise or vibration from the noise or vibration source. Signal processing is performed, and from each output of the adaptive filter at the end of the processing, a contribution rate of a detection signal at the measurement point corresponding to the output to a detection signal at the evaluation point is obtained. A method of detecting a contribution rate of noise or vibration is disclosed.

【0009】また本発明は、上記した適応信号処理方法
を用いて、前記第1の時系列信号を、騒音又は振動発生
手段に印加する騒音又は振動信号とし、且つ前記第2の
時系列信号を、前記騒音又は振動発生手段から発生させ
た騒音又は振動を評価する評価点に於て検出された検出
信号として適応信号処理を行い、該処理終了時の前記各
適応ディジタルフィルタの係数から、前記騒音又は振動
発生手段から前記評価点に到る騒音又は振動の伝達特性
を推定することを特徴とする騒音又は振動の伝達特性推
定方法を開示する。
According to the present invention, by using the above adaptive signal processing method, the first time series signal is a noise or vibration signal to be applied to a noise or vibration generating means, and the second time series signal is , Adaptive signal processing is performed as a detection signal detected at an evaluation point for evaluating noise or vibration generated from the noise or vibration generating means, and the noise is calculated from the coefficient of each adaptive digital filter at the end of the processing. Alternatively, a noise or vibration transfer characteristic estimation method is disclosed, in which the noise or vibration transfer characteristic reaching the evaluation point is estimated from the vibration generating means.

【0010】また本発明は、補償用の騒音又は振動発生
手段から騒音又は振動を発生することにより、騒音又は
振動発生源からの騒音又は振動を評価点に於て抑圧する
ための騒音又は振動の抑圧方法に於て、上記した騒音又
は振動の伝達特性推定方法を用いて、前記補償用騒音又
は振動発生手段から前記評価点に到る騒音又は振動の伝
達特性を推定し、前記騒音又は振動発生源の近傍に於て
検出した基準騒音又は振動を前記推定した伝達特性で処
理した信号と、前記評価点に於て検出した信号とを用い
て、前記基準騒音又は振動を入力とする適応フィルタの
係数を調整し、前記適応フィルタの出力を前記補償用騒
音又は振動発生手段への入力とすることにより前記評価
点に於る騒音又は振動を抑圧することを特徴とする騒音
又は振動の抑圧方法を開示する。
Further, according to the present invention, by generating noise or vibration from a noise or vibration generating means for compensation, noise or vibration for suppressing noise or vibration from a noise or vibration source at an evaluation point. In the suppression method, the transfer characteristic of noise or vibration reaching the evaluation point is estimated from the noise or vibration generating means for compensation by using the above-mentioned method for estimating transfer characteristic of noise or vibration, and the transfer characteristic of noise or vibration is generated. A signal obtained by processing the reference noise or vibration detected in the vicinity of the source with the estimated transfer characteristic and the signal detected at the evaluation point are used for an adaptive filter having the reference noise or vibration as an input. A noise or vibration suppressing method characterized in that the noise or vibration at the evaluation point is suppressed by adjusting a coefficient and using the output of the adaptive filter as an input to the compensating noise or vibration generating means. To disclose.

【0011】また本発明は、騒音又は振動発生源周辺に
設けた複数の測定点に於て検出した検出信号から、評価
点に於る騒音又は振動を模擬的に発生して騒音又は振動
環境を生成するための騒音又は振動環境の模擬方法に於
て、上記した騒音または振動寄与率検出手段を用いて前
記各測定点から検出した検出信号の前記評価点に於る寄
与率を検出し、該検出した寄与率及び寄与率算出時に用
いられた前記適応フィルタの係数とを用いて前記評価点
へ向けて騒音又は振動を発生する騒音又は振動発生手段
の駆動信号を生成することを特徴とする騒音又は振動環
境の模擬方法を開示する。
Further, according to the present invention, the noise or vibration environment at the evaluation point is simulated by generating the noise or vibration at the evaluation point from the detection signals detected at a plurality of measuring points provided around the noise or vibration source. In a method for simulating a noise or vibration environment for generating, detecting the contribution rate at the evaluation point of the detection signal detected from each of the measurement points using the noise or vibration contribution rate detection means described above, Noise that generates noise or vibration toward the evaluation point using the detected contribution rate and the coefficient of the adaptive filter used when calculating the contribution rate, or a drive signal of a vibration generation means is generated. Alternatively, a method of simulating a vibration environment is disclosed.

【0012】[0012]

【作用】適応フィルタを2段設け、さらに、遅延手段に
より1段目と2段目適応フィルタとの間に時間差を設
け、遅延時間を調節することにより、1段目と2段目適
応フィルタをつなぎ合わせた系の特性でより正確に対象
とする系を近似でき、また、直接波と反射波の間隔がか
なり長い系の場合には、直接波に1段目適応フィルタ
を、反射波に2段目適応フィルタを用いれば、より少な
いタップ数でその系を近似できる。
The adaptive filter is provided in two stages, and the delay means adjusts the delay time by providing a time difference between the adaptive filter in the first stage and the adaptive filter in the second stage. The target system can be approximated more accurately by the characteristics of the connected system, and in the case of a system in which the distance between the direct wave and the reflected wave is considerably long, the first-stage adaptive filter for the direct wave and the reflected wave for the reflected wave If a stage adaptive filter is used, the system can be approximated with a smaller number of taps.

【0013】又、上記の方法を用いることで、測定信号
の評価点への寄与率や、与えられた騒音または振動伝達
経路の伝達特性をより容易にかつ正確に検出または推定
できる。又、上記推定した伝達特性を利用して評価点に
於ける騒音または振動を効果的に抑圧することが出来、
更に上記検出した寄与率を用いて評価点に生じる騒音ま
たは振動を模擬的に生成することが出来る。
By using the above method, the contribution rate of the measurement signal to the evaluation point and the transfer characteristic of the given noise or vibration transfer path can be detected or estimated more easily and accurately. Further, it is possible to effectively suppress the noise or vibration at the evaluation point by using the above estimated transfer characteristic,
Further, noise or vibration generated at the evaluation point can be simulated by using the detected contribution rate.

【0014】[0014]

【実施例】以下、本発明の実施例を図面に基づいて説明
する。図2は、本発明を寄与率解析システムに適用した
時の全体構成を示す図である。同図に示した寄与率解析
システムは、機器1から発生された騒音もしくは振動の
情報を得るマイクロフォン等のセンサ2a、2b、2c
と、機器1より一定距離離れた評価点(人間の耳元付
近)に設置したマイクロフォン3と、センサ2a、2
b、2cで検出された信号がマイクロフォン3に対して
どのくらい寄与しているかを解析する寄与率解析装置4
と、表示部5から構成されている。
Embodiments of the present invention will be described below with reference to the drawings. FIG. 2 is a diagram showing an overall configuration when the present invention is applied to a contribution rate analysis system. The contribution rate analysis system shown in FIG. 1 includes sensors 2a, 2b, 2c such as microphones for obtaining information on noise or vibration generated from the device 1.
And a microphone 3 installed at an evaluation point (near human ears) a certain distance away from the device 1, and sensors 2a, 2
Contribution rate analysis device 4 for analyzing how much the signals detected by b and 2c contribute to the microphone 3.
And a display unit 5.

【0015】図1は上記の寄与率解析装置4の一実施例
を示すブロック図で、前記センサ2a、2b、2cで得
られた騒音もしくは振動の信号およびマイクロフォン3
で得られた騒音信号は、それぞれローパスフィルタ40
1〜404で帯域制限を受け、さらにA/D変換器41
1〜414でディジタル信号に変換される。センサ2
a、2b、2cに基づくディジタル信号は、それぞれ1
段目の適応ディジタルフィルタ(ADF)431〜43
3に入力され、ここでフィルタ係数との畳み込み演算が
施されて出力信号y1、y2、y3が生成される。この
とき、必要に応じて遅延回路421〜423によりディ
ジタル信号を遅延する。次に適応ディジタルフィルタ4
31〜433の出力y1〜y3は加算回路451で加算
されてS1とされ、これと評価点に設置されたマイクロ
フォンからの騒音信号をディジタル化したディジタル信
号dとの差信号e1が減算回路46で算出される。適応
アルゴリズム回路441〜443は、上記差信号e1が
最小となるように適応ディジタルフィルタ431〜43
3の係数を更新していく。
FIG. 1 is a block diagram showing an embodiment of the contribution rate analyzing device 4 described above. The noise or vibration signals obtained by the sensors 2a, 2b and 2c and the microphone 3 are shown in FIG.
The noise signals obtained in
1 to 404 limits the band, and further A / D converter 41
At 1 to 414, it is converted into a digital signal. Sensor 2
The digital signals based on a, 2b, and 2c are 1
Stage adaptive digital filters (ADF) 431-43
3 and the output signals y1, y2 and y3 are generated by performing a convolution operation with the filter coefficient. At this time, the delay circuits 421 to 423 delay the digital signal as necessary. Next, the adaptive digital filter 4
The outputs y1 to y3 of 31 to 433 are added by the adding circuit 451 to be S1, and the difference signal e1 between this and the digital signal d obtained by digitizing the noise signal from the microphone installed at the evaluation point is given by the subtracting circuit 46. It is calculated. The adaptive algorithm circuits 441 to 443 set the adaptive digital filters 431 to 43 so that the difference signal e1 is minimized.
The coefficient of 3 will be updated.

【0016】さらに、前記センサ2a、2b、2cから
得られたディジタル信号は、遅延回路424〜426で
必要に応じて一定時間遅延された後、それぞれ2段目の
適応ディジタルフィルタ(ADF)434〜436に入
力され、ここでフィルタ係数との畳み込み演算が施され
て出力信号y4、y5、y6が生成される。ここで遅延
回路424〜426の遅延時間は、遅延設定回路471
により1段目の遅延回路421〜423の遅延時間に連
動して設定される。次に適応ディジタルフィルタ434
〜436の出力y4〜y6は加算回路452で加算され
てs2とされ、これと1段目の適応ディジタルフィルタ
431〜433で処理しきれなかった前記の差信号e1
との差信号e2が減算回路462で算出される。適応ア
ルゴリズム回路444〜446は、上記差信号e2が最
小となるように、適応ディジタルフィルタ434〜43
6の係数を更新していく。なお、残差判定回路72は、
1段目の適応ディジタルフィルタに基づく差信号e1の
大きさに対するしきい値を設定し、差信号e1がそのし
きい値より大きければ、2段目の適応ディジタルフィル
タを用い、小さければ用いないようにするものである。
Further, the digital signals obtained from the sensors 2a, 2b and 2c are delayed by delay circuits 424 to 426 for a predetermined time if necessary, and then the second-stage adaptive digital filters (ADF) 434 to. It is input to 436, and convolution calculation with the filter coefficient is performed here to generate output signals y4, y5, y6. Here, the delay time of the delay circuits 424 to 426 is the delay setting circuit 471.
Is set in association with the delay time of the first-stage delay circuits 421 to 423. Next, the adaptive digital filter 434
The outputs y4 to y6 of ˜436 are added by the adder circuit 452 to be s2, and the difference signal e1 which cannot be processed by the adaptive digital filters 431 to 433 of the first stage is added to s2.
The subtraction circuit 462 calculates the difference signal e2 between The adaptive algorithm circuits 444 to 446 adjust the adaptive digital filters 434 to 43 so that the difference signal e2 is minimized.
The coefficient of 6 is updated. The residual determination circuit 72 is
A threshold value for the magnitude of the difference signal e1 based on the first-stage adaptive digital filter is set, and if the difference signal e1 is larger than the threshold value, the second-stage adaptive digital filter is used, and if it is smaller, not used. It is something to do.

【0017】以上の動作をくり返し、1段目および2段
目の適応ディジタルフィルタ系の係数が一定値に収束し
たとき、各適応ディジタルフィルタ431〜436の出
力y1〜y6から、評価点に対する各センサ設置位置か
らの騒音もしくは振動の寄与率を解析する。また、解析
された寄与率、各適応ディジタルフィルタの係数等は表
示部5に出力し表示する。
The above operation is repeated, and when the coefficients of the first-stage and second-stage adaptive digital filter systems converge to a constant value, from the outputs y1 to y6 of the respective adaptive digital filters 431 to 436, each sensor corresponding to the evaluation point. Analyze the contribution rate of noise or vibration from the installation position. The analyzed contribution rate, the coefficient of each adaptive digital filter, etc. are output to the display unit 5 and displayed.

【0018】寄与率の解析方法としては、前記センサ2
a、2b、2cに対応する適応ディジタルフィルタ出力
信号のパワーを、個々の適応ディジタルフィルタ431
〜436の出力信号y1〜y6及び差信号e2それぞれ
のパワーの和で除した値を、前記各センサに関する寄与
率とする。例えば、前記センサ2aの寄与率は、次式で
表わされる。
As a method of analyzing the contribution rate, the sensor 2 is used.
The powers of the adaptive digital filter output signals corresponding to a, 2b, and 2c are set to the individual adaptive digital filters 431.
The value obtained by dividing the power of each of the output signals y1 to y6 of ˜436 and the difference signal e2 is the contribution rate for each sensor. For example, the contribution rate of the sensor 2a is expressed by the following equation.

【数1】 (Equation 1)

【0019】本実施例によると、1段目の適応ディジタ
ルフィルタ431〜436へ入力する手前に、遅延回路
421〜423を設け、これにより遅延時間を調節する
ことができる。従って直接波の到達に時間のかかるよう
な系の場合には、前記遅延回路にて時間遅れ分だけ遅延
させることにより、適応ディジタルフィルタの同一タッ
プ数で表現できるインパルス応答が長くできる。また、
2段目の適応ディジタルフィルタ434〜436へ入力
する手前に、遅延回路424〜426を設けたことによ
り、1段目の適応ディジタルフィルタと組み合わせて、
様々な音環境でのインパルス応答を表現できる。なお、
これら遅延回路の遅延時間の設定は、適応ディジタルフ
ィルタの遅延時間を種々変化させて適応ディジタルフィ
ルタの係数調整を行い、その誤差がより小さくなるよう
な遅延時間とすることにより行う。
According to this embodiment, the delay circuits 421 to 423 are provided before inputting to the first-stage adaptive digital filters 431 to 436, whereby the delay time can be adjusted. Therefore, in the case of a system in which it takes a long time for a direct wave to reach, the impulse response that can be represented by the same number of taps of the adaptive digital filter can be lengthened by delaying the delay circuit by the time delay. Also,
By providing the delay circuits 424 to 426 before inputting to the second-stage adaptive digital filters 434 to 436, in combination with the first-stage adaptive digital filters,
It can express impulse response in various sound environments. In addition,
The delay time of these delay circuits is set by varying the delay time of the adaptive digital filter to adjust the coefficient of the adaptive digital filter so that the error becomes smaller.

【0020】図3はそのインパルス応答の1例を示して
いる。この図に示すような直接波と反射波の時間間隔が
長いような系では、遅延回路により、直接波は1段目適
応フィルタ、反射波は2段目適応フィルタと分けてイン
パルス応答を求めることができ、タップ数を有効に使え
る。また、1段目、2段目の適応ディジタルフィルタの
タップ数を、それぞれ可変設定可能な構造としておけ
ば、遅延回路と組合せて適応ディジタルフィルタの調整
がより正確且つ容易に行えるようになる。さらに、適応
フィルタは2段だけではなく、状況、周辺環境に応じ
て、3段、4段と並列に複数の適応フィルタを設けるこ
とができる。また、差信号e1及びe2をD/A変換し
て音声として出力することで、適応ディジタルフィルタ
の集束状況を知ることができる。
FIG. 3 shows an example of the impulse response. In a system where the time interval between the direct wave and the reflected wave is long as shown in this figure, the delay circuit should be used to determine the impulse response by dividing the direct wave from the first-stage adaptive filter and the reflected wave from the second-stage adaptive filter. It is possible to effectively use the number of taps. If the number of taps of the first-stage and second-stage adaptive digital filters can be variably set, the adaptive digital filters can be adjusted more accurately and easily in combination with the delay circuit. Furthermore, the number of adaptive filters is not limited to two, and a plurality of adaptive filters can be provided in parallel with three stages and four stages depending on the situation and surrounding environment. Further, the focusing state of the adaptive digital filter can be known by D / A converting the difference signals e1 and e2 and outputting as voice.

【0021】本実施例は、騒音を発生する一般機器、例
えば、掃除機、空調機、車両等にも適用できる。さら
に、騒音もしくは振動の情報を得るセンサ、及びスピー
カ等はその個数を増やすことで測定点を増やすことがで
きる。また、センサとしては、マイクロフォン、振動セ
ンサだけでなく、時間的に変動する時系列信号を検出で
きるセンサ、例えば、アンテナ、光センサ等であっても
本実施例を適用できることはいうまでもない。
The present embodiment can be applied to general equipment that generates noise, such as a vacuum cleaner, an air conditioner, and a vehicle. Furthermore, the number of measurement points can be increased by increasing the number of sensors, speakers, and the like that obtain noise or vibration information. Further, it goes without saying that the present embodiment can be applied not only to the microphone and the vibration sensor but also to a sensor capable of detecting a time-series signal that fluctuates with time, such as an antenna or an optical sensor.

【0022】図4は、本発明になる能動消音装置の一実
施例を示すブロック図である。本装置は、機器6からの
騒音の情報を得るマイクロフォン等の基準センサ7と、
騒音を低減したい目標位置に前記騒音から伝播する音波
と逆位相同振幅の音波(付加音波)を送るための付加音
源であるスピーカ8と、前記目標位置に設置されて機器
6からの音波及びスピーカ8からの付加音波を検出する
消音エラーセンサとしてのマイクロフォン10と、スピ
ーカ8からの付加音波の発生を制御するコントローラ9
から成っている。
FIG. 4 is a block diagram showing an embodiment of the active muffler according to the present invention. This device includes a reference sensor 7 such as a microphone that obtains noise information from the device 6,
A speaker 8 which is an additional sound source for transmitting a sound wave (additional sound wave) having the same phase and an opposite phase as the sound wave propagating from the noise to a target position where noise is desired to be reduced, and a sound wave from the device 6 installed at the target position and a speaker Microphone 10 as a muffling error sensor for detecting the additional sound wave from 8 and controller 9 for controlling the generation of the additional sound wave from speaker 8.
Consists of

【0023】コントローラ9には、ここの制御に必要
な、スピーカ8への入力信号とマイクロフォン10の出
力信号との間の推定伝達関数93が用いられている。図
5は、この推定伝達関数93の推定を行うための測定シ
ステムを示すブロック図で、本発明の適応信号処理方法
を用いたシステムとなっている。そして図5において、
スピーカ8とマイクロフォン10は、図4の実際の減音
装置で用いるときと同じ環境下の同じ位置に配置されて
いるものとする。図5に於て、雑音信号発生源501か
らの基準信号s21はA/D変換器502でディジタル
形態の基準信号s22とされ、該基準信号s22はD/
A変換器503とパワーアンプ504を介してスピーカ
8に与えられ、音波が放射される。また、前記基準信号
s22は、1段目の適応ディジタルフィルタ(ADF)
505に入力され、フィルタ係数との畳み込み演算が行
われて信号s23が生成される。スピーカ8からの放射
音波はマイクロフォン10で検出され、A/D変換器5
09によりディジタル化され、ディジタル信号s25と
される。この信号s25と上記信号s23との差信号e
21が減算回路510で求められ、この差信号e21が
最小となるように、適応アルゴリズム507は適応ディ
ジタルフィルタ505の係数を更新していく。このと
き、必要に応じて遅延回路511により基準信号s22
を遅延する。
The controller 9 uses an estimated transfer function 93 between the input signal to the speaker 8 and the output signal of the microphone 10, which is necessary for the control here. FIG. 5 is a block diagram showing a measurement system for estimating the estimated transfer function 93, which is a system using the adaptive signal processing method of the present invention. And in FIG.
It is assumed that the speaker 8 and the microphone 10 are arranged at the same position in the same environment as that used in the actual sound reduction device of FIG. In FIG. 5, the reference signal s21 from the noise signal source 501 is converted into a digital reference signal s22 by the A / D converter 502, and the reference signal s22 is D /
Sound waves are emitted to the speaker 8 via the A converter 503 and the power amplifier 504 and emitted. The reference signal s22 is the first-stage adaptive digital filter (ADF).
The signal s23 is input to 505 and convolved with the filter coefficient to generate a signal s23. The sound wave radiated from the speaker 8 is detected by the microphone 10, and the A / D converter 5
It is digitized by 09 and becomes a digital signal s25. Difference signal e between the signal s25 and the signal s23
21 is obtained by the subtraction circuit 510, and the adaptive algorithm 507 updates the coefficient of the adaptive digital filter 505 so that the difference signal e21 becomes the minimum. At this time, the reference signal s22 is set by the delay circuit 511 as necessary.
Delay.

【0024】さらに、前記基準信号s22は遅延回路5
12により必要に応じて一定時間遅延され、2段目の適
応ディジタルフィルタ(ADF)506に入力され、フ
ィルタ係数との畳み込み演算が行われて信号s24が生
成される。この出力信号s24と、1段目の適応ディジ
タルフィルタで処理しきれなかった前記差信号e21と
の差e22が減算回路513で求められ、この差信号e
22が最小となるように、適応アルゴリズム508は適
応ディジタルフィルタ506の係数を更新していく。以
上の動作をくり返して各適応ディジタルフィルタの係数
を収束させ、その収束した値から推定伝達関数93が求
められる。
Further, the reference signal s22 is the delay circuit 5
The signal is delayed by a predetermined time by 12 and input to the second-stage adaptive digital filter (ADF) 506, and a convolution operation with the filter coefficient is performed to generate a signal s24. A subtraction circuit 513 obtains a difference e22 between the output signal s24 and the difference signal e21 that cannot be processed by the first-stage adaptive digital filter.
The adaptive algorithm 508 updates the coefficients of the adaptive digital filter 506 so that 22 is minimized. The above operation is repeated to converge the coefficients of each adaptive digital filter, and the estimated transfer function 93 is obtained from the converged values.

【0025】上記のようにして求められた推定伝達関数
93を用いた図4の能動消音装置の動作を次ぎに説明す
る。機器1の騒音情報として前記基準センサ7から出力
された検出信号s11は、前記A/D変換器92でディ
ジタル化されて基準信号s12とされ、さらに適応ディ
ジタルフィルタ90の畳み込み演算が行われて信号s1
6に変換される。この信号s16は、D/A変換器94
でアナログ化され、パワーアンプで増幅されてスピーカ
8から放射される。前記マイクロフォン10はこの放射
音と機器1からの騒音の加わった音波を検出し、その出
力はA/D変換器95でディジタル信号s14とされ
る。一方、前記基準信号s12と前記推定伝達関数93
との畳み込み演算により仮想入力信号s13が生成さ
れ、この信号s13と前記信号s14を用いて、適応ア
ルゴリズム91により信号s14の自乗値が最小となる
ように、適応ディジタルフィルタ90の係数が更新され
る。
The operation of the active silencer of FIG. 4 using the estimated transfer function 93 obtained as described above will be described below. The detection signal s11 output from the reference sensor 7 as the noise information of the device 1 is digitized by the A / D converter 92 to be the reference signal s12, and the signal is obtained by performing the convolution operation of the adaptive digital filter 90. s1
Converted to 6. This signal s16 is the D / A converter 94.
Is converted into analog, is amplified by the power amplifier, and is emitted from the speaker 8. The microphone 10 detects the radiated sound and the sound wave to which the noise from the device 1 is added, and the output thereof is converted into a digital signal s14 by the A / D converter 95. Meanwhile, the reference signal s12 and the estimated transfer function 93
A virtual input signal s13 is generated by a convolution operation with, and the adaptive algorithm 91 updates the coefficient of the adaptive digital filter 90 using the signal s13 and the signal s14 so that the square value of the signal s14 is minimized. .

【0026】本実施例は、騒音を発生する一般機器、例
えば、空調機、車両等にも適用できる。また上記実施例
において、騒音もしくは振動の情報を得るセンサ、及び
スピーカ等は個数を増やすことができる。さらに、図5
の推定伝達関数を求める際の適応フィルタは2段だけで
はなく、状況、周辺環境に応じて、3段、4段と並列に
複数の適応フィルタを設けることができる。
This embodiment can also be applied to general equipment that generates noise, such as air conditioners and vehicles. Further, in the above-described embodiment, the number of sensors, speakers, etc. for obtaining noise or vibration information can be increased. Furthermore, FIG.
The number of adaptive filters for obtaining the estimated transfer function is not only two, but a plurality of adaptive filters can be provided in parallel with three and four stages depending on the situation and the surrounding environment.

【0027】図6は、本発明の適応信号処理方法を用い
た音環境シミュレータの実施例を示すブロック図で、騒
音源である機器11の騒音もしくは振動の情報を得るマ
イクロフォン等のセンサ12a、12b、12cと、機
器11より一定距離離れた評価点(人間の耳元付近)に
設置されたマイクロフォン13と、センサ12a、12
b、12cで検出された信号がマイクロフォン13に対
してどのくらい寄与しているかを解析するための寄与率
解析装置14と、その寄与率解析装置14で得られた寄
与率等をセンサ12a、12b、12cで検出した信号
対応に分配し、評価点での騒音を合成するための音合成
手段15と、合成した音を出力するスピーカ16a、1
6b、及び合成音を再生する部屋17に設置されたダミ
ーヘッド18から成っている。ここで音声合成部15の
構成は図7に示されている。
FIG. 6 is a block diagram showing an embodiment of a sound environment simulator using the adaptive signal processing method of the present invention. Sensors 12a, 12b such as microphones for obtaining information on noise or vibration of the device 11 as a noise source. , 12c, a microphone 13 installed at an evaluation point (near human ears) a certain distance away from the device 11, and sensors 12a, 12
b, 12c, a contribution rate analyzing device 14 for analyzing how much the signals detected by the b and 12c contribute to the microphone 13, and the contribution rates and the like obtained by the contribution rate analyzing device 14 are measured by the sensors 12a, 12b. The sound synthesizing means 15 for synthesizing the noise detected at 12c and synthesizing the noise at the evaluation point, and the speakers 16a, 1 for outputting the synthesized sound.
6b, and a dummy head 18 installed in a room 17 for reproducing a synthetic sound. Here, the configuration of the voice synthesis unit 15 is shown in FIG.

【0028】この構成において、寄与率解析装置14
は、図1の実施例で説明したものであって、センサ12
a〜12cが検出した騒音が、マイクロフォン13で検
出した評価点の騒音にどれくらい寄与しているかの寄与
率を算出して出力する。この各寄与率は、音声合成手段
15の寄与率分配回路151により各センサ12a〜1
2c対応の寄与率に分配され、ディジタルフィルタ15
2a〜152cへそれぞれ入力される。また、前記寄与
率解析手段14の収束時の各適応ディジタルフィルタ係
数がディジタルフィルタ152a〜152cにその係数
としてコピーされる。各ディジタルフィルタ152a〜
152cに於ては、そのフィルタの伝達関数と入力寄与
率との畳み込み演算が行われ、ディジタル信号sa、s
b、scが出力される。これらの信号は加算回路153
で加算され、合成信号sをとして空間補正回路154へ
入力される。
In this configuration, the contribution rate analyzing device 14
Is the same as that described in the embodiment of FIG.
The contribution rate of how much the noise detected by a to 12c contributes to the noise at the evaluation point detected by the microphone 13 is calculated and output. The contribution rate distribution circuit 151 of the voice synthesizing unit 15 determines the contribution rate of each of the sensors 12a-1.
The digital filter 15 is distributed to the contribution rate corresponding to 2c.
2a to 152c, respectively. Further, each adaptive digital filter coefficient at the time of convergence of the contribution rate analyzing means 14 is copied as a coefficient to the digital filters 152a to 152c. Each digital filter 152a-
In 152c, the convolution operation of the transfer function of the filter and the input contribution rate is performed, and the digital signals sa, s
b and sc are output. These signals are added to the addition circuit 153.
Are added and the combined signal s is input to the space correction circuit 154.

【0029】記憶媒体155には、様々な伝達特性(残
響特性等)を持つ空間での伝達特性の測定値がデータベ
ース化して格納されており、空間補正回路154では、
記憶媒体155から呼び出した伝達特性と上記合成信号
sとの畳み込み演算を行う。さらに、合成音を再生する
部屋17内でのスピーカ16a、16bから評価する人
間もしくはダミーヘッド18までの逆伝達特性を表わす
逆フィルタ156と上記空間補正回路154出力との畳
み込み演算を行い、演算結果をD/A変換器でアナログ
信号に変換し、スピーカ16a、16bから出力する。
このようにして様々な特性を有する空間での製品音が再
現できる。
The storage medium 155 stores a database of measured values of transfer characteristics in a space having various transfer characteristics (reverberation characteristics, etc.), and the space correction circuit 154 stores:
The convolution operation of the transfer characteristic called from the storage medium 155 and the composite signal s is performed. Further, a convolution operation of the inverse filter 156 representing the inverse transfer characteristic from the speaker 16a, 16b in the room 17 for reproducing the synthetic sound to the human or the dummy head 18 to be evaluated and the output of the spatial correction circuit 154 is performed, and the operation result Is converted into an analog signal by the D / A converter and output from the speakers 16a and 16b.
In this way, product sounds in a space having various characteristics can be reproduced.

【0030】本実施例は騒音を発生する一般機器、例え
ば、掃除機、空調機、車両等にも適用出来る。さらに、
騒音もしくは振動の情報を得るセンサ、および、スピー
カ等は個数を増やすことが出来る。さらに、寄与率解析
装置14に於る適応フィルタは2段だけではなく、状
況、周辺環境に応じて、3段、4段と並列に複数の適応
フィルタを設けるようにしてもよい。
This embodiment can be applied to general equipment that generates noise, such as a vacuum cleaner, an air conditioner, and a vehicle. further,
It is possible to increase the number of sensors, speakers, and the like for obtaining information on noise or vibration. Furthermore, the number of adaptive filters in the contribution rate analyzing device 14 is not limited to two, and a plurality of adaptive filters may be provided in parallel with three or four stages depending on the situation and the surrounding environment.

【0031】[0031]

【発明の効果】本発明によれば、適応信号処理におい
て、適応フィルタを2段設け、さらに、各適応フィルタ
と遅延時間の組み合わせによって、高残響場等様々な環
境下での音響系のインパルス応答の算出や音場の伝達関
数の推定を高精度に行え、また同様な環境下での消音装
置や音環境シミュレータの実現が可能になる効果があ
る。
According to the present invention, in adaptive signal processing, two stages of adaptive filters are provided, and by combining each adaptive filter and delay time, impulse response of an acoustic system under various environments such as a high reverberant field. Can be calculated and the transfer function of the sound field can be estimated with high accuracy, and a muffling device and a sound environment simulator can be realized under the same environment.

【図面の簡単な説明】[Brief description of the drawings]

【図1】寄与率解析装置の一実施例を示すブロック図で
ある。
FIG. 1 is a block diagram showing an embodiment of a contribution rate analyzing apparatus.

【図2】寄与率解析システムの全体構成を示すブロック
図である。
FIG. 2 is a block diagram showing an overall configuration of a contribution rate analysis system.

【図3】インパルスレスポンスの一例を示す図である。FIG. 3 is a diagram showing an example of an impulse response.

【図4】能動消音装置の全体構成を示すブロック図であ
る。
FIG. 4 is a block diagram showing an overall configuration of an active silencer.

【図5】図4に於る推定伝達関数を求めるシステムの一
実施例を示す図である。
5 is a diagram showing an embodiment of a system for obtaining an estimated transfer function in FIG.

【図6】本発明の音環境シミュレータの一実施例を示す
全体構成図である。
FIG. 6 is an overall configuration diagram showing an embodiment of a sound environment simulator of the present invention.

【図7】図6の音声合成手段の詳細を示すブロック図で
ある。
FIG. 7 is a block diagram showing details of a voice synthesizing unit in FIG.

【符号の説明】[Explanation of symbols]

1、6、11 騒音を発生する機器 2a、2b、2c、7 センサ 3、10、13 マイクロフォン 4、14 寄与率解析装置 421〜426 遅延回路 431〜433 1段目の適応ディジタルフィルタ 434〜436 2段目の適応ディジタルフィルタ 451、452 加算回路 461、462 減算回路 481 残差判定回路 505、506 適応ディジタルフィルタ 510、513 減算回路 511、512 遅延回路 8 消音用スピーカ 9 コントローラ 90 適応ディジタルフィルタ 93 推定伝達関数 16a、16b 音再現用スピーカ 15 音合成手段 18 人間もしくはダミーヘッド 1, 6, 11 Devices that generate noise 2a, 2b, 2c, 7 Sensors 3, 10, 13 Microphones 4, 14 Contribution rate analyzers 421-426 Delay circuits 431-433 First-stage adaptive digital filters 434-436 2 Stage adaptive digital filter 451, 452 Adder circuit 461, 462 Subtractor circuit 481 Residual difference determination circuit 505, 506 Adaptive digital filter 510, 513 Subtractor circuit 511, 512 Delay circuit 8 Silence speaker 9 Controller 90 Adaptive digital filter 93 Estimated transfer Functions 16a, 16b Sound reproduction speaker 15 Sound synthesis means 18 Human or dummy head

───────────────────────────────────────────────────── フロントページの続き (72)発明者 大木 克夫 茨城県土浦市神立町502番地 株式会社日 立製作所機械研究所内 (72)発明者 太田 義注 茨城県土浦市神立町502番地 株式会社日 立製作所機械研究所内 ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Katsuo Oki 502 Jinritsu-cho, Tsuchiura-shi, Ibaraki Machinery Research Institute, Hiritsu Manufacturing Co., Ltd. Machinery Research Laboratory

Claims (14)

【特許請求の範囲】[Claims] 【請求項1】 第1の1個又は複数個の時系列信号を1
個又は複数個の1段目適応フィルタの各々に入力し、該
適応フィルタの出力信号の和と前記第1の時系列信号と
は別の第2の時系列信号との差を表す第1の残差信号が
最小になるように、前記1段目適応フィルタの係数を更
新するとともに、 前記第1の時系列信号を1個又は複数個の2段目適応フ
ィルタの各々に入力し、該適応フィルタの出力信号の和
と前記第1の残差信号との差を表す第2の残差信号が最
小になるように前記2段目適応フィルタを更新すること
を特徴とする適応信号処理方法。
1. A first time series signal or a plurality of time series signals
The first or second plurality of first-stage adaptive filters are input to each of the first-stage adaptive filters, and a first time-sequential signal different from the sum of the output signals of the adaptive filter and the second time-series signal is displayed. The coefficient of the first-stage adaptive filter is updated so that the residual signal is minimized, and the first time-series signal is input to one or more second-stage adaptive filters, respectively, An adaptive signal processing method, wherein the second-stage adaptive filter is updated so that a second residual signal representing a difference between the sum of output signals of the filters and the first residual signal is minimized.
【請求項2】 前記第1の時系列信号を入力とする3段
目あるいは必要に応じてそれ以降の適応フィルタを設
け、該各適応フィルタの係数をその出力の和と1段前の
適応フィルタにより生成された残差信号との差が最小に
なるように更新することを特徴とする請求項1記載の適
応信号処理方法。
2. An adaptive filter provided in a third stage or an optional subsequent stage, which receives the first time-series signal as an input, and the coefficient of each adaptive filter is the sum of its output and the adaptive filter in the previous stage. The adaptive signal processing method according to claim 1, further comprising updating so as to minimize a difference from the residual signal generated by.
【請求項3】 前記1段目及び2段目の適応フィルタに
入力される前記第1の時系列信号を、その遅延時間が可
変な遅延手段により遅延できるようにしたことを特徴と
する請求項1記載の適応信号処理方法。
3. The first time series signal input to the first-stage and second-stage adaptive filters can be delayed by delay means having a variable delay time. 1. The adaptive signal processing method described in 1.
【請求項4】 前記2段目適応フィルタ側の遅延時間
を、前記1段目適応フィルタ側の遅延時間に連動して設
定するようにしたことを特徴とする請求項3記載の適応
信号処理方法。
4. The adaptive signal processing method according to claim 3, wherein the delay time on the second-stage adaptive filter side is set in association with the delay time on the first-stage adaptive filter side. .
【請求項5】 前記1段目適応フィルタ側の遅延時間
を、当該遅延時間を変化させたときに得られる前記第1
の残差信号が最小となる値に設定し、 さらに前記2段目適応フィルタ側の遅延時間を、当該遅
延時間を変化させたときに得られる前記第2の残差信号
が最小となる値に設定することを特徴とする請求項3記
載の適応信号処理方法。
5. The delay time on the side of the first-stage adaptive filter is the first time obtained when the delay time is changed.
Is set to a value that minimizes the residual signal, and the delay time on the second-stage adaptive filter side is set to a value that minimizes the second residual signal obtained when the delay time is changed. The adaptive signal processing method according to claim 3, wherein the adaptive signal processing method is set.
【請求項6】 前記1段目及び2段目適応フィルタの各
々を、そのタップ数が可変な構造としたことを特徴する
請求項1記載の適応信号処理方法。
6. The adaptive signal processing method according to claim 1, wherein each of the first-stage adaptive filter and the second-stage adaptive filter has a structure in which the number of taps is variable.
【請求項7】 前記第1の残差信号に対するしきい値を
設け、前記残差信号が前記しきい値より大きいときだけ
前記2段目の適応フィルタに於る係数の更新処理を行う
ようにしたことを特徴とする請求項1記載の適応信号処
理方法。
7. A threshold value for the first residual signal is provided, and the coefficient updating process in the second-stage adaptive filter is performed only when the residual signal is larger than the threshold value. The adaptive signal processing method according to claim 1, wherein:
【請求項8】 請求項1項ないし7項の内の1つに記載
の適応信号処理方法を用いて、前記第1の時系列信号
を、騒音又は振動発生源周辺に設けた複数の測定点に於
て検出した複数の騒音又は振動の検出信号とし、且つ前
記第2の時系列信号を前記騒音又は振動発生源からの騒
音又は振動を評価するために設けた評価点に於て検出し
た騒音又は振動の検出信号としたときの適応信号処理を
行い、 該処理終了時の前記適応フィルタの各出力から、当該出
力に対応する前記測定点に於る検出信号の、前記評価点
に於る検出信号に対する寄与率を求めることを特徴とす
る騒音又は振動の寄与率の検出方法。
8. A plurality of measurement points at which the first time-series signal is provided around a noise or vibration source by using the adaptive signal processing method according to any one of claims 1 to 7. Noise detected at a plurality of noises or vibrations, and the second time-series signal detected at an evaluation point provided for evaluating noises or vibrations from the noise or vibration source. Alternatively, adaptive signal processing is performed when a detection signal of vibration is detected, and the detection signal at the measurement point corresponding to the output is detected from each output of the adaptive filter at the end of the processing at the evaluation point. A method for detecting a contribution rate of noise or vibration, characterized by obtaining a contribution rate to a signal.
【請求項9】 前記第1の残差信号及び前記第2の残差
信号をそれぞれD/A変換し、音として試聴できるよう
にしたことを特徴とする請求項8記載の騒音又は振動寄
与率の検出方法。
9. The noise or vibration contribution rate according to claim 8, wherein the first residual signal and the second residual signal are respectively D / A converted so that they can be auditioned as sounds. Detection method.
【請求項10】 請求項1項ないし7項の内の1つに記
載の適応信号処理方法を用いて、前記第1の時系列信号
を、騒音又は振動発生手段に印加する騒音又は振動信号
とし、且つ前記第2の時系列信号を、前記騒音又は振動
発生手段から発生させた騒音又は振動を評価する評価点
に於て検出された検出信号として適応信号処理を行い、 該処理終了時の前記各適応ディジタルフィルタの係数か
ら、前記騒音又は振動発生手段から前記評価点に到る騒
音又は振動の伝達特性を推定することを特徴とする騒音
又は振動の伝達特性推定方法。
10. The adaptive signal processing method according to claim 1, wherein the first time-series signal is a noise or vibration signal to be applied to a noise or vibration generating means. And, the second time-series signal is subjected to adaptive signal processing as a detection signal detected at an evaluation point for evaluating noise or vibration generated from the noise or vibration generating means, and at the end of the processing A noise or vibration transfer characteristic estimating method for estimating a transfer characteristic of noise or vibration reaching the evaluation point from the noise or vibration generating means from a coefficient of each adaptive digital filter.
【請求項11】 補償用の騒音又は振動発生手段から騒
音又は振動を発生することにより、騒音又は振動発生源
からの騒音又は振動を評価点に於て抑圧するための騒音
又は振動の抑圧方法に於て、 請求項10記載の騒音又は振動の伝達特性推定方法を用
いて、前記補償用騒音又は振動発生手段から前記評価点
に到る騒音又は振動の伝達特性を推定し、 前記騒音又は振動発生源の近傍に於て検出した基準騒音
又は振動を前記推定した伝達特性で処理した信号と、前
記評価点に於て検出した信号とを用いて、前記基準騒音
又は振動を入力とする適応フィルタの係数を調整し、 前記適応フィルタの出力を前記補償用騒音又は振動発生
手段への入力とすることにより前記評価点に於る騒音又
は振動を抑圧することを特徴とする騒音又は振動の抑圧
方法。
11. A noise or vibration suppressing method for suppressing noise or vibration from a noise or vibration source at an evaluation point by generating noise or vibration from a noise or vibration generating means for compensation. In the noise or vibration transfer characteristic estimating method according to claim 10, the noise or vibration transfer characteristic reaching the evaluation point is estimated from the compensating noise or vibration generating means to generate the noise or vibration. A signal obtained by processing the reference noise or vibration detected in the vicinity of the source with the estimated transfer characteristic and the signal detected at the evaluation point are used for an adaptive filter having the reference noise or vibration as an input. A noise or vibration suppressing method characterized by suppressing a noise or a vibration at the evaluation point by adjusting a coefficient and using an output of the adaptive filter as an input to the compensating noise or vibration generating means. .
【請求項12】 騒音又は振動発生源周辺に設けた複数
の測定点に於て検出した検出信号から、評価点に於る騒
音又は振動を模擬的に発生して騒音又は振動環境を生成
するための騒音又は振動環境の模擬方法に於て、 請求項8記載の騒音または振動寄与率検出手段を用いて
前記各測定点から検出した検出信号の前記評価点に於る
寄与率を検出し、該検出した寄与率及び寄与率算出時に
用いられた前記適応フィルタの係数とを用いて前記評価
点へ向けて騒音又は振動を発生する騒音又は振動発生手
段の駆動信号を生成することを特徴とする騒音又は振動
環境の模擬方法。
12. A noise or vibration environment is generated by simulating noise or vibration at an evaluation point from detection signals detected at a plurality of measurement points provided around a noise or vibration source. In the method for simulating a noise or vibration environment, the noise or vibration contribution rate detecting means according to claim 8 is used to detect a contribution rate of the detection signal detected from each of the measurement points at the evaluation point, Noise that generates noise or vibration toward the evaluation point using the detected contribution rate and the coefficient of the adaptive filter used when calculating the contribution rate, or a drive signal of a vibration generation means is generated. Or a method of simulating a vibration environment.
【請求項13】 第1の1個又は複数個の時系列信号を
処理するための1個又は複数個の1段目適応フィルタ
と、 該適応フィルタの出力信号の和と前記第1の時系列信号
とは別の第2の時系列信号との差を第1の残差信号とし
て算出するための第1の加減算手段と、 前記第1の残差が最小となるように前記1段目適応フィ
ルタの係数を更新するための1段目適応処理手段と、 前記第1の時系列信号を処理するための1個又は複数個
の2段目適応フィルタと、 該適応フィルタの出力信号の和と前記第1の残差信号と
の差を第2の残差信号として算出するための第2の加減
算手段と、 前記第2の残差が最小となるように前記2段目適応フィ
ルタの係数を更新するための2段目適応処理手段と、 を備えたことを特徴とする適応信号処理装置。
13. One or a plurality of first stage adaptive filters for processing a first one or a plurality of time series signals, a sum of output signals of the adaptive filters and the first time series. First addition / subtraction means for calculating a difference from a second time-series signal different from the signal as a first residual signal; and the first stage adaptation so that the first residual is minimized. A first stage adaptive processing means for updating the coefficient of the filter; one or a plurality of second stage adaptive filters for processing the first time series signal; and a sum of output signals of the adaptive filter. Second adding / subtracting means for calculating a difference from the first residual signal as a second residual signal, and coefficients of the second-stage adaptive filter so that the second residual is minimized. An adaptive signal processing device comprising: a second stage adaptive processing means for updating;
【請求項14】 前記1段目及び2段目の適応フィルタ
に入力される前記第1の時系列信号を遅延させるため
の、可変遅延手段を設けたことを特徴とする請求項13
記載の適応信号処理装置。
14. A variable delay means for delaying the first time-series signal input to the first-stage and second-stage adaptive filters.
An adaptive signal processing device as described.
JP16025995A 1995-06-27 1995-06-27 Adaptive signal processing method Expired - Fee Related JP3393738B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP16025995A JP3393738B2 (en) 1995-06-27 1995-06-27 Adaptive signal processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP16025995A JP3393738B2 (en) 1995-06-27 1995-06-27 Adaptive signal processing method

Publications (2)

Publication Number Publication Date
JPH0915035A true JPH0915035A (en) 1997-01-17
JP3393738B2 JP3393738B2 (en) 2003-04-07

Family

ID=15711146

Family Applications (1)

Application Number Title Priority Date Filing Date
JP16025995A Expired - Fee Related JP3393738B2 (en) 1995-06-27 1995-06-27 Adaptive signal processing method

Country Status (1)

Country Link
JP (1) JP3393738B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010236944A (en) * 2009-03-30 2010-10-21 Chubu Electric Power Co Inc Method and system for searching for sound source and vibration source
AT512755A3 (en) * 2013-06-18 2015-01-15 Avl List Gmbh Method for creating a three-dimensional acoustic network

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010236944A (en) * 2009-03-30 2010-10-21 Chubu Electric Power Co Inc Method and system for searching for sound source and vibration source
AT512755A3 (en) * 2013-06-18 2015-01-15 Avl List Gmbh Method for creating a three-dimensional acoustic network
AT512755B1 (en) * 2013-06-18 2015-03-15 Avl List Gmbh Method for creating a three-dimensional acoustic network

Also Published As

Publication number Publication date
JP3393738B2 (en) 2003-04-07

Similar Documents

Publication Publication Date Title
US10607592B2 (en) Noise reducing device, noise reducing method, noise reducing program, and noise reducing audio outputting device
JP6111319B2 (en) Apparatus and method for improving perceived quality of sound reproduction by combining active noise canceling and perceptual noise compensation
US7577262B2 (en) Microphone device and audio player
US5774562A (en) Method and apparatus for dereverberation
KR101149591B1 (en) Audio signal dereverberation
US20150208169A1 (en) Sound field correction apparatus, control method thereof, and computer-readable storage medium
CN105103219A (en) Noise reduction method
US20090169028A1 (en) Sound Field Control Apparatus And Sound Field Control Method
JPH0511771A (en) Noise control device
US20230362542A1 (en) Audio controller for a semi-adaptive active noise reduction device
JP2004078021A (en) Method, device, and program for sound pickup
US20080285768A1 (en) Method and System for Modifying and Audio Signal, and Filter System for Modifying an Electrical Signal
JP3393738B2 (en) Adaptive signal processing method
JP2003250193A (en) Echo elimination method, device for executing the method, program and recording medium therefor
JP3380571B2 (en) Noise canceling device
JP4522509B2 (en) Audio equipment
JP2000099039A (en) Improving method for clearness of amplified voice and device
JP3556427B2 (en) Method for determining control band of audio device
JP4176938B2 (en) Delay time setting method
JPH11167383A (en) Method for checking connection of adaptive equalization system
JP3264141B2 (en) Sound environment simulator
JP2870440B2 (en) Three-dimensional sound field reproduction method
WO2021075108A1 (en) Signal processing device and method, and program
JPH0758223B2 (en) Method of measuring sound deadening / damping effect, measuring device, and signal source search device
WO2023276424A1 (en) Howling suppression device, howling suppression method, and howling suppression program

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees