JPH0769707B2 - Noise control device - Google Patents
Noise control deviceInfo
- Publication number
- JPH0769707B2 JPH0769707B2 JP63058031A JP5803188A JPH0769707B2 JP H0769707 B2 JPH0769707 B2 JP H0769707B2 JP 63058031 A JP63058031 A JP 63058031A JP 5803188 A JP5803188 A JP 5803188A JP H0769707 B2 JPH0769707 B2 JP H0769707B2
- Authority
- JP
- Japan
- Prior art keywords
- noise
- noise control
- signal
- observation
- filter
- 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)
- Circuit For Audible Band Transducer (AREA)
Description
【発明の詳細な説明】 「産業上の利用分野」 この発明は観測マイクロホンで受音した騒音から合成し
た騒音制御信号をスピーカから放射することにより、室
内騒音を抑圧する能動的な騒音制御装置に関するもので
ある。TECHNICAL FIELD The present invention relates to an active noise control device that suppresses indoor noise by radiating a noise control signal synthesized from noise received by an observation microphone from a speaker. It is a thing.
「従来の技術」 第5図は従来の騒音制御装置の構成図である。図におい
て1は騒音源、2は騒音源1が放射する騒音信号を受音
する観測マイクロホン、31,32は観測マイクロホン2が
受音した騒音信号を用いて騒音制御信号を発生するフィ
ルタ、41,42はフィルタ31,32で発生した騒音制御信号を
室内に放射するスピーカ、5は騒音源1から放射される
騒音信号と上記騒音制御信号との差信号で表される誤差
信号を観測する騒音制御点である。"Prior Art" FIG. 5 is a block diagram of a conventional noise control device. In the figure, 1 is a noise source, 2 is an observation microphone that receives a noise signal emitted by the noise source 1, 3 1 and 3 2 are filters that generate a noise control signal using the noise signal received by the observation microphone 2, 4 1 and 4 2 are speakers that radiate the noise control signals generated by the filters 3 1 and 3 2 into the room, and 5 is an error represented by a difference signal between the noise signal emitted from the noise source 1 and the noise control signal. It is a noise control point for observing signals.
第5図で時刻jに騒音源1が発生する騒音信号をxjとす
れば、騒音制御点5で観測される誤差信号yjは次のよう
に表される。When the noise signal generated by the noise source 1 at time j in FIG. 5 is x j , the error signal y j observed at the noise control point 5 is expressed as follows.
yj={gn−gm・(gs1h1+gs2h2)}・xj …… 但し gm:騒音源1と観測マイクロホン2との間の室内伝達関
数 gs1,gs2:スピーカ41,42と騒音制御点5との間の室内伝
達関数 h1,h2:フィルタ31,32の伝達関数 gn:騒音源1と騒音制御点5との間の室内伝達関数 式より誤差信号yjを制御する(零にする)ためには、
次式を満足させる伝達関数h1,h2を有するフィルタ31,32
を実現すれば良いことが分かる。y j = {g n −g m・ (g s1 h 1 + g s2 h 2 )} ・ x j …… where g m : Room transfer function between noise source 1 and observation microphone 2 g s1 , g s2 : Indoor transfer function between speakers 4 1 and 4 2 and noise control point 5 h 1 and h 2 : Transfer function of filter 3 1 and 3 2 g n : Indoor transfer between noise source 1 and noise control point 5 To control (zero) the error signal y j from the function expression,
Filters 3 1 and 3 2 having transfer functions h 1 and h 2 satisfying the following expression
You can see that
gm・(gs1h1+gs2h2)=gn …… 式を変形すれば、次の′式が得られる。g m · (g s1 h 1 + g s2 h 2 ) = g n …… By transforming the formula, the following ′ formula is obtained.
gs1h1+gs2h2=gn・1/gm ……′ ここで伝達関数gs1,gs2が共通の零点を持たないよう、
スピーカ41,42を適当に配置すれば、式左辺gs1h1+g
s2h2を任意の伝達関数に合成するフィルタ31,32の存在
性が、多入力多出力逆フィルタ処理定理MINT〔1〕によ
って保証される。g s1 h 1 + g s2 h 2 = g n · 1 / g m …… ′ Here, the transfer functions g s1 and g s2 have no common zero.
If the speakers 4 1 and 4 2 are properly arranged, the left side of the formula g s1 h 1 + g
The existence of the filters 3 1 , 3 2 that synthesize s2 h 2 into an arbitrary transfer function is guaranteed by the multi-input multi-output inverse filtering theorem MINT [1].
しかしながら観測マイクロホン2を騒音源1の近傍に置
かない場合には、観測マイクロホン2と騒音源1の間の
室内伝達関数gmが非最小位相関数となるので、gn・1/gm
は安定なものとしては存在しない。即ち従来法では観測
マイクロホン2が騒音源1の近傍に設置できず、受音す
る騒音信号に室内反射音が加えられている場合には、騒
音制御点5で観測される室内騒音を制御することができ
ないのである。However, when the observation microphone 2 is not placed near the noise source 1, the room transfer function g m between the observation microphone 2 and the noise source 1 is a non-minimum phase function, so g n · 1 / g m
Does not exist as stable. That is, when the observation microphone 2 cannot be installed in the vicinity of the noise source 1 in the conventional method and the indoor reflected sound is added to the received noise signal, the indoor noise observed at the noise control point 5 should be controlled. You cannot do it.
この発明の目的は観測マイクロホンを騒音源近傍に設置
できず、受音する騒音信号に室内反射音が付加されてい
る場合にも、良好な騒音抑圧性能を示す、能動的な騒音
制御装置を提供することである。An object of the present invention is to provide an active noise control device that exhibits good noise suppression performance even when an observation microphone cannot be installed near a noise source and an indoor reflected sound is added to a received noise signal. It is to be.
「課題を解決するための手段」 この発明は騒音観測マイクロホン・フィルタ・スピーカ
より構成される騒音制御ユニットを複数個用いることに
より、該騒音制御ユニット未満の数の騒音制御点におい
て、任意の音圧を実現するものである。以下、騒音制御
点を1点とした場合について第1図を参照しながら、こ
の発明の説明を行う。図において第5図と重複する部分
には同一の番号を付け説明を省略している。また図中
61,62はこの発明の騒音制御ユニットを表し、71,72,81,
82および91,92はそれぞれ騒音制御ユニット61,62を構成
する騒音観測マイクロホン、フィルタ、スピーカを表し
ている。"Means for Solving the Problem" The present invention uses a plurality of noise control units each consisting of a noise observation microphone, a filter, and a speaker, so that at any noise control point less than the noise control units, an arbitrary sound pressure Is realized. The present invention will be described below with reference to FIG. 1 in the case where the number of noise control points is one. In the figure, the same parts as those in FIG. 5 are designated by the same reference numerals and the description thereof is omitted. Also in the figure
6 1 and 6 2 represent the noise control unit of the present invention, and 7 1 , 7 2 , 8 1 ,
Reference numerals 8 2 and 9 1 and 9 2 represent a noise observation microphone, a filter, and a speaker, which constitute the noise control unit 6 1 and 6 2 , respectively.
第1図の騒音制御系をブロック図を用いて表せば第2図
のようになる。第2図でgm1,gm2は観測騒音源1からマ
イクロホン71,72までの室内伝達関数、gs1,gs2はスピー
カ91,92から騒音制御点5までの間の室内伝達関数、h1,
h2はフィルタ81,82の伝達関数、gnは騒音源1と騒音制
御点5との間の室内伝達関数を表している。さらに室内
伝達関数gs1,gs2とフィルタ81,82の伝達関数h1,h2を入
れ替えると、第2図は第3図のように変形できる。第3
図より時刻jに騒音制御点5で観測される誤差信号y
jは、 yj={gn−(gm1gs1h1+gm2gs2h2)}・xj …… 但し xj:時刻jに騒音源が発生する騒音信号 となる。従って誤差信号yjを制御する(零にする)に
は、次式を満足する伝達関数h1,h2を有するフィルタ81,
82を実現すれば良い。The noise control system of FIG. 1 can be represented by a block diagram in FIG. In Fig. 2, g m1 and g m2 are the indoor transfer functions from the observed noise source 1 to the microphones 7 1 and 7 2 , and g s1 and g s2 are the indoor transfer functions from the speakers 9 1 and 9 2 to the noise control point 5. Function, h 1 ,
h 2 represents a transfer function of the filters 8 1 and 8 2 , and g n represents an indoor transfer function between the noise source 1 and the noise control point 5. Further, when replacing the room transfer function g s1, g s2 and transmission filters 81, 82 function h 1, h 2, FIG. 2 can be modified as Figure 3. Third
From the figure, the error signal y observed at the noise control point 5 at time j
j is y j = {g n − (g m1 g s1 h 1 + g m2 g s2 h 2 )} · x j, where x j is a noise signal generated by the noise source at time j. Therefore, in order to control (set to zero) the error signal y j , the filter 8 1 having the transfer functions h 1 and h 2 satisfying the following expression
8 2 should be realized.
gm1gs1h1+gm2gs2h2=gn …… 式で、 とすると、次式が得られる。g m1 g s1 h 1 + g m2 g s2 h 2 = g n …… Then, the following equation is obtained.
g1h1+g2h2=gn ……′ 従ってg1,g2が共通の零点を持たないよう観測マイクロ
ホン71,72およびスピーカ91,92の配置を定めれば、多入
力多出力逆フィルタ処理定理MINT〔1〕により、誤差信
号yjを零とするフィルタ81,82の存在性が保証され、そ
れらの係数は次のように表される。g 1 h 1 + g 2 h 2 = g n ...... 'Therefore, if the arrangement of observation microphones 7 1 and 7 2 and speakers 9 1 and 9 2 is determined so that g 1 and g 2 do not have a common zero point, many The existence of the filters 8 1 and 8 2 that make the error signal y j zero is guaranteed by the input multi-output inverse filtering theorem MINT [1], and their coefficients are expressed as follows.
z-1:遅れ演算子(unit delay) m:伝達関数gn,g1,g2の次数 i:伝達関数h1,h2の次数 T:転置行列 ところが観測マイクロホン71,72と騒音源1との間の室
内伝達関数gm1,gm2が測定困難となり、式のように該
伝達関数をはっきりと用いてフィルタ81,82を定めるこ
とができない場合も起こり得る。そのような場合を考慮
にいれ、この発明の騒音制御装置では、以下の方法を用
いてフィルタ81,82を求めている。 z -1 : delay operator (unit delay) m: order of transfer function g n , g 1 , g 2 i: order of transfer function h 1 , h 2 T: transposed matrix However, observation microphones 7 1 , 7 2 and noise It may happen that the room transfer functions g m1 and g m2 to and from the source 1 become difficult to measure, and the filters 8 1 and 8 2 cannot be defined using the transfer functions clearly as in the equation. In consideration of such a case, the noise control device of the present invention obtains the filters 8 1 and 8 2 by using the following method.
まず次の条件を仮定する。First, assume the following conditions.
I 騒音源1が放射する騒音信号は定常雑音 II スピーカ91,92から騒音制御点5までの伝達関数
gs1,gs2は既知 ここでIは空調騒音など多くの騒音が定常であると見な
せること、IIはスピーカ91,92からM系列信号等を放射
し、前もって計測しておけば良いこと、から現実的な条
件である。また観測マイクロホン71,72から騒音源1ま
での伝達関数gm1,gm2には特別な拘束条件を与えていな
い。I Noise signal emitted from noise source 1 is steady noise II Transfer function from loudspeakers 9 1 and 9 2 to noise control point 5
g s1 and g s2 are known Here, I can be considered that many noises such as air-conditioning noise are stationary, and II radiates M-sequence signals from speakers 9 1 and 9 2 and should be measured in advance. , Is a realistic condition. Further, no special constraint condition is given to the transfer functions g m1 and g m2 from the observation microphones 7 1 and 7 2 to the noise source 1.
さて第3図で時刻jにおけるh1,h2の入力信号x1j,x
2jは、既知の伝達関数gs1,gs2と室内反射音が付加され
た観測マイクロホン71,72が受音した騒音信号u1,u2との
畳み込み演算で得られる。すると時刻jに騒音制御点5
で観測される誤差信号yjは次式で与えられる。Now, in FIG. 3, input signals x 1j , x of h 1 , h 2 at time j
2j is obtained by the convolution operation of the known transfer functions g s1 and g s2 and the noise signals u 1 and u 2 received by the observation microphones 7 1 and 7 2 to which the indoor reflected sound is added. Then, at time j, noise control point 5
The error signal y j observed at is given by the following equation.
yj=rj−HT X j …… 但し これにより誤差パワーの期待値E{yj 2}を最小にする
フィルタ係数〔H1 TH2 T〕Tは次のように求められる。y j = r j −H T X j …… Accordingly, the filter coefficient [H 1 T H 2 T ] T that minimizes the expected value E {y j 2 } of the error power is obtained as follows.
〔H1 TH2 T〕T=E{X j X j T}-1・E{X jrj} …… フィルタ係数〔H1 TH2 T〕Tの具体的な計算方法には、十
分な時間観測したX j,rjより、式を直接解く方法の
他、逐次近似アルゴリズムに基づく繰り返し演算を用
い、時刻jにおけるフィルタ係数〔H1jTH2jT〕Tを式
の解に収束させられる方法もある。逐次近似アルゴリズ
ムとしては、例えばLMSアルゴリズム〔2〕 Hj+1=Hj+αX jyj ……a 但し α:収束係数 やや計算量は多くなるが、収束の速い再帰最小二乗アル
ゴリズム〔3〕 Hj+1=Hj+μψj -1 X jyj ……b ψj+1=αψj+X j X T ……c 但し μ,α:係数 等を利用することができる。[H 1 T H 2 T ] T = E { X j X j T } −1 · E { X j r j } ... The specific calculation method of the filter coefficient [H 1 T H 2 T ] T is as follows. In addition to the method of directly solving the expression from X j , r j observed for a sufficient time, the filter coefficient [H 1 j T H 2 j T ] T at time j is calculated using the iterative calculation based on the successive approximation algorithm. There is also a method that can be converged to. As the successive approximation algorithm, for example, the LMS algorithm [2] H j + 1 = H j + α X j y j. alpha: high computational convergence factor is increased, fast recursive least squares algorithm convergent [3] H j + 1 = H j + μψ j -1 X j y j ...... b ψ j + 1 = αψ j + X j X T ...... c however μ, α: Coefficients can be used.
次に式で与えられる解〔H1 TH2 T〕Tが式の解と一致
することを確かめよう。X j =〔G1G2〕TXj …… であることを用いて式を変形すると、次の式が得ら
れる。Next, make sure that the solution [H 1 T H 2 T ] T given by the equation agrees with the solution of the equation. When the equation is modified by using the fact that X j = [G 1 G 2 ] T X j ..., The following equation is obtained.
〔H1 TH2 T〕T=E{〔G1G2〕TXjXj T〔G1G2〕}-1 ・E{〔G1G2〕TXjXj TR} …… 更に仮定Iを考慮して上式を変形すると、 〔H1 TH2 T〕T=〔G1G2〕-1R ……′ となり、式で表される解に一致する。従って騒音制御
点5で観測される誤差信号yjを零にするフィルタ81,82
は、式の示すフィルタ係数を持つことが確かめられ
た。[H 1 T H 2 T ] T = E {[G 1 G 2 ] T X j X j T [G 1 G 2 ]} -1 · E {[G 1 G 2 ] T X j X j T R} ...... Further, if the above formula is modified in consideration of the assumption I, [H 1 T H 2 T ] T = [G 1 G 2 ] -1 R ...... ', which agrees with the solution represented by the formula. Therefore, the filters 8 1 , 8 2 that make the error signal y j observed at the noise control point 5 zero
Was confirmed to have the filter coefficient shown by the formula.
即ちこの発明の騒音制御装置によれば、観測マイクロホ
ン71,72が騒音源近傍に設置できず、このため受音する
騒音信号に室内反射音が付加されてる場合にも、騒音制
御点5で観測される室内騒音を制御することができるの
である。That is, according to the noise control device of the present invention, even when the observation microphones 7 1 and 7 2 cannot be installed in the vicinity of the noise source, and therefore the indoor reflected sound is added to the received noise signal, the noise control point 5 It is possible to control the indoor noise observed at.
以上、室内の1点で観測される騒音を制御する場合を例
に説明したが、室内の複数点を対象とする騒音制御の場
合でも、「従来の技術」及び「課題を解決するための手
段」で述べたのと全く同様の説明ができる。In the above, the case of controlling the noise observed at one point in the room has been described as an example, but even in the case of noise control targeting a plurality of points in the room, the “prior art” and the “means for solving the problem” It can be explained exactly as described in ".
「実施例」 この発明の実施例について第4図を参照しながら説明す
る。図において騒音源1から放射される騒音信号は室内
反射音が付加された後、騒音制御ユニット61〜6k(k=
2,3……)の構成要素である観測マイクロホン71〜7kお
よび騒音ユニット数未満の数の騒音制御点51〜5Lに設置
された騒音制御点マイクロホン1001〜100Lで受音され
る。観測マイクロホン71〜7kは受音信号をディジタル信
号u1〜ukに変換し、それぞれの信号をフィルタ81〜8kお
よび畳み込み演算器101へ出力する。フィルタ81〜8kで
はu1〜ukとフィルタ係数H1〜Hkとの畳み込み演算によっ
て騒音制御信号s1〜skが合成される。フィルタ81〜8kに
おいて合成された騒音制御信号s1〜skはスピーカ91〜9k
へ入力され、アナログ信号に変換された後、これらスピ
ーカから室内へ放射され、騒音源1から放射される騒音
信号との差信号で表される誤差信号となり、騒音制御点
マイクロホン1001〜100Lで観測される。畳み込み演算器
101では予め測定したスピーカ91〜9kから騒音制御点マ
イクロホン1001〜100Lまでの伝達関数gVW(v=1,2……
k;w=1,2……L)と時刻jにおける観測マイクロホン71
〜7kの出力信号u1〜ukとから、 xVW=gVW・uV …… 但し uV:第v番目の騒音制御ユニットを構成する観測マイク
ロホンの出力信号 が合成される。そして上記畳み込み演算器101の出力信
号XXjと騒音制御点マイクロホン1001〜100Lで観測さ
れ、ディジタル化された誤差信号Yjはフィルタ係数演算
器102に入力される。フィルタ係数演算器102では「課題
を解決するための手段」で述べた逐次近似アルゴリズム
aを多点制御用に拡張した次のアルゴリズム〔4〕 Hj+1=Hj+αXXjYj …… 但し、 Hj:時刻jにおけるフィルタ81〜8kの係数 i:フィルタ81〜8kの次数 y1(j)…yL(j):時刻jに観測された誤差信号 α:収束係数 を実行し、時刻j+1に対するフィルタ81〜8kの係数を
定めた後、各フィルタの係数設定入力に係数を送り出
す。フィルタ81〜8kの係数は、上記,式の処理を繰
り返すことにより、騒音制御点マイクロホン1001〜100L
で観測される誤差信号パワーの期待値を零とする係数へ
と収束する。"Embodiment" An embodiment of the present invention will be described with reference to FIG. After the noise signal emitted from a noise source 1, the room reflected sound are added in the figure, the noise control unit 6 1 ~6 k (k =
2, 3 ...) Observation microphones 7 1 to 7 k and noise control points less than the number of noise units 5 1 to 5 L installed at noise control points microphones 100 1 to 100 L To be done. The observation microphones 7 1 to 7 k convert the received signals into digital signals u 1 to u k and output the respective signals to the filters 8 1 to 8 k and the convolution calculator 101. In the filters 8 1 to 8 k , the noise control signals s 1 to s k are combined by the convolution operation of u 1 to u k and the filter coefficients H 1 to H k . Noise control signal s 1 ~s k synthesized in the filter 8 1 to 8 k is the speaker 9 1 to 9 k
After being input to the analog signal and converted into an analog signal, it is radiated into the room from these speakers and becomes an error signal represented by the difference signal with the noise signal radiated from the noise source 1. The noise control point microphones 100 1 to 100 L Observed in. Convolution calculator
The transfer function g VW from the speaker 9 1 to 9 k measured in advance in 101 to the noise control point microphone 100 1 ~100 L (v = 1,2 ......
k; w = 1,2 ... L) and observation microphone 7 1 at time j
And an output signal u 1 ~u k of ~7 k, x VW = g VW · u V ...... where u V: output signal of the observation microphone constituting the first v-th noise control unit is synthesized. The output signal XX j of the convolution calculator 101 and the error signal Y j observed and digitized by the noise control point microphones 100 1 to 100 L are input to the filter coefficient calculator 102. In the filter coefficient calculator 102, the following algorithm [4] obtained by expanding the successive approximation algorithm a described in “Means for solving the problem” for multipoint control [4] H j + 1 = H j + α XX j Y j , H j : Coefficients of filters 8 1 to 8 k at time j i: order of the filter 8 1 ~8 k y 1 (j) ... y L (j): Error signal observed at time j α: Convergence coefficient is executed, coefficients of filters 8 1 to 8 k for time j + 1 are determined, and then coefficient setting input of each filter Send the coefficient to. The coefficients of the filters 8 1 to 8 k are set to the noise control point microphones 100 1 to 100 L by repeating the above process.
It converges to a coefficient that makes the expected value of the error signal power observed at 0 be zero.
即ちこの発明の騒音制御装置を用いれば、騒音源近傍に
騒音信号を観測するマイクロホンが設置できず、受音信
号に室内反射音が加わる場合にも騒音制御点において観
測される室内騒音を制御することが可能となる。That is, if the noise control device of the present invention is used, a microphone for observing a noise signal cannot be installed in the vicinity of the noise source, and the indoor noise observed at the noise control point is controlled even when the indoor reflected sound is added to the received signal. It becomes possible.
なお本実施例においては、スピーカ91〜9kから放射され
る騒音制御信号が観測マイクロホン71〜7kで受音される
ことによって生じる音響結合(ハウリング)の除去手段
を組み込んでいない。しかしスピーカ91〜9kと観測マイ
クロホン71〜7kとの間の室内伝達関数が測定可能である
ことを考えれば、上記音響結合が音響用エコーキャンセ
ラを用いることにより除去できることは明らかであろ
う。It should be noted that the present embodiment does not incorporate means for removing acoustic coupling (howling) that occurs when the noise control signals radiated from the speakers 9 1 to 9 k are received by the observation microphones 7 1 to 7 k . However, considering that the room transfer function between the speakers 9 1 to 9 k and the observation microphones 7 1 to 7 k can be measured, it is clear that the above acoustic coupling can be eliminated by using the acoustic echo canceller. Let's do it.
「発明の効果」 以上のようにこの発明の騒音制御装置によれば、騒音源
近傍に騒音信号を受音する観測マイクロホンの設置でき
る場合は勿論、騒音源の近傍に騒音信号を受音する観測
マイクロホンが設置できず、受音信号に室内反射音が付
加される場合でも、騒音制御点において観測される室内
騒音を制御することができる。[Advantages of the Invention] As described above, according to the noise control device of the present invention, not only when an observation microphone that receives a noise signal can be installed near the noise source, but also when an observation microphone that receives a noise signal is detected near the noise source. Even if the microphone cannot be installed and the indoor reflected sound is added to the received signal, the indoor noise observed at the noise control point can be controlled.
「文 献」 〔1〕三好、金田:音場の逆フィルタ問題、第2回ディ
ジタル信号処理シンポジウム講演論文集B−1−1(19
87) 〔2〕ヘイキン著、武部訳:適応フィルタ入門、現代工
学社、p.108(1987) 〔3〕ヘイキン著、武部訳:適応フィルタ入門、現代工
学社、p.140(1987) 〔4〕S.J.Elliott et al.:A Multiple Error LMS Algo
rithm and Its Application to the Active Control of
Sound and Vibration,IEEE trans.ASSP,vol.35,pp1423
−1434(1987)"Documents" [1] Miyoshi, Kanada: Inverse filter problem of sound field, Proc. Of the 2nd Digital Signal Processing Symposium B-1-1 (19
87) [2] Haykin, Translated by Takebe: Introduction to Adaptive Filters, Contemporary Engineering, p.108 (1987) [3] Haykin, Translated by Takebe: Introduction to Adaptive Filters, Contemporary Engineering, p.140 (1987) [4] ] SJ Elliott et al .: A Multiple Error LMS Algo
rithm and Its Application to the Active Control of
Sound and Vibration, IEEE trans.ASSP, vol.35, pp1423
-1434 (1987)
第1図は室内の1点において騒音制御を行う場合に用い
られるこの発明の騒音制御装置の構成図、第2図は第1
図に示すこの発明の騒音制御装置のブロック図、第3図
は第2図において伝達関数gs1,gs2とh1,h2を入れ換えて
得られるブロック図、第4図は室内の複数点において騒
音制御を行う場合に用いられるこの発明の騒音制御装置
の一実施例を示す図、第5図は室内の1点において騒音
制御を行う場合に用いられる従来の騒音制御装置の構成
図である。 1:騒音源、2:従来の騒音制御装置で用いられている観測
マイクロホン、31,32:従来の騒音制御装置で用いられて
いるフィルタ、41,42:従来の騒音制御装置で用いられて
いるスピーカ、5,51〜5k:室内に設けられる騒音制御
点、61〜6k:本発明の騒音制御ユニット、71〜7k:騒音制
御ユニットを構成する観測マイクロホン、81〜8k:騒音
制御ユニットを構成するフィルタ、91〜9k:騒音制御ユ
ニットを構成するスピーカ、1001〜100L:騒音制御点に
設置する騒音制御点マイクロホン、101:畳み込み演算
器、102:フィルタ係数演算器、gn:騒音源1から騒音制
御点5までの伝達関数、gm1,gm2:騒音源1から観測マイ
クロホン71,72までの室内伝達関数、gs1,gs2:スピーカ9
1,92から騒音制御点5までの室内伝達関数、h1,h2:フィ
ルタ81,82の伝達関数、N:騒音源、u1〜uk:観測マイクロ
ホン71〜7kが受音し、ディジタル化した騒音信号、s1〜
sk:フィルタ81〜8kで合成した騒音制御信号、yj:室内の
一点において騒音制御を行う場合に騒音制御点で時刻j
に観測される誤差信号、xj:時刻jに騒音源1から放射
される騒音信号、rj:時刻jに騒音源1から騒音制御点
5に放射される騒音信号、x1j,x2j:それぞれ時刻jにお
ける観測マイクロホンの出力信号u1,u2と伝達関数g1,g2
の畳み込みで得られる信号、y1〜yL:この発明の実施例
において、時刻jに騒音制御点マイクロホン1001〜100L
で観測される誤差信号、XXj:畳み込み演算器101の出力
信号、Hj+1:時刻j+1に対するフィルタ81〜8kの係数
(フィルタ係数演算器102の出力信号)。FIG. 1 is a block diagram of a noise control device of the present invention used for noise control at one point in a room, and FIG.
FIG. 3 is a block diagram of the noise control device of the present invention, FIG. 3 is a block diagram obtained by exchanging transfer functions g s1 , g s2 and h 1 , h 2 in FIG. 2, and FIG. FIG. 5 is a diagram showing an embodiment of a noise control device of the present invention used when performing noise control in a room, and FIG. 5 is a configuration diagram of a conventional noise control device used when performing noise control at one point in a room. . 1: Noise source, 2: Observation microphone used in conventional noise control device, 3 1 , 3 2 : Filter used in conventional noise control device, 4 1 , 4 2 : In conventional noise control device Speakers used, 5, 5 1 to 5 k : noise control points provided in the room, 6 1 to 6 k : noise control unit of the present invention, 7 1 to 7 k : observation microphones constituting the noise control unit, 8 1 to 8 k : Filter configuring noise control unit, 9 1 to 9 k : Speaker configuring noise control unit, 100 1 to 100 L : Noise control point microphone installed at noise control point, 101: Convolution calculator , 102: filter coefficient calculator, g n : transfer function from noise source 1 to noise control point 5, g m1 , g m2 : indoor transfer function from noise source 1 to observation microphones 7 1 , 7 2 , g s1 , g s2 : Speaker 9
Room transfer function from 1 , 9 2 to noise control point 5, h 1 , h 2 : transfer function of filter 8 1 , 8 2 , N: noise source, u 1 ~ u k : observation microphone 7 1 ~ 7 k Sound signal received and digitized, s 1 ~
s k : noise control signal synthesized by the filters 8 1 to 8 k , y j : time j at the noise control point when noise control is performed at one point in the room
Error signal observed at x j : noise signal emitted from noise source 1 at time j, r j : noise signal emitted from noise source 1 to noise control point 5 at time j, x 1j , x 2j : Output signals u 1 and u 2 of the observation microphone and transfer functions g 1 and g 2 at time j, respectively.
, Y 1 to y L : noise control point microphones 100 1 to 100 L at time j in the embodiment of the present invention.
XX j : output signal of convolutional computing unit 101, H j + 1 : coefficients of filters 8 1 to 8 k with respect to time j + 1 (output signal of filter coefficient computing unit 102).
Claims (1)
マイクロホンと、 その観測マイクロホンが受音した騒音信号を用いて上記
室内騒音を抑圧する騒音制御信号を合成するフィルタ
と、 上記騒音信号を室内に放射するスピーカとから成る騒音
制御ユニットの複数個と、 これら騒音ユニットの数未満の数の騒音制御点それぞれ
に設けられた騒音制御点マイクロホンと、 上記各騒音観測マイクロホンの出力信号と、上記各スピ
ーカから上記各騒音制御点までの伝達関数とを畳み込み
演算する畳み込み演算手段と、 その畳み込み演算手段の演算結果と、上記各騒音制御点
マイクロホンの出力とを入力して、これら騒音制御点マ
イクロホンの出力の二乗和の最小化(零化)を行う逐次
近似アルゴリズムにより上記各フィルタの係数を演算し
上記各フィルタに設定するフィルタ係数演算手段とを具
備することを特徴とする騒音制御装置。1. An observation microphone for observing a noise signal generated by a noise source, a filter for synthesizing a noise control signal for suppressing the indoor noise by using the noise signal received by the observation microphone, and the noise signal A plurality of noise control units each consisting of a speaker radiating into the room, noise control point microphones provided at each of the noise control points less than the number of these noise units, output signals of the noise observation microphones, and Convolution calculation means for performing convolution calculation of the transfer function from each speaker to each noise control point, the calculation result of the convolution calculation means, and the output of each noise control point microphone are input, and these noise control point microphones are input. The coefficient of each filter is calculated by the successive approximation algorithm that minimizes (zeroizes) the sum of squares of the output of Noise control apparatus characterized by comprising a filter coefficient calculating means for setting the filter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP63058031A JPH0769707B2 (en) | 1988-03-11 | 1988-03-11 | Noise control device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP63058031A JPH0769707B2 (en) | 1988-03-11 | 1988-03-11 | Noise control device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH01231597A JPH01231597A (en) | 1989-09-14 |
JPH0769707B2 true JPH0769707B2 (en) | 1995-07-31 |
Family
ID=13072579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP63058031A Expired - Fee Related JPH0769707B2 (en) | 1988-03-11 | 1988-03-11 | Noise control device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH0769707B2 (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2674252B2 (en) * | 1989-12-29 | 1997-11-12 | 日産自動車株式会社 | Active noise control device |
JP2541020B2 (en) * | 1991-01-25 | 1996-10-09 | 日立プラント建設株式会社 | Muffler for radiated noise in openings |
JP2796001B2 (en) * | 1991-10-05 | 1998-09-10 | 三菱電機株式会社 | Silencer |
JPH0830278A (en) * | 1994-07-14 | 1996-02-02 | Honda Motor Co Ltd | Active vibration control device |
JPH0830902A (en) * | 1994-07-14 | 1996-02-02 | Sony Corp | Noise reduction device for camcorder |
DE10362073A1 (en) * | 2003-11-06 | 2005-11-24 | Herbert Buchner | Apparatus and method for processing an input signal |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5926039A (en) * | 1982-08-03 | 1984-02-10 | Komatsu Ltd | Method and apparatus for deciding material quality of cv graphite cast iron |
JPS62209295A (en) * | 1986-03-10 | 1987-09-14 | 日立プラント建設株式会社 | Multipoint electronic silencing system |
-
1988
- 1988-03-11 JP JP63058031A patent/JPH0769707B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH01231597A (en) | 1989-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0195416B1 (en) | Inverse control system | |
Kuo et al. | Nonlinear adaptive bilinear filters for active noise control systems | |
US8462962B2 (en) | Sound processor, sound processing method and recording medium storing sound processing program | |
JPH1152988A (en) | Control method of adaptive array and adaptive array device | |
US10037755B2 (en) | Method and system for active noise reduction | |
US7181026B2 (en) | Post-processing scheme for adaptive directional microphone system with noise/interference suppression | |
JP3359460B2 (en) | Adaptive filter and echo canceller | |
JPWO2018163810A1 (en) | Signal processing apparatus and method, and program | |
US6381272B1 (en) | Multi-channel adaptive filtering | |
WO2007123051A1 (en) | Adaptive array controlling device, method, program, and adaptive array processing device, method, program | |
Shi et al. | A partial-update minimax algorithm for practical implementation of multi-channel feedforward active noise control | |
JP2007096389A (en) | Recursive noise canceller | |
CN113114865A (en) | Combined function linkage type kernel self-response nonlinear echo cancellation method | |
Wang et al. | Residue-to-binary number converters for three moduli sets | |
JPH0769707B2 (en) | Noise control device | |
JPH0511771A (en) | Noise control device | |
JP3864914B2 (en) | Echo suppression device | |
JP4041770B2 (en) | Acoustic echo cancellation method, apparatus, program, and recording medium | |
US20220272450A1 (en) | Design method for feedforward active noise control system using analog filter | |
Sachau et al. | Real-time implementation of the frequency-domain FxLMS algorithm without block delay for an adaptive noise blocker | |
JP3616341B2 (en) | Multi-channel echo cancellation method, apparatus thereof, program thereof, and recording medium | |
JP2529464B2 (en) | Noise canceller | |
EP2257082A1 (en) | Background noise estimation in a loudspeaker-room-microphone system | |
JPS63311396A (en) | Electronic muffling system | |
Lu et al. | Power-of-two quantizer FLANN filter for nonlinear active noise control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |