JP3293922B2 - Active noise control device - Google Patents
Active noise control deviceInfo
- Publication number
- JP3293922B2 JP3293922B2 JP01413293A JP1413293A JP3293922B2 JP 3293922 B2 JP3293922 B2 JP 3293922B2 JP 01413293 A JP01413293 A JP 01413293A JP 1413293 A JP1413293 A JP 1413293A JP 3293922 B2 JP3293922 B2 JP 3293922B2
- Authority
- JP
- Japan
- Prior art keywords
- filter
- noise
- signal
- identification
- transfer function
- 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.)
- Expired - Fee Related
Links
Landscapes
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)
- Exhaust Silencers (AREA)
Description
【0001】[0001]
【産業上の利用分野】この発明は、騒音源から伝達され
る騒音に制御音源から発生される制御音を干渉させるこ
とにより騒音の低減を図る能動型騒音制御装置に関し、
特に、制御音を発生するラウドスピーカと残留騒音を検
出するマイクロフォンとの間の伝達関数を、少ない演算
量で測定できるようにしたものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an active noise control device for reducing noise by causing a control sound generated from a control sound source to interfere with noise transmitted from a noise source.
In particular, a transfer function between a loudspeaker that generates a control sound and a microphone that detects residual noise can be measured with a small amount of calculation.
【0002】[0002]
【従来の技術】従来の能動型騒音制御装置として、英国
特許第2149614号や特表平1−501344号公
報等に記載のものがある。これら従来の装置は、例えば
航空機の客室等の閉空間に適用される騒音低減装置であ
って、そのような閉空間内の複数の位置に設置され音圧
を検出するマイクロフォンと、その閉空間に制御音を発
生する複数のラウドスピーカとを備え、騒音源の騒音発
生状態に基づいて、閉空間に伝達される騒音と逆位相の
制御音をラウドスピーカから発生させて騒音を打ち消し
ている。2. Description of the Related Art Conventional active noise control devices include those described in British Patent No. 2,149,614 and Japanese Patent Publication No. 1-501344. These conventional devices are, for example, noise reduction devices applied to a closed space such as a cabin of an aircraft, and microphones installed at a plurality of positions in such a closed space to detect sound pressure, and a microphone installed in the closed space. A plurality of loudspeakers that generate control sounds are provided, and based on the noise generation state of the noise source, the loudspeakers generate control sounds having a phase opposite to the noise transmitted to the closed space to cancel the noises.
【0003】そして、ラウドスピーカから発せられる制
御音の生成方法として、PROCEEDINGS OF THE IEEE,VOL.
63 PAGE 1692,1975,“ADAPTIVE NOISE CANSELLATION :
PRINCIPLES AND APPLICATIONS ”で述べられている‘WI
DROW LMS’アルゴリズムを多チャンネルに展開したアル
ゴリズムを適用している。その内容は、上記特許の発明
者による論文、“A MULTIPLE ERROR LMS ALGORITHM AND
ITS APPLICATION TOTHE ACTIVE CONTROL OF SOUND AND
VIBRATION ”,IEEE TRANS.ACOUST.,SPEECH,SIGNAL PRO
CESSING,VOL.ASSP −35,PP.1423−1434,1987 にも述べ
られている。[0003] As a method of generating a control sound emitted from a loudspeaker, PROCEEDINGS OF THE IEEE, VOL.
63 PAGE 1692, 1975, “ADAPTIVE NOISE CANSELLATION:
PRINCIPLES AND APPLICATIONS ”
DROW LMS 'algorithm is applied to multiple channels. The contents are described in a paper by the inventor of the above patent, “A MULTIPLE ERROR LMS ALGORITHM AND
ITS APPLICATION TOTHE ACTIVE CONTROL OF SOUND AND
VIBRATION ”, IEEE TRANS.ACOUST., SPEECH, SIGNAL PRO
CESSING, VOL. ASSP-35, PP. 1423-1434, 1987.
【0004】即ち、LMSアルゴリズムは、適応型ディ
ジタルフィルタのフィルタ係数を更新するのに好適なア
ルゴリズムの一つであって、例えば、いわゆるFilt
ered−X LMSアルゴリズムにあっては、ラウド
スピーカからマイクロフォンまでの伝達関数を表すフィ
ルタを、全てのラウドスピーカとマイクロフォンとの組
み合わせについて設定し、騒音源の騒音発生状態を表す
基準信号をそのフィルタで処理した値と、各マイクロフ
ォンが検出した残留騒音とに基づいて、各ラウドスピー
カ毎に設けられた適応型ディジタルフィルタのフィルタ
係数を更新している。That is, the LMS algorithm is one of the suitable algorithms for updating the filter coefficient of the adaptive digital filter.
In the ered-X LMS algorithm, a filter representing a transfer function from a loudspeaker to a microphone is set for all combinations of loudspeakers and microphones, and a reference signal representing a noise generation state of a noise source is set by the filter. Based on the processed value and the residual noise detected by each microphone, the filter coefficient of the adaptive digital filter provided for each loudspeaker is updated.
【0005】ここで、このような能動型騒音制御装置で
は、ラウドスピーカからマイクロフォンまでの伝達関数
を表すフィルタが正確にその伝達関数を表しているとい
うのが前提となっており、フィルタが表す伝達関数と、
実際の物理的な空間の伝達関数との間のズレが大きい
と、騒音低減効果が低下するばかりか、周波数領域で9
0度近い位相差が生じると逆に発散してしまうこともあ
る。Here, in such an active noise control device, it is premised that a filter representing a transfer function from a loudspeaker to a microphone accurately represents the transfer function. Function and
When the deviation from the transfer function of the actual physical space is large, not only does the noise reduction effect decrease, but also the noise reduction effect is reduced by 9% in the frequency domain.
Conversely, when a phase difference close to 0 degree occurs, the light may diverge.
【0006】このような不具合の解決を図る従来の技術
として、特開平3−259722号公報に記載されたも
のがあり、これは、冷蔵庫のコンプレッサで発生し機械
室ダクトを通じて外部に放射される騒音を、その機械室
ダクトから放射される前に打ち消す装置であって、機械
室ダクト内に騒音制御を行うラウドスピーカ及びマイク
ロフォンを備えていて、コンプレッサの駆動状態に応じ
てラウドスピーカから制御音を発生して騒音低減を図る
一方、騒音制御特性が劣化しないように、コンプレッサ
が停止する度に、同定音としてホワイトノイズを発生し
てラウドスピーカ及びマイクロフォン間の伝達関数を測
定し、フィルタの同定を行っている。As a conventional technique for solving such a problem, there is a technique disclosed in Japanese Patent Application Laid-Open No. 3-259722, which discloses a noise generated in a compressor of a refrigerator and radiated outside through a machine room duct. Loudspeaker and a microphone for noise control in the machine room duct, and generate a control sound from the loudspeaker according to the driving state of the compressor. Each time the compressor is stopped, white noise is generated as an identification sound to measure the transfer function between the loudspeaker and the microphone, and the filter is identified, so that the noise control characteristics are not degraded while the compressor is stopped. ing.
【0007】[0007]
【発明が解決しようとする課題】しかしながら、上述し
た冷蔵庫に関する技術を例えば車両等にそのまま適用す
ることを考えた場合、車室内の伝達関数は、温度,湿
度,窓の開閉,乗員数等の種々の要因によって短時間に
且つ大きく変化するため、例えばエンジンを停止する度
にフィルタ係数を更新しても、フィルタが表す伝達関数
と、実際の物理的な空間の伝達関数との間のズレが時間
の経過とともに大きくなってしまい、良好な騒音制御が
行えない。つまり、車両等のように伝達関数の変動が激
しい場合には、騒音制御と並行して伝達関数を測定する
ことにより、フィルタを常時同定することが望ましい。However, when the technology relating to the refrigerator described above is considered to be applied to a vehicle or the like as it is, the transfer function in the passenger compartment is varied depending on the temperature, humidity, opening and closing of windows, number of occupants, and the like. For example, even if the filter coefficient is updated every time the engine is stopped, the difference between the transfer function represented by the filter and the transfer function in the actual physical space is time-dependent. The noise increases with the lapse of time, and good noise control cannot be performed. That is, when the transfer function varies greatly, such as in a vehicle, it is desirable to always identify the filter by measuring the transfer function in parallel with the noise control.
【0008】また、上記従来の技術では、伝達関数の測
定にホワイトノイズを用いているため、ホワイトノイズ
を発生させるための演算処理若しくはホワイトノイズを
記憶しておく記憶領域が必要であり、しかもホワイトノ
イズに基づいて伝達関数を測定する場合には、高速フー
リエ変換(FFT)や逆フーリエ変換(IFFT)を実
行する必要があるばかりか、同定用フィルタとホワイト
ノイズとの畳み込み演算等が必須となるため、伝達関数
の常時測定を可能とするのに高速演算が可能な高価な演
算素子が必要となり、コストを引き上げる要因となって
しまう。Further, in the above-mentioned conventional technique, since white noise is used for measuring a transfer function, an arithmetic processing for generating white noise or a storage area for storing white noise is required. When a transfer function is measured based on noise, it is necessary not only to execute a fast Fourier transform (FFT) or an inverse Fourier transform (IFFT), but also to perform a convolution operation between an identification filter and white noise. Therefore, an expensive operation element capable of high-speed operation is required to enable constant measurement of the transfer function, which causes a cost increase.
【0009】本発明は、このような従来の技術が有する
未解決の課題に着目してなされたものであって、制御音
を発生するラウドスピーカと残留騒音を検出するマイク
ロフォンとの間の伝達関数を少ない演算量で測定可能と
することにより、伝達関数が短時間に且つ大幅に変動す
る場合であっても制御特性の劣化を防止できる能動型騒
音制御装置を提供することを目的としている。The present invention has been made in view of such an unsolved problem of the prior art, and has a transfer function between a loudspeaker for generating a control sound and a microphone for detecting a residual noise. It is an object of the present invention to provide an active noise control device capable of preventing deterioration of control characteristics even when the transfer function fluctuates in a short time and greatly by making it possible to measure the noise with a small amount of calculation.
【0010】[0010]
【課題を解決するための手段】上記目的を達成するため
に、請求項1記載の発明は、騒音源から騒音が伝達され
る空間に制御音を発生可能な制御音源と、前記騒音源の
騒音発生状態を検出し基準信号として出力する騒音発生
状態検出手段と、前記空間内の所定位置における残留騒
音を検出し残留騒音信号として出力する残留騒音検出手
段と、フィルタ係数可変の適応ディジタルフィルタと、
前記基準信号を前記適応ディジタルフィルタでフィルタ
処理して前記制御音源の駆動信号を生成する駆動信号生
成手段と、前記制御音源及び前記残留騒音検出手段間の
伝達関数をモデル化した伝達関数フィルタと、前記基準
信号を前記伝達関数フィルタでフィルタ処理した結果と
前記残留騒音信号とに基づいて前記空間内の騒音が低減
するように前記適応ディジタルフィルタのフィルタ係数
を更新する適応処理手段と、を備えた能動型騒音制御装
置において、前記基準信号に同期したパルス信号を一時
的に前記駆動信号に加算可能な信号加算手段と、フィル
タ係数可変で且つ前記伝達関数一つに対して二つの同定
用フィルタと、前記信号加算手段が前記加算を行った際
の前記残留騒音信号に基づいて前記一方の同定用フィル
タを設定し且つ前記信号加算手段が前記加算を行わなか
った際の前記残留騒音信号に基づいて前記他方の同定用
フィルタを設定する同定用フィルタ設定手段と、前記一
方の同定用フィルタのフィルタ係数及び前記他方の同定
用フィルタのフィルタ係数の差に基づいて前記伝達関数
フィルタを設定する伝達関数フィルタ設定手段と、を備
えた。According to one aspect of the present invention, there is provided a control sound source capable of generating a control sound in a space where noise is transmitted from a noise source, and a noise source of the noise source. A noise generation state detection means for detecting a generation state and outputting it as a reference signal, a residual noise detection means for detecting residual noise at a predetermined position in the space and outputting it as a residual noise signal, and an adaptive digital filter having a variable filter coefficient;
A drive signal generating unit that generates a drive signal of the control sound source by filtering the reference signal with the adaptive digital filter, a transfer function filter that models a transfer function between the control sound source and the residual noise detection unit, Adaptive processing means for updating a filter coefficient of the adaptive digital filter so as to reduce noise in the space based on a result of filtering the reference signal by the transfer function filter and the residual noise signal. In the active noise control device, a signal addition unit that can temporarily add a pulse signal synchronized with the reference signal to the drive signal, and two identification filters that have variable filter coefficients and one transfer function. And setting the one identification filter based on the residual noise signal when the signal addition means has performed the addition, and Identification filter setting means for setting the other identification filter based on the residual noise signal when the signal addition means does not perform the addition, and a filter coefficient of the one identification filter and a filter coefficient of the other identification filter Transfer function filter setting means for setting the transfer function filter based on a difference between filter coefficients of the filter.
【0011】また、上記目的を達成するために、請求項
2記載の発明は、騒音源から騒音が伝達される空間に制
御音を発生可能な制御音源と、前記騒音源の騒音発生状
態を検出し基準信号として出力する騒音発生状態検出手
段と、前記空間内の所定位置における残留騒音を検出し
残留騒音信号として出力する残留騒音検出手段と、フィ
ルタ係数可変の適応ディジタルフィルタと、前記基準信
号を前記適応ディジタルフィルタでフィルタ処理して前
記制御音源の駆動信号を生成する駆動信号生成手段と、
前記制御音源及び前記残留騒音検出手段間の伝達関数を
モデル化した伝達関数フィルタと、前記基準信号を前記
伝達関数フィルタでフィルタ処理した結果と前記残留騒
音信号とに基づいて前記空間内の騒音が低減するように
前記適応ディジタルフィルタのフィルタ係数を更新する
適応処理手段と、を備えた能動型騒音制御装置におい
て、前記基準信号に同期したパルス信号を一時的に前記
駆動信号に加算可能な信号加算手段と、フィルタ係数可
変で且つ前記伝達関数一つに対して二つの同定用フィル
タと、前記信号加算手段が前記加算を行った際の前記残
留騒音信号の同期をとった平均値を前記一方の同定用フ
ィルタのフィルタ係数とし且つ前記信号加算手段が前記
加算を行わなかった際の前記残留騒音信号の同期をとっ
た平均値を前記他方の同定用フィルタのフィルタ係数と
する同定用フィルタ設定手段と、前記一方の同定用フィ
ルタのフィルタ係数及び前記他方の同定用フィルタのフ
ィルタ係数の差に基づいて前記伝達関数フィルタを設定
する伝達関数フィルタ設定手段と、を設けた。According to another aspect of the present invention, a control sound source capable of generating a control sound in a space to which noise is transmitted from a noise source and a noise generation state of the noise source are detected. A noise generation state detecting means for outputting a reference signal, a residual noise detecting means for detecting a residual noise at a predetermined position in the space and outputting the residual noise signal, an adaptive digital filter having a variable filter coefficient, and the reference signal. Drive signal generating means for generating a drive signal of the control sound source by performing a filtering process with the adaptive digital filter;
A transfer function filter that models a transfer function between the control sound source and the residual noise detection unit; and a noise in the space based on a result of filtering the reference signal with the transfer function filter and the residual noise signal. An adaptive processing means for updating a filter coefficient of the adaptive digital filter so as to reduce the number of the adaptive digital filters, wherein a signal addition capable of temporarily adding a pulse signal synchronized with the reference signal to the drive signal is provided. Means, a filter coefficient variable and two identification filters for one of the transfer functions, and an average value obtained by synchronizing the residual noise signal when the signal adding means performs the addition with the one The filter average of the remaining noise signal when the signal addition means does not perform the addition is used as the filter coefficient of the identification filter of Identification filter setting means as a filter coefficient of the identification filter; and transfer function filter setting for setting the transfer function filter based on a difference between a filter coefficient of the one identification filter and a filter coefficient of the other identification filter. Means.
【0012】[0012]
【作用】請求項1又は請求項2記載の発明にあっては、
駆動信号生成手段が、騒音源の騒音発生状態を表す基準
信号を適応ディジタルフィルタでフィルタ処理して制御
音源の駆動信号を生成するから、この駆動信号によって
制御音源が駆動されると、制御音源からは、空間内に伝
達される騒音に相関のある制御音が発生するが、制御開
始直後は、適応ディジタルフィルタのフィルタ係数が最
適値に収束しているとは限らないので、必ずしも制御音
と騒音とが干渉し合わないから空間内の騒音が低減され
るとはいえない。According to the invention described in claim 1 or 2,
The drive signal generating means generates a drive signal for the control sound source by filtering the reference signal representing the noise generation state of the noise source with an adaptive digital filter. Generates a control sound that is correlated with the noise transmitted to the space.However, immediately after the control starts, the filter coefficient of the adaptive digital filter does not always converge to the optimum value. Do not interfere with each other, it cannot be said that noise in the space is reduced.
【0013】しかし、適応処理手段が、基準信号を伝達
関数フィルタでフィルタ処理した結果と残留騒音信号と
に基づいて空間内の騒音が低減するように適応ディジタ
ルフィルタのフィルタ係数を更新すると、制御が進行す
るに従って適応ディジタルフィルタのフィルタ係数は最
適値に向かって収束していくから、制御音と騒音とが干
渉し合って空間内の騒音が低減されるようになる。However, when the adaptive processing means updates the filter coefficients of the adaptive digital filter based on the result of filtering the reference signal with the transfer function filter and the residual noise signal so as to reduce the noise in the space, the control is performed. As the filter coefficient of the adaptive digital filter converges toward the optimum value as it progresses, the control sound and the noise interfere with each other to reduce the noise in the space.
【0014】また、信号加算手段が、基準信号に同期し
たパルス信号を一時的に駆動信号に加算するため、その
加算が行われた際には、制御音源の駆動信号は、基準信
号を適応ディジタルフィルタでフィルタ処理した結果と
基準信号に同期したパルス信号とが重畳された信号とな
り、その加算が行われなかった際には、制御音源の駆動
信号は基準信号を適応ディジタルフィルタでフィルタ処
理した結果そのものとなる。Further, since the signal adding means temporarily adds a pulse signal synchronized with the reference signal to the drive signal, when the addition is performed, the drive signal of the control sound source is obtained by converting the reference signal to an adaptive digital signal. When the result of filtering by the filter and the pulse signal synchronized with the reference signal are superimposed, and the addition is not performed, the drive signal of the control sound source is the result of filtering the reference signal by the adaptive digital filter. It becomes itself.
【0015】そして、請求項1記載の発明にあっては、
同定用フィルタ設定手段は、上記加算が行われた際に
は、一方の同定用フィルタC^1 を残留騒音信号に基づ
いて設定し、上記加算が行われなかった際には、他方の
同定用フィルタC^0 を残留騒音信号に基づいて設定す
る。すると、一方の同定用フィルタC^1 と他方の同定
用フィルタC^0 とでは、パルス信号の発生の分だけ差
が生じることになる。そして、残留騒音信号に含まれて
いるパルス信号に起因する成分は、そのパルス信号の大
きさを“1”と考えれば、制御音源及び残留騒音検出手
段間の伝達関数のインパルス応答である。[0015] According to the first aspect of the present invention,
The identification filter setting means sets one identification filter C ^ 1 based on the residual noise signal when the addition is performed, and sets the other identification filter C ^ 1 when the addition is not performed. The filter C ^ 0 is set based on the residual noise signal. Then, a difference is generated between the one identification filter C # 1 and the other identification filter C # 0 by an amount corresponding to the generation of the pulse signal. The component caused by the pulse signal included in the residual noise signal is an impulse response of a transfer function between the control sound source and the residual noise detecting means, assuming that the magnitude of the pulse signal is “1”.
【0016】従って、伝達関数フィルタ設定手段におい
て、同定用フィルタC^1 のフィルタ係数及び同定用フ
ィルタC^0 のフィルタ係数の差に基づけば、伝達関数
フィルタが設定されることになる。一方、請求項2記載
の発明にあっては、同定用フィルタ設定手段は、上記加
算が行われた際には、一方の同定用フィルタC^1 のフ
ィルタ係数を残留騒音信号の同期をとった平均値(時間
的な起点をそろえた平均値)とし、上記加算が行われな
かった際には、他方の同定用フィルタC^0 のフィルタ
係数を残留騒音信号の同期をとった平均値とする。な
お、これら平均値の演算方法は特に限定されるものでは
なく、所定時間内の累計から平均値を計算することとし
てもよいし、移動平均を求めるようにしてもよい。ま
た、LMSアルゴリズムを利用して同期をとった平均値
に収束するようなアルゴリズムとしてもよい。Therefore, the transfer function filter is set by the transfer function filter setting means based on the difference between the filter coefficient of the identification filter C # 1 and the filter coefficient of the identification filter C # 0 . On the other hand, in the invention according to claim 2, the identification filter setting means synchronizes the filter coefficient of one identification filter C ^ 1 with the residual noise signal when the addition is performed. When the above addition is not performed, the filter coefficient of the other identification filter C ^ 0 is set to an average value synchronized with the residual noise signal when the addition is not performed. . The method of calculating these average values is not particularly limited, and the average value may be calculated from the cumulative total within a predetermined time, or the moving average may be obtained. Alternatively, an algorithm that converges to a synchronized average value using the LMS algorithm may be used.
【0017】ここで、騒音源で発生した騒音をX、基準
信号をx、基準信号に同期したパルス信号をx' 、残留
騒音をE、騒音源と残留騒音検出手段との間の伝達関数
をG、適応ディジタルフィルタをW、制御音源と残留騒
音検出手段との間の伝達関数をCとすれば、上記加算が
行われ制御音源の駆動信号が上記重畳された信号である
場合には、残留騒音Eは下記の(1)式のように表され
る。Here, the noise generated by the noise source is X, the reference signal is x, the pulse signal synchronized with the reference signal is x ', the residual noise is E, and the transfer function between the noise source and the residual noise detecting means is X. G, the adaptive digital filter is W, and the transfer function between the control sound source and the residual noise detection means is C. If the addition is performed and the drive signal of the control sound source is the superimposed signal, the residual The noise E is represented by the following equation (1).
【0018】 E=GX+CWx+Cx' ……(1) 一方、上記加算が行われず駆動信号が基準信号を適応デ
ィジタルフィルタでフィルタ処理した結果そのものであ
る場合には、パルス信号x' が発生していないから、残
留騒音Eは下記の(2)式のように表される。 E=GX+CWx ……(2) そして、同定用フィルタC^1 の各フィルタ係数は上記
(1)式で表される残留騒音Eの同期をとった平均値で
あり、同定用フィルタC^0 の各フィルタ係数は上記
(2)式で表される残留騒音Eの同期をとった平均値で
あるが、いずれも瞬時値ではなく時間的な起点をそろえ
た上での平均値であることから、同定用フィルタC^1
を求めるための残留騒音信号と、同定用フィルタC^0
を求めるための残留騒音信号とが異なる時刻の残留騒音
に基づくものであっても、同定用フィルタC^1 のフィ
ルタ係数と同定用フィルタC^0 のフィルタ係数との差
を求めると、上記(1)式の右辺第1項,第2項と上記
(2)式の右辺第1項,第2項とが相殺され、下記の
(3)式のような結果が得られる。E = GX + CWx + Cx ′ (1) On the other hand, if the addition is not performed and the drive signal is the result of filtering the reference signal with the adaptive digital filter, the pulse signal x ′ is not generated. , The residual noise E is expressed by the following equation (2). E = GX + CWx ...... (2 ) Then, the filter coefficient of the identifying filter C ^ 1 is the average value taken synchronization residual noise E represented by the above formula (1), of the identifying filter C ^ 0 Each filter coefficient is an average value obtained by synchronizing the residual noise E expressed by the above equation (2). However, since each of the filter coefficients is not an instantaneous value but an average value obtained by adjusting the starting points in time, Identification filter C ^ 1
And a noise filter C た め0
If the difference between the filter coefficient of the identification filter C # 1 and the filter coefficient of the identification filter C # 0 is determined even if the residual noise signal for calculating The first and second terms on the right side of the equation (1) and the first and second terms on the right side of the equation (2) are canceled, and a result as shown in the following equation (3) is obtained.
【0019】 C^1 −C^0 =Cx' ……(3) つまり、伝達関数フィルタ設定手段が、同定用フィルタ
設定手段によって設定された同定用フィルタC^1 のフ
ィルタ係数及び同定用フィルタC^0 のフィルタ係数の
差を求めると、上記(3)式からも明らかなように、実
際の物理的空間の制御音源及び残留騒音検出手段間の伝
達関数Cに、基準信号xに同期したパルス信号x' を印
加した応答Cx' が得られることになる。C ^ 1 −C ^ 0 = Cx ′ (3) That is, the transfer function filter setting means sets the filter coefficient of the identification filter C ^ 1 and the identification filter C set by the identification filter setting means. When the difference between the filter coefficients of ^ 0 is obtained, as is apparent from the above equation (3), the transfer function C between the control sound source and the residual noise detection means in the actual physical space is added to the pulse synchronized with the reference signal x. A response Cx 'to which the signal x' has been applied is obtained.
【0020】従って、例えばパルス信号x' を大きさ1
のパルス信号(インパルス信号)とすれば、上記(3)
式の結果は伝達関数Cのインパルス応答を表すことにな
るから、伝達関数フィルタ設定手段によって簡易且つ正
確に伝達関数フィルタが設定されるようになる。Therefore, for example, when the pulse signal x ′ is
If the pulse signal (impulse signal) is (3)
Since the result of the expression represents the impulse response of the transfer function C, the transfer function filter can be easily and accurately set by the transfer function filter setting means.
【0021】[0021]
【実施例】以下、この発明の実施例を図面に基づいて説
明する。図1は、本発明の第1実施例の全体構成を示す
図であり、この実施例は、騒音源としてのエンジン4か
ら空間としての車室6内に伝達されるこもり音の低減を
図る車両用能動型騒音制御装置1に本発明を適用したも
のである。Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a diagram showing an overall configuration of a first embodiment of the present invention. This embodiment is directed to a vehicle for reducing a muffled sound transmitted from an engine 4 as a noise source to a vehicle interior 6 as a space. The present invention is applied to an active noise control device 1 for use.
【0022】先ず、構成を説明すると、車体3は、前輪
2a,2b,後輪2c,2d及び各車輪2a〜2dと車
体3との間に介在するサスペンションによって支持され
ている。なお、図1に示す車両は、前輪2a及び2bが
車体3前部に配置されたエンジン4によって回転駆動さ
れるいわゆる前置きエンジン前輪駆動車である。エンジ
ン4には、クランク角センサ5が取り付けられていて、
このクランク角センサ5は、エンジン4のクランク角の
回転に同期した(例えばクランクが1°回転する度に一
つのパルス信号でなる)クランク角信号Xをコントロー
ラ10に供給する。First, the construction will be described. The vehicle body 3 is supported by suspensions interposed between the vehicle body 3 and the front wheels 2a, 2b, the rear wheels 2c, 2d and the wheels 2a to 2d. The vehicle shown in FIG. 1 is a so-called front-engine front-wheel drive vehicle in which front wheels 2a and 2b are rotationally driven by an engine 4 arranged in the front part of the vehicle body 3. A crank angle sensor 5 is attached to the engine 4,
The crank angle sensor 5 supplies the controller 10 with a crank angle signal X synchronized with the rotation of the crank angle of the engine 4 (for example, one pulse signal every time the crank rotates 1 °).
【0023】また、車体3の車室6内には、制御音源と
してのラウドスピーカ7a,7b,7c及び7dが、前
部座席S1 ,S2 及び後部座席S3 ,S4 のそれぞれに
対向するドア部に配置されている。さらに、各座席S1
〜S4 のヘッドレスト位置には、残留騒音検出手段とし
てのマイクロフォン8a〜8hが、それぞれ二つずつ配
設されていて、これらマイクロフォン8a〜8hが音圧
として測定した残留騒音信号e1 〜e8 が、コントロー
ラ10に供給される。In the cabin 6 of the vehicle body 3, loudspeakers 7a, 7b, 7c and 7d as control sound sources face the front seats S 1 and S 2 and the rear seats S 3 and S 4 respectively. It is arranged in the door part. In addition, each seat S 1
The head restraint position of the to S 4 is a microphone 8a~8h as residual noise detecting means, have been respectively disposed two by two, the residual noise signal e 1 to e these microphones 8a~8h was measured as a sound pressure 8 Is supplied to the controller 10.
【0024】そして、コントローラ10は、マイクロコ
ンピュータや必要なインタフェース回路等を含んで構成
されていて、クランク角センサ5から供給されるクラン
ク角信号Xと、マイクロフォン8a〜8hから供給され
る残留騒音信号e1 〜e8 とに基づいて、後述する演算
処理を実行し、車室6内に伝達されるこもり音を打ち消
すような制御音がラウドスピーカ7a〜7dから発せら
れるように、それらラウドスピーカ7a〜7dに駆動信
号y1 〜y4 を出力する。The controller 10 includes a microcomputer, necessary interface circuits, and the like, and includes a crank angle signal X supplied from the crank angle sensor 5 and a residual noise signal supplied from the microphones 8a to 8h. Based on e 1 to e 8 , the loudspeakers 7a perform arithmetic processing described later so that control sounds that cancel the muffled sound transmitted into the vehicle interior 6 are emitted from the loudspeakers 7a to 7d. and outputs a drive signal y 1 ~y 4 to ~7D.
【0025】図2は、コントローラ10の機能構成を示
すブロック図であって、このコントローラ10は、クラ
ンク角信号Xに基づき、こもり音の原因となるエンジン
4で発生する振動と同じ周期の(例えば、レシプロ4気
筒の場合は、クランクの180度を一周期とした)正弦
波でなる基準信号xを生成し出力する基準信号生成部1
1と、基準信号xに基づきその基準信号xと同じ周期の
パルス信号x' を生成し出力するパルス信号生成部12
と、を有している。FIG. 2 is a block diagram showing a functional configuration of the controller 10. The controller 10 has a cycle (for example, the same cycle as that of vibration generated in the engine 4 causing a muffled sound, based on the crank angle signal X). In the case of a four-cylinder reciprocating cylinder, a reference signal generating unit 1 for generating and outputting a reference signal x formed of a sine wave (180 degrees of the crank is defined as one cycle)
1 and a pulse signal generator 12 for generating and outputting a pulse signal x 'having the same cycle as the reference signal x based on the reference signal x.
And
【0026】また、コントローラ10は、ラウドスピー
カ7a〜7dに対応した個数(M個:本実施例では、M
=4)の適応ディジタルフィルタWm (m=1〜M)
と、その適応ディジタルフィルタWm の各フィルタ係数
Wmi(i=0〜I−1:Iは適応ディジタルフィルタW
m のタップ数(フィルタ係数の個数))と基準信号xと
を畳み込んで駆動信号y1 〜y4 を生成し出力する駆動
信号生成手段としての駆動信号生成部13と、各ラウド
スピーカ7a〜7d及びマイクロフォン8a〜8h間の
伝達関数を有限インパルス応答関数の形でモデル化した
伝達関数フィルタC^lm(l=1〜L:Lはマイクロフ
ォン8a〜8hの個数であり、本実施例ではL=8)
と、その伝達関数フィルタC^lmの各フィルタ係数C^
lmj (j=0〜J−1:Jは伝達関数フィルタC^lmの
タップ数)と基準信号xとを畳み込んで基準処理信号r
lmを生成し出力する基準処理信号生成部14と、基準処
理信号rlm及び残留騒音信号e1 〜e8 に基づいて車室
6内のこもり音が低減するように駆動信号生成部13内
の適応ディジタルフィルタWm の各フィルタ係数Wmiを
更新する適応処理手段としての適応処理部15と、を有
している。The controller 10 has a number corresponding to the loudspeakers 7a to 7d (M: in this embodiment, M
= 4) adaptive digital filter W m (m = 1 to M)
If, each filter coefficient of the adaptive digital filter W m W mi (i = 0~I -1: I adaptive digital filter W
A drive signal generation unit 13 as drive signal generation means for generating and outputting drive signals y 1 to y 4 by convolving the number of taps of m (the number of filter coefficients) and the reference signal x; A transfer function filter C 関 数lm (l = 1 to L: L is the number of the microphones 8 a to 8 h, in which the transfer function between the microphones 7 a and 8 a to 8 h is modeled in the form of a finite impulse response function. = 8)
And each filter coefficient C 係数 of the transfer function filter C ^ lm
lmj (j = 0 to J-1: J is the number of taps of the transfer function filter C ^ lm ) and the reference signal x to convolve the reference processed signal r
a reference processed signal generating unit 14 for generating and lm output, the reference processed signal r lm and muffled sound in the passenger compartment 6 based on the residual noise signal e 1 to e 8 is in the drive signal generating unit 13 so as to reduce An adaptive processing unit 15 as an adaptive processing means for updating each filter coefficient W mi of the adaptive digital filter W m .
【0027】さらに、コントローラ10は、パルス信号
生成部12によって生成されたパルス信号x' を駆動信
号ym に加算可能な加算部16と、パルス信号生成部1
2及び加算部16間に介在し且つパルス信号x' を断続
的に加算部16に供給する断続スイッチ部17と、伝達
関数フィルタC^lm一つに対して二つ設けられ且つパル
ス信号x' が入力された時点から後述するクロックパル
スCP に同期してフィルタ係数を順に一つずつ出力する
フィルタ係数可変の同定用フィルタC^1lm ,C^0lm
と、残留騒音信号el と同定用フィルタC^1lm 若しく
は同定用フィルタC^0lm の出力との差を演算する減算
部18と、同定用フィルタC^1lm ,C^0lm と減算部
18との間に介在し且つ同定用フィルタC^1lm ,C^
0lm の何れか一方の出力を減算部18に供給する切換ス
イッチ部19と、同定用フィルタC^1lm ,C^0lm の
各フィルタ係数C^1lmj,C^0lmjが減算部18の出力
に一致するように各フィルタ係数C^1lmj,C^0lmjを
更新するフィルタ係数更新部20と、フィルタ係数更新
部20及び同定用フィルタC^1lm ,C^0lm 間に介在
し且つ同定用フィルタC^1lm ,C^0lm の何れか一方
のフィルタ係数が更新されるように切り換わる切換スイ
ッチ部21と、同定用フィルタC^1lm ,C^0lm の各
フィルタ係数C^1lmj,C^0lmjの差に基づいて伝達関
数フィルタC^lmを設定する減算部22と、を有してい
る。Furthermore, the controller 10 includes an adder capable of adding unit 16 into a pulse signal x 'the drive signal y m that is generated by the pulse signal generator 12, the pulse signal generator 1
2 and an intermittent switch section 17 interposed between the adder section 16 and intermittently supplying the pulse signal x 'to the adder section 16, and two pulse signal x' provided for one transfer function filter C ^ lm. There filter identifying the filter coefficient variable one by one output in synchronization with the clock pulses C P to be described later, from the time it is entered the filter coefficients in the order C ^ 1lm, C ^ 0lm
When, a subtracting unit 18 for calculating a difference between the output of the residual noise signal e l and identifying filter C ^ 1 lm or identifying filter C ^ 0lm, identifying filter C ^ 1 lm, and a C ^ 0lm subtraction section 18 And an identification filter C {1lm , C }
A change-over switch 19 for supply to the subtraction unit 18 to either the output of 0Lm, identifying filter C ^ 1 lm, each filter coefficient of the C ^ 0lm C ^ 1lmj, C ^ 0lmj matches the output of the subtraction section 18 A filter coefficient updating unit 20 for updating the filter coefficients C ^ 1lmj and C ^ 0lmj as described above , and a filter coefficient updating unit 20 interposed between the filter coefficient updating unit 20 and the identification filters C ^ 1lm , C ^ 0lm and the identification filters C ^ 1lm , The changeover switch unit 21 that switches so that one of the filter coefficients of C ^ 0lm is updated, and the difference between the filter coefficients C ^ 1lmj and C ^ 0lmj of the identification filters C ^ 1lm and C ^ 0lm. And a subtraction unit 22 for setting a transfer function filter C ^ lm .
【0028】ただし、断続スイッチ部17,切換スイッ
チ部19及び切換スイッチ部21は互いに連動してい
て、断続スイッチ部17が接続状態の際には、切換スイ
ッチ19及び21は同定用フィルタC^1lm 側に切り換
わり、断続スイッチ17が遮断状態の際には、切換スイ
ッチ19及び21は同定用フィルタC^0lm 側に切り換
わるようになっている。However, the intermittent switch section 17, the changeover switch section 19 and the changeover switch section 21 are interlocked with each other, and when the intermittent switch section 17 is in the connected state, the changeover switches 19 and 21 are connected to the identification filter C ^ 1lm. Side, and when the intermittent switch 17 is in the cut-off state, the changeover switches 19 and 21 are switched to the identification filter C # 0lm side.
【0029】またさらに、コントローラ10は、クラン
ク角信号Xに基づき、基準信号xのN(Nは整数)倍の
周波数の(例えば、レシプロ4気筒の場合、N=10で
あれば、クランクが18度回転する度に一つのパルスで
なる)クロックパルスCP を生成し出力するクロックパ
ルス生成部23を有しており、コントローラ10内の各
処理は、基本的にはそのクロックパルスCP に同期して
実行される。従って、適応ディジタルフィルタWm ,同
定用フィルタC^1lm ,C^0lm は、こもり音の周期と
同じ周期である基準信号xの一周期当たり、N回の出力
を行うことになる。Further, based on the crank angle signal X, the controller 10 controls the frequency of the reference signal x to be N (N is an integer) times the frequency (for example, in the case of a reciprocating four-cylinder engine, if N = 10, the crank is 18 times). degrees consisting of rotating every a single pulse) has a clock pulse generator 23 to the clock pulse C P to produce an output, each processing in the controller 10 is basically synchronized with the clock pulses C P And executed. Therefore, the adaptive digital filter W m and the identification filters C ^ 1lm and C ^ 0lm output N times per cycle of the reference signal x having the same cycle as the muffled sound.
【0030】ここで、本実施例にあっては、フィルタ係
数更新部15は、適応ディジタルフィルタのフィルタ係
数を更新するのに好適なアルゴリズムの一つであるLM
Sアルゴリズムに基づいて、適応ディジタルフィルタW
m の各フィルタ係数Wmiを更新するが、特に基準信号x
を伝達関数フィルタC^lmでフィルタ処理した値rlmを
用いていることから、Filtered−X LMSア
ルゴリズムが実行されることになり、適応ディジタルフ
ィルタWm のフィルタ係数Wmiの更新式は下記の(4)
式又は(5)式のようになる。Here, in the present embodiment, the filter coefficient updating unit 15 is an LM which is one of the algorithms suitable for updating the filter coefficient of the adaptive digital filter.
Based on the S algorithm, the adaptive digital filter W
m , each filter coefficient W mi is updated, but in particular the reference signal x
Is filtered using a transfer function filter C ^ lm , the filtered r- lm algorithm is executed, and the updating equation of the filter coefficient W mi of the adaptive digital filter W m is (4)
Equation or equation (5) is obtained.
【0031】 λW は発散抑制係数であって、1以下の値を採る。ま
た、αW は収束係数と呼ばれる係数であって、フィルタ
が最適に収束する速度やその安定性に関与する。なお、
(n)が付されている項は、サンプリング時刻nにおけ
る値であることを示し、添字kは、一周期(N回)内の
何番目の処理に対応するかを表す変数であって、0から
N−1までの値を採る。[0031] λ W is a divergence suppression coefficient and takes a value of 1 or less. Α W is a coefficient called a convergence coefficient, and is related to the speed at which the filter converges optimally and its stability. In addition,
The term to which (n) is attached indicates that it is a value at the sampling time n, and the subscript k is a variable representing the number of processing in one cycle (N times), To N-1.
【0032】一方、同定用フィルタC^1lm ,C^0lm
のフィルタ係数を更新するフィルタ係数更新部20も、
本実施例では、フィルタ係数更新部15と同様にLMS
アルゴリズムに基づいて減算部18の出力が零となるよ
うに同定用フィルタC^1lm,C^0lm の各フィルタ係
数C^1lmj,C^0lmjを更新するのであるが、パルス信
号x' がフィルタ係数更新部20に入力されていること
から、いわゆる同期式LMSアルゴリズムが実行される
ことになる。なお、同期式LMSアルゴリズムとは、騒
音の基本周波数に同期したパルス信号を基準信号として
用いたLMSアルゴリズムのことである(“日本音響学
会講演論文集 平成4年3月”515〜516頁に詳し
い。)。On the other hand, the identification filters C ^ 1lm , C ^ 0lm
The filter coefficient updating unit 20 for updating the filter coefficient of
In the present embodiment, the LMS
Identifying filter so that the output of the subtracting unit 18 based on the algorithm becomes zero C ^ 1 lm, each filter coefficient of the C ^ 0lm C ^ 1lmj, although to update the C ^ 0lmj, pulse signal x 'is the filter coefficient Since the data is input to the updating unit 20, a so-called synchronous LMS algorithm is executed. Note that the synchronous LMS algorithm is an LMS algorithm using a pulse signal synchronized with the fundamental frequency of noise as a reference signal (see “Acoustic Society of Japan, March 1994, March 199”, pp. 515-516). .).
【0033】従って、同定用フィルタC^1lm ,C^
0lm の各フィルタ係数C^1lmj,C^ 0lmjの更新式は、
下記の(6)式のようになる。 C^*lmk(n+1)=λC C^*lmk(n)−αC el ……(6) ただし、λC は発散抑制係数、αC は収束係数であり、
*は1又は0である。なお、この(6)式は、同定用フ
ィルタC^1lm ,C^0lm のタップ数Jが、N以下であ
ることを前提としていて、タップ数JがNよりも大きい
場合には、この(6)式に代えて下記の(7)式を用い
ることになる。Therefore, the identification filter C ^1lm, C ^
0lmEach filter coefficient C ^1lmj, C ^ 0lmjThe update formula for
The following equation (6) is obtained. C ^* lmk(N + 1) = λCC ^* lmk(N) -αCel …… (6) where λCIs the divergence suppression coefficient, αCIs the convergence coefficient,
* Is 1 or 0. Note that this equation (6) is
Filter C ^1lm, C ^0lmIs less than or equal to N
The tap number J is larger than N
In this case, the following equation (7) is used instead of equation (6).
Will be.
【0034】 C^*lm(k+Np) (n+1)=λC C^*lm(k+Np) (n)−αC el ……(7) この(7)式を用いる場合、p,kは、加算部16にお
いて駆動信号ym にパルス信号x' が加算された時点か
ら計数したクロックパルスCP の個数をNで割ったとき
の商と余りである。即ち、上記(6)式は、上記(7)
式においてp=0の場合と考えることもできる。[0034] C ^ * lm (k + Np ) (n + 1) = λ C C ^ * lm (k + Np) (n) -α C e l ...... (7) the case of using the equation (7), p , k is the quotient and remainder when the number of clock pulses C P of the pulse signal x 'has counted from the time it is added to the drive signal y m at the addition section 16 is divided by N. That is, the above equation (6) is equivalent to the above (7)
In the equation, it can be considered that p = 0.
【0035】図3はコントローラ10内で実行される処
理の概要を示したフローチャートであり、以下、本実施
例の動作を説明する。即ち、コントローラ10内におけ
る処理は上述したようにクロックパルスCPに同期して
1サンプリングの処理が実行されるようになっていて、
先ず、そのステップ101において、残留騒音信号el
を読み込み、次いでステップ102では、基準信号xを
伝達関数フィルタC^lmでフィルタ処理して基準処理信
号rlmを演算する。FIG. 3 is a flowchart showing the outline of the processing executed in the controller 10. The operation of this embodiment will be described below. That is, processing in the controller 10 in is not so the process of synchronization to one sampling clock pulses C P as described above is executed,
First, in step 101, the residual noise signal e l
Reading, then in step 102, and the reference signal x filtered by a transfer function filter C ^ lm calculates a reference processed signal r lm.
【0036】そして、ステップ103に移行し、上記
(4)又は(5)式に従って、適応ディジタルフィルタ
Wm の各フィルタ係数Wmiを更新する。ステップ103
において全てのフィルタ係数Wmiを更新したら、ステッ
プ104に移行し、フラグF1 が1であるか否かを判定
する。このフラグF1 は、ラウドスピーカ7a〜7dの
駆動信号y1 〜y4 の何れかに、同定信号たるパルス信
号x' を加算する場合(F1 =1)か加算しない場合
(F1 =0)かを表すフラグであり、後述する他の処理
で適宜設定される。[0036] Then, the process proceeds to step 103, in accordance with the above (4) or (5), to update the filter coefficients W mi of the adaptive digital filter W m. Step 103
After updating all the filter coefficients W mi in step, the process proceeds to step 104, and it is determined whether or not the flag F 1 is 1. The flag F 1, if either of the drive signal y 1 ~y 4 loudspeaker 7a to 7d, when adding the serving identification signal pulse signals x 'to (F 1 = 1) or not added (F 1 = 0 ), And is appropriately set in other processing described later.
【0037】ステップ104の判定が「NO」の場合
は、何れの駆動信号ym にも同定信号は加算されていな
いと判断できるから、ステップ105に移行し、同定用
フィルタC^0lm のフィルタ係数C^0lm(k+Np) を上記
(7)式に従って更新する。なお、この場合、フィルタ
係数更新部20に基準信号として入力されるのがパルス
信号x' であるため、全てのフィルタ係数C^0lmjを更
新する必要はなく、k+Np番目の係数のみである。[0037] If the determination in step 104 is "NO", since any identification signal to the drive signal y m can be judged not to be added, the process proceeds to step 105, the filter coefficient of the identifying filter C ^ 0lm C ^ 0lm (k + Np) is updated according to the above equation (7). In this case, since the pulse signal x 'is input to the filter coefficient updating unit 20 as the reference signal, it is not necessary to update all the filter coefficients C ^ 0lmj, and only the k + Np-th coefficient.
【0038】一方、ステップ104の判定が「YES」
の場合は、何れかの駆動信号ym にパルス信号x' が加
算されていると判断し、ステップ106に移行する。こ
のステップ106では、フラグF2 の値が1であるか否
かを判定するが、このフラグF2 は、何れの駆動信号y
1 〜y4 にパルス信号x' が加算されているかを表すフ
ラグであって、本実施例では、四つのラウドスピーカ7
a〜7dを有することから、フラグF2 は1〜4の値を
採る。On the other hand, the determination in step 104 is "YES".
For the either it determines that pulse signal x 'is added to the drive signal y m, the process proceeds to step 106. In step 106, the value of the flag F 2 is determines whether or not 1, the flag F 2 is one of the drive signal y
A flag 1 ~y 4 to the pulse signal x 'represents whether the summed, in this embodiment, four loudspeakers 7
since it has A~7d, flag F 2 takes a value of 1-4.
【0039】そこで、ステップ106,107,108
においてフラグF2 が1,2,3であるか否かを判定
し、ステップ106の判定が「YES」の場合(F2 =
1)にはステップ109に移行し、ステップ107の判
定が「YES」の場合(F2 =2)にはステップ110
に移行し、ステップ108の判定が「YES」の場合
(F2 =3)にはステップ111に移行し、ステップ1
08の判定が「NO」の場合(F2 =4)にはステップ
112に移行する。Therefore, steps 106, 107, 108
If the flag F 2 is equal to or 1, 2, 3, the determination in step 106 of "YES" in (F 2 =
In step 1), the process proceeds to step 109. If the determination in step 107 is “YES” (F 2 = 2), step 110 is executed.
When the determination in step 108 is “YES” (F 2 = 3), the process proceeds to step 111,
If the determination at 08 is “NO” (F 2 = 4), the process proceeds to step 112.
【0040】そして、ステップ109〜112の各処理
では、上記(7)式に従って、同定用フィルタC
^1lm ,C^0lm の各フィルタ係数C^1lm(k+Np) を更
新するが、ステップ109ではm=1、ステップ110
ではm=2、ステップ111ではm=3,ステップ11
2ではm=4として上記(7)式の更新を行う。これ
は、同定用フィルタC^1lm ,C^0lm は、伝達関数フ
ィルタC^lmを同定するために更新しているのである
が、ラウドスピーカ7a〜7d及びマイクロフォン8a
〜8h間の伝達関数を測定するには、複数のラウドスピ
ーカ7a〜7dから同定音を発生させてしまうと一つの
マイクロフォン8a〜8hに複数のラウドスピーカ7a
〜7dから発せられた同定音が同時に到達してしまい正
確な伝達関数が測定できなくなるため、後述のように一
度には一つのラウドスピーカ7a〜7dからのみ同定音
を発生させるようにしていることと対応させるためであ
る。In each of the processes of steps 109 to 112, the identification filter C is calculated according to the above equation (7).
The filter coefficients C ^ 1lm (k + Np) of ^ 1lm and C ^ 0lm are updated. In step 109, m = 1 and step 110
M = 2, step 111, m = 3, step 11
In 2, the above equation (7) is updated with m = 4. This is because the identification filters C ^ 1lm and C ^ 0lm are updated in order to identify the transfer function filter C 、 lm , but the loudspeakers 7a to 7d and the microphone 8a
In order to measure the transfer function between the loudspeakers 7a to 7h, if the identification sounds are generated from the loudspeakers 7a to 7d, the loudspeakers 7a
Since the identification sounds emitted from the speakers 7 to 7d arrive at the same time and an accurate transfer function cannot be measured, the identification sound is generated only from one loudspeaker 7a to 7d at a time as described later. This is for correspondence.
【0041】ステップ105の更新処理又はステップ1
09〜112の何れかの更新処理を終えたら、ステップ
113に移行し、同定用フィルタC^1lm の各フィルタ
係数C^1lmjから同定用フィルタC^0lm の各フィルタ
係数C^0lmjを差し引いた結果に基づいて、伝達関数フ
ィルタC^lmを設定する。なお、このような演算に基づ
いて伝達関数フィルタC^lmを設定することができる理
由については、後述する。Update processing of step 105 or step 1
After completing any of the update processing of 09 to 112, and proceeds to step 113, by subtracting the respective filter coefficients C ^ 0lmj of identifying filters C ^ 0lm from each filter coefficient C ^ 1lmj of identifying filters C ^ 1 lm results , A transfer function filter C ^ lm is set. The reason why the transfer function filter C ^ lm can be set based on such an operation will be described later.
【0042】ステップ113で伝達関数フィルタC^lm
の設定を行ったら、ステップ114に移行し、全ての駆
動信号y1 〜y4 を零クリアする。そして、ステップ1
15に移行し、変数kがNに達しているか否かを判定
し、未だ達していないと判定された場合は、ステップ1
16に移行して変数kをインクリメントする。一方、ス
テップ115の判定が「YES」の場合は、クロックパ
ルスCP が発生する度にステップ116で一つずつイン
クリメントされる変数がNに達したということであり、
換言すれば、現時点は、基準信号xの一周期が経過し新
たな周期に入る時点であるということである。In step 113, the transfer function filter C ^ lm
After making the settings, the process proceeds to step 114, to zero clear all of the drive signal y 1 ~y 4. And step 1
Then, it is determined whether or not the variable k has reached N. If it is determined that the variable k has not yet reached N, step 1 is executed.
The process proceeds to 16 to increment the variable k. On the other hand, if the determination in step 115 is “YES”, it means that the variable incremented by one in step 116 reaches N each time the clock pulse CP is generated.
In other words, the current time point is a point in time when one cycle of the reference signal x has passed and a new cycle has started.
【0043】そこで、ステップ117に移行し、一周期
内の処理の回数を計数する変数kを零クリアし、ステッ
プ118に移行して、変数pが定数P(Pは整数)に達
しているか否かを判定する。ここで、変数pは、上述し
たように駆動信号ym にパルス信号x' が加算された時
点から起算したクロックパルスCP の個数をNで割った
ときの商であり、より具体的には、パルス信号x' が駆
動信号ym に加算された時点から、基準信号xの周期が
いくつ過ぎたかを表している。そして、定数Pは、伝達
関数フィルタC^lmのタップ数Jを、周期Nの何倍にす
るかを表す整数であって、車室内の残響時間等を考慮し
予め設定される。例えば、タップ数Jを周期Nの2.5倍
にするには、定数Pを3に設定する。Then, the process proceeds to step 117 to clear the variable k for counting the number of processes in one cycle to zero, and then proceeds to step 118 to determine whether the variable p has reached a constant P (P is an integer). Is determined. Here, the variable p is the quotient when the number of clock pulses C P which is counted from the time when the pulse signal x 'to the drive signal y m is added as described above divided by N, and more specifically, from the time the pulse signal x 'is added to the drive signal y m, which indicates whether the period of the reference signal x has passed many. The constant P is an integer indicating how many times the cycle number N of the tap J of the transfer function filter C ^ lm is set, and is set in advance in consideration of the reverberation time in the vehicle interior. For example, to make the number of taps J 2.5 times the period N, the constant P is set to 3.
【0044】ステップ118の判定が「NO」の場合
は、パルス信号x' を何れかの駆動信号ym に加算し同
定音を発生させた時点から未だタップ数Jに対応する時
間が経過していないと判断できる、より具体的には、ス
テップ105,109〜112における同定用フィルタ
の更新処理が完了していないと判断できるから、新たな
同定音を発生させてしまうと正確な伝達関数の同定が妨
げられると考え、ステップ119に移行して変数pをイ
ンクリメントする。[0044] If step 118 is "NO", yet elapsed time corresponding to the pulse signal x 'to one of the drive signal y m still number of taps from the time that caused the identified sound is added to J is More specifically, since it can be determined that the updating process of the identification filter in steps 105 and 109 to 112 has not been completed, if a new identification sound is generated, accurate identification of the transfer function will be performed. Is determined to be disturbed, the process proceeds to step 119, and the variable p is incremented.
【0045】しかし、ステップ118の判定が「YE
S」の場合は、ステップ105,109〜112におけ
る同定用フィルタの更新処理が完了したと判断できるか
ら、ステップ120に移行して変数pを零クリアした後
に、ステップ121以降の処理を実行する。ステップ1
21では、フラグF1 が1であるか否かを判定する。こ
こで、フラグF1 が1でないと判定された場合は、前回
はステップ105の処理が実行されて同定用フィルタC
^0lm のフィルタ係数C^0lmjが更新されたと判断でき
るから、次回の処理では、同定用フィルタC^1lm のう
ち、m=1の同定用フィルタC^1l1 のフィルタ係数C
^1l1jを更新するようにフラグ設定等を行う。具体的に
は、ステップ122に移行してフラグF1 を1に設定す
るとともに、ステップ123に移行してフラグF2 を1
に設定し、そして、ステップ124に移行してm=1の
駆動信号y1 を所定値bに設定する。この所定値bは、
駆動信号ym に大きさ1のインパルスが重畳されたとみ
なされる定数であって、図2では、パルス信号x' に対
応するものである。However, the determination at step 118 is "YE
In the case of "S", it can be determined that the update processing of the identification filter in steps 105 and 109 to 112 has been completed, so the processing shifts to step 120 to clear the variable p to zero, and then executes the processing of step 121 and thereafter. Step 1
In 21, the flag F 1 is determined whether or not 1. Here, if it is determined that the flag F 1 is not 1, the process of step 105 is executed last time, and the identification filter C
Since ^ filter coefficient C ^ 0Lmj of 0lm it can be determined to have been updated, in the next processing, of identifying filters C ^ 1 lm, filter coefficients of the identification filter C ^ 1L1 of m = 1 C
フ ラ グ Set a flag to update 1l1j . Specifically, sets the flag F 1 to 1 proceeds to step 122, the flag F 2 proceeds to step 123 1
Then, the process proceeds to step 124 to set the drive signal y 1 of m = 1 to a predetermined value b. This predetermined value b is
This is a constant that is regarded as a superimposition of an impulse of magnitude 1 on the drive signal y m , and in FIG. 2 corresponds to the pulse signal x ′.
【0046】ステップ121でフラグF1 が1であると
判定された場合は、前回はステップ109〜112の何
れかの更新処理が実行されたと判断できるから、ステッ
プ125に移行し、フラグF2 が1であるか否かを判定
する。そして、ステップ125の判定が「NO」の場合
にはステップ126に移行してフラグF2 が2であるか
否かを判定し、ステップ126の判定が「NO」の場合
にはステップ127に移行してフラグF2 が3であるか
否かを判定する。[0046] When the flag F 1 in step 121 is determined to be 1, since the last time can be determined with any of the update process of step 109 to 112 is executed, the process proceeds to step 125, the flag F 2 is It is determined whether it is 1 or not. When the determination in step 125 is “NO”, the process proceeds to step 126 to determine whether or not the flag F 2 is 2, and when the determination in step 126 is “NO”, the process proceeds to step 127. it is determined whether the flag F 2 is 3.
【0047】即ち、ステップ125の判定が「YES」
の場合には、前回はステップ109の処理が実行されて
同定用フィルタC^1l1 のフィルタ係数C^1l1jが更新
されたと判断できるから、次回の処理では、m=2に対
応する更新処理が実行されるように、ステップ128に
移行してフラグF2 を2に設定し、駆動信号y2 を所定
値bに設定する。That is, the determination in step 125 is “YES”
In the case of, since the last time can be judged to processing in step 109 the filter coefficients C ^ 1l1j of identifying filters C ^ 1L1 is executed is updated, in the next process, the update process corresponding to m = 2 is performed as is, sets the flag F 2 to 2 proceeds to step 128, sets the driving signals y 2 to a predetermined value b.
【0048】そして、ステップ126の判定が「YE
S」の場合には、前回はステップ110においてm=2
に対応する更新処理が行われ、ステップ127の判定が
「YES」の場合には、前回はステップ111において
m=3に対応する更新処理が行われたと判断できるか
ら、それぞれ、次回はm=3,m=4に対応する更新処
理が実行されるように、ステップ130,131若しく
はステップ132,133の処理を実行する。Then, the determination in step 126 is "YE
In the case of "S", m = 2 in the previous step in step 110
Is performed, and if the determination in step 127 is "YES", it can be determined that the update processing corresponding to m = 3 was performed in the previous step 111, so that the next time m = 3 , M = 4, the processes of steps 130 and 131 or steps 132 and 133 are executed.
【0049】なお、ステップ127の判定が「NO」の
場合は、前回はステップ112においてm=4に対応す
る更新処理が行われたと判断できるから、m=1〜4に
対応する一巡の更新処理が完了したと判断し、次回はス
テップ105における更新処理が行われるように、ステ
ップ134に移行してフラグF1 を0に設定する。そし
て、ステップ116,119,124,129,13
1,133又は134の処理を実行したら、ステップ1
35に移行し、下記の(8)式に従って駆動信号ym を
演算する。If the determination in step 127 is "NO", it can be determined that the update processing corresponding to m = 4 was performed in step 112 last time. There was judged complete, the next time, as update processing at the step 105 is performed, the flag F 1 proceeds to step 134 to set to zero. Then, steps 116, 119, 124, 129, 13
After executing the processing of 1, 133 or 134, step 1
The process proceeds to 35, where the drive signal y m is calculated according to the following equation (8).
【0050】 ym =ym +x*Wm ……(8) なお、上記(8)式の右辺の「*」は、畳み込み積分を
表す。ここで、上記(8)式の右辺第1項のym は、ス
テップ114が実行される結果、ステップ115,11
8又は127の判定が「NO」の場合には、全ての駆動
信号ym において0である。従って、各駆動信号y1 〜
y4 は、何れも基準信号xと適応ディジタルフィルタW
m の各フィルタ係数Wmiとを畳み込んだ値であり、図2
で説明すれば、駆動信号生成部13の出力が、そのまま
ラウドスピーカ7a〜7dの駆動信号ym となる。Y m = y m + x * W m (8) Note that “*” on the right side of the above equation (8) represents a convolution integral. Here, the value y m of the first term on the right side of the above equation (8) is obtained by executing steps 114 and 115 and 11
8 or 127 is determined in the case of "NO" is 0 for all of the drive signal y m. Therefore, each of the drive signals y 1 to y 1
y 4 is the reference signal x and the adaptive digital filter W
m is a value obtained by convolving each filter coefficient W mi with m , and FIG.
In will be described, the output of the drive signal generator 13, as a driving signal y m loudspeakers 7a to 7d.
【0051】しかし、ステップ124,129,131
又は133の何れかが実行された場合には、何れかの駆
動信号y1 〜y4 の値が所定値bに設定されているた
め、所定値bに設定されていない駆動信号ym について
は、基準信号xと適応ディジタルフィルタWm の各フィ
ルタ係数Wmiとを畳み込んだ値がラウドスピーカ7a〜
7dの駆動信号となるが、所定値bに設定されている駆
動信号ym については、基準信号xと適応ディジタルフ
ィルタWm の各フィルタ係数Wmiとを畳み込んだ値に、
さらに所定値bが加算された値(図2で説明すれば、駆
動信号生成部13の出力ym に加算部16においてパル
ス信号x' を加算した値)が、ラウドスピーカ7a〜7
dの駆動信号ym となる。However, steps 124, 129, 131
Or if the 133 either is executed, either the value of the drive signal y 1 ~y 4 is set to a predetermined value b, for the drive signal y m that is not set to a predetermined value b , convolving the filter coefficients W mi of the reference signal x and an adaptive digital filter W m values loudspeaker 7a~
Although the 7d driving signal, for driving signals y m which is set to a predetermined value b, and the respective filter coefficients W mi of the reference signal x and an adaptive digital filter W m to the value convolution,
(Will be described in FIG. 2, a value obtained by adding a pulse signal x 'in the addition unit 16 to an output y m of the drive signal generating unit 13) further predetermined value b is added value is, loudspeakers 7a~7
the drive signals y m of d.
【0052】ステップ135の処理を終えたら、ステッ
プ136に移行し、各駆動信号ymを各ラウドスピーカ
7a〜7dに出力する。すると、ラウドスピーカ7a〜
7dから車室6内に制御音が発生するが、制御開始直後
は適応ディジタルフィルタWm の各フィルタ係数Wmiが
最適な値に収束しているとは限らないので、必ずしも車
室6内に伝達されたこもり音が低減されるとはいえな
い。[0052] After completing the process of step 135, the process proceeds to step 136, and outputs each of the driving signals y m to each loudspeaker 7a to 7d. Then, the loudspeakers 7a to
Although control sound into the passenger compartment 6 from 7d is generated, since after start of control is not necessarily the filter coefficient W mi of the adaptive digital filter W m is converged to an optimal value, necessarily the passenger compartment 6 It cannot be said that the transmitted muffled sound is reduced.
【0053】しかし、図3に示す処理が繰り返し実行さ
れると、ステップ103でフィルタ係数WmiがLMSア
ルゴリズムに応じた上記(4)式又は(5)式に従って
逐次更新され、最適値に向かって収束していくから、車
室6内に伝達されるこもり音がラウドスピーカ7a〜7
dから発せられる制御音によって打ち消されるようにな
り、車室6内の騒音の低減が図られる。However, when the processing shown in FIG. 3 is repeatedly executed, the filter coefficient W mi is sequentially updated in step 103 according to the above equation (4) or (5) according to the LMS algorithm, and toward the optimum value. Since the sound converges, the muffled sound transmitted into the cabin 6 is transmitted to the loudspeakers 7 a to 7.
The sound is canceled by the control sound emitted from d, and the noise in the vehicle interior 6 is reduced.
【0054】そして、本実施例にあっては、駆動信号y
m をそのまま出力する場合と、駆動信号xに同期したタ
イミングで何れか一つの駆動信号ym に所定値bを加算
して出力する場合とがあり、加算せずに出力する場合に
は、ステップ105において同定用フィルタC^0lm を
上記(7)式に従って更新する一方、加算して出力する
場合には、ステップ109〜112の何れかにおいて、
その加算された駆動信号ym に対応する同定用フィルタ
C^1lm を上記(7)式に従って更新することとしてい
る。In this embodiment, the drive signal y
the case of directly outputting m, there is a case of outputting by adding a predetermined value b to any one of the drive signal y m at a timing synchronized with the drive signal x, when outputting without the addition, the step In step 105, while updating the identification filter C ^ 0lm in accordance with the above equation (7), when adding and outputting, in one of steps 109 to 112,
The identifying filter C ^ 1 lm corresponding to the summed drive signal y m is set to be updated in accordance with equation (7).
【0055】ここで、上記(7)式は、LMSアルゴリ
ズムに基づいた更新式であるが、同定用フィルタC^
1lm ,C^0lm への入力はパルス信号x' であり、その
大きさを1と考えれば、これら同定用フィルタC
^1lm ,C^0lm の出力は、そのフィルタ係数C
^1lmk,C^0lmkそのものであるため、上記(7)式に
従ってそれらフィルタ係数C^1lmk,C^0lmkを逐次更
新すると、フィルタ係数C^1lmk,C^0lmkは、変数k
毎の残留騒音信号el の平均値(同期をとった平均値)
に収束することになる。Here, the above equation (7) is an updating equation based on the LMS algorithm.
The input to 1lm , C ^ 0lm is a pulse signal x ′, and if the magnitude is considered to be 1, these identification filters C ′
The output of ^ 1lm , C ^ 0lm is the filter coefficient C
Since ^ 1lmk and C ^ 0lmk are themselves , if the filter coefficients C ^ 1lmk and C ^ 0lmk are sequentially updated according to the above equation (7), the filter coefficients C ^ 1lmk and C ^ 0lmk become variable k
Average value of the residual noise signal e l (average value synchronized)
Converges to
【0056】つまり、パルス信号x' を加算しない場合
に上記(7)式に従って更新される同定用フィルタC^
0lm は、例えば図4(a)に示すように、その時の残留
騒音信号el の波形を表すようになる。同様に、パルス
信号x' を加算した場合に更新される同定用フィルタC
^1lmもその時の残留騒音信号el の波形を表すように
なるが、この場合には、パルス信号x' に応じた同定音
が制御音とともにラウドスピーカから発せられているた
め、残留騒音信号el には、こもり音等に起因する残留
騒音の他、ラウドスピーカから発せられてマイクロフォ
ン8a〜8hに到達した同定音も含まれることになる。That is, when the pulse signal x ′ is not added, the identification filter C ^ updated according to the above equation (7)
0lm, for example, as shown in FIG. 4 (a), it will represent the waveform of the residual noise signal e l at that time. Similarly, the identification filter C updated when the pulse signal x 'is added
^ 1 lm also the residual will represent the waveforms of the noise signal e l at that time, in this case, since the identification sound corresponding to the pulse signal x 'is emitted from the loudspeaker with the control sound, residual noise signal e l includes not only residual noise due to muffled sound and the like, but also identification sound emitted from the loudspeaker and reaching the microphones 8a to 8h.
【0057】従って、同定用フィルタC^1lm は、例え
ば図4(b)に示すように、同定用フィルタC^0lm の
形状を表す図4(a)に比べてやや歪んだ波形となる
が、その歪んだ分は、同定音に起因するものである。こ
のため、同定用フィルタC^1lm の各フィルタ係数C^
1lmjから同定用フィルタC^0lm の各フィルタ係数C^
0lmjを差し引くと、図4(c)に示すように、同定音に
起因する波形のみが残ることになり、この波形は、パル
ス信号x' に応じた同定音をラウドスピーカから発生さ
せ、それをマイクロフォン8a〜8hで測定した音に相
当するから、それらラウドスピーカ及びマイクロフォン
8a〜8h間の伝達関数のインパルス応答に一致するは
ずである。Accordingly, the identification filter C ^ 1lm has a slightly distorted waveform as compared with FIG. 4A showing the shape of the identification filter C ^ 0lm , for example, as shown in FIG. 4B. The distorted portion is caused by the identification sound. Therefore, each filter coefficient C 係数 of the identification filter C ^ 1lm
From 1lmj , identification filter C ^ 0lm each filter coefficient C ^
When 0lmj is subtracted, as shown in FIG. 4 (c), only a waveform due to the identification sound remains, and this waveform causes the loudspeaker to generate an identification sound corresponding to the pulse signal x '. Since this corresponds to the sound measured by the microphones 8a to 8h, it should match the impulse response of the transfer function between the loudspeakers and the microphones 8a to 8h.
【0058】なお、同定用フィルタC^1lm に含まれる
同定音以外の残留騒音成分と、同定用フィルタC^0lm
に含まれる同定音以外の残留騒音成分とは、それらの同
定には時間的なズレがあることから瞬間的には一致しな
いが、本実施例では、残留騒音信号el の瞬時値ではな
く同期をとった平均値が同定用フィルタC^1lm ,C^
0lm のフィルタ係数C^1lmj,C^0lmjとなるようなア
ルゴリズムを採用しており、しかも、同定用フィルタC
^1lm を更新する処理と同定用フィルタC^0l m を更新
する処理とを頻繁に切り換えるようにしているため、図
3に示す処理を適宜繰り返した後には、同定音以外の残
留騒音成分は一致したものとみなすことができる。[0058] Incidentally, the residual noise component other than the identified sound included in the identification filter C ^ 1 lm, identifying filter C ^ 0lm
Are not instantaneously coincident with the residual noise components other than the identification sound included in the residual noise signal e l because they are temporally shifted in their identification. Is the average value of the filters for identification C {1lm , C }
The algorithm is such that the filter coefficients C ^ 1lmj and C ^ 0lmj of 0 lm are obtained.
^ For 1lm a process for updating and the identifying filter C ^ 0l m to frequently switch between process of updating, after appropriately repeating the process shown in FIG. 3, the residual noise component other than the identified sound matching It can be regarded as having done.
【0059】よって、ステップ113において、同定用
フィルタC^1lm の各フィルタ係数C^1lmjから同定用
フィルタC^0lm の各フィルタ係数C^0lmjを差し引い
た結果に基づいて伝達関数フィルタC^lmを設定する
と、例えば車室6内の温度,湿度,窓の開閉状況,乗員
数の変化等に応じてラウドスピーカ7a〜7d及びマイ
クロフォン8a〜8h間の伝達関数が当初の状態から変
化しても、それを追従するように伝達関数フィルタC^
lmが変更されるため、実際の物理的空間における伝達関
数と、伝達関数フィルタC^lmが表す伝達関数との間の
ズレを極力小さくすることができ、良好な騒音低減効果
を得ることができる。[0059] Therefore, in step 113, the transfer function filter C ^ lm on the basis of the result of subtracting the respective filter coefficients C ^ 0lmj of identifying filters C ^ 0lm from each filter coefficient C ^ 1lmj of identifying filters C ^ 1 lm If set, even if the transfer function between the loudspeakers 7a to 7d and the microphones 8a to 8h changes from the initial state according to, for example, changes in the temperature, humidity, window opening / closing status, and the number of occupants in the passenger compartment 6, The transfer function filter C ^
Since lm is changed, the difference between the transfer function in the actual physical space and the transfer function represented by the transfer function filter C ^ lm can be minimized, and a good noise reduction effect can be obtained. .
【0060】しかも、同定用フィルタC^1lm ,C^
0lm も、結局は残留騒音信号el の同期をとった平均値
に収束することからその同定精度も高く、且つ、それら
同定用フィルタC^1lm ,C^0lm のフィルタ係数C^
1lmj,C^0lmjの減算という簡易な演算で伝達関数フィ
ルタC^lmが同定されるから、高い精度で伝達関数フィ
ルタC^lmを設定することができる。このことは、同定
音のパワーをそれほど大きくしなくても済むということ
であるから、同定音を発生することにより乗員に不快感
を与えてしまうようなこともない。Moreover, the identification filters C {1lm , C }
0Lm also eventually higher the identification accuracy since it converges to the average value taken synchronization residual noise signal e l, and, for those identified filter C ^ 1 lm, filter coefficients C ^ 0lm C ^
Since the transfer function filter C ^ lm is identified by a simple operation of subtraction of 1lmj and C ^ 0lmj , the transfer function filter C ^ lm can be set with high accuracy. This means that it is not necessary to increase the power of the identification sound so much, so that the generation of the identification sound does not give an occupant a discomfort.
【0061】そして、本実施例にあっては、ホワイトノ
イズではなく、パルス信号x' に基づいて伝達関数フィ
ルタC^lmを同定する構成であるため、下記のような利
点がある。第1に、ホワイトノイズを発生させるために
は、全周波数帯域に渡って一定のパワーを有するランダ
ム信号を生成しなければならず、演算によって生成する
には計算量の増大を招き、ランダム信号を記憶しておく
構成とするとメモリ容量の増大を招いてしまうのに対
し、パルス信号x' の生成は非常に容易であるから、演
算負荷の軽減又はメモリ容量の削減が図られる。In this embodiment, since the transfer function filter C ^ lm is identified based on the pulse signal x ′ instead of the white noise, the following advantages are provided. First, in order to generate white noise, it is necessary to generate a random signal having constant power over the entire frequency band. While the configuration in which the data is stored causes an increase in the memory capacity, the generation of the pulse signal x 'is very easy, so that the calculation load or the memory capacity can be reduced.
【0062】第2に、同定用フィルタC^1lm ,C^
0lm への入力がパルス信号x' であるため、本来ならば
畳み込み積分を行う必要があるのに対し、クロックパル
スCPに同期してフィルタ係数C^1lmj,C^0lmjを出
力するだけで済むから、これによっても大幅な演算負荷
の軽減が図られる。第3に、本実施例では、同定用フィ
ルタC^1lm ,C^0lm のフィルタ係数C^1lmj,C^
0lmjの更新にLMSアルゴリズムを適用しているが、基
準信号としてパルス信号x' を適用しているため、上記
(7)式からも明らかなように、k+Np番目のフィル
タ係数C^*lm(k+Np) のみを更新すればよいので、これ
によっても演算負荷の軽減が図られる。Second, an identification filter C {1lm , C }
Since the input to 0lm is a pulse signal x ', whereas it is necessary to perform convolution would otherwise need only outputs filter coefficient C ^ 1lmj, a C ^ 0lmj in synchronism with the clock pulses C P Therefore, the calculation load can be greatly reduced. Third, in the present embodiment, the filter coefficients C ^ 1lmj , C ^ of the identification filters C ^ 1lm , C ^ 0lm.
Although the LMS algorithm is applied to update 0lmj, since the pulse signal x 'is applied as the reference signal, the k + Np-th filter coefficient C ^ * lm (k Since only + Np) needs to be updated, the calculation load can be reduced.
【0063】以上から、本実施例の構成であれば、従来
に比べて非常に少ない演算量で高精度に伝達関数フィル
タC^lmを同定することができるので、高速演算可能な
高価な演算素子を導入することなく、且つ、乗員に不快
感を与えることなく、車両のように頻繁に伝達関数が変
化する空間であっても、良好な騒音低減効果が得られる
ものである。As described above, according to the configuration of the present embodiment, the transfer function filter Clmlm can be identified with a very small amount of operation and a high accuracy compared with the conventional one, so that an expensive operation element capable of high-speed operation is provided. , And a good noise reduction effect can be obtained even in a space where the transfer function changes frequently, such as a vehicle, without introducing discomfort to the occupant.
【0064】ここで、本実施例にあっては、クランク角
センサ5及び基準信号生成部11によって騒音発生状態
検出手段が構成され、加算部16,断続スイッチ部17
及びステップ124,129,131,133,135
の処理によって信号加算手段が構成され、減算部18,
切換スイッチ部19,21,フィルタ係数更新部20及
びステップ105〜112の処理によって同定用フィル
タ設定手段が構成され、減算部22及びステップ113
の処理によって伝達関数フィルタ設定手段が構成され
る。Here, in this embodiment, the noise generation state detecting means is constituted by the crank angle sensor 5 and the reference signal generating section 11, and the adding section 16 and the intermittent switch section 17 are provided.
And steps 124,129,131,133,135
Constitutes a signal adding means, and the subtracting section 18,
The selection switch units 19 and 21, the filter coefficient update unit 20, and the processing of steps 105 to 112 constitute an identification filter setting unit, and the subtraction unit 22 and step 113
The transfer function filter setting means is constituted by the above processing.
【0065】図5は、本発明の第2実施例の構成を示す
ブロック図であって、上記第1実施例の図2と同様に、
コントローラ10の機能構成を示している。なお、本実
施例も、上記第1実施例と同様に、エンジンから車室内
に伝達されるこもり音の低減を図る車両用能動型騒音制
御装置に本発明を適用したものであり、本実施例のその
他の構成は、上記第1実施例と同様である。FIG. 5 is a block diagram showing the configuration of the second embodiment of the present invention. As in FIG. 2 of the first embodiment, FIG.
2 shows a functional configuration of the controller 10. In this embodiment, similarly to the first embodiment, the present invention is applied to an active noise control device for a vehicle for reducing a muffled sound transmitted from an engine into a vehicle interior. The other configuration is the same as that of the first embodiment.
【0066】即ち、本実施例は、駆動信号ym を生成す
る適応ディジタルフィルタWm に関しても同期式LMS
アルゴリズム(より具体的には、同期式Filtere
d−X LMSアルゴリズム)を適用したものであり、
従って、基準信号生成部11は、上記第1実施例とは異
なり、クランク角信号Xに基づき、こもり音に同期した
インパルス列でなる基準信号xを生成し出力する。[0066] That is, this embodiment, synchronous LMS regard adaptive digital filter W m to generate a drive signal y m
Algorithm (more specifically, synchronous Filtere
d-X LMS algorithm).
Therefore, unlike the first embodiment, the reference signal generator 11 generates and outputs a reference signal x composed of an impulse train synchronized with the muffled sound, based on the crank angle signal X.
【0067】そして、この基準信号xを駆動信号生成部
13に供給する一方、基準信号xに同期したパルス信号
として加算部16,同定用フィルタC^1lm ,C^0lm
及びフィルタ係数更新部20に供給している。さらに、
上記第1実施例では存在した基準処理信号生成部14を
省略し、適応ディジタルフィルタWm のフィルタ係数W
miを更新するフィルタ係数更新部15には、同定用フィ
ルタC^1lm の出力と同定用フィルタC^0lm の出力と
の差を演算する減算部30の出力をそのまま供給してい
る。つまり、本実施例では、同定用フィルタC^1lm 及
びC^0lm が、上記第1実施例における伝達関数フィル
タC^lmの機能を兼ね備えているのである。Then, the reference signal x is supplied to the drive signal generation unit 13, while the addition unit 16, the identification filters C # 1lm and C # 0lm are supplied as pulse signals synchronized with the reference signal x.
And a filter coefficient updating unit 20. further,
In the first embodiment to omit the reference processed signal generating unit 14 that existed, the filter coefficient W of the adaptive digital filter W m
The output of the subtraction unit 30 that calculates the difference between the output of the identification filter C # 1lm and the output of the identification filter C # 1lm is supplied to the filter coefficient update unit 15 that updates mi as it is. That is, in the present embodiment, the identification filters C # 1lm and C # 0lm also have the function of the transfer function filter C # lm in the first embodiment.
【0068】図6は、本実施例のコントローラ10内に
おける処理の概要を示すフローチャートであり、上記第
1実施例で説明した図3と同等の処理には同じ符号を付
し、その重複する説明は省略する。即ち、本実施例にあ
っては、ステップ105で同定用フィルタC^0lm のフ
ィルタ係数C^0lm(k+Np) の更新を行った後に、ステッ
プ201に移行し、下記の(9)式に従って、フィルタ
係数C^0lm(k+Np) の負の変化量Δrlmk を演算する。FIG. 6 is a flowchart showing the outline of the processing in the controller 10 of the present embodiment. The same reference numerals are given to the same processes as those of FIG. 3 described in the first embodiment, and the same explanation will be given. Is omitted. That is, in this embodiment, after the filter coefficient C 係数 0lm (k + Np) of the identification filter C ^ 0lm is updated in step 105, the process proceeds to step 201, and according to the following equation (9). , The negative change amount Δr lmk of the filter coefficient C ^ 0lm (k + Np) is calculated.
【0069】 Δrlmk =−{C^0lm(k+Np) (n+1)−C^0lm(k+Np) (n)} ……(9) 次いで、ステップ202に移行し、下記の(10)式に従
って、ステップ103の適応ディジタルフィルタWm の
フィルタ係数Wmiの更新に使用する基準処理信号rlmk
を累算する。Δr lmk = − {C ^ 0lm (k + Np) (n + 1) −C ^ 0lm (k + Np) (n)} (9) Then, the process proceeds to step 202 and the following (10) According to the equation, the reference processing signal r lmk used for updating the filter coefficient W mi of the adaptive digital filter W m in step 103
Is accumulated.
【0070】 rlmk =rlmk +Δrlmk ……(10) 一方、ステップ109,110,111又は112で同
定用フィルタC^1lmのフィルタ係数C^1lm(k+Np) の
更新を行った後に、ステップ203,205,207又
は209に移行し、下記の(11)式に従って、フィルタ
係数C^1lm(k+ Np) の変化量Δrlmk を演算する。R lmk = r lmk + Δr lmk (10) On the other hand, after updating the filter coefficient C ^ 1lm (k + Np) of the identification filter C ^ 1lm in step 109, 110, 111 or 112, The process proceeds to steps 203, 205, 207 or 209, and calculates a change amount Δr lmk of the filter coefficient C ^ 1lm (k + Np) according to the following equation (11).
【0071】 Δrlmk =C^0lm(k+Np) (n+1)−C^0lm(k+Np) (n) ……(11) 次いで、ステップ204,206,208又は210に
移行し、下記の(12)式に従って、基準処理信号rlmk
を累算する。 rlmk =rlmk +Δrlmk ……(12) このような演算によって基準処理信号rlmk を求められ
るのは、基準信号xをインパルス列とした結果、伝達関
数フィルタC^lmへの入力と、同定用フィルタC
^1lm ,C^0lm への入力とが等しく、何れの出力もフ
ィルタ係数C^lmk ,C^1lmk,C^0lmkであり、しか
も、上記第1実施例で説明したように、伝達関数フィル
タC^lmのフィルタ係数C^lmj が同定用フィルタC^
1lm ,C^0lmのフィルタ係数C^1lmj,C^0lmjの差
から求められるのであるから、それらフィルタ係数C^
1lmk,C^0lmkの変化量を極性を異ならせて累積すれ
ば、伝達関数フィルタC^lmの各フィルタ係数C^lmj
と基準信号xとを畳み込んで基準処理信号rlmk を演算
したのと同じ結果が得られるからである。Δr lmk = C ^ 0lm (k + Np) (n + 1) −C ^ 0lm (k + Np) (n) (11) Then, the process proceeds to step 204, 206, 208 or 210, and According to equation (12), the reference processing signal r lmk
Is accumulated. r lmk = r lmk + Δr lmk (12) The reference processing signal r lmk can be obtained by such an operation because the reference signal x is an impulse train, the input to the transfer function filter C ^ lm and the identification. Filter C
^ 1lm and C ^ 0 lm are equal to the input, and both outputs are filter coefficients C ^ lmk , C ^ 1lmk , and Clm 0 lmk . Further, as described in the first embodiment, the transfer function filter C フ ィ ル タ lm filter coefficient Clm lmj is the identification filter C 同 定
1 lm, C ^ filter coefficient 0lm C ^ 1lmj, because of being determined from the difference between the C ^ 0lmj, they filter coefficient C ^
If the changes of 1lmk and Clm0lmk are accumulated with different polarities, each filter coefficient C ^ lmj of the transfer function filter C ^ lm
This is because the same result as that of calculating the reference processing signal r lmk by convolving the reference processing signal x with the reference signal x is obtained.
【0072】さらに、基準信号xをインパルス列とした
結果、ステップ135における駆動信号ym の演算は、
下記の(13)式に従って行われる。 ym =ym +Wmk ……(13) このように、本実施例の構成であれば、上記第1実施例
よりもさらに演算負荷の軽減が図られるという利点があ
る。その他の作用効果は、上記第1実施例と同様であ
る。[0072] Further, as a result of the reference signal x and an impulse train, the operation of the drive signal y m at step 135,
This is performed according to the following equation (13). y m = y m + W mk ...... (13) Thus, with the configuration of this embodiment, there is an advantage that the further reduction of calculation load than the first embodiment can be achieved. Other functions and effects are the same as those of the first embodiment.
【0073】ここで、本実施例では、演算負荷の軽減を
図った結果、形の上では伝達関数フィルタC^lmは存在
せず、同定用フィルタC^1lm ,C^0lm が伝達関数フ
ィルタC^lmの機能を兼ね備えていることから、同定用
フィルタC^1lm ,C^0lm,減算部30及びステップ
201〜210における処理が伝達関数フィルタ設定手
段に対応するものである。In this embodiment, as a result of reducing the calculation load, the transfer function filter C ^ lm does not exist in the form, and the identification filters C ^ 1lm and C ^ 0lm are used in the transfer function filter CC. Since it also has the function of ^ lm , the processing in the identification filters C ^ 1lm , C ^ 0lm , the subtraction unit 30, and steps 201 to 210 corresponds to the transfer function filter setting means.
【0074】なお、上記各実施例では、本発明に係る能
動型騒音制御装置を、周期的な騒音である車両のこもり
音の低減を図る装置に適用した場合について説明した
が、本発明は、例えばロード・ノイズのようなランダム
・ノイズの低減を図る装置であって適用可能である。即
ち、例えばロード・ノイズの低減を図る装置であれば、
基準信号xはサスペンションに生じる上下加速度となる
が、その基準信号xを30〜40Hz程度の低周波を抽
出するローパス・フィルタで処理し、その低周波成分に
同期してパルス信号x' を生成すれば、上記第1実施例
と同等の作用効果が得られる。In each of the above embodiments, the case where the active noise control device according to the present invention is applied to a device for reducing the muffled sound of a vehicle, which is a periodic noise, has been described. For example, the present invention is applicable to a device for reducing random noise such as road noise. That is, for example, if the device aims to reduce road noise,
The reference signal x is a vertical acceleration generated in the suspension. The reference signal x is processed by a low-pass filter that extracts a low frequency of about 30 to 40 Hz, and a pulse signal x 'is generated in synchronization with the low-frequency component. In this case, the same operation and effect as in the first embodiment can be obtained.
【0075】また、上記各実施例では、本発明に係る能
動型騒音制御装置を車両に適用した場合について説明し
たが、本発明の適用対象はこれに限定されるものではな
く、車室以外の空間であっても当然に適用可能である。
さらに、上記各実施例では、同期式LMSアルゴリズム
を利用して、同定用フィルタC^1lm ,C^0lm の各フ
ィルタ係数が、残留騒音信号el の同期をとった平均値
となるようにしているが、残留騒音信号el の同期をと
った平均値の求め方はこれに限定されるものではなく、
例えば、所定時間内の残留騒音信号elを同期をとって
累算しておき、その累算結果を累算回数で割って求める
ようにしてもよいし、或いは、重みを付けて移動平均を
求める手法を利用してもよい。In each of the above embodiments, the case where the active noise control device according to the present invention is applied to a vehicle has been described. However, the application object of the present invention is not limited to this. Naturally, it can be applied to a space.
Furthermore, in the above embodiments, by utilizing the synchronous LMS algorithm, identifying filter C ^ 1 lm, each filter coefficient of the C ^ 0lm is, as a mean value synchronizing the residual noise signal e l However, the method of obtaining the synchronized average value of the residual noise signal el is not limited to this.
For example, advance accumulates synchronize the residual noise signal e l within a predetermined time may be calculated by dividing the accumulated result in accumulation times, or a moving average with a weighting You may use the technique of finding.
【0076】またさらに、上記各実施例では、同定用フ
ィルタC^1lm ,C^0lm の各フィルタ係数を残留騒音
信号el の同期をとった平均値に設定する場合について
説明しているが、空間内に残留する騒音が変化の少ない
安定した状態にあることが明らかであれば、残留騒音信
号el の瞬時値を同定用フィルタC^1lm ,C^0lmの
各フィルタ係数としても、伝達関数の同定は可能であ
る。[0076] Furthermore, in the above embodiments, identifying filter C ^ 1 lm, it has been described for the case of setting the average value synchronized for each filter coefficient residual noise signal e l of C ^ 0lm, if it is clear that the noise remaining in the space is in the less stable state of change, the remaining noise signal e identifying filter the instantaneous value of l C ^ 1 lm, even the filter coefficients of C ^ 0lm, transfer function Is possible.
【0077】[0077]
【発明の効果】以上説明したように、本発明によれば、
伝達関数一つに対して二つの同定用フィルタを設け、そ
れら同定用フィルタを適宜更新し、そして、それら二つ
の同定用フィルタのフィルタ係数の差に基づいて伝達関
数フィルタを設定する構成としたため、従来に比べて非
常に少ない演算量で高精度に伝達関数フィルタを同定す
ることができ、高速演算可能な高価な演算素子を導入す
ることなく、車両のように頻繁に伝達関数が変化する空
間であっても、良好な騒音低減効果を得ることができる
という効果がある。As described above, according to the present invention,
Because two identification filters are provided for one transfer function, the identification filters are appropriately updated, and the transfer function filter is set based on the difference between the filter coefficients of the two identification filters. A space in which the transfer function changes frequently like a vehicle without introducing expensive operation elements capable of performing high-speed operation, which makes it possible to identify a transfer function filter with a very small amount of operation compared to the conventional art. However, there is an effect that a good noise reduction effect can be obtained.
【図1】本発明の第1実施例の全体構成を示す図であ
る。FIG. 1 is a diagram showing an overall configuration of a first embodiment of the present invention.
【図2】第1実施例のコントローラの機能構成を示すブ
ロック図である。FIG. 2 is a block diagram illustrating a functional configuration of a controller according to the first embodiment.
【図3】第1実施例のコントローラ内で実行される処理
の概要を示すフローチャートである。FIG. 3 is a flowchart illustrating an outline of processing executed in a controller according to the first embodiment;
【図4】伝達関数を設定できる理由を説明する波形図で
ある。FIG. 4 is a waveform chart for explaining the reason why a transfer function can be set.
【図5】第2実施例のコントローラの機能構成を示すブ
ロック図である。FIG. 5 is a block diagram illustrating a functional configuration of a controller according to a second embodiment.
【図6】第2実施例のコントローラ内で実行される処理
の概要を示すフローチャートである。FIG. 6 is a flowchart illustrating an outline of processing executed in a controller according to a second embodiment.
【符号の説明】 1 車両用能動型騒音制御装置 4 エンジン(騒音源) 5 クランク角センサ 6 車室(空間) 7a〜7d ラウドスピーカ(制御音源) 8a〜8h マイクロフォン(残留騒音検出手
段) 10 コントローラ 11 基準信号生成部 12 パルス信号生成部 13 駆動信号生成部(駆動信号生成手
段) 14 基準処理信号生成部 15 フィルタ係数更新部(適応処理手
段) 16 加算部 17 断続スイッチ部 18,22 減算部 19,21 切換スイッチ部 20 フィルタ係数更新部 Wm 適応ディジタルフィルタ C^lm 伝達関数フィルタ C^1lm ,C^0lm 同定用フィルタ[Description of Signs] 1 Active noise control device for vehicle 4 Engine (noise source) 5 Crank angle sensor 6 Cab (space) 7a to 7d Loudspeaker (control sound source) 8a to 8h Microphone (residual noise detection means) 10 Controller Reference Signs List 11 Reference signal generation unit 12 Pulse signal generation unit 13 Drive signal generation unit (drive signal generation unit) 14 Reference processing signal generation unit 15 Filter coefficient update unit (adaptive processing unit) 16 Addition unit 17 Intermittent switch unit 18, 22 Subtraction unit 19 , 21 changeover switch unit 20 the filter coefficient updating unit W m adaptive digital filter C ^ lm transfer function filter C ^ 1lm, C ^ 0lm identifying filters
フロントページの続き (56)参考文献 特開 昭62−193310(JP,A) 特開 平4−267298(JP,A) (58)調査した分野(Int.Cl.7,DB名) G10K 11/178 B60R 11/02 F01N 1/00 F01N 1/06 H03H 17/00 601 H03H 17/02 601 H03H 21/00 Continuation of the front page (56) References JP-A-62-193310 (JP, A) JP-A-4-267298 (JP, A) (58) Fields investigated (Int. Cl. 7 , DB name) G10K 11 / 178 B60R 11/02 F01N 1/00 F01N 1/06 H03H 17/00 601 H03H 17/02 601 H03H 21/00
Claims (2)
音を発生可能な制御音源と、前記騒音源の騒音発生状態
を検出し基準信号として出力する騒音発生状態検出手段
と、前記空間内の所定位置における残留騒音を検出し残
留騒音信号として出力する残留騒音検出手段と、フィル
タ係数可変の適応ディジタルフィルタと、前記基準信号
を前記適応ディジタルフィルタでフィルタ処理して前記
制御音源の駆動信号を生成する駆動信号生成手段と、前
記制御音源及び前記残留騒音検出手段間の伝達関数をモ
デル化した伝達関数フィルタと、前記基準信号を前記伝
達関数フィルタでフィルタ処理した結果と前記残留騒音
信号とに基づいて前記空間内の騒音が低減するように前
記適応ディジタルフィルタのフィルタ係数を更新する適
応処理手段と、を備えた能動型騒音制御装置において、 前記基準信号に同期したパルス信号を一時的に前記駆動
信号に加算可能な信号加算手段と、フィルタ係数可変で
且つ前記伝達関数一つに対して二つの同定用フィルタ
と、前記信号加算手段が前記加算を行った際の前記残留
騒音信号に基づいて前記一方の同定用フィルタを設定し
且つ前記信号加算手段が前記加算を行わなかった際の前
記残留騒音信号に基づいて前記他方の同定用フィルタを
設定する同定用フィルタ設定手段と、前記一方の同定用
フィルタのフィルタ係数及び前記他方の同定用フィルタ
のフィルタ係数の差に基づいて前記伝達関数フィルタを
設定する伝達関数フィルタ設定手段と、を備えたことを
特徴とする能動型騒音制御装置。1. A control sound source capable of generating a control sound in a space to which noise is transmitted from a noise source, a noise generation state detecting means for detecting a noise generation state of the noise source and outputting as a reference signal, Residual noise detecting means for detecting residual noise at a predetermined position and outputting the same as a residual noise signal, an adaptive digital filter having a variable filter coefficient, and filtering the reference signal with the adaptive digital filter to generate a drive signal of the control sound source. A drive signal generating means for generating, a transfer function filter that models a transfer function between the control sound source and the residual noise detecting means, and a result of filtering the reference signal with the transfer function filter and the residual noise signal. Adaptive processing means for updating the filter coefficient of the adaptive digital filter based on the noise in the space based on the noise. An active noise control device, comprising: a signal adding means capable of temporarily adding a pulse signal synchronized with the reference signal to the drive signal; and a filter coefficient variable and two identification functions for one transfer function. A filter, and setting the one filter for identification based on the residual noise signal when the signal adding means has performed the addition, and applying the residual noise signal when the signal adding means has not performed the addition. An identification filter setting means for setting the other identification filter based on the difference, and a transmission setting the transfer function filter based on a difference between a filter coefficient of the one identification filter and a filter coefficient of the other identification filter. An active noise control device comprising: a function filter setting unit.
音を発生可能な制御音源と、前記騒音源の騒音発生状態
を検出し基準信号として出力する騒音発生状態検出手段
と、前記空間内の所定位置における残留騒音を検出し残
留騒音信号として出力する残留騒音検出手段と、フィル
タ係数可変の適応ディジタルフィルタと、前記基準信号
を前記適応ディジタルフィルタでフィルタ処理して前記
制御音源の駆動信号を生成する駆動信号生成手段と、前
記制御音源及び前記残留騒音検出手段間の伝達関数をモ
デル化した伝達関数フィルタと、前記基準信号を前記伝
達関数フィルタでフィルタ処理した結果と前記残留騒音
信号とに基づいて前記空間内の騒音が低減するように前
記適応ディジタルフィルタのフィルタ係数を更新する適
応処理手段と、を備えた能動型騒音制御装置において、 前記基準信号に同期したパルス信号を一時的に前記駆動
信号に加算可能な信号加算手段と、フィルタ係数可変で
且つ前記伝達関数一つに対して二つの同定用フィルタ
と、前記信号加算手段が前記加算を行った際の前記残留
騒音信号の同期をとった平均値を前記一方の同定用フィ
ルタのフィルタ係数とし且つ前記信号加算手段が前記加
算を行わなかった際の前記残留騒音信号の同期をとった
平均値を前記他方の同定用フィルタのフィルタ係数とす
る同定用フィルタ設定手段と、前記一方の同定用フィル
タのフィルタ係数及び前記他方の同定用フィルタのフィ
ルタ係数の差に基づいて前記伝達関数フィルタを設定す
る伝達関数フィルタ設定手段と、を備えたことを特徴と
する能動型騒音制御装置。2. A control sound source capable of generating a control sound in a space to which noise is transmitted from a noise source; a noise generation state detection means for detecting a noise generation state of the noise source and outputting the noise generation state as a reference signal; Residual noise detecting means for detecting residual noise at a predetermined position and outputting the same as a residual noise signal, an adaptive digital filter having a variable filter coefficient, and filtering the reference signal with the adaptive digital filter to generate a drive signal of the control sound source. A drive signal generating means for generating, a transfer function filter that models a transfer function between the control sound source and the residual noise detecting means, and a result of filtering the reference signal with the transfer function filter and the residual noise signal. Adaptive processing means for updating the filter coefficient of the adaptive digital filter based on the noise in the space based on the noise. An active noise control device, comprising: a signal adding means capable of temporarily adding a pulse signal synchronized with the reference signal to the drive signal; and a filter coefficient variable and two identification functions for one transfer function. A filter and a synchronized average value of the residual noise signal when the signal addition means performs the addition, as a filter coefficient of the one identification filter, and when the signal addition means does not perform the addition. Identification filter setting means for setting an average value synchronized with the residual noise signal as a filter coefficient of the other identification filter, and a filter coefficient of the one identification filter and a filter coefficient of the other identification filter And a transfer function filter setting means for setting the transfer function filter based on the difference between the two.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP01413293A JP3293922B2 (en) | 1993-01-29 | 1993-01-29 | Active noise control device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP01413293A JP3293922B2 (en) | 1993-01-29 | 1993-01-29 | Active noise control device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH06230786A JPH06230786A (en) | 1994-08-19 |
JP3293922B2 true JP3293922B2 (en) | 2002-06-17 |
Family
ID=11852614
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP01413293A Expired - Fee Related JP3293922B2 (en) | 1993-01-29 | 1993-01-29 | Active noise control device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3293922B2 (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3728837B2 (en) * | 1996-12-12 | 2005-12-21 | 住友電気工業株式会社 | Active noise control device |
WO1999053476A1 (en) * | 1998-04-15 | 1999-10-21 | Fujitsu Limited | Active noise controller |
JP3804275B2 (en) * | 1998-06-22 | 2006-08-02 | 日産自動車株式会社 | Active noise vibration control device |
CN113421541B (en) * | 2021-08-23 | 2021-10-29 | 南京南大电子智慧型服务机器人研究院有限公司 | Frequency domain active noise control system based on coefficient smoothing |
-
1993
- 1993-01-29 JP JP01413293A patent/JP3293922B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH06230786A (en) | 1994-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3094517B2 (en) | Active noise control device | |
JP4077383B2 (en) | Active vibration noise control device | |
JPH05249984A (en) | Active noise control device for vehicle | |
KR960011152B1 (en) | Active vibration control system | |
JP3293922B2 (en) | Active noise control device | |
JP3579898B2 (en) | Vehicle vibration control device and vibration control method | |
JP2003195951A (en) | Active type vibration control device and same for vehicle | |
JP2876896B2 (en) | Active noise control system for vehicles | |
JPH07219560A (en) | Active noise controller | |
JP3367133B2 (en) | Active noise control device for vehicle and active vibration control device for vehicle | |
JP3517887B2 (en) | Active noise control system for vehicles | |
JP3617079B2 (en) | Active noise control device and active vibration control device | |
JPH04342296A (en) | Active type noise controller | |
JP3674963B2 (en) | Active noise control device and active vibration control device | |
JP3371477B2 (en) | Active vibration control device for vehicle and active noise control device for vehicle | |
JPH06230787A (en) | Active noise controller | |
JP3503155B2 (en) | Active noise control device and active vibration control device | |
JPH07239690A (en) | On-vehicle active noise controller and on-vehicle active vibration controller | |
JP3467830B2 (en) | Method for determining convergence coefficient of adaptive control device, active noise control device, and active vibration control device | |
JPH06250674A (en) | Active noise controller | |
JP2743639B2 (en) | Active noise control device | |
JPH06314097A (en) | Active noise controller | |
JPH08190390A (en) | Active noise controller | |
JPH0728477A (en) | Active vibration control device and active noise control device | |
JPH0588683A (en) | Noise controller for vehicle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090405 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090405 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100405 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110405 Year of fee payment: 9 |
|
LAPS | Cancellation because of no payment of annual fees |