JPH10215194A - 雑音消去装置 - Google Patents
雑音消去装置Info
- Publication number
- JPH10215194A JPH10215194A JP9014410A JP1441097A JPH10215194A JP H10215194 A JPH10215194 A JP H10215194A JP 9014410 A JP9014410 A JP 9014410A JP 1441097 A JP1441097 A JP 1441097A JP H10215194 A JPH10215194 A JP H10215194A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- noise
- value
- estimated
- power
- 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
Links
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/10—Applications
- G10K2210/108—Communication systems, e.g. where useful sound is kept and noise is cancelled
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3012—Algorithms
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3023—Estimation of noise, e.g. on error signals
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3045—Multiple acoustic inputs, single acoustic output
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3054—Stepsize variation
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
- G10L2021/02161—Number of inputs available containing the signal or the noise to be suppressed
- G10L2021/02165—Two microphones, one receiving mainly the noise signal and the other one mainly the speech signal
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
Abstract
(57)【要約】
【課題】 雑音消去装置において、十分な雑音信号の消
去と所望信号の参照マイクロフォンへの混入による信号
歪みの低減を実現する。 【解決手段】 主信号と参照信号を受けて、主信号にお
ける所望信号の平均電力と雑音信号の平均電力から主信
号における信号対雑音電力比を推定すると共に、参照信
号における所望信号の平均電力と雑音信号の平均電力か
ら参照信号における信号対雑音電力比を推定する信号対
雑音電力比推定回路１０を有する。主信号中の雑音信号
を推定する適応フィルタ３の係数更新のステップサイズ
は推定された主信号における信号対雑音電力比によって
制御される。参照信号中の所望信号を推定する適応フィ
ルタ６の係数更新のステップサイズは推定された参照信
号における信号対雑音電力比によって制御される。遅延
回路８，９は、信号対雑音電力比推定回路１０において
信号対雑音電力比の推定値が算出される際に電力平均操
作によって生じる時間遅延を補償する。
去と所望信号の参照マイクロフォンへの混入による信号
歪みの低減を実現する。 【解決手段】 主信号と参照信号を受けて、主信号にお
ける所望信号の平均電力と雑音信号の平均電力から主信
号における信号対雑音電力比を推定すると共に、参照信
号における所望信号の平均電力と雑音信号の平均電力か
ら参照信号における信号対雑音電力比を推定する信号対
雑音電力比推定回路１０を有する。主信号中の雑音信号
を推定する適応フィルタ３の係数更新のステップサイズ
は推定された主信号における信号対雑音電力比によって
制御される。参照信号中の所望信号を推定する適応フィ
ルタ６の係数更新のステップサイズは推定された参照信
号における信号対雑音電力比によって制御される。遅延
回路８，９は、信号対雑音電力比推定回路１０において
信号対雑音電力比の推定値が算出される際に電力平均操
作によって生じる時間遅延を補償する。
Description
【０００１】
【発明の属する技術分野】本発明は雑音消去装置に関
し、特マイクロフォンやハンドセット等から入力された
音声信号に対して背景雑音信号を適応フィルタを使用し
て除去するようにした雑音消去装置に関するものであ
る。
し、特マイクロフォンやハンドセット等から入力された
音声信号に対して背景雑音信号を適応フィルタを使用し
て除去するようにした雑音消去装置に関するものであ
る。
【０００２】
【従来の技術】マイクロフォンやハンドセット等から入
力された音声信号に混入した背景雑音信号は、情報圧縮
度の高い狭帯域音声符号化装置や音声認識装置等におい
ては大きな問題となる。この様な音響的に重畳した雑音
成分の消去を目的とした雑音消去装置として、適応フィ
ルタを用いた２入力型雑音消去装置が１９７５年発行の
刊行物「プロシーディングズ・オブ・アイ・イー・イー
・イー」の６３巻１２号の１６９２〜１７１６頁に記載
されている（Ｂ．Ｗｉｄｒｏｗ ｅｔ．ａｌ．，“Ａｄ
ａｐｔｉｖｅ Ｎｏｉｓｅ Ｃａｎｃｅｌｌｉｎｇ：Ｐ
ｒｉｎｃｉｐｌｅａｎｄ Ａｐｐｌｉｃａｔｉｏｎ
ｓ，”ＰＲＯＣＥＥＤＩＮＧＳ ＯＦ ＩＥＥＥ，ＶＯ
Ｌ．６３，ＮＯ．１２，１９７５，ｐｐ．１６９２〜１
７１６；以下、文献１という）。
力された音声信号に混入した背景雑音信号は、情報圧縮
度の高い狭帯域音声符号化装置や音声認識装置等におい
ては大きな問題となる。この様な音響的に重畳した雑音
成分の消去を目的とした雑音消去装置として、適応フィ
ルタを用いた２入力型雑音消去装置が１９７５年発行の
刊行物「プロシーディングズ・オブ・アイ・イー・イー
・イー」の６３巻１２号の１６９２〜１７１６頁に記載
されている（Ｂ．Ｗｉｄｒｏｗ ｅｔ．ａｌ．，“Ａｄ
ａｐｔｉｖｅ Ｎｏｉｓｅ Ｃａｎｃｅｌｌｉｎｇ：Ｐ
ｒｉｎｃｉｐｌｅａｎｄ Ａｐｐｌｉｃａｔｉｏｎ
ｓ，”ＰＲＯＣＥＥＤＩＮＧＳ ＯＦ ＩＥＥＥ，ＶＯ
Ｌ．６３，ＮＯ．１２，１９７５，ｐｐ．１６９２〜１
７１６；以下、文献１という）。
【０００３】この２入力型雑音消去装置は、参照信号用
マイクロフォンに入力された雑音信号が主信号マイクロ
フォンに到達するまでに通る音響経路（ノイズパス）の
インパルス応答を近似する適応フィルタを用いて、主信
号用マイクロフォンに混入する雑音を推定し、主信号用
マイクロフォンに入力された主信号（所望信号と雑音の
混在信号）から推定雑音信号を減算することによって、
雑音信号を抑圧する様に動作する。
マイクロフォンに入力された雑音信号が主信号マイクロ
フォンに到達するまでに通る音響経路（ノイズパス）の
インパルス応答を近似する適応フィルタを用いて、主信
号用マイクロフォンに混入する雑音を推定し、主信号用
マイクロフォンに入力された主信号（所望信号と雑音の
混在信号）から推定雑音信号を減算することによって、
雑音信号を抑圧する様に動作する。
【０００４】この時、適応フィルタのフィルタ係数は、
主信号から推定雑音信号を減算して得られる誤差信号と
参照用マイクロフォンにて得られる参照信号との相関を
取ることにより修正される。この様な適応フィルタの係
数修正、すなわち、収束アルゴリズムの代表的なものと
して文献１に記載されている「エル・エム・エス・アル
ゴリズム（ＬＭＳ ＡＬＧＯＲＩＴＨＭ）」や「アイ・
イー・イー・イー・トランザクションズ・オン・オート
マチック・コントロール，１２巻，３号，１９６７年，
２８２〜２８７頁（ＩＥＥＥ ＴＲＡＮＳＡＣＴＩＯＮ
Ｓ ＯＮ ＡＵＴＯＭＡＴＩＣ ＣＯＮＴＲＯＬ，ＶＯ
Ｌ．１２，ＮＯ．３，１９６７，ｐｐ．２８２〜２８
７）」（以下、文献２という）に記載されている「ラー
ニング・アイデンティフィケーション・メソッド（ＬＥ
ＡＲＮＩＮＧ ＩＤＥＮＴＩＦＩＣＡＴＩＯＮ ＭＥＴ
ＨＯＤ；ＬＩＭ）（「学習同定法」ともいう）が知られ
ている。
主信号から推定雑音信号を減算して得られる誤差信号と
参照用マイクロフォンにて得られる参照信号との相関を
取ることにより修正される。この様な適応フィルタの係
数修正、すなわち、収束アルゴリズムの代表的なものと
して文献１に記載されている「エル・エム・エス・アル
ゴリズム（ＬＭＳ ＡＬＧＯＲＩＴＨＭ）」や「アイ・
イー・イー・イー・トランザクションズ・オン・オート
マチック・コントロール，１２巻，３号，１９６７年，
２８２〜２８７頁（ＩＥＥＥ ＴＲＡＮＳＡＣＴＩＯＮ
Ｓ ＯＮ ＡＵＴＯＭＡＴＩＣ ＣＯＮＴＲＯＬ，ＶＯ
Ｌ．１２，ＮＯ．３，１９６７，ｐｐ．２８２〜２８
７）」（以下、文献２という）に記載されている「ラー
ニング・アイデンティフィケーション・メソッド（ＬＥ
ＡＲＮＩＮＧ ＩＤＥＮＴＩＦＩＣＡＴＩＯＮ ＭＥＴ
ＨＯＤ；ＬＩＭ）（「学習同定法」ともいう）が知られ
ている。
【０００５】図５は従来の雑音消去装置の原理を説明す
るブロック図である。雑音消去装置は、主信号用マイク
ロフォン１，参照信号用マイクロフォン２，適応フィル
タ３，減算器４，出力端子５から構成されている。話者
の発した所望信号Ｓ（ｚ）は、音響伝達特性ＨＡ（ｚ）
の経路を通って話者の口元に置かれた主信号用マイクロ
フォン１に入力される。但し、 ｚ＝exp （２πｊ／ＦＳ）…（１） である。尚、ＦＳはサンプリング周波数である。
るブロック図である。雑音消去装置は、主信号用マイク
ロフォン１，参照信号用マイクロフォン２，適応フィル
タ３，減算器４，出力端子５から構成されている。話者
の発した所望信号Ｓ（ｚ）は、音響伝達特性ＨＡ（ｚ）
の経路を通って話者の口元に置かれた主信号用マイクロ
フォン１に入力される。但し、 ｚ＝exp （２πｊ／ＦＳ）…（１） である。尚、ＦＳはサンプリング周波数である。
【０００６】一方、雑音源から発生している雑音Ｎ
（ｚ）は、音響伝達特性ＧＡ（ｚ）の経路を通って主信
号マイクロフォンに入力されると共に、音響伝達特性Ｇ
Ｂ（ｚ）の経路を通って話者から離れた位置に置かれた
参照信号用マイクロフォン２に入力される。適応フィル
タ３は、例えば前述した「ＬＩＭ」に従って、上記主信
号ＸＡ（ｚ）及び参照信号ＸＢ（ｚ）を基に、所望信号
Ｓ（ｚ）が入力されていない時の雑音源Ｎ（ｚ）から、
参照信号用マイクロフォン２に入力された雑音が主信号
用マイクロフォン１に到達するまでに通る音響伝達特性
（ノイズパス）Ｐ（ｚ）を推定する。
（ｚ）は、音響伝達特性ＧＡ（ｚ）の経路を通って主信
号マイクロフォンに入力されると共に、音響伝達特性Ｇ
Ｂ（ｚ）の経路を通って話者から離れた位置に置かれた
参照信号用マイクロフォン２に入力される。適応フィル
タ３は、例えば前述した「ＬＩＭ」に従って、上記主信
号ＸＡ（ｚ）及び参照信号ＸＢ（ｚ）を基に、所望信号
Ｓ（ｚ）が入力されていない時の雑音源Ｎ（ｚ）から、
参照信号用マイクロフォン２に入力された雑音が主信号
用マイクロフォン１に到達するまでに通る音響伝達特性
（ノイズパス）Ｐ（ｚ）を推定する。
【０００７】推定すべき音響伝達特性Ｐ（ｚ）は、 Ｐ（ｚ）＝ＧＡ（ｚ）／ＧＢ（ｚ）…（２） の様に表される。すなわち、適応フィルタ３は音響伝達
特性Ｐ（ｚ）と同じ伝達関数Ｗｌ（ｚ）を有するフィル
タを構成して、主信号に含まれる雑音成分と同じ推定雑
音信号Ｆ1 （ｚ）を生成する様動作する。減算器４は主
信号ＸＡ（ｚ）から適応フィルタ３が出力する推定雑音
信号Ｆ1 （ｚ）を減算し、結果として出力信号Ｅ1
（ｚ）を出力する。所望信号Ｓ（ｚ）が入力されていな
い時の出力信号Ｅ1 （ｚ）は、 Ｅ1 （ｚ）＝ＸＡ（ｚ）−Ｆ1 （ｚ） ＝ＸＡ（ｚ）−Ｗ1 （ｚ）ＸＢ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）−Ｗ1 （ｚ）ＧＢ（ｚ）Ｎ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）−｛ＧＡ（ｚ）／ＧＢ（ｚ）｝ ＧＢ（ｚ）Ｎ（ｚ） ＝０…（３） と表される。
特性Ｐ（ｚ）と同じ伝達関数Ｗｌ（ｚ）を有するフィル
タを構成して、主信号に含まれる雑音成分と同じ推定雑
音信号Ｆ1 （ｚ）を生成する様動作する。減算器４は主
信号ＸＡ（ｚ）から適応フィルタ３が出力する推定雑音
信号Ｆ1 （ｚ）を減算し、結果として出力信号Ｅ1
（ｚ）を出力する。所望信号Ｓ（ｚ）が入力されていな
い時の出力信号Ｅ1 （ｚ）は、 Ｅ1 （ｚ）＝ＸＡ（ｚ）−Ｆ1 （ｚ） ＝ＸＡ（ｚ）−Ｗ1 （ｚ）ＸＢ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）−Ｗ1 （ｚ）ＧＢ（ｚ）Ｎ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）−｛ＧＡ（ｚ）／ＧＢ（ｚ）｝ ＧＢ（ｚ）Ｎ（ｚ） ＝０…（３） と表される。
【０００８】この様に、適応フィルタ３は、所望信号Ｓ
（ｚ）が含まれていない時に出力信号Ｅ1 （ｚ）が零に
なる様に係数更新を行うことによって音響伝達特性Ｐ
（ｚ）を推定することができる。この出力信号Ｅ1
（ｚ）は、適応フィルタの学習動作において、誤差を表
すことから、誤差信号と呼ばれている。
（ｚ）が含まれていない時に出力信号Ｅ1 （ｚ）が零に
なる様に係数更新を行うことによって音響伝達特性Ｐ
（ｚ）を推定することができる。この出力信号Ｅ1
（ｚ）は、適応フィルタの学習動作において、誤差を表
すことから、誤差信号と呼ばれている。
【０００９】適応フィルタ３が収束した後の出力信号Ｅ
1 （ｚ）は、 Ｅ1 （ｚ）＝ＸＡ（ｚ）−Ｆ1 （ｚ） ＝ＸＡ（ｚ）−Ｗ1 （ｚ）ＸＢ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ） −Ｗ1 （ｚ）ＧＢ（ｚ）Ｎ（ｚ）＋ＨＢ（ｚ）Ｓ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ） −Ｗ1 （ｚ）ＧＢ（ｚ）Ｎ（ｚ）−Ｗ1 （ｚ）ＨＢ（ｚ）Ｓ（ｚ） ＝ＨＡ（ｚ）Ｓ（ｚ）−Ｗ1 （ｚ）ＨＢ（ｚ）Ｓ（ｚ） ＝ＨＡ（ｚ）Ｓ（ｚ） ［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］…（４） と表される。
1 （ｚ）は、 Ｅ1 （ｚ）＝ＸＡ（ｚ）−Ｆ1 （ｚ） ＝ＸＡ（ｚ）−Ｗ1 （ｚ）ＸＢ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ） −Ｗ1 （ｚ）ＧＢ（ｚ）Ｎ（ｚ）＋ＨＢ（ｚ）Ｓ（ｚ） ＝ＧＡ（ｚ）Ｎ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ） −Ｗ1 （ｚ）ＧＢ（ｚ）Ｎ（ｚ）−Ｗ1 （ｚ）ＨＢ（ｚ）Ｓ（ｚ） ＝ＨＡ（ｚ）Ｓ（ｚ）−Ｗ1 （ｚ）ＨＢ（ｚ）Ｓ（ｚ） ＝ＨＡ（ｚ）Ｓ（ｚ） ［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］…（４） と表される。
【００１０】上式から分かる様に、出力信号Ｅ1 （ｚ）
には、雑音信号Ｎ（ｚ）は含まれておらず雑音が消去さ
れていることが分かる。しかしながら、参照用マイクロ
フォン２に所望信号成分Ｓ（ｚ）が含まれている場合、
すなわち、所望信号Ｓ（ｚ）から参照用マイクロフォン
２への音響伝達特性ＨＢ（ｚ）が零でない場合には、信
号歪み［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］
が発生するという問題があった。
には、雑音信号Ｎ（ｚ）は含まれておらず雑音が消去さ
れていることが分かる。しかしながら、参照用マイクロ
フォン２に所望信号成分Ｓ（ｚ）が含まれている場合、
すなわち、所望信号Ｓ（ｚ）から参照用マイクロフォン
２への音響伝達特性ＨＢ（ｚ）が零でない場合には、信
号歪み［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］
が発生するという問題があった。
【００１１】この問題を解決する方法として、出力信号
Ｅ1 （ｚ）に含まれる信号歪みを補正する適応フィルタ
を付加する方法が特開平８−５６１８０号公報で提案さ
れている。図６はこの特開平８−５６１８０号公報で提
案されている従来の雑音消去装置の構成の一例を示すブ
ロック図である。
Ｅ1 （ｚ）に含まれる信号歪みを補正する適応フィルタ
を付加する方法が特開平８−５６１８０号公報で提案さ
れている。図６はこの特開平８−５６１８０号公報で提
案されている従来の雑音消去装置の構成の一例を示すブ
ロック図である。
【００１２】この雑音消去装置では、図５の雑音消去装
置に信号歪みを補正するための適応フィルタ６と減算器
７が付加されている。適応フィルタ６は、主信号ＸＡ
（ｚ）に所望信号Ｓ（ｚ）が含まれ、かつ、雑音がない
かまたは一定レベル以下である時に減算器７の出力Ｅ2
（ｚ）が小さくなる様に学習動作を行う。適応フィルタ
６の有する伝達特性をＷ2 （ｚ）とすると、この学習
は、Ｎ（ｚ）が零、あるいは無視できる時に、Ｅ2
（ｚ）が、 Ｅ2 （ｚ）＝ＸＡ（ｚ）−Ｆ2 （ｚ） ＝ＸＡ（ｚ）−Ｗ2 （ｚ）Ｅ1 （ｚ） ＝ＨＡ（ｚ）Ｓ（ｚ）−Ｗ2 （ｚ）ＨＡ（ｚ）Ｓ（ｚ） ［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］ ＝ＨＡ（ｚ）Ｓ（ｚ）−Ｗ2 （ｚ）ＨＡ（ｚ）Ｓ（ｚ） ＋Ｗ2 （ｚ）Ｓ（ｚ）ＨＢ（ｚ）Ｗ1 （ｚ） ＝Ｗ2 （ｚ）Ｓ（ｚ）ＨＢ（ｚ）Ｗ1 （ｚ） −ＨＡ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ） ＝０…（５） となる様に、例えば「学習同定法」などによる適応アル
ゴリズムにより行われる。
置に信号歪みを補正するための適応フィルタ６と減算器
７が付加されている。適応フィルタ６は、主信号ＸＡ
（ｚ）に所望信号Ｓ（ｚ）が含まれ、かつ、雑音がない
かまたは一定レベル以下である時に減算器７の出力Ｅ2
（ｚ）が小さくなる様に学習動作を行う。適応フィルタ
６の有する伝達特性をＷ2 （ｚ）とすると、この学習
は、Ｎ（ｚ）が零、あるいは無視できる時に、Ｅ2
（ｚ）が、 Ｅ2 （ｚ）＝ＸＡ（ｚ）−Ｆ2 （ｚ） ＝ＸＡ（ｚ）−Ｗ2 （ｚ）Ｅ1 （ｚ） ＝ＨＡ（ｚ）Ｓ（ｚ）−Ｗ2 （ｚ）ＨＡ（ｚ）Ｓ（ｚ） ［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］ ＝ＨＡ（ｚ）Ｓ（ｚ）−Ｗ2 （ｚ）ＨＡ（ｚ）Ｓ（ｚ） ＋Ｗ2 （ｚ）Ｓ（ｚ）ＨＢ（ｚ）Ｗ1 （ｚ） ＝Ｗ2 （ｚ）Ｓ（ｚ）ＨＢ（ｚ）Ｗ1 （ｚ） −ＨＡ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ） ＝０…（５） となる様に、例えば「学習同定法」などによる適応アル
ゴリズムにより行われる。
【００１３】従って、適応フィルタ６の有する伝達特性
Ｗ2 （ｚ）は、 Ｗ2 （ｚ）＝ＨＡ（ｚ）／｛ＨＡ（ｚ）−ＨＢ（ｚ）Ｗ1 （ｚ）｝ ＝１／［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］…（６） となる。
Ｗ2 （ｚ）は、 Ｗ2 （ｚ）＝ＨＡ（ｚ）／｛ＨＡ（ｚ）−ＨＢ（ｚ）Ｗ1 （ｚ）｝ ＝１／［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］…（６） となる。
【００１４】この様にして、学習された適応フィルタ６
の出力する信号Ｆ2 （ｚ）は、 Ｆ2 （ｚ）＝Ｗ2 （ｚ）Ｅ1 （ｚ） ＝｛１／［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］｝ ＨＡ（ｚ）Ｓ（ｚ）［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］ ＝ＨＡ（ｚ）Ｓ（ｚ）…（７） となり、信号歪みの補正された所望信号ＨＡ（ｚ）Ｓ
（ｚ）を出力する。
の出力する信号Ｆ2 （ｚ）は、 Ｆ2 （ｚ）＝Ｗ2 （ｚ）Ｅ1 （ｚ） ＝｛１／［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］｝ ＨＡ（ｚ）Ｓ（ｚ）［１−｛ＨＢ（ｚ）／ＨＡ（ｚ）｝Ｗ1 （ｚ）］ ＝ＨＡ（ｚ）Ｓ（ｚ）…（７） となり、信号歪みの補正された所望信号ＨＡ（ｚ）Ｓ
（ｚ）を出力する。
【００１５】
【発明が解決しようとする課題】上述した様に、従来の
雑音消去装置では、雑音信号Ｎ（ｚ）が存在し、所望信
号成分Ｓ（ｚ）がないか無視できる程度に小さい区間
で、適応フィルタ３の係数更新動作を行い、雑音の音響
伝達特性を学習する。更に、所望信号成分Ｓ（ｚ）が存
在し、雑音成分Ｎ（ｚ）がないか無視できる程度に小さ
い区間で、適応フィルタ４の係数更新動作を行い、信号
歪みの補正フィルタを学習する。従って、上記学習を行
うためには、所望信号成分Ｓ（ｚ）がない（あるいは少
ない）区間と雑音信号成分Ｎ（ｚ）がない（あるいは少
ない）区間を検出し、外部から適応フィルタに学習動作
を指示する必要がある。
雑音消去装置では、雑音信号Ｎ（ｚ）が存在し、所望信
号成分Ｓ（ｚ）がないか無視できる程度に小さい区間
で、適応フィルタ３の係数更新動作を行い、雑音の音響
伝達特性を学習する。更に、所望信号成分Ｓ（ｚ）が存
在し、雑音成分Ｎ（ｚ）がないか無視できる程度に小さ
い区間で、適応フィルタ４の係数更新動作を行い、信号
歪みの補正フィルタを学習する。従って、上記学習を行
うためには、所望信号成分Ｓ（ｚ）がない（あるいは少
ない）区間と雑音信号成分Ｎ（ｚ）がない（あるいは少
ない）区間を検出し、外部から適応フィルタに学習動作
を指示する必要がある。
【００１６】しかしながら、雑音消去装置が設置される
状況によっては、一般に所望信号と雑音信号の各レベル
に応じて適応フィルタの学習区間を外部から指示するこ
とは困難な場合が多い。従来の雑音消去装置では、適切
な学習区間を指示して各適応フィルタを学習させない限
り、十分な雑音消去性能と歪み補正特性を得ることがで
きないという問題がある。
状況によっては、一般に所望信号と雑音信号の各レベル
に応じて適応フィルタの学習区間を外部から指示するこ
とは困難な場合が多い。従来の雑音消去装置では、適切
な学習区間を指示して各適応フィルタを学習させない限
り、十分な雑音消去性能と歪み補正特性を得ることがで
きないという問題がある。
【００１７】本発明の目的は、外部から適切な学習区間
を指示することができない場合においても、十分な雑音
消去性能と信号歪みの低減を実現できる雑音消去装置を
提供することである。
を指示することができない場合においても、十分な雑音
消去性能と信号歪みの低減を実現できる雑音消去装置を
提供することである。
【００１８】
【課題を解決するための手段】本発明によれば、所望信
号に雑音信号が混在した主信号に予め定められた遅延時
間を与えて遅延主信号として出力する第１の遅延手段
と、前記雑音信号を参照信号として受け予め定められた
遅延時間を与えて遅延参照信号として出力する第２の遅
延手段と、前記遅延主信号から第１の推定雑音信号を減
算して第１の所望信号出力を生成する第１の減算手段
と、前記遅延参照信号から第１の推定所望信号を減算し
て第１の雑音信号出力を生成する第２の減算手段と、前
記第１の雑音信号出力を入力として前記遅延主信号に含
まれる雑音信号を適応的に推定して前記第１の推定雑音
信号として出力する第１の適応フィルタと、前記第１の
所望信号出力を入力として前記遅延参照信号に含まれる
所望信号を適応的に推定して前記第１の推定所望信号と
して出力する第２の適応フィルタと、前記主信号及び参
照信号を入力として前記主信号における所望信号電力及
び雑音信号電力、更に前記参照信号における所望信号電
力及び雑音信号電力を夫々算出して主信号対雑音信号電
力比の推定値及び参照信号対雑音信号電力比の推定値を
出力する信号対雑音電比推定手段と、前記主信号対雑音
信号電力比の推定値及び前記参照信号対雑音信号電力比
の推定値を夫々入力として前記第１及び第２の適応フィ
ルタのフィルタ係数の修正量を夫々決定するための第１
及び第２のステップサイズを出力するステップサイズ出
力手段とを含むことを特徴とする雑音消去装置が得られ
る。
号に雑音信号が混在した主信号に予め定められた遅延時
間を与えて遅延主信号として出力する第１の遅延手段
と、前記雑音信号を参照信号として受け予め定められた
遅延時間を与えて遅延参照信号として出力する第２の遅
延手段と、前記遅延主信号から第１の推定雑音信号を減
算して第１の所望信号出力を生成する第１の減算手段
と、前記遅延参照信号から第１の推定所望信号を減算し
て第１の雑音信号出力を生成する第２の減算手段と、前
記第１の雑音信号出力を入力として前記遅延主信号に含
まれる雑音信号を適応的に推定して前記第１の推定雑音
信号として出力する第１の適応フィルタと、前記第１の
所望信号出力を入力として前記遅延参照信号に含まれる
所望信号を適応的に推定して前記第１の推定所望信号と
して出力する第２の適応フィルタと、前記主信号及び参
照信号を入力として前記主信号における所望信号電力及
び雑音信号電力、更に前記参照信号における所望信号電
力及び雑音信号電力を夫々算出して主信号対雑音信号電
力比の推定値及び参照信号対雑音信号電力比の推定値を
出力する信号対雑音電比推定手段と、前記主信号対雑音
信号電力比の推定値及び前記参照信号対雑音信号電力比
の推定値を夫々入力として前記第１及び第２の適応フィ
ルタのフィルタ係数の修正量を夫々決定するための第１
及び第２のステップサイズを出力するステップサイズ出
力手段とを含むことを特徴とする雑音消去装置が得られ
る。
【００１９】そして、前記信号対雑音電力比推定手段
は、前記主信号から第２の推定雑音信号を減算して第２
の所望信号出力を生成する第３の減算手段と、前記参照
信号から第２の推定所望信号を減算して第２の雑音信号
出力を生成する第４の減算手段と、前記第２の雑音信号
出力を入力として前記主信号に含まれる雑音信号を適応
的に推定して前記第２の推定雑音信号として出力する第
３の適応フィルタと、前記第２の所望信号出力を入力と
し前記参照信号に含まれる所望信号を適応的に推定して
前記第２の推定所望信号として出力する第４の適応フィ
ルタと、前記第２の所望信号出力を入力としてその自乗
平均値を算出して前記主信号における所望信号電力電力
として出力する第１の電力平均手段と、前記第２の推定
雑音信号を入力としてその自乗平均値を算出して前記主
信号における雑音信号電力として出力する第２の電力平
均手段と、前記第２の推定所望信号を入力としてその自
乗平均値を算出して前記参照信号における所望信号電力
として出力する第３の電力平均手段と、前記第２の雑音
信号出力を入力としてその自乗平均値を算出して前記参
照信号における雑音信号電力として出力する第４の電力
平均手段と、前記主信号における所望信号電力を前記主
信号における雑音信号電力で除算して前記主信号対雑音
信号電力比の推定値を出力する第１の除算手段と、前記
参照信号における所望信号電力を前記参照信号における
雑音信号電力で除算して前記参照信号対雑音信号電力比
の推定値を出力する第２の除算手段とを有することを特
徴としている。
は、前記主信号から第２の推定雑音信号を減算して第２
の所望信号出力を生成する第３の減算手段と、前記参照
信号から第２の推定所望信号を減算して第２の雑音信号
出力を生成する第４の減算手段と、前記第２の雑音信号
出力を入力として前記主信号に含まれる雑音信号を適応
的に推定して前記第２の推定雑音信号として出力する第
３の適応フィルタと、前記第２の所望信号出力を入力と
し前記参照信号に含まれる所望信号を適応的に推定して
前記第２の推定所望信号として出力する第４の適応フィ
ルタと、前記第２の所望信号出力を入力としてその自乗
平均値を算出して前記主信号における所望信号電力電力
として出力する第１の電力平均手段と、前記第２の推定
雑音信号を入力としてその自乗平均値を算出して前記主
信号における雑音信号電力として出力する第２の電力平
均手段と、前記第２の推定所望信号を入力としてその自
乗平均値を算出して前記参照信号における所望信号電力
として出力する第３の電力平均手段と、前記第２の雑音
信号出力を入力としてその自乗平均値を算出して前記参
照信号における雑音信号電力として出力する第４の電力
平均手段と、前記主信号における所望信号電力を前記主
信号における雑音信号電力で除算して前記主信号対雑音
信号電力比の推定値を出力する第１の除算手段と、前記
参照信号における所望信号電力を前記参照信号における
雑音信号電力で除算して前記参照信号対雑音信号電力比
の推定値を出力する第２の除算手段とを有することを特
徴としている。
【００２０】また、前記ステップサイズ出力手段は、前
記主信号対雑音信号電力比の推定値を予め定められた単
調減少する関数に入力して第１の関数値を算出する手段
と、前記第１の関数値が予め定められた第１の最大値と
第１の最小値との間の値にあるときには前記第１の関数
値を、前記第１の最大値より大なるときには前記第１の
最大値を、前記第１の最小値より小なるときには前記第
１の最小値を、夫々前記第１のステップサイズとして出
力する手段と、前記参照信号対雑音信号電力比の推定値
を予め定められた単調増加する関数に入力して第２の関
数値を算出する手段と、前記第２の関数値が予め定めら
れた第２の最大値と第２の最小値との間の値にあるとき
には前記第２の関数値を、前記第２の最大値より大なる
ときには前記第２の最大値を、前記第２の最小値より小
なるときには前記第２の最小値を、夫々前記第２のステ
ップサイズとして出力する手段とを有することを特徴と
すしている。
記主信号対雑音信号電力比の推定値を予め定められた単
調減少する関数に入力して第１の関数値を算出する手段
と、前記第１の関数値が予め定められた第１の最大値と
第１の最小値との間の値にあるときには前記第１の関数
値を、前記第１の最大値より大なるときには前記第１の
最大値を、前記第１の最小値より小なるときには前記第
１の最小値を、夫々前記第１のステップサイズとして出
力する手段と、前記参照信号対雑音信号電力比の推定値
を予め定められた単調増加する関数に入力して第２の関
数値を算出する手段と、前記第２の関数値が予め定めら
れた第２の最大値と第２の最小値との間の値にあるとき
には前記第２の関数値を、前記第２の最大値より大なる
ときには前記第２の最大値を、前記第２の最小値より小
なるときには前記第２の最小値を、夫々前記第２のステ
ップサイズとして出力する手段とを有することを特徴と
すしている。
【００２１】更に、前記第１及び第２のステップサイズ
の値を比較して小なる方のステップサイズをゼロとする
ようにし、また前記第１及び第２の遅延手段の時間遅延
は、前記信号対雑音電力比推定手段において前記主信号
対雑音信号電力比及び参照信号対雑音信号電力比の各推
定値が算出される際に生ずる時間遅延に等しく設定され
ていることを特徴としている。
の値を比較して小なる方のステップサイズをゼロとする
ようにし、また前記第１及び第２の遅延手段の時間遅延
は、前記信号対雑音電力比推定手段において前記主信号
対雑音信号電力比及び参照信号対雑音信号電力比の各推
定値が算出される際に生ずる時間遅延に等しく設定され
ていることを特徴としている。
【００２２】本装置では、主信号に混入する雑音信号を
消去するために主信号中の雑音信号を推定する第１の適
応フィルタと、参照信号に混入する所望信号を消去する
ために参照信号中の所望信号を推定する第２の適応フィ
ルタと、主信号に混入する雑音信号を消去するために主
信号中の雑音信号を推定する第３の適応フィルタと、参
照信号に混入する所望信号を消去するために参照信号中
の所望信号を推定する第４の適応フィルタとを設ける。
消去するために主信号中の雑音信号を推定する第１の適
応フィルタと、参照信号に混入する所望信号を消去する
ために参照信号中の所望信号を推定する第２の適応フィ
ルタと、主信号に混入する雑音信号を消去するために主
信号中の雑音信号を推定する第３の適応フィルタと、参
照信号に混入する所望信号を消去するために参照信号中
の所望信号を推定する第４の適応フィルタとを設ける。
【００２３】第３の適応フィルタ及び第４の適応フィル
タを用いて主信号における信号対雑音電力比と参照信号
における信号対雑音電力比を推定する。そして、主信号
における信号対雑音信号電力比が小さい時には、第１の
適応フィルタが推定すべき主信号中の雑音信号成分が所
望信号に比べ多いと判断し、第１の適応フィルタには大
きなステップサイズを供給し、収束を速める。逆に、主
信号における信号対雑音信号電力比が大きい時には、第
１の適応フィルタが推定すべき雑音信号成分より所望信
号が多いと判断して第１の適応フィルタには小さなステ
ップサイズを供給し、誤った収束へ向かうことを防止す
る。
タを用いて主信号における信号対雑音電力比と参照信号
における信号対雑音電力比を推定する。そして、主信号
における信号対雑音信号電力比が小さい時には、第１の
適応フィルタが推定すべき主信号中の雑音信号成分が所
望信号に比べ多いと判断し、第１の適応フィルタには大
きなステップサイズを供給し、収束を速める。逆に、主
信号における信号対雑音信号電力比が大きい時には、第
１の適応フィルタが推定すべき雑音信号成分より所望信
号が多いと判断して第１の適応フィルタには小さなステ
ップサイズを供給し、誤った収束へ向かうことを防止す
る。
【００２４】同様に、参照信号における信号対雑音電力
比が大きい時には、第２の適応フィルタが推定すべき参
照信号中の所望信号成分が雑音信号に比べ多いと判断
し、第２の適応フィルタには大きなステップサイズを供
給し、収束を速める。逆に、参照信号における信号対雑
音信号電力比が小さい時には、第２の適応フィルタが推
定すべき所望信号成分より雑音信号が多いと判断して第
２の適応フィルタには小さなステップサイズを供給し、
誤った収束へ向かうことを防止する。
比が大きい時には、第２の適応フィルタが推定すべき参
照信号中の所望信号成分が雑音信号に比べ多いと判断
し、第２の適応フィルタには大きなステップサイズを供
給し、収束を速める。逆に、参照信号における信号対雑
音信号電力比が小さい時には、第２の適応フィルタが推
定すべき所望信号成分より雑音信号が多いと判断して第
２の適応フィルタには小さなステップサイズを供給し、
誤った収束へ向かうことを防止する。
【００２５】
【発明の実施の形態】次に、本発明の実施の形態につい
て図面と共に説明する。
て図面と共に説明する。
【００２６】図１は本発明の一実施の形態のブロック図
を示す。図１に示す様に、本実施の形態は、主信号用マ
イクロフォン１，参照信号用マイクロフォン２，出力端
子５，適応フィルタ３，適応フィルタ６，減算器４，減
算器７，遅延回路８，遅延回路９，信号対雑音電力比推
定回路１０，ステップサイズ出力回路１１を備えてい
る。
を示す。図１に示す様に、本実施の形態は、主信号用マ
イクロフォン１，参照信号用マイクロフォン２，出力端
子５，適応フィルタ３，適応フィルタ６，減算器４，減
算器７，遅延回路８，遅延回路９，信号対雑音電力比推
定回路１０，ステップサイズ出力回路１１を備えてい
る。
【００２７】先ず、適応フィルタ３及び適応フィルタ６
の動作について説明する。主信号ＸＡ（ｚ）は、遅延回
路８によってＤサンプルの遅延を与えられ、ＸＡ（ｚ）
Ｚ-Dとなって、減算器４に入力される。但し、Ｚ-DはＤ
サンプルの遅延を表す。
の動作について説明する。主信号ＸＡ（ｚ）は、遅延回
路８によってＤサンプルの遅延を与えられ、ＸＡ（ｚ）
Ｚ-Dとなって、減算器４に入力される。但し、Ｚ-DはＤ
サンプルの遅延を表す。
【００２８】一方、参照信号ＸＢ（ｚ）は遅延回路９に
よってＤサンプルの遅延を与えられ、ＸＢ（ｚ）Ｚ-Dと
なって減算器７に入力される。このＤサンプルの遅延は
後述する信号対雑音電力推定回路１０で算出される信号
対雑音電力比の算出のために生じる時間的遅延を補償す
るためのものであり、適応フィルタ３と適応フィルタ６
の主信号側と参照信号側の両方に同じ遅延量が挿入され
ているため、主信号と参照信号の相対関係には影響を与
えない。従って、以下の説明では、Ｄを零と仮定する。
よってＤサンプルの遅延を与えられ、ＸＢ（ｚ）Ｚ-Dと
なって減算器７に入力される。このＤサンプルの遅延は
後述する信号対雑音電力推定回路１０で算出される信号
対雑音電力比の算出のために生じる時間的遅延を補償す
るためのものであり、適応フィルタ３と適応フィルタ６
の主信号側と参照信号側の両方に同じ遅延量が挿入され
ているため、主信号と参照信号の相対関係には影響を与
えない。従って、以下の説明では、Ｄを零と仮定する。
【００２９】適応フィルタ３は主信号ＸＡ（ｚ）に含ま
れている雑音信号を推定する様に動作し、適応フィルタ
６は参照信号ＸＢ（ｚ）に含まれる所望信号を推定する
様に動作する。この時、主信号中の雑音信号を推定する
適応フィルタ３の入力としては、参照信号から適応フィ
ルタ６が推定した所望信号を減算して得られた雑音信号
が用いられる。
れている雑音信号を推定する様に動作し、適応フィルタ
６は参照信号ＸＢ（ｚ）に含まれる所望信号を推定する
様に動作する。この時、主信号中の雑音信号を推定する
適応フィルタ３の入力としては、参照信号から適応フィ
ルタ６が推定した所望信号を減算して得られた雑音信号
が用いられる。
【００３０】同様に、参照信号中の所望信号を推定する
適応フィルタ６の入力としては、主信号から、適応フィ
ルタ３が推定した雑音信号を減算して得られた所望信号
が用いられる。従って、適応フィルタ３と適応フィルタ
６は図１に示す通り、交差接続（クロスカップル）の形
態となっている。
適応フィルタ６の入力としては、主信号から、適応フィ
ルタ３が推定した雑音信号を減算して得られた所望信号
が用いられる。従って、適応フィルタ３と適応フィルタ
６は図１に示す通り、交差接続（クロスカップル）の形
態となっている。
【００３１】今、減算器４の出力信号をＥ1 （ｚ）、減
算器７の出力信号をＥ2 （ｚ）とし、適応フィルタ３の
伝達特性とフィルタ出力を夫々、Ｗ1 （ｚ），Ｆ1
（ｚ）とし、適応フィルタ６の伝達特性とフィルタ出力
を夫々、Ｗ2 （ｚ），Ｆ2 （ｚ）とすると、 Ｅ1 （ｚ）＝ＸＡ（ｚ）−Ｆ1 （ｚ） ＝ＸＡ（ｚ）−Ｗ1 （ｚ）Ｅ2 （ｚ）…（８） Ｅ2 （ｚ）＝ＸＢ（ｚ）−Ｆ2 （ｚ） ＝ＸＢ（ｚ）−Ｗ2 （ｚ）Ｅ1 （ｚ）…（９） が得られる。
算器７の出力信号をＥ2 （ｚ）とし、適応フィルタ３の
伝達特性とフィルタ出力を夫々、Ｗ1 （ｚ），Ｆ1
（ｚ）とし、適応フィルタ６の伝達特性とフィルタ出力
を夫々、Ｗ2 （ｚ），Ｆ2 （ｚ）とすると、 Ｅ1 （ｚ）＝ＸＡ（ｚ）−Ｆ1 （ｚ） ＝ＸＡ（ｚ）−Ｗ1 （ｚ）Ｅ2 （ｚ）…（８） Ｅ2 （ｚ）＝ＸＢ（ｚ）−Ｆ2 （ｚ） ＝ＸＢ（ｚ）−Ｗ2 （ｚ）Ｅ1 （ｚ）…（９） が得られる。
【００３２】また、主信号ＸＡ（ｚ）と参照信号ＸＢ
（ｚ）は、図５で用いた所望信号Ｓ（ｚ），雑音Ｎ
（ｚ），音響伝達特性ＨＡ（ｚ），ＨＢ（ｚ），ＧＡ
（ｚ），ＧＢ（ｚ）を使って、 ＸＡ（ｚ）＝ＧＡ（ｚ）Ｎ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ）…（１０） ＸＢ（ｚ）＝ＧＢ（ｚ）Ｎ（ｚ）＋ＨＢ（ｚ）Ｓ（ｚ）…（１１） と表される。
（ｚ）は、図５で用いた所望信号Ｓ（ｚ），雑音Ｎ
（ｚ），音響伝達特性ＨＡ（ｚ），ＨＢ（ｚ），ＧＡ
（ｚ），ＧＢ（ｚ）を使って、 ＸＡ（ｚ）＝ＧＡ（ｚ）Ｎ（ｚ）＋ＨＡ（ｚ）Ｓ（ｚ）…（１０） ＸＢ（ｚ）＝ＧＢ（ｚ）Ｎ（ｚ）＋ＨＢ（ｚ）Ｓ（ｚ）…（１１） と表される。
【００３３】これ等の式から、Ｅ1 （ｚ）及びＥ2
（ｚ）を求めると、 Ｅ1 （ｚ）＝［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＨＡ（ｚ）−Ｗ1 （ｚ）ＨＢ（ｚ）｝Ｓ（ｚ） ＋［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＧＡ（ｚ）−Ｗ1 （ｚ）ＧＢ（ｚ）｝Ｎ（ｚ）…（１２） Ｅ2 （ｚ）＝［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＨＢ（ｚ）−Ｗ1 （ｚ）ＨＡ（ｚ）｝Ｓ（ｚ） ＋［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＧＢ（ｚ）−Ｗ1 （ｚ）ＧＡ（ｚ）｝Ｎ（ｚ）…（１３） となる。
（ｚ）を求めると、 Ｅ1 （ｚ）＝［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＨＡ（ｚ）−Ｗ1 （ｚ）ＨＢ（ｚ）｝Ｓ（ｚ） ＋［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＧＡ（ｚ）−Ｗ1 （ｚ）ＧＢ（ｚ）｝Ｎ（ｚ）…（１２） Ｅ2 （ｚ）＝［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＨＢ（ｚ）−Ｗ1 （ｚ）ＨＡ（ｚ）｝Ｓ（ｚ） ＋［１／｛１−Ｗ1 （ｚ）Ｗ2 （ｚ）｝］ ｛ＧＢ（ｚ）−Ｗ1 （ｚ）ＧＡ（ｚ）｝Ｎ（ｚ）…（１３） となる。
【００３４】従って、下記の式 Ｗ1 （ｚ）＝ＧＡ（ｚ）／ＧＢ（ｚ）…（１５） Ｗ2 （ｚ）＝ＨＢ（ｚ）／ＨＡ（ｚ）…（１６） が満足されれば、 Ｅ1 （ｚ）＝Ｓ（ｚ）…（１７） Ｅ2 （ｚ）＝Ｎ（ｚ）…（１８） となって、Ｅ1 （ｚ）には、雑音の消去された所望信号
が得られる。
が得られる。
【００３５】ところで、適応フィルタ３が主信号中の雑
音信号を正確に推定するためには、推定の妨害になる主
信号中の所望信号が推定すべき雑音信号より相対的に小
さい時に適応フィルタの係数の更新量を大きくする必要
がある。逆に、主信号中の所望信号が雑音信号より大き
い時には、推定の妨害となる信号の方が多いので適応フ
ィルタの係数の更新量を下げる必要がある。
音信号を正確に推定するためには、推定の妨害になる主
信号中の所望信号が推定すべき雑音信号より相対的に小
さい時に適応フィルタの係数の更新量を大きくする必要
がある。逆に、主信号中の所望信号が雑音信号より大き
い時には、推定の妨害となる信号の方が多いので適応フ
ィルタの係数の更新量を下げる必要がある。
【００３６】一方、適応フィルタ６が参照信号中の所望
信号を正確に推定するためには、推定の妨害になる参照
信号中の雑音信号が推定すべき所望信号により相対的に
小さい時に適応フィルタの係数の更新量を大きくする必
要がある。逆に、参照信号中の雑音信号が所望信号より
大きい時には、推定の妨害となる信号の方が多いので適
応フィルタの係数の更新量を下げる必要がある。
信号を正確に推定するためには、推定の妨害になる参照
信号中の雑音信号が推定すべき所望信号により相対的に
小さい時に適応フィルタの係数の更新量を大きくする必
要がある。逆に、参照信号中の雑音信号が所望信号より
大きい時には、推定の妨害となる信号の方が多いので適
応フィルタの係数の更新量を下げる必要がある。
【００３７】こうした、適応フィルタの係数の更新量の
制約は、適応フィルタの学習アルゴリズムにおけるステ
ップサイズを制御することで可能となる。
制約は、適応フィルタの学習アルゴリズムにおけるステ
ップサイズを制御することで可能となる。
【００３８】ここで、適応フィルタの学習アルゴリズム
として、「学習同定法」を用いるものと仮定し、適応フ
ィルタ３を例にとり、係数更新の方法を説明する。今、
主信号ＸＡ（ｚ）の時間軸表現xa（ｋ），適応フィルタ
３の入力であるＥ2 （ｚ）の時間軸表現をe2（ｋ），適
応フィルタ３の出力であるＦ1 （ｚ）の時間軸表現をf1
（ｋ），減算器４の出力であるＥ1 （ｚ）の時間軸表現
をe1（ｋ）とする。但し、ｋは時刻を表す指標である。
として、「学習同定法」を用いるものと仮定し、適応フ
ィルタ３を例にとり、係数更新の方法を説明する。今、
主信号ＸＡ（ｚ）の時間軸表現xa（ｋ），適応フィルタ
３の入力であるＥ2 （ｚ）の時間軸表現をe2（ｋ），適
応フィルタ３の出力であるＦ1 （ｚ）の時間軸表現をf1
（ｋ），減算器４の出力であるＥ1 （ｚ）の時間軸表現
をe1（ｋ）とする。但し、ｋは時刻を表す指標である。
【００３９】時刻ｋにおける適応フィルタ３のｊ番目の
係数をｗ1j（ｋ）とすると、適応フィルタ３の出力する
推定雑音信号f1（ｋ）は、次式で表現される。ここで、
Ｎは適応フィルタ３のタップ数を表す。
係数をｗ1j（ｋ）とすると、適応フィルタ３の出力する
推定雑音信号f1（ｋ）は、次式で表現される。ここで、
Ｎは適応フィルタ３のタップ数を表す。
【００４０】
【数１】 減算器４で求められた誤差信号e1（ｋ）を用いて、時刻
（ｋ＋１）における係数ｗ1j（ｋ＋１）は次式で計算さ
れる。
（ｋ＋１）における係数ｗ1j（ｋ＋１）は次式で計算さ
れる。
【００４１】
【数２】 上式において、μ1 （ｋ）は、適応フィルタ３における
係数更新のステップサイズである。ステップサイズμ1
（ｋ）は、その値が大きい場合には、係数の修正量が多
くなるため収束が速くなる一方、修正量が大きい分だけ
係数更新の妨害となる成分が存在する場合にはその影響
を強く受けて残留誤差量が多くなる。反対に、ステップ
サイズの値が小さい場合には、収束時間が増加するが、
妨害信号成分の影響が少なく残留誤差量は小さくなる。
従って、ステップサイズの設定には、「収束時間」と
「残留誤差」にトレードオフが存在する。
係数更新のステップサイズである。ステップサイズμ1
（ｋ）は、その値が大きい場合には、係数の修正量が多
くなるため収束が速くなる一方、修正量が大きい分だけ
係数更新の妨害となる成分が存在する場合にはその影響
を強く受けて残留誤差量が多くなる。反対に、ステップ
サイズの値が小さい場合には、収束時間が増加するが、
妨害信号成分の影響が少なく残留誤差量は小さくなる。
従って、ステップサイズの設定には、「収束時間」と
「残留誤差」にトレードオフが存在する。
【００４２】同様に、適応フィルタ６についても、今、
参照信号ＸＢ（ｚ）の時間軸表現をxb（ｋ），適応フィ
ルタ６の入力であるＥ1 （ｚ）の時間軸表現をe1
（ｋ），適応フィルタ６の出力であるＦ2 （ｚ）の時間
軸表現をf2（ｋ），減算器７の出力であるＥ2 （ｚ）の
時間軸表現をe2（ｋ）とし、適応フィルタ６のｊ番目の
係数をｗ2j（ｋ）とすると、適応フィルタ６の出力する
推定所望信号f2（ｋ）は、次式で表現される。
参照信号ＸＢ（ｚ）の時間軸表現をxb（ｋ），適応フィ
ルタ６の入力であるＥ1 （ｚ）の時間軸表現をe1
（ｋ），適応フィルタ６の出力であるＦ2 （ｚ）の時間
軸表現をf2（ｋ），減算器７の出力であるＥ2 （ｚ）の
時間軸表現をe2（ｋ）とし、適応フィルタ６のｊ番目の
係数をｗ2j（ｋ）とすると、適応フィルタ６の出力する
推定所望信号f2（ｋ）は、次式で表現される。
【００４３】
【数３】 減算器７で求められた誤差信号e2（ｋ）を用いて、時刻
（ｋ＋１）における係数ｗ2j（ｋ＋１）は次式で計算さ
れる。
（ｋ＋１）における係数ｗ2j（ｋ＋１）は次式で計算さ
れる。
【００４４】
【数４】 上式において、μ2 （ｋ）は、適応フィルタ６における
係数更新のステップサイズである。
係数更新のステップサイズである。
【００４５】以上説明した様に、適応フィルタのステッ
プサイズを制御することによって、係数の更新量は制御
できる。
プサイズを制御することによって、係数の更新量は制御
できる。
【００４６】次に、信号対雑音電力比推定回路１０の動
作について説明する。図２は信号対雑音電力比推定回路
１０を示すブロック図である。信号対雑音電力比推定回
路１０は、適応フィルタ１２，適応フィルタ１３，減算
器１４，減算器１５，電力平均回路１６，電力平均回路
１７，電力平均回路１８，電力平均回路１９，除算回路
２０，除算回路２１から構成されている。
作について説明する。図２は信号対雑音電力比推定回路
１０を示すブロック図である。信号対雑音電力比推定回
路１０は、適応フィルタ１２，適応フィルタ１３，減算
器１４，減算器１５，電力平均回路１６，電力平均回路
１７，電力平均回路１８，電力平均回路１９，除算回路
２０，除算回路２１から構成されている。
【００４７】信号対雑音電力比推定回路１０の中の２つ
の適応フィルタと２つの減算器の接続は、前述した適応
フィルタのクロスカップル形態と全く同じ形になってい
る。但し、適応フィルタ１２のステップサイズμ3 と適
応フィルタ１３のステップサイズμ4 は収束速度を速め
るために大きめの値の固定値を設定するものとし、係数
更新のアルゴリズムとして、「学習同定法」を仮定する
場合には、ステップサイズμとして例えば、「０．２」
から「０．５」程度の値に設定される。この様に設定さ
れた適応フィルタ１２と適応フィルタ１３は、大きい値
に設定されたステップサイズのために、残留誤差が大き
くなるが、高速収束するという利点がある。
の適応フィルタと２つの減算器の接続は、前述した適応
フィルタのクロスカップル形態と全く同じ形になってい
る。但し、適応フィルタ１２のステップサイズμ3 と適
応フィルタ１３のステップサイズμ4 は収束速度を速め
るために大きめの値の固定値を設定するものとし、係数
更新のアルゴリズムとして、「学習同定法」を仮定する
場合には、ステップサイズμとして例えば、「０．２」
から「０．５」程度の値に設定される。この様に設定さ
れた適応フィルタ１２と適応フィルタ１３は、大きい値
に設定されたステップサイズのために、残留誤差が大き
くなるが、高速収束するという利点がある。
【００４８】今、適応フィルタ１２と適応フィルタ１３
が共に収束していると仮定すると、適応フィルタ１２の
出力f3（ｋ）は、主信号中の雑音信号を出力している。
また、減算器１４の出力e3（ｋ）は、主信号中の所望信
号を出力している。電力平均回路１６は、減算器１４の
出力を自乗して、その時間平均を求め、主信号における
所望信号電力ＰＳＡ（ｋ）を出力する。また、電力平均
回路１７は、適応フィルタ１２の出力を自乗して、その
時間平均を求め、主信号における雑音信号電力ＰＮＡ
（ｋ）を出力する。
が共に収束していると仮定すると、適応フィルタ１２の
出力f3（ｋ）は、主信号中の雑音信号を出力している。
また、減算器１４の出力e3（ｋ）は、主信号中の所望信
号を出力している。電力平均回路１６は、減算器１４の
出力を自乗して、その時間平均を求め、主信号における
所望信号電力ＰＳＡ（ｋ）を出力する。また、電力平均
回路１７は、適応フィルタ１２の出力を自乗して、その
時間平均を求め、主信号における雑音信号電力ＰＮＡ
（ｋ）を出力する。
【００４９】同様に、適応フィルタ１３の出力f4（ｋ）
は、参照信号中の所望信号を出力している。また、減算
器１５の出力e4（ｋ）は、参照信号中の雑音信号を出力
している。電力平均回路１９は、減算器１５の出力を自
乗して、その時間平均を求め、参照信号における雑音信
号電力ＰＮＢ（ｋ）を出力する。また、電力平均回路１
８は、適応フィルタ１３の出力を自乗して、その時間平
均を求め、参照信号における所望信号電力ＰＳＢ（ｋ）
を出力する。
は、参照信号中の所望信号を出力している。また、減算
器１５の出力e4（ｋ）は、参照信号中の雑音信号を出力
している。電力平均回路１９は、減算器１５の出力を自
乗して、その時間平均を求め、参照信号における雑音信
号電力ＰＮＢ（ｋ）を出力する。また、電力平均回路１
８は、適応フィルタ１３の出力を自乗して、その時間平
均を求め、参照信号における所望信号電力ＰＳＢ（ｋ）
を出力する。
【００５０】除算回路２０は、電力平均回路１６の出力
する所望信号電力を電力平均回路１７の出力する雑音信
号電力で除算し、結果として主信号における信号対雑音
電力比の推定値ＳＮＲＡ（ｋ）を出力する。
する所望信号電力を電力平均回路１７の出力する雑音信
号電力で除算し、結果として主信号における信号対雑音
電力比の推定値ＳＮＲＡ（ｋ）を出力する。
【００５１】同様に、除算回路２１は、電力平均回路１
８の出力する所望信号電力を電力平均回路１９の出力す
る雑音信号電力で除算し、結果として参照信号における
信号対雑音電力比の推定値ＳＮＲＢ（ｋ）を出力する。
８の出力する所望信号電力を電力平均回路１９の出力す
る雑音信号電力で除算し、結果として参照信号における
信号対雑音電力比の推定値ＳＮＲＢ（ｋ）を出力する。
【００５２】ところで、電力平均回路１６，１７，１
８，１９における平均処理を移動平均等で行う場合を仮
定すると、算出された電力平均値は実際の電力変化に対
して平均回数に依存する遅延ΔAVを生じてしまう。そこ
で、この実施の形態では、この遅延ΔAVを補償するた
め、遅延回路８，９を備えている。従って、遅延回路
８，９の遅延量Ｚ-Dは、ΔAVに等しく設定することが望
ましい。
８，１９における平均処理を移動平均等で行う場合を仮
定すると、算出された電力平均値は実際の電力変化に対
して平均回数に依存する遅延ΔAVを生じてしまう。そこ
で、この実施の形態では、この遅延ΔAVを補償するた
め、遅延回路８，９を備えている。従って、遅延回路
８，９の遅延量Ｚ-Dは、ΔAVに等しく設定することが望
ましい。
【００５３】以上、説明した様に信号対雑音電力比推定
回路１０は、クロスカップル接続された適応フィルタの
ステップサイズを大きい値に設定することにより高速収
束させ、収束後の各信号を用いて主信号における信号対
雑音電力比の推定値ＳＮＲＡ（ｋ）と参照信号における
信号対雑音電力比の推定値ＳＮＲＢ（ｋ）を出力する。
回路１０は、クロスカップル接続された適応フィルタの
ステップサイズを大きい値に設定することにより高速収
束させ、収束後の各信号を用いて主信号における信号対
雑音電力比の推定値ＳＮＲＡ（ｋ）と参照信号における
信号対雑音電力比の推定値ＳＮＲＢ（ｋ）を出力する。
【００５４】次に、ステップサイズ出力回路１１の動作
について図３，４のフローチャートを用いて説明する。
ステップサイズ出力回路１１では、先ず、信号対雑音電
力比推定回路１０から出力された主信号における信号対
雑音電力比の推定値ＳＮＲＡ（ｋ）が単調減少関数へ入
力される（ステップ３１）。今、ｆ（・）をＳＮＲＡ
（ｋ）に対する単調減少関数とすると、単調減少関数の
出力ＯＵＴ1 （ｋ）は、以下の式で表される（ステップ
３２）。
について図３，４のフローチャートを用いて説明する。
ステップサイズ出力回路１１では、先ず、信号対雑音電
力比推定回路１０から出力された主信号における信号対
雑音電力比の推定値ＳＮＲＡ（ｋ）が単調減少関数へ入
力される（ステップ３１）。今、ｆ（・）をＳＮＲＡ
（ｋ）に対する単調減少関数とすると、単調減少関数の
出力ＯＵＴ1 （ｋ）は、以下の式で表される（ステップ
３２）。
【００５５】 ＯＵＴ1 （ｋ）＝ｆ（ＳＮＲＡ（ｋ））…（２３） 関数値ＯＵＴ1 （ｋ）を用いて、適応フィルタ３のステ
ップサイズμ1 （ｋ）は以下の式により算出される。
ップサイズμ1 （ｋ）は以下の式により算出される。
【００５６】 μ1 （ｋ）＝clip［ＯＵＴ1 （ｋ），μ1max，μ1min］…（２４） ここで、clip［ａ，ｂ，ｃ］は、最小値，最大値を設定
するための関数で、次式の様に定義される。
するための関数で、次式の様に定義される。
【００５７】 clip［ａ，ｂ，ｃ］＝ａ（ｃ≦ａ≦ｂ） clip［ａ，ｂ，ｃ］＝ｂ（ａ＞ｂ） clip［ａ，ｂ，ｃ］＝ｃ（ａ＜ｃ） …（２５） この様に、最大値μ1maxと最小値μ1minの設定によるス
テップサイズの制限は、適応フィルタの安定した動作の
ために有効である。適応フィルタ３に関しては、主信号
における信号対雑音電力比の推定値ＳＮＲＡ（ｋ）を単
調減少関数へ入力して得られた関数値をステップサイズ
に設定することにより、信号対雑音電力比が大きい時は
小さいステップサイズを設定し、逆に、信号対雑音電力
比が小さい時は大きなステップサイズを設定することに
なる（ステップ３３〜３６）。
テップサイズの制限は、適応フィルタの安定した動作の
ために有効である。適応フィルタ３に関しては、主信号
における信号対雑音電力比の推定値ＳＮＲＡ（ｋ）を単
調減少関数へ入力して得られた関数値をステップサイズ
に設定することにより、信号対雑音電力比が大きい時は
小さいステップサイズを設定し、逆に、信号対雑音電力
比が小さい時は大きなステップサイズを設定することに
なる（ステップ３３〜３６）。
【００５８】同様に、参照信号における信号対雑音電力
比の推定値ＳＮＲＢ（ｋ）は、単調増加関数へ入力され
る（ステップ４１）。今、ｇ（・）をＳＮＲＢ（ｋ）に
対する単調増加関数とすると、単調増加関数の出力ＯＵ
Ｔ2 （ｋ）は、以下の式で表される（ステップ４２）。
比の推定値ＳＮＲＢ（ｋ）は、単調増加関数へ入力され
る（ステップ４１）。今、ｇ（・）をＳＮＲＢ（ｋ）に
対する単調増加関数とすると、単調増加関数の出力ＯＵ
Ｔ2 （ｋ）は、以下の式で表される（ステップ４２）。
【００５９】 ＯＵＴ2 （ｋ）＝ｇ（ＳＮＲＢ（ｋ））…（２６） 関数値ＯＵＴ2 （ｋ）を用いて、適応フィルタ６のステ
ップサイズμ2 （ｋ）は以下の式により算出される。
ップサイズμ2 （ｋ）は以下の式により算出される。
【００６０】 μ2 （ｋ）＝clip［ＯＵＴ2 （ｋ），μ2max，μ2min］…（２７） 適応フィルタ６に関しては、参照信号における信号対雑
音電力比の推定値ＳＮＲＢ（ｋ）を単調増加関数へ入力
して得られた関数値をステップサイズに設定することに
より、信号対雑音電力比が大きい時は大きなステップサ
イズを設定し、逆に、信号対雑音電力比が小さい時は小
さなステップサイズを設定することになる（ステップ４
３〜４６）。
音電力比の推定値ＳＮＲＢ（ｋ）を単調増加関数へ入力
して得られた関数値をステップサイズに設定することに
より、信号対雑音電力比が大きい時は大きなステップサ
イズを設定し、逆に、信号対雑音電力比が小さい時は小
さなステップサイズを設定することになる（ステップ４
３〜４６）。
【００６１】以上、説明した様にステップサイズ出力回
路１１は、主信号における信号対雑音電力比の推定値Ｓ
ＮＲＡ（ｋ）の値に応じて適応フィルタ３に供給するス
テップサイズを制御する。同様に、参照信号における信
号対雑音電力比の推定値ＳＮＲＢ（ｋ）の値に応じて適
応フィルタ６に供給するステップサイズを制御する。
路１１は、主信号における信号対雑音電力比の推定値Ｓ
ＮＲＡ（ｋ）の値に応じて適応フィルタ３に供給するス
テップサイズを制御する。同様に、参照信号における信
号対雑音電力比の推定値ＳＮＲＢ（ｋ）の値に応じて適
応フィルタ６に供給するステップサイズを制御する。
【００６２】また、ステップサイズ出力回路１１で求め
られた適応フィルタ３のステップサイズμ1 （ｋ）と、
適応フィルタ６のステップサイズμ2 （ｋ）の値を比較
し、小さい方のステップサイズを零に設定することで、
一方の適応フィルタの学習動作を停止することができ
る。この制御は、同時に二つの適応フィルタが学習動作
を行うことによる干渉を低減することができるため、よ
り正確な学習のために有効である。
られた適応フィルタ３のステップサイズμ1 （ｋ）と、
適応フィルタ６のステップサイズμ2 （ｋ）の値を比較
し、小さい方のステップサイズを零に設定することで、
一方の適応フィルタの学習動作を停止することができ
る。この制御は、同時に二つの適応フィルタが学習動作
を行うことによる干渉を低減することができるため、よ
り正確な学習のために有効である。
【００６３】
【発明の効果】以上説明した様に、本発明によれば主信
号における信号対雑音電力比の推定値から主信号中の雑
音信号を推定する適応フィルタにとって妨害信号となる
所望信号と消去すべき雑音信号の大小関係を求め、適応
フィルタに供給するステップサイズを制御する様にした
ため、主信号中の雑音信号を正確に推定できる。同様
に、参照信号における信号対雑音電力比の推定値から参
照信号中の所望信号を推定する適応フィルタにとって妨
害信号となる雑音信号と消去すべき所望信号の大小関係
を求め、適応フィルタに供給するステップサイズを制御
する様にしたため、参照信号中の所望信号を正確に推定
できる。結果として、参照信号に所望信号が洩れ込んで
いる場合において、外部から適応フィルタの学習動作を
指示することなく、高速収束と信号歪みの少ない雑音消
去装置を実現できる。
号における信号対雑音電力比の推定値から主信号中の雑
音信号を推定する適応フィルタにとって妨害信号となる
所望信号と消去すべき雑音信号の大小関係を求め、適応
フィルタに供給するステップサイズを制御する様にした
ため、主信号中の雑音信号を正確に推定できる。同様
に、参照信号における信号対雑音電力比の推定値から参
照信号中の所望信号を推定する適応フィルタにとって妨
害信号となる雑音信号と消去すべき所望信号の大小関係
を求め、適応フィルタに供給するステップサイズを制御
する様にしたため、参照信号中の所望信号を正確に推定
できる。結果として、参照信号に所望信号が洩れ込んで
いる場合において、外部から適応フィルタの学習動作を
指示することなく、高速収束と信号歪みの少ない雑音消
去装置を実現できる。
【図１】本発明の一実施の形態のブロック図である。
【図２】図１における信号対雑音電力比推定回路のブロ
ック図である。
ック図である。
【図３】図１のステップサイズ出力回路１１の動作を示
すフローチャートである。
すフローチャートである。
【図４】図１のステップサイズ出力回路１１の動作を示
すフローチャートである。
すフローチャートである。
【図５】従来の雑音消去装置の構成を示す図である。
【図６】雑音消去装置の原理を示すブロック図である。
１ 主信号用マイクロフォン ２ 参照信号用マイクロフォン ８，９ 遅延回路 ３，６，１２，１３ 適応フィルタ ４，７，１４，１５ 減算器 ５ 出力端子 １０ 信号対雑音電力比推定回路 １１ ステップサイズ出力回路 １６，１７，１８，１９ 電力平均回路 ２０，２１ 除算回路
Claims (6)
- 【請求項１】 所望信号に雑音信号が混在した主信号に
予め定められた遅延時間を与えて遅延主信号として出力
する第１の遅延手段と、 前記雑音信号を参照信号として受け予め定められた遅延
時間を与えて遅延参照信号として出力する第２の遅延手
段と、 前記遅延主信号から第１の推定雑音信号を減算して第１
の所望信号出力を生成する第１の減算手段と、 前記遅延参照信号から第１の推定所望信号を減算して第
１の雑音信号出力を生成する第２の減算手段と、 前記第１の雑音信号出力を入力として前記遅延主信号に
含まれる雑音信号を適応的に推定して前記第１の推定雑
音信号として出力する第１の適応フィルタと、 前記第１の所望信号出力を入力として前記遅延参照信号
に含まれる所望信号を適応的に推定して前記第１の推定
所望信号として出力する第２の適応フィルタと、 前記主信号及び参照信号を入力として前記主信号におけ
る所望信号電力及び雑音信号電力、更に前記参照信号に
おける所望信号電力及び雑音信号電力を夫々算出して主
信号対雑音信号電力比の推定値及び参照信号対雑音信号
電力比の推定値を出力する信号対雑音電比推定手段と、 前記主信号対雑音信号電力比の推定値及び前記参照信号
対雑音信号電力比の推定値を夫々入力として前記第１及
び第２の適応フィルタのフィルタ係数の修正量を夫々決
定するための第１及び第２のステップサイズを出力する
ステップサイズ出力手段と、を含むことを特徴とする雑
音消去装置。 - 【請求項２】 前記信号対雑音電力比推定手段は、 前記主信号から第２の推定雑音信号を減算して第２の所
望信号出力を生成する第３の減算手段と、 前記参照信号から第２の推定所望信号を減算して第２の
雑音信号出力を生成する第４の減算手段と、 前記第２の雑音信号出力を入力として前記主信号に含ま
れる雑音信号を適応的に推定して前記第２の推定雑音信
号として出力する第３の適応フィルタと、 前記第２の所望信号出力を入力とし前記参照信号に含ま
れる所望信号を適応的に推定して前記第２の推定所望信
号として出力する第４の適応フィルタと、 前記第２の所望信号出力を入力としてその自乗平均値を
算出して前記主信号における所望信号電力電力として出
力する第１の電力平均手段と、 前記第２の推定雑音信号を入力としてその自乗平均値を
算出して前記主信号における雑音信号電力として出力す
る第２の電力平均手段と、 前記第２の推定所望信号を入力としてその自乗平均値を
算出して前記参照信号における所望信号電力として出力
する第３の電力平均手段と、 前記第２の雑音信号出力を入力としてその自乗平均値を
算出して前記参照信号における雑音信号電力として出力
する第４の電力平均手段と、 前記主信号における所望信号電力を前記主信号における
雑音信号電力で除算して前記主信号対雑音信号電力比の
推定値を出力する第１の除算手段と、 前記参照信号における所望信号電力を前記参照信号にお
ける雑音信号電力で除算して前記参照信号対雑音信号電
力比の推定値を出力する第２の除算手段と、を有するこ
とを特徴とする請求項１記載の雑音除去装置。 - 【請求項３】 前記ステップサイズ出力手段は、 前記主信号対雑音信号電力比の推定値を予め定められた
単調減少する関数に入力して第１の関数値を算出する手
段と、 前記第１の関数値が予め定められた第１の最大値と第１
の最小値との間の値にあるときには前記第１の関数値
を、前記第１の最大値より大なるときには前記第１の最
大値を、前記第１の最小値より小なるときには前記第１
の最小値を、夫々前記第１のステップサイズとして出力
する手段と、を有することを特徴とする請求項１または
２記載の雑音消去装置。 - 【請求項４】 前記ステップサイズ出力手段は、 前記参照信号対雑音信号電力比の推定値を予め定められ
た単調増加する関数に入力して第２の関数値を算出する
手段と、 前記第２の関数値が予め定められた第２の最大値と第２
の最小値との間の値にあるときには前記第２の関数値
を、前記第２の最大値より大なるときには前記第２の最
大値を、前記第２の最小値より小なるときには前記第２
の最小値を、夫々前記第２のステップサイズとして出力
する手段と、を有することを特徴とする請求項１〜３い
ずれか記載の雑音消去装置。 - 【請求項５】 前記第１及び第２のステップサイズの値
を比較して小なる方のステップサイズをゼロとするよう
にしたことを特徴とする請求項４記載の雑音消去装置。 - 【請求項６】 前記第１及び第２の遅延手段の時間遅延
は、前記信号対雑音電力比推定手段において前記主信号
対雑音信号電力比及び参照信号対雑音信号電力比の各推
定値が算出される際に生ずる時間遅延に等しく設定され
ていることを特徴とする請求項１〜５いずれか記載の雑
音消去装置。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP9014410A JP2930101B2 (ja) | 1997-01-29 | 1997-01-29 | 雑音消去装置 |
| CA002228121A CA2228121C (en) | 1997-01-29 | 1998-01-28 | Noise canceller |
| AU52787/98A AU727476B2 (en) | 1997-01-29 | 1998-01-28 | Noise canceler |
| EP98101469A EP0856834A3 (en) | 1997-01-29 | 1998-01-28 | Noise canceler |
| US09/015,525 US5978824A (en) | 1997-01-29 | 1998-01-29 | Noise canceler |
| KR1019980002544A KR100307882B1 (ko) | 1997-01-29 | 1998-01-30 | 잡음제거장치 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP9014410A JP2930101B2 (ja) | 1997-01-29 | 1997-01-29 | 雑音消去装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10215194A true JPH10215194A (ja) | 1998-08-11 |
| JP2930101B2 JP2930101B2 (ja) | 1999-08-03 |
Family
ID=11860284
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP9014410A Expired - Lifetime JP2930101B2 (ja) | 1997-01-29 | 1997-01-29 | 雑音消去装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US5978824A (ja) |
| EP (1) | EP0856834A3 (ja) |
| JP (1) | JP2930101B2 (ja) |
| KR (1) | KR100307882B1 (ja) |
| AU (1) | AU727476B2 (ja) |
| CA (1) | CA2228121C (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100334911B1 (ko) * | 2000-01-28 | 2002-05-04 | 오길록 | 적응형 전송시스템에서의 수신 신호 레벨 예측방법 |
| US7366309B2 (en) | 2004-03-08 | 2008-04-29 | Nec Corporation | Robot |
| US7720233B2 (en) | 2003-09-02 | 2010-05-18 | Nec Corporation | Signal processing method and apparatus |
Families Citing this family (59)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU739217C (en) * | 1997-06-11 | 2002-06-06 | Nec Corporation | Adaptive filter, step size control method thereof, and record medium therefor |
| JPH11345000A (ja) * | 1998-06-03 | 1999-12-14 | Nec Corp | 雑音消去方法及び雑音消去装置 |
| US6707912B2 (en) * | 1999-03-11 | 2004-03-16 | Motorola, Inc. | Method and apparatus for setting a step size for an adaptive filter coefficient of an echo canceller |
| GB9910448D0 (en) * | 1999-05-07 | 1999-07-07 | Ensigma Ltd | Cancellation of non-stationary interfering signals for speech recognition |
| GB9922654D0 (en) | 1999-09-27 | 1999-11-24 | Jaber Marwan | Noise suppression system |
| SG97885A1 (en) * | 2000-05-05 | 2003-08-20 | Univ Nanyang | Noise canceler system with adaptive cross-talk filters |
| US6963649B2 (en) * | 2000-10-24 | 2005-11-08 | Adaptive Technologies, Inc. | Noise cancelling microphone |
| FI110564B (fi) * | 2001-03-29 | 2003-02-14 | Nokia Corp | Järjestelmä automaattisen kohinanvaimennuksen (ANC) kytkemiseksi päälle ja poiskytkemiseksi matkapuhelimessa |
| US7245726B2 (en) * | 2001-10-03 | 2007-07-17 | Adaptive Technologies, Inc. | Noise canceling microphone system and method for designing the same |
| US20060116874A1 (en) * | 2003-10-24 | 2006-06-01 | Jonas Samuelsson | Noise-dependent postfiltering |
| JP2005249816A (ja) * | 2004-03-01 | 2005-09-15 | Internatl Business Mach Corp <Ibm> | 信号強調装置、方法及びプログラム、並びに音声認識装置、方法及びプログラム |
| KR100698630B1 (ko) * | 2004-06-28 | 2007-03-21 | 삼성전자주식회사 | 스텝사이즈 조정기능을 구비한 등화기 및 등화방법 |
| US20060035593A1 (en) * | 2004-08-12 | 2006-02-16 | Motorola, Inc. | Noise and interference reduction in digitized signals |
| US7602844B2 (en) * | 2005-12-05 | 2009-10-13 | National Taiwan University | On-line step-size calculation using signal power estimation and tone grouping of the frequency-domain equalizer for DMT-based transceiver |
| US8345890B2 (en) | 2006-01-05 | 2013-01-01 | Audience, Inc. | System and method for utilizing inter-microphone level differences for speech enhancement |
| CN1809105B (zh) * | 2006-01-13 | 2010-05-12 | 北京中星微电子有限公司 | 适用于小型移动通信设备的双麦克语音增强方法及系统 |
| US8744844B2 (en) * | 2007-07-06 | 2014-06-03 | Audience, Inc. | System and method for adaptive intelligent noise suppression |
| US8204252B1 (en) | 2006-10-10 | 2012-06-19 | Audience, Inc. | System and method for providing close microphone adaptive array processing |
| US9185487B2 (en) * | 2006-01-30 | 2015-11-10 | Audience, Inc. | System and method for providing noise suppression utilizing null processing noise subtraction |
| US8194880B2 (en) | 2006-01-30 | 2012-06-05 | Audience, Inc. | System and method for utilizing omni-directional microphones for speech enhancement |
| US8934641B2 (en) | 2006-05-25 | 2015-01-13 | Audience, Inc. | Systems and methods for reconstructing decomposed audio signals |
| US8849231B1 (en) | 2007-08-08 | 2014-09-30 | Audience, Inc. | System and method for adaptive power control |
| US8204253B1 (en) | 2008-06-30 | 2012-06-19 | Audience, Inc. | Self calibration of audio device |
| US8150065B2 (en) * | 2006-05-25 | 2012-04-03 | Audience, Inc. | System and method for processing an audio signal |
| US8949120B1 (en) | 2006-05-25 | 2015-02-03 | Audience, Inc. | Adaptive noise cancelation |
| EP1914726A1 (en) * | 2006-10-16 | 2008-04-23 | SiTel Semiconductor B.V. | base station of a telephone system and telephone system comprising such as a base station |
| US8259926B1 (en) | 2007-02-23 | 2012-09-04 | Audience, Inc. | System and method for 2-channel and 3-channel acoustic echo cancellation |
| US8189766B1 (en) | 2007-07-26 | 2012-05-29 | Audience, Inc. | System and method for blind subband acoustic echo cancellation postfiltering |
| US8326617B2 (en) | 2007-10-24 | 2012-12-04 | Qnx Software Systems Limited | Speech enhancement with minimum gating |
| US8015002B2 (en) | 2007-10-24 | 2011-09-06 | Qnx Software Systems Co. | Dynamic noise reduction using linear model fitting |
| US8606566B2 (en) * | 2007-10-24 | 2013-12-10 | Qnx Software Systems Limited | Speech enhancement through partial speech reconstruction |
| EP2058803B1 (en) | 2007-10-29 | 2010-01-20 | Harman/Becker Automotive Systems GmbH | Partial speech reconstruction |
| US8180064B1 (en) | 2007-12-21 | 2012-05-15 | Audience, Inc. | System and method for providing voice equalization |
| US8143620B1 (en) | 2007-12-21 | 2012-03-27 | Audience, Inc. | System and method for adaptive classification of audio sources |
| US8194882B2 (en) | 2008-02-29 | 2012-06-05 | Audience, Inc. | System and method for providing single microphone noise suppression fallback |
| US8355511B2 (en) | 2008-03-18 | 2013-01-15 | Audience, Inc. | System and method for envelope-based acoustic echo cancellation |
| US8521530B1 (en) | 2008-06-30 | 2013-08-27 | Audience, Inc. | System and method for enhancing a monaural audio signal |
| US8774423B1 (en) * | 2008-06-30 | 2014-07-08 | Audience, Inc. | System and method for controlling adaptivity of signal modification using a phantom coefficient |
| US9838784B2 (en) | 2009-12-02 | 2017-12-05 | Knowles Electronics, Llc | Directional audio capture |
| US9008329B1 (en) | 2010-01-26 | 2015-04-14 | Audience, Inc. | Noise reduction using multi-feature cluster tracker |
| US8700394B2 (en) * | 2010-03-24 | 2014-04-15 | Microsoft Corporation | Acoustic model adaptation using splines |
| CN102859591B (zh) * | 2010-04-12 | 2015-02-18 | 瑞典爱立信有限公司 | 用于语音编码器中的噪声消除的方法和装置 |
| US9558755B1 (en) | 2010-05-20 | 2017-01-31 | Knowles Electronics, Llc | Noise suppression assisted automatic speech recognition |
| CN102456341A (zh) * | 2010-10-28 | 2012-05-16 | 美商富迪科技股份有限公司 | 麦克风装置与声音处理方法 |
| US8838445B1 (en) * | 2011-10-10 | 2014-09-16 | The Boeing Company | Method of removing contamination in acoustic noise measurements |
| US8615394B1 (en) * | 2012-01-27 | 2013-12-24 | Audience, Inc. | Restoration of noise-reduced speech |
| WO2013187932A1 (en) | 2012-06-10 | 2013-12-19 | Nuance Communications, Inc. | Noise dependent signal processing for in-car communication systems with multiple acoustic zones |
| US9805738B2 (en) | 2012-09-04 | 2017-10-31 | Nuance Communications, Inc. | Formant dependent speech signal enhancement |
| US9640194B1 (en) | 2012-10-04 | 2017-05-02 | Knowles Electronics, Llc | Noise suppression for speech processing based on machine-learning mask estimation |
| WO2014070139A2 (en) | 2012-10-30 | 2014-05-08 | Nuance Communications, Inc. | Speech enhancement |
| KR20140111480A (ko) * | 2013-03-11 | 2014-09-19 | 삼성전자주식회사 | 보코더 잡음 억제 방법 및 장치 |
| US9484044B1 (en) | 2013-07-17 | 2016-11-01 | Knuedge Incorporated | Voice enhancement and/or speech features extraction on noisy audio signals using successively refined transforms |
| US9530434B1 (en) | 2013-07-18 | 2016-12-27 | Knuedge Incorporated | Reducing octave errors during pitch determination for noisy audio signals |
| US9536540B2 (en) | 2013-07-19 | 2017-01-03 | Knowles Electronics, Llc | Speech signal separation and synthesis based on auditory scene analysis and speech modeling |
| US9208794B1 (en) * | 2013-08-07 | 2015-12-08 | The Intellisis Corporation | Providing sound models of an input signal using continuous and/or linear fitting |
| US9510096B2 (en) * | 2014-05-04 | 2016-11-29 | Yang Gao | Noise energy controlling in noise reduction system with two microphones |
| DE112015003945T5 (de) | 2014-08-28 | 2017-05-11 | Knowles Electronics, Llc | Mehrquellen-Rauschunterdrückung |
| US9978388B2 (en) | 2014-09-12 | 2018-05-22 | Knowles Electronics, Llc | Systems and methods for restoration of speech components |
| US9820042B1 (en) | 2016-05-02 | 2017-11-14 | Knowles Electronics, Llc | Stereo separation and directional suppression with omni-directional microphones |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4649505A (en) * | 1984-07-02 | 1987-03-10 | General Electric Company | Two-input crosstalk-resistant adaptive noise canceller |
| JP3089082B2 (ja) * | 1991-07-10 | 2000-09-18 | シャープ株式会社 | 適応型ディジタルフィルタ |
| JP3008763B2 (ja) * | 1993-12-28 | 2000-02-14 | 日本電気株式会社 | 適応フィルタによるシステム同定の方法および装置 |
| US5644596A (en) * | 1994-02-01 | 1997-07-01 | Qualcomm Incorporated | Method and apparatus for frequency selective adaptive filtering |
| JPH0856180A (ja) * | 1994-08-10 | 1996-02-27 | Toshiba Corp | ノイズキャンセラ |
| JP2760373B2 (ja) * | 1995-03-03 | 1998-05-28 | 日本電気株式会社 | 雑音消去装置 |
| JP2685031B2 (ja) * | 1995-06-30 | 1997-12-03 | 日本電気株式会社 | 雑音消去方法及び雑音消去装置 |
-
1997
- 1997-01-29 JP JP9014410A patent/JP2930101B2/ja not_active Expired - Lifetime
-
1998
- 1998-01-28 EP EP98101469A patent/EP0856834A3/en not_active Ceased
- 1998-01-28 CA CA002228121A patent/CA2228121C/en not_active Expired - Fee Related
- 1998-01-28 AU AU52787/98A patent/AU727476B2/en not_active Ceased
- 1998-01-29 US US09/015,525 patent/US5978824A/en not_active Expired - Fee Related
- 1998-01-30 KR KR1019980002544A patent/KR100307882B1/ko not_active IP Right Cessation
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100334911B1 (ko) * | 2000-01-28 | 2002-05-04 | 오길록 | 적응형 전송시스템에서의 수신 신호 레벨 예측방법 |
| US7720233B2 (en) | 2003-09-02 | 2010-05-18 | Nec Corporation | Signal processing method and apparatus |
| US9543926B2 (en) | 2003-09-02 | 2017-01-10 | Nec Corporation | Signal processing method and device |
| US7366309B2 (en) | 2004-03-08 | 2008-04-29 | Nec Corporation | Robot |
Also Published As
| Publication number | Publication date |
|---|---|
| KR100307882B1 (ko) | 2001-10-19 |
| CA2228121C (en) | 2001-05-01 |
| CA2228121A1 (en) | 1998-07-29 |
| JP2930101B2 (ja) | 1999-08-03 |
| AU727476B2 (en) | 2000-12-14 |
| US5978824A (en) | 1999-11-02 |
| AU5278798A (en) | 1998-08-06 |
| EP0856834A3 (en) | 1999-02-17 |
| EP0856834A2 (en) | 1998-08-05 |
| KR19980070927A (ko) | 1998-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2930101B2 (ja) | 雑音消去装置 | |
| JP2874679B2 (ja) | 雑音消去方法及びその装置 | |
| JP2882364B2 (ja) | 雑音消去方法及び雑音消去装置 | |
| JP2685031B2 (ja) | 雑音消去方法及び雑音消去装置 | |
| JP4632047B2 (ja) | 信号処理方法および装置 | |
| JP2760373B2 (ja) | 雑音消去装置 | |
| JP2643852B2 (ja) | エコーキャンセラ | |
| JP3566158B2 (ja) | エコーキャンセラ装置 | |
| JPH11345000A (ja) | 雑音消去方法及び雑音消去装置 | |
| JP2947093B2 (ja) | 適応フィルタによるシステム同定の方法および装置 | |
| JPWO2006049260A1 (ja) | 信号処理の方法、信号処理の装置および信号処理用プログラム | |
| JPH10304489A (ja) | エコー・雑音成分除去装置 | |
| JP5321914B2 (ja) | システム同定装置及びシステム同定方法 | |
| JP3293575B2 (ja) | 雑音消去方法及びそれを用いた雑音消去装置 | |
| US20050025273A1 (en) | Method and apparatus for echo compensation | |
| JP2002076999A (ja) | システム同定方法および装置 | |
| JP5058052B2 (ja) | エコー消去装置 | |
| JP2023102181A (ja) | 信号処理装置、信号処理方法、および信号処理プログラム | |
| JPH05134678A (ja) | 適応型雑音抑圧装置 | |
| JP4247158B2 (ja) | 多チャネル音響エコー消去方法、多チャネル音響エコー消去装置、多チャネル音響エコー消去プログラム、記録媒体 | |
| JP2023100220A (ja) | 信号処理装置、信号処理方法、および信号処理プログラム | |
| JP2023100564A (ja) | 信号処理装置、信号処理方法、および信号処理プログラム | |
| WO2019092798A1 (ja) | 信号処理装置、信号処理方法、および信号処理プログラム | |
| JP2004357053A (ja) | エコーキャンセラ装置およびエコーキャンセラ方法 | |
| JP2000278186A (ja) | 適応型ノイズキャンセラ |