JPH08328569A - Fuzzy active noise erase device - Google Patents

Fuzzy active noise erase device

Info

Publication number
JPH08328569A
JPH08328569A JP7131648A JP13164895A JPH08328569A JP H08328569 A JPH08328569 A JP H08328569A JP 7131648 A JP7131648 A JP 7131648A JP 13164895 A JP13164895 A JP 13164895A JP H08328569 A JPH08328569 A JP H08328569A
Authority
JP
Japan
Prior art keywords
fuzzy
noise
error
signal
membership 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.)
Granted
Application number
JP7131648A
Other languages
Japanese (ja)
Other versions
JP3445872B2 (en
Inventor
Tadashi Ohashi
正 大橋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP13164895A priority Critical patent/JP3445872B2/en
Publication of JPH08328569A publication Critical patent/JPH08328569A/en
Application granted granted Critical
Publication of JP3445872B2 publication Critical patent/JP3445872B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Soundproofing, Sound Blocking, And Sound Damping (AREA)
  • Filters That Use Time-Delay Elements (AREA)
  • Image Processing (AREA)

Abstract

PURPOSE: To erase a noise flexibly by estimating error scatter without a C filter and by fuzzy control, related to a noise erase device actively erasing the noise such as an active noise controller, etc., applying a fuzzy theory. CONSTITUTION: This device is the noise erase device applying a noise erase signal to the noise and erasing the noise. Then, the noise erase device is the device provided with an error scatter system that an error erasing the noise is scattered timewisely, and is provided with a fuzzy operation part 8 performing a fuzzy inference for the error scattered by the error scatter system and an input signal and an adaptive filter 2 inputting the operation result of the fuzzy operation part 8 and the input signal and outputting the noise erase signal, and is constituted so that the optimum coefficient of the adaptive filter 2 is obtained based on the operation result of the fuzzy operation part 8, and the noise erase signal is generated.

Description

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

【0001】[0001]

【産業上の利用分野】本発明は,ファジィ理論を応用し
た能動騒音制御装置等の能動ノイズ消去装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an active noise canceling device such as an active noise control device to which a fuzzy theory is applied.

【0002】特に,誤差検出部から適応フィルタに到る
伝達系に誤差を時間軸方向に散乱させる伝達系(誤差散
乱系)をもつ能動騒音制御装置,及び電話機のエコーキ
ャンセラー,画像のゴースト,制振制御,耐震制御等の
能動ノイズ消去装置に関するものである。
In particular, an active noise control device having a transmission system (error scattering system) that scatters errors in the time axis direction from the error detector to the adaptive filter, and an echo canceller for telephones, image ghosting, and control The present invention relates to an active noise canceller such as vibration control and seismic control.

【0003】[0003]

【従来の技術】図12は従来のノイズ消去装置の構成を
示す。図12において,210は未知系である。
2. Description of the Related Art FIG. 12 shows the structure of a conventional noise canceller. In FIG. 12, 210 is an unknown system.

【0004】211は適応フィルタである。212はノ
イズ消去フィルタであって,ノイズを消去する信号を生
成するものである。
Reference numeral 211 is an adaptive filter. A noise elimination filter 212 generates a signal for eliminating noise.

【0005】213は重み制御部であって,正規化最小
自乗法(NLMS)によりノイズ消去フィルタ212の
最適フィルタ係数を求めるものである。214はCフィ
ルタであって,誤差散乱系の伝達関数を推定するフィル
タである。
A weight control unit 213 obtains an optimum filter coefficient of the noise elimination filter 212 by the normalized least squares method (NLMS). A C filter 214 is a filter for estimating the transfer function of the error scattering system.

【0006】215は加算部であって,未知系210を
通過する信号h(i)に外乱(ノイズ)Nj を加算する
ものである。216は加算部であって,未知系210を
通過してノイズを加算された信号にノイズ消去信号を加
算するものである。
Reference numeral 215 denotes an adder which adds a disturbance (noise) N j to the signal h (i) passing through the unknown system 210. Reference numeral 216 denotes an adder, which adds a noise elimination signal to the signal that has passed through the unknown system 210 and added with noise.

【0007】217は誤差散乱系であって,誤差検出部
(図示せず)から適応フィルタ211にいたる伝達系で
あって,時間軸に対して散乱をもつものである。Xj
時点jの入力信号である。
Reference numeral 217 denotes an error scattering system, which is a transmission system from an error detection section (not shown) to the adaptive filter 211, and has scattering on the time axis. X j is the input signal at time j.

【0008】h(i)は未知系210を伝わる信号であ
る。Hj は時点jの消去フィルタの伝達関数である。H
j+1 は時点j+1の消去フィルタの伝達関数である。
H (i) is a signal transmitted through the unknown system 210. H j is the transfer function of the cancellation filter at time j. H
j + 1 is the transfer function of the cancellation filter at time j + 1.

【0009】Yj はCフィルタ214の出力である。G
j はノイズ消去フィルタHj (212)の出力である。
j はシステムの出力であり,誤差散乱系の入力となる
ものである。
Y j is the output of the C filter 214. G
j is the output of the noise elimination filter H j (212).
E j is the output of the system and is the input of the error scattering system.

【0010】ej は誤差散乱系の出力である。Nj はノ
イズ信号である。図12のシステムの伝達関数は次のよ
うに表される。
E j is the output of the error scattering system. N j is a noise signal. The transfer function of the system of FIG. 12 is expressed as:

【0011】Hj+1 =Hj +Kej j /ΣYj -t(Σ
はt=0〜N−1の和)となる。但し,Nはノイズ消去
フィルタ212のタップ数,Kはステップゲインであ
る。
H j + 1 = H j + Ke j Y j / ΣY j -t
Is the sum of t = 0 to N−1). However, N is the number of taps of the noise elimination filter 212, and K is a step gain.

【0012】図12の構成の動作を説明する。入力信号
j は,信号伝達系h(i)を伝わって加算部215で
ノイズNj を加算する。一方,Cフィルタ(散乱推定フ
ィルタ)214は入力信号Xj に対して誤差散乱系を推
定し,Yj を出力する。重み制御部213はYj とej
とにより最適フィルタ係数を求める。時刻jのノイズ消
去フィルタ212の伝達関数をHj としたとき,Hj+1
=Hj +Kej j /ΣYj -t(Σはt=0〜N−1の
和)で計算される。
The operation of the configuration of FIG. 12 will be described. The input signal X j is transmitted through the signal transmission system h (i) and the noise N j is added by the addition unit 215. Meanwhile, the C filter (scatter estimation filter) 214 estimates the error scattering based on the input signals X j, outputs a Y j. The weight control unit 213 uses Y j and e j.
The optimum filter coefficient is obtained by and. When the transfer function of the noise elimination filter 212 at time j is H j , H j + 1
= H j + Ke j Y j / ΣY j −t (Σ is the sum of t = 0 to N−1).

【0013】ノイズ消去フィルタ212はXj (j=0
〜N−1)とHj+1 によりたたみこみの計算を行いGj
を求める。加算部216により,ノイズを含む信号にノ
イズを消去する信号Gj を加算し,ノイズを除去した信
号Ej を生成する。
The noise elimination filter 212 uses X j (j = 0).
~ N-1) and H j + 1 to calculate the convolution, G j
Ask for. The adder 216 adds the noise-eliminating signal G j to the noise-containing signal to generate a noise- free signal E j .

【0014】さらに次の時点j+1で同様の処理を行
い,Gj+1 を求める。
Further, at the next time point j + 1, similar processing is performed to obtain G j + 1 .

【0015】[0015]

【発明が解決しようとする課題】従来のノイズ消去装置
は,誤差散乱系を推定するCフィルタを設け,学習同定
することにより適応フィルタの最適フィルタ係数を求め
るようにしていたので,多数のDSPを必要とし,処理
を長時間必要とした。
In the conventional noise canceller, a C filter for estimating an error scattering system is provided, and the optimum filter coefficient of the adaptive filter is obtained by learning identification. Therefore, a large number of DSPs are used. It took a long time to process.

【0016】本発明は,Cフィルタを必要としないで,
誤差散乱を推定し,ファジィ制御により柔軟にノイズ消
去を行うことのできるファジィ能動ノイズ消去装置を提
供することを目的とする。
The present invention does not require a C filter,
It is an object of the present invention to provide a fuzzy active noise canceller capable of estimating error scattering and flexibly canceling noise by fuzzy control.

【0017】[0017]

【課題を解決するための手段】本発明は,ノイズに対し
て能動的にノイズ消去信号を印加してノイズを消去する
能動ノイズ消去装置において,ノイズ消去装置はノイズ
を消去した誤差が時間的に散乱する誤差散乱系を持つも
のであって,該誤差散乱系により散乱された誤差と入力
信号に対してファジィ推論をするファジィ演算部とファ
ジィ演算部の演算結果と入力信号を入力してノイズ消去
信号を出力する適応フィルタとを備え,ファジィ演算部
の演算結果に基づいて適応フィルタの最適係数を求め,
能動的にノイズ消去信号を生成するようにした。
SUMMARY OF THE INVENTION The present invention is an active noise canceller for actively applying a noise canceling signal to noise to cancel the noise. A fuzzy arithmetic unit having an error scattering system that scatters and performing fuzzy inference on the error and the input signal scattered by the error scattering system, and the operation result of the fuzzy arithmetic unit and the input signal are input to cancel noise. An adaptive filter that outputs a signal, and obtains the optimum coefficient of the adaptive filter based on the calculation result of the fuzzy calculation unit.
A noise canceling signal is actively generated.

【0018】図1は本発明の基本構成を示す図である。
図1において,1は未知系である。
FIG. 1 is a diagram showing the basic configuration of the present invention.
In FIG. 1, 1 is an unknown system.

【0019】2は適応フィルタである。5はファジィ演
算部であって,ファジィ演算部5の入力値,出力値のメ
ンバシップ関数,メンバシップ関数を選択するファジィ
制御規則に基づいてファジィ推論をし,出力値を決定す
るものである。
Reference numeral 2 is an adaptive filter. Reference numeral 5 denotes a fuzzy arithmetic unit which performs fuzzy inference based on a fuzzy control rule for selecting a membership function of an input value and an output value of the fuzzy arithmetic unit 5 and a membership function to determine an output value.

【0020】6は誤差散乱系である。ファジィ演算部5
において,11は誤差散乱演算部であって,誤差散乱系
のインパルス応答の時間軸上の標本値a(t)によりY
j を求めるものである。例えば,Yj =Σa(t)Xj
と表すことができる。Σはtについての和であり,t=
0〜T−1(Tはタップ数に相当する)である。
Reference numeral 6 is an error scattering system. Fuzzy operation unit 5
In the figure, 11 is an error scattering calculation unit, and Y is determined by a sample value a (t) on the time axis of the impulse response of the error scattering system
It seeks j . For example, Y j = Σa (t) X j
It can be expressed as. Σ is the sum of t, and t =
0 to T-1 (T corresponds to the number of taps).

【0021】16は誤差散乱推論部であって,誤差散乱
系のインパルス応答の時間軸上の標本値a(t)のファ
ジィ推論をするものである。17は評価部であって,フ
ァジィ演算によりノイズ消去信号Gj が最適に求められ
たかを評価し,適切な出力が得られていない場合にファ
ジィ制御規則,メンバシップ関数の選択等を行うもので
ある。
Reference numeral 16 denotes an error scattering inference unit, which performs fuzzy inference on the sample value a (t) on the time axis of the impulse response of the error scattering system. An evaluation unit 17 evaluates whether or not the noise elimination signal G j is optimally obtained by fuzzy operation, and selects a fuzzy control rule, membership function, etc. when an appropriate output is not obtained. is there.

【0022】21はメンバシップ関数保持部である。2
2はファジィ制御規則保持部である。Xj は入力信号で
ある。
Reference numeral 21 is a membership function holding unit. Two
Reference numeral 2 is a fuzzy control rule holding unit. X j is an input signal.

【0023】Gj は適応フィルタの出力であって,ノイ
ズ消去信号である。Ej は出力である。ej は誤差散乱
系6の出力である。
G j is the output of the adaptive filter and is a noise elimination signal. E j is the output. e j is the output of the error scattering system 6.

【0024】Hj ,Hj+1 は適用フィルタの伝達関数で
ある。
H j and H j + 1 are transfer functions of the applied filter.

【0025】[0025]

【作用】図1の本発明の基本構成の動作を説明する。未
知系1を通過したノイズを含む入力信号Xj は加算部8
でノイズ消去信号G j と加算され,Ej を出力する。時
刻jで誤差散乱系6はej を出力する。
The operation of the basic configuration of the present invention shown in FIG. 1 will be described. Not yet
Input signal X containing noise that has passed through intelligence system 1jIs the addition unit 8
Noise elimination signal G jIs added and EjIs output. Time
At time j, the error scattering system 6 is ejIs output.

【0026】ファジィ演算部5は入力値Xj ,ej を入
力し,Xj ,ej ,Kおよび誤差散乱の標本値a(t)
についてもそれぞれのメンバシップ関数により,例え
ば,MIN−MAX法等によりファジィ推論をする。そ
して,ファジィ推論の結果に基づいて重心法等により非
ファジィ化し,適応フィルタのフィルタ係数の制御デー
タ(ステップゲインK,誤差散乱の推定値Yj ,入力値
j ,誤差ej 等)を求める。
The fuzzy operation section 5 inputs the input values X j , e j, and inputs the sample values a (t) of X j , e j , K and error scattering.
For each of these, fuzzy inference is performed by the respective membership functions, for example, the MIN-MAX method. Then, based on the result of the fuzzy inference, the control data (step gain K, estimated value Y j of error scattering, input value X j , error e j, etc.) of the filter coefficient of the adaptive filter is defuzzified by the centroid method or the like. .

【0027】評価部17はej を評価し,ノイズ消去信
号Gj が適切でない場合にファジィ規則およびメンバシ
ップ関数の変更等を行う。
The evaluation unit 17 evaluates e j and changes the fuzzy rule and the membership function if the noise elimination signal G j is not appropriate.

【0028】[0028]

【実施例】図2は本発明の実施例1を示す図である。図
2において,1は未知系である。
EXAMPLE 1 FIG. 2 is a diagram showing Example 1 of the present invention. In FIG. 2, 1 is an unknown system.

【0029】2は適応フィルタである。3はノイズ消去
フィルタである。4は重み制御部であって,NLMS法
によりフィルタ係数を決定するものである。
Reference numeral 2 is an adaptive filter. 3 is a noise elimination filter. A weight control unit 4 determines the filter coefficient by the NLMS method.

【0030】5はファジィ演算部であって,ファジィ演
算部5の入力値,出力値に対するメンバシップ関数,メ
ンバシップ関数を選択するファジィ規則に基づいてファ
ジィ推論をし,非ファジィ化をするものである。
Reference numeral 5 denotes a fuzzy arithmetic unit, which performs fuzzy inference based on a fuzzy rule for selecting a membership function and a membership function with respect to an input value and an output value of the fuzzy arithmetic unit 5 for defuzzification. is there.

【0031】6は誤差散乱系である。ファジィ演算部5
において,10は非ファジィ化部であって,ファジィ推
論の結果得られた統合結果をもとに重み制御部の演算デ
ータを算出するものである。
Reference numeral 6 is an error scattering system. Fuzzy operation unit 5
In the above, reference numeral 10 denotes a defuzzification unit, which calculates the operation data of the weight control unit based on the integrated result obtained as a result of the fuzzy inference.

【0032】11は誤差散乱演算部であって,誤差散乱
系のインパルス応答の時間軸上の標本値a(t)により
j を求めるものである。Yj =Σa(t)Xj であ
る。12はEd 演算部であって,Ed =Yj /ΣYj 2
を算出するものである。誤差散乱を表す関数は正規分
布,ポアソン分布,二項分布等である。
Reference numeral 11 denotes an error scattering calculation unit, which calculates Y j from the sample value a (t) on the time axis of the impulse response of the error scattering system. Y j = Σa (t) X j . Reference numeral 12 is an E d calculation unit, and E d = Y j / ΣY j 2
Is calculated. Functions that represent error scattering are normal distribution, Poisson distribution, binomial distribution, and the like.

【0033】13はその他の数値演算部であって,その
他の値(Xj ,ej ,K等)を算出するものである。1
5はファジィ推論部である。
Reference numeral 13 denotes another numerical operation unit for calculating other values (X j , e j , K, etc.). 1
5 is a fuzzy inference unit.

【0034】16は誤差散乱推論部であって,誤差散乱
系のインパルス応答の時間軸上の標本値a(t)のファ
ジィ推論をするものである。17は評価部であって,フ
ァジィ演算によりノイズ消去信号Gj が最適に求められ
たかを評価し,適切な出力が得られていない場合にファ
ジィ制御規則,メンバシップ関数の選択等を行うもので
ある。
Reference numeral 16 denotes an error scattering inference unit, which performs fuzzy inference on the sample value a (t) on the time axis of the impulse response of the error scattering system. An evaluation unit 17 evaluates whether or not the noise elimination signal G j is optimally obtained by fuzzy operation, and selects a fuzzy control rule, membership function, etc. when an appropriate output is not obtained. is there.

【0035】21はメンバシップ関数保持部である。2
2はファジィ制御規則保持部である。Xj は入力信号で
ある。
Reference numeral 21 is a membership function holding unit. Two
Reference numeral 2 is a fuzzy control rule holding unit. X j is an input signal.

【0036】Nj はノイズである。Gj はノイズ消去フ
ィルタの出力である。Ej は出力である。
N j is noise. G j is the output of the noise elimination filter. E j is the output.

【0037】ej は誤差散乱系6の出力である。Hj
j+1 はノイズ消去フィルタの伝達関数である。Kはス
テップゲインであって,ファジィ演算の結果の統合結果
に基づいて非ファジィ化して得られるものである。
E j is the output of the error scattering system 6. H j ,
H j + 1 is the transfer function of the noise elimination filter. K is a step gain, which is obtained by defuzzification based on the integrated result of the fuzzy calculation results.

【0038】Ed =Yj /ΣYj 2 である。図2の本発
明の実施例1の動作を説明する。図2の系は,入力信号
j が未知系1を伝わり,その過程でノイズNj の影響
を受ける。そして,そのノイズを含む信号に対してノイ
ズNj を消去するノイズ消去信号が加えられその理想的
な誤差信号はEj であるが,誤差散乱の影響を受け,誤
差がej で観測される系を表している。
E d = Y j / ΣY j 2 . The operation of the first embodiment of the present invention shown in FIG. 2 will be described. In the system of FIG. 2, the input signal X j is transmitted through the unknown system 1 and is affected by noise N j in the process. Then, a noise elimination signal for eliminating the noise N j is added to the signal including the noise, and the ideal error signal is E j , but the error is affected by error scattering and the error is observed at e j. It represents the system.

【0039】図2の系において,入力信号Xj は未知系
1を伝わり,加算部7でノイズNjを加算され,加算部
8でノイズ消去信号Gj を加算され,誤差Ej となる。
jは誤差散乱系6を伝播し誤差ej として観測され
る。
In the system of FIG. 2, the input signal X j is transmitted through the unknown system 1, the noise N j is added by the adder 7, and the noise canceling signal G j is added by the adder 8, resulting in an error E j .
E j propagates through the error scattering system 6 and is observed as an error e j .

【0040】ファジィ演算部5は入力値Xj ,ej を入
力し,Xj ,ej ,K,a(t)のメンバシップ関数に
よりそれぞれのファジィ推論をする。例えば,MIN−
MAX法等により推論結果を求める。非ファジィ化部1
0は,推論結果を統合し,例えば,重心法等により
j ,ej ,K,a(t)を求める。そして,誤差散乱
演算部11は,Xj とa(t)により,Yj =Σa
(t)Xj を算出する。Ed 演算部12はその値により
d =Yj /ΣYj 2 を算出する。
The fuzzy operation unit 5 inputs the input values X j and e j , and performs fuzzy inference on the basis of the membership functions of X j , e j , K and a (t). For example, MIN-
The inference result is obtained by the MAX method or the like. Defuzzification unit 1
0 integrates the inference results and obtains X j , e j , K, and a (t) by, for example, the centroid method. Then, the error scattering calculator 11 calculates Y j = Σa from X j and a (t).
(T) Calculate X j . The E d calculation unit 12 calculates E d = Y j / ΣY j 2 from the value.

【0041】そして,重み制御部4は,時刻jのノイズ
消去フィルタ3の伝達関数をHj としたとき,ノイズ消
去フィルタ3の伝達関数Hj+1 をHj+1 =Hj +Kej
dで計算する。さらに,ノイズ消去フィルタ3はXj
(j=0〜N−1)とHj+1によりたたみこみ演算を行
いGj を求める。加算部8はノイズを含む信号にノイズ
消去信号Gj を加算し,ノイズを消去した信号Ej を生
成する。さらに次の時点j+1で同様の処理を行い,E
j+1 を求める。
Then, when the transfer function of the noise canceling filter 3 at time j is H j , the weight controller 4 sets the transfer function H j + 1 of the noise canceling filter 3 to H j + 1 = H j + K e j.
Calculate with E d . Furthermore, the noise elimination filter 3 has X j
Request (j = 0~N-1) and H j + 1 performs a convolution operation by G j. The addition unit 8 adds the noise elimination signal G j to the signal containing noise to generate a noise elimination signal E j . Further, at the next time point j + 1, similar processing is performed, and E
Find j + 1 .

【0042】評価部17は,ej の収束の度合い(収束
の大きさ),有効に消去された周波数の範囲(周波数依
存性),収束の安定性(収束の時間変化)をそれぞれの
メンバシップ関数を基にファジィ推論し,非ファジィ化
して評価する。そして,評価値がよりよくなるように,
メンバシップ制御規則,メンバシップ関数を変更する。
The evaluator 17 determines the degree of convergence of e j (the size of the convergence), the range of the frequency that is effectively erased (frequency dependence), and the stability of the convergence (time change of the convergence) for each membership. Fuzzy inference is performed based on the function, and defuzzification is performed for evaluation. And so that the evaluation value becomes better,
Change membership control rules and membership functions.

【0043】図3は本発明の実施例2であって,ダクト
の騒音を消去する能動騒音制御装置に適用したものであ
る。図3において,2は適応フィルタである。
FIG. 3 shows a second embodiment of the present invention, which is applied to an active noise control system for eliminating noise in a duct. In FIG. 3, reference numeral 2 is an adaptive filter.

【0044】3はノイズ消去フィルタである。4は重み
制御部である。6は誤差散乱系である。
Reference numeral 3 is a noise elimination filter. Reference numeral 4 is a weight control unit. 6 is an error scattering system.

【0045】5はファジィ演算部である。41’はダク
トである。42’はセンサマイクであって,ダクトに入
る騒音を検出するものである。
Reference numeral 5 is a fuzzy operation section. 41 'is a duct. 42 'is a sensor microphone for detecting noise entering the duct.

【0046】43’はエラーマイクであって,ダクトか
ら外に出る騒音(誤差信号)を検出するものである。4
4’は二次音源であって,騒音を消去する音を出力する
スピーカである。
Reference numeral 43 'is an error microphone for detecting noise (error signal) output from the duct. Four
Reference numeral 4'denotes a secondary sound source, which is a speaker that outputs a sound that eliminates noise.

【0047】45’は加算部であって,センサマイク4
2の検出信号と回り込み防止フィルタ46の出力を加算
するものである。46は回り込み防止フィルタであっ
て,二次音源44’から発生する音がセンサマイク42
に伝達する系を模擬し,二次音源44’から発生する音
がセンサマイク42に検出されることによる影響を防止
するものである。
Reference numeral 45 'is an addition unit, which is used for the sensor microphone 4
The detection signal of 2 and the output of the wraparound prevention filter 46 are added. Reference numeral 46 denotes a sneak-in prevention filter, which allows the sound generated from the secondary sound source 44 'to be transmitted to the sensor microphone 42.
The system for transmitting the sound is simulated to prevent the sound generated from the secondary sound source 44 'from being detected by the sensor microphone 42.

【0048】Xsjはセンサマイク42の出力である。X
S は回り込み防止フィルタ46の出力である。Xj は加
算部45の加算結果である。
X sj is the output of the sensor microphone 42. X
S is the output of the wraparound prevention filter 46. X j is the addition result of the addition unit 45.

【0049】Ej はエラーマイク43の出力であって,
誤差散乱系の出力である(前述のe j に相当する)。E
j は誤差散乱系の出力である。
EjIs the output of the error microphone 43,
This is the output of the error scattering system. jEquivalent to). E
jIs the output of the error scattering system.

【0050】K,a(t)はファジィ演算部5の出力で
ある。従来のCフィルタの出力Yj は,ファジィ演算部
において,正規分布等の散乱を表す分布のメンバシップ
関数によりファジィ推論し,非ファジィ化して得られる
a(t)を基に演算により算出する。散乱を表す時間関
数a(t)によりYj=Σa(t)Xj-t (Σはt=0
〜T−1の累積(TはCフィルタのタップ数)))によ
り求める。誤差散乱を表す関数は正規分布,二項分布,
ポアソン分布等である。
K, a (t) are outputs of the fuzzy operation section 5. The output Y j of the conventional C filter is calculated by a fuzzy operation unit based on a (t) obtained by performing fuzzy inference using a membership function of a distribution representing scattering such as a normal distribution and defuzzification. By the time function a (t) representing the scattering, Y j = Σa (t) X jt (Σ is t = 0.
~ Accumulation of T-1 (T is the number of taps of the C filter))). The function that represents the error scattering is the normal distribution, the binomial distribution,
Poisson distribution, etc.

【0051】図3の構成の動作を説明する。ダクト41
に入る騒音の時刻jにおける騒音Xsjはセンサマイク4
2で検出される。加算部45はその信号Xsjと回り込み
防止フィルタ46の出力Xsjとを加算し,信号Xj を生
成する。
The operation of the configuration of FIG. 3 will be described. Duct 41
The noise X sj at the time j of the incoming noise is the sensor microphone 4
Detected in 2. The adder 45 adds the signal X sj and the output X sj of the wraparound prevention filter 46 to generate a signal X j .

【0052】一方,二次音源44’は適応フィルタ2の
出力するGj に基づいて騒音を消去するための音を出力
する。二次音源で発生した音は騒音に加えられ,誤差散
乱系6で散乱されて誤差がエラーマイク43で検出され
る。エラーマイク43は誤差Ej を出力する(前述のe
j に対応する)。
On the other hand, the secondary sound source 44 'outputs a sound for eliminating noise based on G j output by the adaptive filter 2. The sound generated by the secondary sound source is added to the noise, scattered by the error scattering system 6, and the error is detected by the error microphone 43. The error microphone 43 outputs the error E j (e described above).
corresponding to j ).

【0053】ファジィ演算部5はXj ,Ej ,K,a
(t)のメンバシップ関数に基づいてファジィ推論をす
る。そして推論結果を非ファジィ化してK,a(t),
Edを求め,重み制御部4に入力する。Yj = Σa
(t)Yj-t であり,Ed =Yj/ΣYj -i 2 である。
The fuzzy operation unit 5 operates at X j , E j , K, a
Fuzzy inference is performed based on the membership function of (t). Then, the inference result is defuzzified so that K, a (t),
Ed is calculated and input to the weight control unit 4. Y j = Σa
(T) Y jt , and E d = Y j / ΣY j -i 2 .

【0054】重み制御部4はファジィ演算部5の算出し
たEj ,ステップゲインKとEdとにより,時刻j+1
のノイズ消去フィルタ3の伝達関数Hj+1 を求める。H
j+1=Hj +Kej Eである。
The weight controller 4 uses the E j calculated by the fuzzy calculator 5 and the step gains K and Ed to calculate the time j + 1.
Then, the transfer function H j + 1 of the noise elimination filter 3 is obtained. H
j + 1 = H j + Ke j E.

【0055】ノイズ消去フィルタ3は入力Xj (j=0
〜N−1)とHj とにより,たたみ込みの計算を行い,
ノイズ消去信号Gj を出力する。二次音源44’はGj
により消去音を発生する騒音に加えることにより騒音を
消去する。
The noise canceling filter 3 receives the input X j (j = 0).
~ N-1) and H j are used to calculate the convolution,
The noise elimination signal G j is output. The secondary sound source 44 'is G j
The noise is deleted by adding the noise to the generated noise.

【0056】また,ファジィ演算部5において,出力E
j がファジィ推論により評価され,ファジィ規則,メン
バシップを適切なものに変更する。図4は本発明のファ
ジィ演算部の実施例である。
In the fuzzy operation section 5, the output E
j is evaluated by fuzzy reasoning, and the fuzzy rules and membership are changed to appropriate ones. FIG. 4 shows an embodiment of the fuzzy arithmetic unit of the present invention.

【0057】図4において,5はファジィ演算部であ
る。10は非ファジィ化部である。
In FIG. 4, reference numeral 5 is a fuzzy arithmetic unit. Reference numeral 10 is a defuzzification section.

【0058】15はファジィ推論部である。ファジィ推
論部15において,31はXj 推論部であって,Xj
メンバシップ関数とXj の入力信号とによりMIN−M
AX法によりXj のファジィ推論をするものである。
Reference numeral 15 is a fuzzy inference unit. In the fuzzy inference section 15, 31 is X a j inference unit, X j MIN-M by the input signal of the membership function and X j of
Fuzzy inference of X j is performed by the AX method.

【0059】32はEj 推論部であって,Ej のメンバ
シップ関数とEj の入力信号とによりMIN−MAX法
によりEj をファジィ推論をするものである。33はK
推論部であって,Kのメンバシップ関数およびXj とE
j の推論結果に基づいてMIN−MAX法によりK(ス
テップゲイン)をファジィ推論するものである。
[0059] 32 is to be E j inference unit, the fuzzy inference to E j by MIN-MAX method by the input signal of the membership function and E j of E j. 33 is K
An inference unit, which has a membership function of K and X j and E
Fuzzy inference of K (step gain) is performed by the MIN-MAX method based on the inference result of j .

【0060】34は誤差散乱推論部であって,散乱誤差
のメンバシップ関数とEj ,Xj の推論結果に基づい
て,MIN−MAX法により誤差散乱の係数a(t)を
ファジィ推論をするものである。
An error scattering inference unit 34 performs fuzzy inference on the error scattering coefficient a (t) by the MIN-MAX method based on the scattering error membership function and the inference results of E j and X j. It is a thing.

【0061】非ファジィ化部10において,41はXj
演算部であって,Xj のファジィ推論結果を統合し,重
心法によりXj を非ファジィ化した値を求めるものであ
る。
In the defuzzification section 10, 41 is X j
An arithmetic unit integrates the fuzzy inference results of X j, and requests the defuzzification values of X j by gravity method.

【0062】42はEj 演算部であって,Ej のファジ
ィ推論結果を統合し,重心法によりEj を非ファジィ化
した値を求めるものである。43はK演算部であって,
Kのファジィ推論結果を統合し,重心法によりKを非フ
ァジィ化した値を求めるものである。
[0062] 42 is a E j arithmetic unit integrates the fuzzy inference results of E j, and requests a value obtained by defuzzification the E j by gravity method. 43 is a K operation unit,
The fuzzy inference result of K is integrated, and the value obtained by defuzzifying K by the centroid method is obtained.

【0063】44は誤差散乱演算部であって,誤差散乱
推論部のファジィ推論結果およびX j 推論部のファジィ
推論結果を統合し,誤差散乱の標本値a(t)を演算す
るものである。Yj =Σa(t)Xj で求める。
Reference numeral 44 denotes an error scattering calculation unit
Fuzzy reasoning result of the reasoning section and X jFuzzy inference section
Integrate inference results and calculate sample value a (t) of error scattering
Things. Yj= Σa (t) XjAsk in.

【0064】45はEd 演算部であって,Yj によりE
d を算出するものである。図5は本発明のメンバシップ
関数の実施例である。図5 (a)は入力信号Xj に対する
メンバシップ関数であり,VS(非常に小さい),S
(小さい),M(中ぐらいの大きさ),B(大きい),
VB(非常に大きい)の5種類である。横軸Suppo
rtの単位は電圧である。
[0064] 45 is a E d arithmetic unit, E by Y j
It calculates d . FIG. 5 is an embodiment of the membership function of the present invention. FIG. 5 (a) shows the membership function for the input signal X j , VS (very small), S
(Small), M (medium size), B (large),
There are five types of VB (very large). Horizontal axis Suppo
The unit of rt is voltage.

【0065】図5 (b)は出力信号Ej に対するメンバシ
ップ関数であり,VS(非常に小さい),S(小さ
い),M(中ぐらいの大きさ),B(大きい),VB
(非常に大きい)の5種類である。横軸Support
の単位は電圧である。
FIG. 5B is a membership function for the output signal E j , which is VS (very small), S (small), M (medium size), B (large), VB.
There are five types (very large). Horizontal axis Support
The unit of is voltage.

【0066】図5 (c)はステップゲインKに対するメン
バシップ関数であり,VS(非常に小さい),S(小さ
い),M(中ぐらいの大きさ),B(大きい),VB
(非常に大きい)の5種類である。
FIG. 5C shows the membership function with respect to the step gain K. VS (very small), S (small), M (medium size), B (large), VB
There are five types (very large).

【0067】図5 (d)はa(t)に対するメンバシップ
関数であり,VS(非常に小さい),S(小さい),M
(中ぐらいの大きさ),B(大きい),VB(非常に大
きい)の5種類である。
FIG. 5 (d) is a membership function for a (t), VS (very small), S (small), M
There are five types: (medium size), B (large), and VB (very large).

【0068】図6は本発明のファジィ制御規則の例であ
る。図6 (a)はステップゲインKの制御規則の例であっ
て,Xj とEj のメンバシップ関数を前件部として後件
部Kを選択する規則である。
FIG. 6 shows an example of the fuzzy control rule of the present invention. FIG. 6A is an example of the control rule of the step gain K, which is a rule for selecting the consequent part K with the membership functions of X j and E j as the antecedent part.

【0069】例えば,前件部Xj がB,Ej がMのと
き,後件部KはBである。図6 (b)は前件部Xj とEj
に対して後件部である誤差散乱系のインパルス応答の標
本値a(t)を選択するファジィ制御規則である。
For example, when the antecedent part X j is B and E j is M, the consequent part K is B. Figure 6 (b) shows the antecedents X j and E j.
Is a fuzzy control rule for selecting the sample value a (t) of the impulse response of the error scattering system which is the consequent part.

【0070】例えば,前件部Xj がM,Ej がBのと
き,後件部a(t)のファジィ制御規則はMである。図
7は本発明の誤算散乱の標本値a(t)の算出方法の例
であり,メンバシップ関数をVS,M,VBの3種類に
設定した場合について説明する。
For example, when the antecedent part X j is M and E j is B, the fuzzy control rule of the antecedent part a (t) is M. FIG. 7 is an example of a method of calculating the sample value a (t) of the miscalculation scattering of the present invention, and a case where the membership functions are set to three types of VS, M, and VB will be described.

【0071】図7 (a)はXj のメンバシップ関数であ
る。図7 (b)はEj のメンバシップ関数である。図7
(c)はXj のメンバシップ関数がVSでそのサポートの
値がXj0,Ej のメンバシップ関数がMでそのサポート
の値がEj0のときに,a(t)をMIN−MAX法によ
りファジィ推論して得られる結果を示す。
FIG. 7A shows the membership function of X j . FIG. 7B shows the membership function of E j . Figure 7
(c) shows that when the membership function of X j is VS, its support value is X j0 , the membership function of E j is M, and its support value is E j0 , a (t) is MIN-MAX method. The results obtained by fuzzy inference are shown by.

【0072】図7 (d)はXj のメンバシップ関数がM
で,そのサポートの値がXj0,Ej のメンバシップ関数
がMでそのサポートの値がEj0のときに,a(t)をM
IN−MAX法によりファジィ推論して得られる結果を
示す。
In FIG. 7D, the membership function of X j is M.
Then, when the support value is X j0 , the membership function of E j is M, and the support value is E j0 , a (t) is M
The result obtained by fuzzy inference by the IN-MAX method is shown.

【0073】図7 (e)は,図7 (c)と図7 (d)のファジ
ィ推論の結果を統合したものである。入力信号Xj がX
j0,出力信号Ej がEj0であるとする。誤算散乱の標本
値a(t)のメンバシップ関数のファジィ制御規則は,
次のとおりである(図6のa(t)のメンバシップ関数
のファジィ制御規則に従うとする)。
FIG. 7 (e) is an integration of the results of the fuzzy inference of FIGS. 7 (c) and 7 (d). Input signal X j is X
It is assumed that j0 and the output signal E j are E j0 . The fuzzy control rule of the membership function of sample value a (t) of miscalculation is
It is as follows (assuming that the fuzzy control rule of the membership function of a (t) in FIG. 6 is followed).

【0074】 前件部Xj がVS,Ej がMのとき,後件部a(t)はVSである・・(1) 前件部Xj がM,Ej がMのとき,後件部a(t)はMである。・・・(2) Xj のメンバシップ関数VSにおいてXj0に対するグレ
ードμ(Xj0)は0.7である。また,メンバシップ関
数MにおいてXj0に対するグレードμ(Xj0)は0.3
である。
When the antecedent part X j is VS and the E j is M, the antecedent part a (t) is VS ... (1) When the antecedent part X j is M and E j is M, the The subject part a (t) is M. (2) In the membership function VS of X j grade for X j0 μ (X j0) is 0.7. Further, Grade mu (X j0) for X j0 in the membership function M is 0.3
Is.

【0075】Ej のメンバシップ関数MにおいてEj0
対するグレードμ(Ej0)は0.6である。そこで,X
j のメンバシップ関数VSを前件部とした場合には,上
記条件(1)により,a( t)のメンバシップ関数はVS
である。そして,MIN−MAX法のファジィ演算規則
に従い,μ(Xj0)=0.3とμ(Ej0)=0.6の小
さい方を採用して,a(t)のグレードμ(a)=0.
3とし,a(t)のメンバシップ関数とグレードμ
(a)=0.3とのMAX演算を行い図7 (c)のファジ
ィ推論の結果を得る。
[0075] grade for E j0 in the membership function M of E j μ (E j0) is 0.6. So X
When the membership function VS of j is the antecedent part, the membership function of a (t) is VS according to the above condition (1).
Is. Then, according to the fuzzy operation rule of the MIN-MAX method, the smaller one of μ (X j0 ) = 0.3 and μ (E j0 ) = 0.6 is adopted, and the grade μ (a) = a (t) = 0.
3, the membership function of a (t) and the grade μ
The MAX operation with (a) = 0.3 is performed to obtain the result of the fuzzy inference of FIG. 7 (c).

【0076】同様に,Xj のメンバシップ関数Mを前件
部とした場合には,上記条件(2) により,a(t)のメ
ンバシップ関数はMである。そして,MIN−MAX法
のファジィ演算規則に従い,μ(Xj0)=0・7とμ
(Ej0)=0・6の小さい方を採用して,a(t)のグ
レードμ(a)=0.6とし,a(t)のメンバシップ
関数とグレードμ(a)=0.6とのMAX演算を行い
ファジィ推論の結果として図7 (d)を得る。
Similarly, when the membership function M of X j is the antecedent part, the membership function of a (t) is M according to the above condition (2). Then, according to the fuzzy operation rule of the MIN-MAX method, μ (X j0 ) = 0.7 and μ
Adopting the smaller one of (E j0 ) = 0.6, the grade μ (a) of a (t) = 0.6, and the membership function of a (t) and the grade μ (a) = 0.6. MAX calculation with and is performed, and FIG. 7 (d) is obtained as a result of fuzzy inference.

【0077】さらに,MIN−MAX法の演算規則に従
い,図7 (c)と図7 (d)のファジィ推論の結果を統合し
てMAXをとりa(t)に対するファジィ推論の結果と
して図7 (e)を得る。
Further, according to the operation rule of the MIN-MAX method, the fuzzy inference results of FIGS. 7 (c) and 7 (d) are integrated to obtain MAX and the fuzzy inference result for a (t) is obtained as shown in FIG. get e).

【0078】そのときの標本値a0 は,例えば図7 (e)
の重心を求め,その値を標本値a0とする。本発明は,
出力Ej の変化を観察し,Ej が高速に目標値に近づく
ようにメンバシップ関数のファジィ制御規則の変更をす
る。
The sample value a 0 at that time is, for example, as shown in FIG.
The center of gravity of is calculated, and the value is set as the sample value a 0 . The present invention
Observing the change of the output E j , the fuzzy control rule of the membership function is changed so that E j quickly approaches the target value.

【0079】図8は本発明のファジィ制御規則の変更方
法の例である。図8 (a)は,収束の度合いが小さい場合
である。(a)−1はステップゲインKの制御規則の変更
方法である。左下がりの方向に1つ後退させる。
FIG. 8 shows an example of the method of changing the fuzzy control rule of the present invention. FIG. 8A shows the case where the degree of convergence is small. (a) -1 is a method of changing the control rule of the step gain K. Move back one direction to the left.

【0080】例えば,Xj がVS,Ej がBでステップ
ゲインKがMのとき,Kは左斜めに1つ下のSになり,
j はS,Ej はBに変更をする。(a)−2はインパル
ス応答の標本値a(t)の変更規則である。左上がりの
方向に1つ進行させる。
For example, when X j is VS, E j is B, and the step gain K is M, K is one diagonally lower left S,
X j is changed to S and E j is changed to B. (a) -2 is a rule for changing the sample value a (t) of the impulse response. Advance one to the left.

【0081】例えば,Xj がVS,Ej がVSでインパ
ルス応答の標本値a(t)がVSのとき,a(t)は左
に斜め1つ上のVSを選択し,Xj はS,Ej はSに変
更する。
For example, when X j is VS, E j is VS, and the sample value a (t) of the impulse response is VS, a (t) selects a VS one diagonally higher to the left, and X j is S , E j is changed to S.

【0082】図8 (b)は消去周波数帯域幅が狭い場合で
ある。(a)−3はステップゲインKの制御規則の変更方
法である。y軸方向(上)に進行させる。
FIG. 8B shows a case where the erase frequency bandwidth is narrow. (a) -3 is a method of changing the control rule of the step gain K. Proceed in the y-axis direction (up).

【0083】(a)−4はインパルス応答の標本値a
(t)の制御規則である。右下がりの方向に後退させ
る。図8 (c)は収束が安定していない場合である。
(A) -4 is a sample value a of the impulse response
It is a control rule of (t). Move backward in the downward direction. FIG. 8 (c) shows the case where the convergence is not stable.

【0084】(a)−5はステップゲインKの制御規則の
変更方法である。y軸方向(下)に進行させる。(a)−
6はインパルス応答の標本値a(t)の制御規則の変更
方法である。右上がりの方向に後退させる。
(A) -5 is a method of changing the control rule of the step gain K. Proceed in the y-axis direction (down). (a)-
6 is a method of changing the control rule of the sample value a (t) of the impulse response. Move backward in the upward direction.

【0085】図9は本発明の出力のファジィ評価方法の
説明図である。図9 (a)は消去前と消去後のEj の変化
の説明図である。縦軸はEj であり,横軸は周波数fで
ある。消去前と消去後とでの差を観察し,どの程度減っ
たかにより評価する。
FIG. 9 is an explanatory diagram of an output fuzzy evaluation method of the present invention. FIG. 9A is an explanatory diagram of changes in E j before and after erasing. The vertical axis represents E j and the horizontal axis represents frequency f. Observe the difference between before and after erasing, and evaluate the degree of reduction.

【0086】評価は,収束性の度合いの判断,周波
数領域における判断,収束性の安定性について判断す
る。 収束性の度合いの判断 図9(a) により収束性の度合いの判断について説明す
る。
In the evaluation, the degree of convergence is judged, the judgment is made in the frequency domain, and the stability of convergence is judged. Judgment of Convergence Degree The judgment of the degree of convergence will be described with reference to FIG.

【0087】収束の度合の判断はEj のサンプリング値
を数回,例えば16回ないし64回採集し,それを平均
させることにより,突発的に大きいEj を除き,平均化
させて収束性の判断を行う。
The determination of the degree of convergence is performed by collecting the sampled values of E j several times, for example, 16 to 64 times, and averaging them to eliminate the suddenly large E j and average them. Make a decision.

【0088】 周波数領域における判断 図9 (b)により周波数領域についての判断を説明する。
図9 (b)の縦軸はEj (f)であって,Ej の周波数に
対する関数である。横軸は周波数である。
Judgment in Frequency Domain The judgment in the frequency domain will be described with reference to FIG.
The vertical axis of FIG. 9B is E j (f), which is a function of E j with respect to frequency. The horizontal axis is frequency.

【0089】fL は消去できた周波数の最小値であり,
U は消去できた周波数の最大値である。FFT(高速
フーリエ変換)により,消去量と周波数の関係E
j (f)を求め,消去量を評価する。
F L is the minimum frequency that can be erased ,
f U is the maximum frequency that can be erased. By FFT (Fast Fourier Transform) E
j (f) is obtained and the erase amount is evaluated.

【0090】 収束性の安定性 図9 (c)を参照して収束性の判断をする。図9 (c)にお
いて,縦軸は出力の時間関数Ej (t)である。横軸は
時間である。
Convergence stability is determined with reference to FIG. 9C. In FIG. 9C, the vertical axis is the output time function E j (t). The horizontal axis is time.

【0091】Ej の最小値を検出したらその値を基準と
して何回かサンプリングし,それでも最小値近傍であれ
ば安定した収束状態であると判断する。図10は本発明
の誤差の評価方法の説明図である。
When the minimum value of E j is detected, it is sampled several times with that value as a reference, and if it is still near the minimum value, it is determined that the convergence state is stable. FIG. 10 is an explanatory diagram of the error evaluation method of the present invention.

【0092】図10 (a)は本発明を総合評価する方法を
示す。本発明は,収束の度合いEj ,収束の安定性Ej
(t),収束の周波数帯域E j (f)について総合的に
評価し,それぞれが最小となる状態を求める。
FIG. 10A shows a method for comprehensively evaluating the present invention.
Show. In the present invention, the degree of convergence Ej, Convergence stability Ej
(T), convergence frequency band E jAbout (f)
Evaluate and find the state that minimizes each.

【0093】図10 (b)は本発明の誤差の評価をするた
めのメンバシップ関数の例である。横軸は誤差の評価
(Ej ,Ej (f),Ej (t))である。縦軸はグレ
イドを表す。
FIG. 10 (b) shows an example of the membership function for evaluating the error of the present invention. The horizontal axis is the error evaluation (E j , E j (f), E j (t)). The vertical axis represents the grade.

【0094】Very Bad(非常に悪い),Sos
o(容認できる程度),Excelent(非常によ
い)のメンバシップ関数をもとにMIN−MAX法等で
ファジィ推論をし,求められた推論結果に対して重心法
等で非ファジィ化をし,評価値を求める。
Very Bad (very bad), Sos
Based on the membership function of o (acceptable) and Excellent (very good), fuzzy inference is performed by the MIN-MAX method, etc., and the obtained inference result is defuzzified by the centroid method, etc. Obtain an evaluation value.

【0095】図11は本発明のメンバシップ関数の変更
方法の例を示す。図11 (a)は,VSのメンバシップ関
数とその演算式の例である。S2 ,Smの値を変更する
ことにより横方向に平行移動させることができる。
FIG. 11 shows an example of a membership function changing method according to the present invention. FIG. 11A shows an example of the VS membership function and its arithmetic expression. By changing the values of S 2 and Sm, it is possible to make a parallel translation in the lateral direction.

【0096】図11 (b)は,Mのメンバシップ関数とそ
の演算式の例である。S1 ,S3 ,Sm の値を変更する
ことにより横方向に平行移動させることができる。図1
1 (c)は,VBのメンバシップ関数とその演算式の例で
ある。S4 ,Smの値を変更することにより横方向に平
行移動させることができる。
FIG. 11B shows an example of the membership function of M and its arithmetic expression. By changing the values of S 1 , S 3 and Sm, it is possible to perform parallel translation in the lateral direction. FIG.
1 (c) is an example of the membership function of VB and its arithmetic expression. It can be moved in parallel in the lateral direction by changing the value of S 4, Sm.

【0097】[0097]

【発明の効果】本発明によれば,Cフィルタを使用する
ことなく誤差散乱系をプログラムにより推定でき,誤差
散乱系の誤差を時間軸方向に散乱させた場合のインパル
ス応答の標本値a(t)の算出を短縮できる。そのた
め,適応フィルタの学習を高速化することができる。ま
た,システムの稼動中でもメンバシップ関数のファジィ
制御規則,メンバシップ関数を補正することができ学習
による消去効果を大幅に向上させることができる。
According to the present invention, the error scattering system can be estimated by a program without using the C filter, and the sample value a (t) of the impulse response when the error of the error scattering system is scattered in the time axis direction. ) Can be shortened. Therefore, learning of the adaptive filter can be speeded up. In addition, the fuzzy control rule of the membership function and the membership function can be corrected even while the system is operating, and the erasing effect by learning can be greatly improved.

【0098】また,ノイズ消去の誤差を誤差散乱系の確
率過程として捉えファジィ写像により制御しているの
で,柔軟にノイズ消去を行うことができる。
Further, since the error of noise elimination is regarded as the stochastic process of the error scattering system and controlled by the fuzzy mapping, the noise elimination can be flexibly performed.

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

【図1】本発明の基本構成を示す図である。FIG. 1 is a diagram showing a basic configuration of the present invention.

【図2】本発明の実施例1を示す図である。FIG. 2 is a diagram showing a first embodiment of the present invention.

【図3】本発明の実施例2を示す図である。FIG. 3 is a diagram showing a second embodiment of the present invention.

【図4】本発明のファジィ演算部の実施例を示す図であ
る。
FIG. 4 is a diagram showing an embodiment of a fuzzy operation unit of the present invention.

【図5】メンバシップ関数の実施例を示す図である。FIG. 5 is a diagram showing an example of a membership function.

【図6】ファジィ制御規則の例を示す図である。FIG. 6 is a diagram showing an example of a fuzzy control rule.

【図7】誤差散乱の標本値a(t)の算出方法の例を示
す図である。
FIG. 7 is a diagram illustrating an example of a method of calculating a sample value a (t) of error scattering.

【図8】ファジィ規則の変更方法の例を示す図である。FIG. 8 is a diagram showing an example of a method for changing a fuzzy rule.

【図9】ファジィ評価の方法の説明図である。FIG. 9 is an explanatory diagram of a fuzzy evaluation method.

【図10】本発明の誤差の評価方法を示す図である。FIG. 10 is a diagram showing an error evaluation method of the present invention.

【図11】本発明のメンバシップ関数の変更方法の例を
示す図である。
FIG. 11 is a diagram showing an example of a membership function changing method of the present invention.

【図12】従来の技術を示す図である。FIG. 12 is a diagram showing a conventional technique.

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

1:未知系 2:適応フィルタ 5:ファジィ演算部 6:誤差散乱系 11:誤差散乱演算部 16:誤差散乱推論部 17:評価部 21:メンバシップ関数保持部 22:ファジィ制御規則保持部 1: Unknown system 2: Adaptive filter 5: Fuzzy operation unit 6: Error scattering system 11: Error scattering operation unit 16: Error scattering inference unit 17: Evaluation unit 21: Membership function holding unit 22: Fuzzy control rule holding unit

Claims (4)

【特許請求の範囲】[Claims] 【請求項1】 ノイズに対して能動的にノイズ消去信号
を印加してノイズを消去する能動ノイズ消去装置におい
て,能動ノイズ消去装置はノイズを消去した誤差が時間
的に散乱する誤差散乱系を持つものであって,該誤差散
乱系により散乱された誤差と入力信号に対してファジィ
推論をするファジィ演算部とファジィ演算部の演算結果
と入力信号を入力してノイズ消去信号を出力する適応フ
ィルタとを備え,ファジィ演算部の演算結果に基づいて
適応フィルタの最適係数を求め,ノイズ消去信号を生成
することを特徴とするファジィ能動ノイズ消去装置。
1. In an active noise canceller for canceling noise by actively applying a noise canceling signal to noise, the active noise canceller has an error scattering system in which an error canceling noise is scattered over time. A fuzzy operation unit for performing fuzzy inference on the error scattered by the error scattering system and the input signal, and an adaptive filter for inputting the operation result of the fuzzy operation unit and the input signal and outputting a noise cancellation signal. A fuzzy active noise eliminator, which is characterized in that it obtains an optimum coefficient of an adaptive filter based on a calculation result of a fuzzy calculator and generates a noise canceling signal.
【請求項2】 ファジィ演算部は誤差散乱の標本値a
(t)をファジィ演算し,Yj =Σa(t)Xj (但
し,Xj は入力信号であり,Σはt=0〜T−1につい
ての和であって,Tはタップ数に相当する値)で誤差散
乱系を推定することを特徴とする請求項1に記載のファ
ジィ能動ノイズ消去装置。
2. The fuzzy operation unit is a sample value a of error scattering.
(T) is fuzzy calculated and Y j = Σa (t) X j (where X j is the input signal, Σ is the sum for t = 0 to T−1, and T corresponds to the number of taps. The fuzzy active noise canceller according to claim 1, wherein the error scattering system is estimated by the value of
【請求項3】 ファジィ演算部はノイズ消去の誤差を評
価する評価部と誤差評価のためのメンバシップ関数を備
え,誤差をファジィ評価することを特徴とするファジィ
能動ノイズ消去装置。
3. A fuzzy active noise canceller, characterized in that the fuzzy arithmetic unit comprises an evaluation unit for evaluating an error of noise cancellation and a membership function for error evaluation, and fuzzy evaluates the error.
【請求項4】 ダクトの騒音を消去するものであること
を特徴とする請求項1,2もしくは3に記載のファジィ
能動ノイズ消去装置。
4. The fuzzy active noise canceller according to claim 1, wherein the fuzzy active noise canceller is for canceling noise in the duct.
JP13164895A 1995-05-30 1995-05-30 Fuzzy active noise canceller Expired - Fee Related JP3445872B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP13164895A JP3445872B2 (en) 1995-05-30 1995-05-30 Fuzzy active noise canceller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP13164895A JP3445872B2 (en) 1995-05-30 1995-05-30 Fuzzy active noise canceller

Publications (2)

Publication Number Publication Date
JPH08328569A true JPH08328569A (en) 1996-12-13
JP3445872B2 JP3445872B2 (en) 2003-09-08

Family

ID=15062978

Family Applications (1)

Application Number Title Priority Date Filing Date
JP13164895A Expired - Fee Related JP3445872B2 (en) 1995-05-30 1995-05-30 Fuzzy active noise canceller

Country Status (1)

Country Link
JP (1) JP3445872B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539229A (en) * 2021-07-30 2021-10-22 北京安声浩朗科技有限公司 Noise reduction parameter determination method and device, active noise reduction method and device
CN113539228A (en) * 2021-07-30 2021-10-22 北京安声浩朗科技有限公司 Noise reduction parameter determination method and device, active noise reduction method and device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539229A (en) * 2021-07-30 2021-10-22 北京安声浩朗科技有限公司 Noise reduction parameter determination method and device, active noise reduction method and device
CN113539228A (en) * 2021-07-30 2021-10-22 北京安声浩朗科技有限公司 Noise reduction parameter determination method and device, active noise reduction method and device
CN113539229B (en) * 2021-07-30 2023-10-31 北京安声浩朗科技有限公司 Noise reduction parameter determining method and device, active noise reduction method and device
CN113539228B (en) * 2021-07-30 2023-10-31 北京安声浩朗科技有限公司 Noise reduction parameter determining method and device, active noise reduction method and device

Also Published As

Publication number Publication date
JP3445872B2 (en) 2003-09-08

Similar Documents

Publication Publication Date Title
US8199924B2 (en) System for active noise control with an infinite impulse response filter
JP5496418B2 (en) Adaptive equalizer, acoustic echo canceller device and active noise control device
US5701349A (en) Active vibration controller
Bergamasco et al. Active noise control with on-line estimation of non-Gaussian noise characteristics
US11087735B2 (en) Active noise control method and system
JPH05265468A (en) Active type noise controller
Yin et al. Hermite functional link artificial-neural-network-assisted adaptive algorithms for IoV nonlinear active noise control
JP2023502076A (en) Active denoising system with convergence detection
JP3445872B2 (en) Fuzzy active noise canceller
Kawamura et al. A noise reduction method based on linear prediction with variable step-size
Ellison et al. An Acoustic Scene Perspective on Spatial, Temporal, and Spectral Aspects of Marine Mammal Behavioral Responses to Noise.
JPH0511771A (en) Noise control device
JP2004078021A (en) Method, device, and program for sound pickup
CN116547989A (en) System and method for adaptive beamforming
JPH07160508A (en) Filter coefficient deciding method for adaptive filter
US5949890A (en) Active noise control apparatus and waveform transforming apparatus through neural network
JP3550751B2 (en) Active control device
JP3089794B2 (en) Method and apparatus for identifying unknown system using adaptive filter
JP3636772B2 (en) Fuzzy noise canceler
JPH0844375A (en) Noise eliminating device and noise eliminating method
JPH0588684A (en) Adaptive signal processing method, adaptive signal processor, and active noise controller
Polo Development and testing of a number of MATLAB based fuzzy system applications
JPH0758223B2 (en) Method of measuring sound deadening / damping effect, measuring device, and signal source search device
Wahbia et al. A Robust Embedded Non-Linear Acoustic Noise Cancellation (ANC) Using Artificial Neural Network (ANN) for Improving the Quality of Voice Communications
Sakai An Improved structure of Acoustic Echo Cancellation Using Blind Source Separation

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20030617

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090627

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100627

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110627

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120627

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees