JP2010211190A - バックグラウンドノイズ推定 - Google Patents
バックグラウンドノイズ推定 Download PDFInfo
- Publication number
- JP2010211190A JP2010211190A JP2010012611A JP2010012611A JP2010211190A JP 2010211190 A JP2010211190 A JP 2010211190A JP 2010012611 A JP2010012611 A JP 2010012611A JP 2010012611 A JP2010012611 A JP 2010012611A JP 2010211190 A JP2010211190 A JP 2010211190A
- Authority
- JP
- Japan
- Prior art keywords
- spectral density
- power spectral
- signal
- value
- background noise
- 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
- 230000003595 spectral effect Effects 0.000 claims abstract description 289
- 238000009499 grossing Methods 0.000 claims abstract description 109
- 238000004364 calculation method Methods 0.000 claims description 143
- 238000000034 method Methods 0.000 claims description 53
- 230000003044 adaptive effect Effects 0.000 claims description 29
- 230000001965 increasing effect Effects 0.000 claims description 23
- 230000002829 reductive effect Effects 0.000 claims description 19
- 230000009467 reduction Effects 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 15
- 230000007423 decrease Effects 0.000 claims description 10
- 230000001419 dependent effect Effects 0.000 claims description 10
- 230000001953 sensory effect Effects 0.000 claims description 9
- 230000000630 rising effect Effects 0.000 claims description 7
- 230000021317 sensory perception Effects 0.000 claims description 4
- 230000000670 limiting effect Effects 0.000 claims description 2
- 230000001174 ascending effect Effects 0.000 claims 1
- 230000000873 masking effect Effects 0.000 description 79
- 238000004422 calculation algorithm Methods 0.000 description 64
- 238000012360 testing method Methods 0.000 description 47
- 230000000694 effects Effects 0.000 description 27
- 230000005236 sound signal Effects 0.000 description 22
- 230000005284 excitation Effects 0.000 description 12
- 230000004044 response Effects 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 238000001228 spectrum Methods 0.000 description 9
- 230000008901 benefit Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 8
- 230000008447 perception Effects 0.000 description 8
- 230000002123 temporal effect Effects 0.000 description 8
- 238000001914 filtration Methods 0.000 description 7
- 210000003027 ear inner Anatomy 0.000 description 6
- 210000002469 basement membrane Anatomy 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 210000003477 cochlea Anatomy 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000282412 Homo Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000009849 deactivation Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000001976 improved effect Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 210000002751 lymph Anatomy 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 210000004379 membrane Anatomy 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 210000002985 organ of corti Anatomy 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 210000003582 temporal bone Anatomy 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques 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
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques 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
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Circuit For Audible Band Transducer (AREA)
- Noise Elimination (AREA)
Abstract
【解決手段】音響的なバックグラウンドノイズのパワースペクトル密度を推定するシステムであって、バックグラウンドノイズを表すノイズ信号を生成するセンサユニットと、パワースペクトル密度計算ユニットと、時間領域信号平滑化ユニットと、周波数領域信号平滑化ユニットと、増分計算ユニットと、減分計算ユニットと、推定信号平滑化ユニットとを備える。
【選択図】図3
Description
本発明は、バックグラウンドノイズを推定するシステムおよび方法に関し、特に、同時のスピーチアクティビティの間にバックグラウンドノイズを推定するシステムおよび方法に関する。
受信器の情報内容に寄与しない、そしてそのために擾乱とはみなされない音波は、一般的に、バックグラウンドノイズといわれる。バックグラウンドノイズの発展の過程は、典型的に3つのさまざまな段階に分類され得る。これらは、1つ以上の供給源によるノイズの放出、ノイズの移送、そしてノイズの受容である。第1にノイズ信号(例えば、バックグラウンドノイズ)をノイズ自体の供給源において抑制し、その後信号の移送を抑制することによる試行がなされることが明らかである。しかしながら、ノイズ信号の放出は、多くの場合において所望のレベルまで低減され得ない。なぜなら、例えば、時間および位置に対して自然発生的に発生する周囲ノイズは、不十分に制御され得るのみであるか、またはまったく制御されないからである。
音響的なバックグラウンドノイズのパワースペクトル密度を推定するシステムが提供され、該システムは、バックグラウンドノイズを表すノイズ信号を生成するセンサユニットと、パワースペクトル密度計算ユニットであって、該パワースペクトル密度計算ユニットは、連続的な計算サイクルを配置することによってノイズ信号から現在のパワースペクトル密度を連続的に決定するように適合され、対応するパワースペクトル密度出力信号を提供するように適合された、計算ユニットと、時間領域信号平滑化ユニットであって、該時間領域信号平滑化ユニットは、パワースペクトル密度出力信号を時間領域内で平滑化するように適合され、結果として生じる時間的に平滑化された信号を提供するように適合された、時間領域信号平滑化ユニットと、周波数領域信号平滑化ユニットであって、該周波数領域信号平滑化ユニットは、時間領域信号平滑化ユニットから受信された時間的に平滑化された信号を周波数領域内で平滑化するように適合され、結果として生じる平滑化パワースペクトル密度信号を提供する用に適合された、周波数領域信号平滑化ユニットと、増分計算ユニットであって、バックグラウンドノイズのパワースペクトル密度の推定値に依存した増分の計算のために適合された、増分計算ユニットと、減分計算ユニットであって、バックグラウンドノイズのパワースペクトル密度の推定値に依存した減分の計算のために適合された、減分計算ユニットと、推定信号平滑化ユニットであって、増分および減分から、バックグラウンドノイズのパワースペクトル密度の推定値を計算するように適合されている、推定信号平滑化ユニットとを備える。平滑化パワースペクトル密度信号のレベルが増加する場合において、最大増分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも大きくなる場合に達成するまで、最小増分値から開始して、所定量だけ該増分が増加する。平滑化パワースペクトル密度信号のレベルが減少する場合において、最大減分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも小さくなる場合に達成するまで、最小減分値から開始して、所定量だけ該減分が増加する。
(項目1)
音響的なバックグラウンドノイズのパワースペクトル密度を推定するシステムであって、該システムは、
該バックグラウンドノイズを表すノイズ信号を生成するセンサユニットと、
連続計算サイクルを配置することによって、該ノイズ信号から現在のパワースペクトル密度を連続的に決定するように適合され、対応するパワースペクトル密度出力信号を提供するように適合されたパワースペクトル密度計算ユニットと、
該パワースペクトル密度出力信号を時間領域で平滑化するように適合され、結果として生じる時間的に平滑化された信号を提供するように適合された時間領域信号平滑化ユニットと、
該時間領域信号平滑化ユニットから受信された該時間的に平滑化された信号を周波数領域で平滑化するように適合され、結果として生じる平滑化パワースペクトル密度信号を提供するように適合された周波数領域信号平滑化ユニットと、
該バックグラウンドノイズの該パワースペクトル密度の推定値に依存した増分の計算のために適合された増分計算ユニットと、
該バックグラウンドノイズの該パワースペクトル密度の該推定値に依存した減分の計算のために適合された減分計算ユニットと、
該増分および減分から、該バックグラウンドノイズの該パワースペクトル密度の該推定値を計算するように適合された推定信号平滑化ユニットと
を備え、
該平滑化パワースペクトル密度信号のレベルが増加する場合において、最大増分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも大きくなる場合に達成するまで、最小増分値から開始して、所定量だけ該増分が増加し、
該平滑化パワースペクトル密度信号のレベルが減少する場合において、最大減分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも小さくなる場合に達成するまで、最小減分値から開始して、所定量だけ該減分が増加する、システム。
(項目2)
誤差信号を提供する適応フィルタをさらに備え、上記パワースペクトル密度計算ユニットは、連続計算サイクルを配置する該適応フィルタの該誤差信号から現在のパワースペクトル密度を決定するように適合され、上記システムは、対応するパワースペクトル密度出力信号および対応する平滑化パワースペクトル密度信号を提供するように適合される、上記項目に記載のシステム。
(項目3)
上記システムは、
新しい計算サイクルにおいて決定された上記パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された上記バックグラウンドノイズの該パワースペクトル密度の推定値よりも小さい場合に、上記増分値の計算のモードから、上記減分値の計算のモードまで該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化することであって、該システムは該増分値の現在の値を最小増分値に再設定するように適合されている、ことと、
新しい計算サイクルにおいて決定された該パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された該バックグラウンドノイズの該パワースペクトル密度の推定値よりも大きい場合に、該減分値の計算のモードから、該増分値の計算のモードまで該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化することであって、該システムは該減分値の現在の値を最小減分値に再設定するように適合されている、ことと
を行うように適合されている、上記項目のいずれかに記載のシステム。
(項目4)
上記システムは、上記バックグラウンドノイズの上記パワースペクトル密度の上記推定値を減分する場合には、該推定値の低減を一定の規定値まで制限するように適合され、その結果、該バックグラウンドノイズの該パワースペクトル密度の該推定値は、上記現在計算された値にかかわらず最小値未満に減少する、上記項目のいずれかに記載のシステム。
(項目5)
上記時間領域信号平滑化ユニットは、2つの異なる時間定数を利用して時間に対する現在測定されたパワースペクトル密度の平滑化のために適合され、該2つの異なる時間定数のうちの1つは上昇信号の場合のためのものであり、該2つの異なる時間定数のうちの1つは減少信号の場合のためのものである、上記項目のいずれかに記載のシステム。
(項目6)
上記周波数領域信号平滑化ユニットは、周波数平滑化第3係数を用いて最小周波数から上方に開始する、そして/または周波数平滑化第4係数を用いて最大周波数から下方に開始する、上記時間領域信号ユニットからの上記時間的に平滑化された信号の平滑化のために適合されている、上記項目のいずれかに記載のシステム。
(項目7)
現在測定されたパワースペクトル密度の時間に対する平滑化のための第1および第2の係数は、ヒトの耳の心理音響学的感覚特性を表し、そして/または、
該現在測定されたパワースペクトル密度の周波数に対する平滑化のための第3および第4の係数は、該ヒトの耳の心理音響学的感覚特性を表す、上記項目のいずれかに記載のシステム。
(項目8)
上記増分値の増加の値は、現在測定されたパワースペクトル密度の上記平滑化パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて個別に選択され、上記減分値の増加の値は、現在測定されたパワースペクトル密度の該平滑化パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて選択される、上記項目のいずれかに記載のシステム。
(項目9)
上記システムは、心理音響学的感覚知覚に対応する周波数群内の平滑化パワースペクトル密度または非平滑化パワースペクトル密度のスペクトル成分を、さらなる処理の前に、各周波数群に対する単一の結合信号に統合するように適合されている、上記項目のいずれかに記載のシステム。
(項目10)
音響的なバックグラウンドノイズのパワースペクトル密度の推定方法であって、該方法は、
パワースペクトル密度計算ユニットによって、マイクロフォン信号からの現在のパワースペクトル密度を決定し、対応するパワースペクトル密度出力信号を提供するステップと、
該提供されたパワースペクトル密度出力信号を時間領域内で平滑化して、結果として生じる時間的に平滑化された信号を提供するステップと、
該時間的に平滑化された信号を周波数領域で平滑化して、結果として生じる平滑化パワースペクトル密度信号を提供するステップと、
該バックグラウンドノイズのパワースペクトル密度の推定値に依存して増分を計算するステップと、
該バックグラウンドノイズの該パワースペクトル密度の推定値に依存して減分を計算するステップと、
該増分および減分から、該バックグラウンドノイズの該パワースペクトル密度の該推定値を計算するステップと
を包含し、
該平滑化パワースペクトル密度信号のレベルが増加する場合において、最大増分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズの該パワースペクトル密度の該推定値よりも大きくなる場合に達成するまで、最小増分値から開始して、所定量だけ該増分が増加し、
該平滑化パワースペクトル密度信号のレベルが減少する場合において、最大減分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値よりも小さくなる場合に達成するまで、最小減分値から開始して、所定量だけ該減分が増加する、方法。
(項目11)
連続的な計算サイクルを配置することによって、適応フィルタから導出された誤差信号から現在のパワースペクトル密度を決定するステップと、
対応するパワースペクトル密度出力信号および対応する平滑化パワースペクトル密度信号を提供するステップと、
をさらに包含する、上記項目に記載の方法。
(項目12)
新しい計算サイクルにおいて決定された上記パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された上記バックグラウンドノイズの該パワースペクトル密度の推定値よりも小さい場合に、上記増分値の計算のモードから、上記減分値の計算のモードに、該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化するステップであって、該増分値の現在の値は最小増分値に再設定される、ステップと、
新しい計算サイクルにおいて決定された該パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された該バックグラウンドノイズの該パワースペクトル密度の推定値よりも大きい場合に、該減分値の計算のモードから、該増分値の計算のモードに、該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化するステップであって、該減分値の現在の値は最小減分値に再設定される、ステップと
をさらに包含する、上記項目のいずれかに記載の方法。
(項目13)
上記バックグラウンドノイズの上記パワースペクトル密度の上記推定値を減分する場合には、該推定値の低減を一定の規定値まで制限するステップをさらに包含し、その結果、該バックグラウンドノイズの該パワースペクトル密度の該推定値は、上記現在計算された値にかかわらず最小値未満に減少する、上記項目のいずれかに記載の方法。
(項目14)
2つの異なる時間定数を利用して時間に対する現在測定されたパワースペクトル密度を上記時間領域内で平滑化するステップをさらに包含し、該2つの異なる時間定数のうちの1つは上昇信号の場合のためのものであり、該2つの異なる時間定数のうちの1つは減少信号の場合のためのものである、上記項目のいずれかに記載の方法。
(項目15)
周波数平滑化第3係数を用いて最小周波数から上方に開始する、そして/または周波数平滑化第4係数を用いて最大周波数から下方に開始する、上記時間領域信号平滑化ユニットからの時間的に平滑化された信号を周波数領域内で平滑化するステップを包含する、上記項目のいずれかに記載の方法。
(項目16)
現在測定されたパワースペクトル密度の時間に対する平滑化のための第1および第2の係数は、ヒトの耳の心理音響学的感覚特性を表し、そして/または、
該現在測定されたパワースペクトル密度の周波数に対する平滑化のための第3および第4の係数は、該ヒトの耳の心理音響学的感覚特性を表す、上記項目のいずれかに記載の方法。
(項目17)
上記増分値の増加の値は、現在測定されたパワースペクトル密度の上記(平滑化)パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて個別に選択され、上記減分値の増加の値は、現在測定されたパワースペクトル密度の該(平滑化)パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて個別に選択される、上記項目のいずれかに記載の方法。
(項目18)
心理音響学的感覚知覚に対応する周波数群内の(平滑化)パワースペクトル密度のスペクトル成分を、さらなる処理の前に、各周波数群に対する単一の結合信号に統合する、上記項目のいずれかに記載の方法。
音響的なバックグラウンドノイズのパワースペクトル密度を推定するシステムおよび方法が提示され、ここでは、平滑化パワースペクトル密度信号のレベルが増加する場合において、最大増分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも大きくなる場合に達成するまで、最小増分値から開始して、所定量だけ該増分が増加する。平滑化パワースペクトル密度信号のレベルが減少する場合において、最大減分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも小さくなる場合に達成するまで、最小減分値から開始して、所定量だけ該減分が増加する。
以下に開示される例において、バックグラウンドノイズのパワースペクトル密度は、マイクロフォン信号から直接的に、または適応フィルタの誤差信号から推定される。適応方法およびシステムは、アルゴリズムが、周囲状況の変化(例えば、時間に対するレベルおよびスペクトル成分の変化を受けるノイズ信号に対する変化)への、それらのフィルタ係数の一定の修正に対して自動的に適合されることである。この能力は、パラメータを連続的に最適化するシステム構造によって提供される。このようなシステムにおいて、入力センサ(例えば、マイクロフォン)は、1つ以上のノイズ供給源によって生成された望ましくないノイズ(例えば、バックグラウンドノイズ)を表す信号を得るために用いられる。次いで、信号は、適応フィルタの入力にルーティングされ、フィルタによって処理されて出力信号を生成し、この出力信号は、望ましくないノイズ信号が与えられると有用な信号(例えば、音声信号)から差し引かれ、適応フィルタの入力信号と、望ましくないノイズ信号との相関が、有用な信号と共に発生する。この差し引きから得られる出力信号は、また、適応フィルタに対する誤差信号ともいわれる。望ましくないノイズを表す入力センサの信号と共に、誤差信号は、観測されるエコーのレベル全体を適応して最小化させるために、パラメータの修正および適応フィルタの特性の基本を形成する。
MSE=E{e2[n]}
である。
k=0,...N−1に対して、bk[n+1]=bk[n]=2・μ・e[n]・x[n−k]
のように説明される。
0<μ<μmax=1/[(N+1)・E{x2[n]}]
ここで、Nは、FIRフィルタの次数を表し、E{x2[n]}は、基準信号x[n]の信号出力を表す。実際には、使用されるステップサイズおよび収束パラメータμは、それぞれ、しばしばμ=μmax/10であるように選ばれる。LMSフィルタの最小二乗平均アルゴリズムは、こうして、以下で概略されるように実現され得る。
1 制御変数をn=0に設定することと、アルゴリズムの実行の始まりにおいてk=0,…,N−1に対する開始係数bk[n=0]を選択することと(例えば、k=0...N−1においてbk[0]=0およびe[0]=d[0])、増幅因子μ<μmax(例えば、μ=μmax/10)を選択することとによるアルゴリズムの初期化。
2 基準信号x[n]および信号d[n]の格納。
3
4 誤差の決定e[n]=d[n]−y[n]。
5 k=0,...,Nに対して、bk[n+1]=bk[n]+2・μ・e[n]・x[n−k]
に従う係数の更新。
6 次の繰り返しステップn=n+1の実行およびステップ2〜6を繰り返す。
Inc(ω)=IncMin+100*ΔInc
として計算されるからである。
Dec(ω)=DecMin+100*ΔDec
として計算されるからである。
7 時間領域信号平滑化ユニット
8 周波数領域信号平滑化ユニット
9 増分計算ユニット
10 減分計算ユニット
11 推定信号平滑化ユニット
Claims (18)
- 音響的なバックグラウンドノイズのパワースペクトル密度を推定するシステムであって、該システムは、
該バックグラウンドノイズを表すノイズ信号を生成するセンサユニットと、
連続計算サイクルを配置することによって、該ノイズ信号から現在のパワースペクトル密度を連続的に決定するように適合され、対応するパワースペクトル密度出力信号を提供するように適合されたパワースペクトル密度計算ユニットと、
該パワースペクトル密度出力信号を時間領域で平滑化するように適合され、結果として生じる時間的に平滑化された信号を提供するように適合された時間領域信号平滑化ユニットと、
該時間領域信号平滑化ユニットから受信された該時間的に平滑化された信号を周波数領域で平滑化するように適合され、結果として生じる平滑化パワースペクトル密度信号を提供するように適合された周波数領域信号平滑化ユニットと、
該バックグラウンドノイズの該パワースペクトル密度の推定値に依存した増分の計算のために適合された増分計算ユニットと、
該バックグラウンドノイズの該パワースペクトル密度の該推定値に依存した減分の計算のために適合された減分計算ユニットと、
該増分および減分から、該バックグラウンドノイズの該パワースペクトル密度の該推定値を計算するように適合された推定信号平滑化ユニットと
を備え、
該平滑化パワースペクトル密度信号のレベルが増加する場合において、最大増分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも大きくなる場合に達成するまで、最小増分値から開始して、所定量だけ該増分が増加し、
該平滑化パワースペクトル密度信号のレベルが減少する場合において、最大減分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値該増分値よりも小さくなる場合に達成するまで、最小減分値から開始して、所定量だけ該減分が増加する、システム。 - 誤差信号を提供する適応フィルタをさらに備え、前記パワースペクトル密度計算ユニットは、連続計算サイクルを配置する該適応フィルタの該誤差信号から現在のパワースペクトル密度を決定するように適合され、前記システムは、対応するパワースペクトル密度出力信号および対応する平滑化パワースペクトル密度信号を提供するように適合される、請求項1に記載のシステム。
- 前記システムは、
新しい計算サイクルにおいて決定された前記パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された前記バックグラウンドノイズの該パワースペクトル密度の推定値よりも小さい場合に、前記増分値の計算のモードから、前記減分値の計算のモードまで該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化することであって、該システムは該増分値の現在の値を最小増分値に再設定するように適合されている、ことと、
新しい計算サイクルにおいて決定された該パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された該バックグラウンドノイズの該パワースペクトル密度の推定値よりも大きい場合に、該減分値の計算のモードから、該増分値の計算のモードまで該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化することであって、該システムは該減分値の現在の値を最小減分値に再設定するように適合されている、ことと
を行うように適合されている、請求項1または請求項2に記載のシステム。 - 前記システムは、前記バックグラウンドノイズの前記パワースペクトル密度の前記推定値を減分する場合には、該推定値の低減を一定の規定値まで制限するように適合され、その結果、該バックグラウンドノイズの該パワースペクトル密度の該推定値は、前記現在計算された値にかかわらず最小値未満に減少する、請求項1〜請求項3のいずれか1項に記載のシステム。
- 前記時間領域信号平滑化ユニットは、2つの異なる時間定数を利用して時間に対する現在測定されたパワースペクトル密度の平滑化のために適合され、該2つの異なる時間定数のうちの1つは上昇信号の場合のためのものであり、該2つの異なる時間定数のうちの1つは減少信号の場合のためのものである、請求項1〜請求項4のいずれか1項に記載のシステム。
- 前記周波数領域信号平滑化ユニットは、周波数平滑化第3係数を用いて最小周波数から上方に開始する、そして/または周波数平滑化第4係数を用いて最大周波数から下方に開始する、前記時間領域信号ユニットからの前記時間的に平滑化された信号の平滑化のために適合されている、請求項1〜請求項5のいずれか1項に記載のシステム。
- 現在測定されたパワースペクトル密度の時間に対する平滑化のための第1および第2の係数は、ヒトの耳の心理音響学的感覚特性を表し、そして/または、
該現在測定されたパワースペクトル密度の周波数に対する平滑化のための第3および第4の係数は、該ヒトの耳の心理音響学的感覚特性を表す、請求項1〜請求項6のいずれか1項に記載のシステム。 - 前記増分値の増加の値は、現在測定されたパワースペクトル密度の前記平滑化パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて個別に選択され、前記減分値の増加の値は、現在測定されたパワースペクトル密度の該平滑化パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて選択される、請求項1〜請求項7のいずれか1項に記載のシステム。
- 前記システムは、心理音響学的感覚知覚に対応する周波数群内の平滑化パワースペクトル密度または非平滑化パワースペクトル密度のスペクトル成分を、さらなる処理の前に、各周波数群に対する単一の結合信号に統合するように適合されている、請求項1〜請求項8のいずれか1項に記載のシステム。
- 音響的なバックグラウンドノイズのパワースペクトル密度の推定方法であって、該方法は、
パワースペクトル密度計算ユニットによって、マイクロフォン信号からの現在のパワースペクトル密度を決定し、対応するパワースペクトル密度出力信号を提供するステップと、
該提供されたパワースペクトル密度出力信号を時間領域内で平滑化して、結果として生じる時間的に平滑化された信号を提供するステップと、
該時間的に平滑化された信号を周波数領域で平滑化して、結果として生じる平滑化パワースペクトル密度信号を提供するステップと、
該バックグラウンドノイズのパワースペクトル密度の推定値に依存して増分を計算するステップと、
該バックグラウンドノイズの該パワースペクトル密度の推定値に依存して減分を計算するステップと、
該増分および減分から、該バックグラウンドノイズの該パワースペクトル密度の該推定値を計算するステップと
を包含し、
該平滑化パワースペクトル密度信号のレベルが増加する場合において、最大増分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズの該パワースペクトル密度の該推定値よりも大きくなる場合に達成するまで、最小増分値から開始して、所定量だけ該増分が増加し、
該平滑化パワースペクトル密度信号のレベルが減少する場合において、最大減分値が、同時に、新しい計算サイクルにおいて現在決定された該パワースペクトル密度の値が、事前の計算サイクルにおいて決定された該バックグラウンドノイズのパワースペクトル密度の該推定値よりも小さくなる場合に達成するまで、最小減分値から開始して、所定量だけ該減分が増加する、方法。 - 連続的な計算サイクルを配置することによって、適応フィルタから導出された誤差信号から現在のパワースペクトル密度を決定するステップと、
対応するパワースペクトル密度出力信号および対応する平滑化パワースペクトル密度信号を提供するステップと、
をさらに包含する、請求項10に記載の方法。 - 新しい計算サイクルにおいて決定された前記パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された前記バックグラウンドノイズの該パワースペクトル密度の推定値よりも小さい場合に、前記増分値の計算のモードから、前記減分値の計算のモードに、該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化するステップであって、該増分値の現在の値は最小増分値に再設定される、ステップと、
新しい計算サイクルにおいて決定された該パワースペクトル密度の現在の値が、事前の計算サイクルにおいて計算された該バックグラウンドノイズの該パワースペクトル密度の推定値よりも大きい場合に、該減分値の計算のモードから、該増分値の計算のモードに、該バックグラウンドノイズの該パワースペクトル密度を推定する計算を変化するステップであって、該減分値の現在の値は最小減分値に再設定される、ステップと
をさらに包含する、請求項10または請求項11に記載の方法。 - 前記バックグラウンドノイズの前記パワースペクトル密度の前記推定値を減分する場合には、該推定値の低減を一定の規定値まで制限するステップをさらに包含し、その結果、該バックグラウンドノイズの該パワースペクトル密度の該推定値は、前記現在計算された値にかかわらず最小値未満に減少する、請求項10〜請求項12のいずれか1項に記載の方法。
- 2つの異なる時間定数を利用して時間に対する現在測定されたパワースペクトル密度を前記時間領域内で平滑化するステップをさらに包含し、該2つの異なる時間定数のうちの1つは上昇信号の場合のためのものであり、該2つの異なる時間定数のうちの1つは減少信号の場合のためのものである、請求項10〜請求項13のいずれか1項に記載の方法。
- 周波数平滑化第3係数を用いて最小周波数から上方に開始する、そして/または周波数平滑化第4係数を用いて最大周波数から下方に開始する、前記時間領域信号平滑化ユニットからの時間的に平滑化された信号を周波数領域内で平滑化するステップを包含する、請求項10〜請求項14のいずれか1項に記載の方法。
- 現在測定されたパワースペクトル密度の時間に対する平滑化のための第1および第2の係数は、ヒトの耳の心理音響学的感覚特性を表し、そして/または、
該現在測定されたパワースペクトル密度の周波数に対する平滑化のための第3および第4の係数は、該ヒトの耳の心理音響学的感覚特性を表す、請求項10〜請求項15のいずれか1項に記載の方法。 - 前記増分値の増加の値は、現在測定されたパワースペクトル密度の前記(平滑化)パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて個別に選択され、前記減分値の増加の値は、現在測定されたパワースペクトル密度の該(平滑化)パワースペクトル密度信号内の各スペクトル位置に対して異なる値を用いて個別に選択される、請求項10〜請求項16のいずれか1項に記載の方法。
- 心理音響学的感覚知覚に対応する周波数群内の(平滑化)パワースペクトル密度のスペクトル成分を、さらなる処理の前に、各周波数群に対する単一の結合信号に統合する、請求項10〜請求項17のいずれか1項に記載の方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP09154541.8 | 2009-03-06 | ||
EP09154541.8A EP2226794B1 (en) | 2009-03-06 | 2009-03-06 | Background noise estimation |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010211190A true JP2010211190A (ja) | 2010-09-24 |
JP5439200B2 JP5439200B2 (ja) | 2014-03-12 |
Family
ID=40910903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010012611A Active JP5439200B2 (ja) | 2009-03-06 | 2010-01-22 | バックグラウンドノイズ推定 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8422697B2 (ja) |
EP (1) | EP2226794B1 (ja) |
JP (1) | JP5439200B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015129760A1 (ja) * | 2014-02-28 | 2015-09-03 | 日本電信電話株式会社 | 信号処理装置、方法及びプログラム |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9589580B2 (en) | 2011-03-14 | 2017-03-07 | Cochlear Limited | Sound processing based on a confidence measure |
US9516371B2 (en) | 2012-04-24 | 2016-12-06 | Mobitv, Inc. | Closed captioning management system |
US9135920B2 (en) * | 2012-11-26 | 2015-09-15 | Harman International Industries, Incorporated | System for perceived enhancement and restoration of compressed audio signals |
DE102013111784B4 (de) * | 2013-10-25 | 2019-11-14 | Intel IP Corporation | Audioverarbeitungsvorrichtungen und audioverarbeitungsverfahren |
EP2980798A1 (en) | 2014-07-28 | 2016-02-03 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Harmonicity-dependent controlling of a harmonic filter tool |
AT516046B1 (de) * | 2014-12-30 | 2016-02-15 | Audio Lab Swiss Ag | Verfahren und vorrichtung zur bestimmung der qualität eines übertragungssystems |
EP3252769B8 (en) * | 2016-06-03 | 2020-04-01 | Sony Corporation | Adding background sound to speech-containing audio data |
SE541331C2 (en) * | 2017-11-30 | 2019-07-09 | Creo Dynamics Ab | Active noise control method and system |
SE1850077A1 (en) | 2018-01-24 | 2019-07-25 | Creo Dynamics Ab | Active noise control method and system using variable actuator and sensor participation |
CN108922556B (zh) * | 2018-07-16 | 2019-08-27 | 百度在线网络技术(北京)有限公司 | 声音处理方法、装置及设备 |
CN109643554B (zh) * | 2018-11-28 | 2023-07-21 | 深圳市汇顶科技股份有限公司 | 自适应语音增强方法和电子设备 |
CN109727605B (zh) * | 2018-12-29 | 2020-06-12 | 苏州思必驰信息科技有限公司 | 处理声音信号的方法及*** |
CN110031083A (zh) * | 2018-12-31 | 2019-07-19 | 瑞声科技(新加坡)有限公司 | 一种噪音总声压级测量方法、***及计算机可读存储介质 |
CN110021305B (zh) * | 2019-01-16 | 2021-08-20 | 上海惠芽信息技术有限公司 | 一种音频滤波方法、音频滤波装置及可穿戴设备 |
US11011182B2 (en) * | 2019-03-25 | 2021-05-18 | Nxp B.V. | Audio processing system for speech enhancement |
CN112133322A (zh) * | 2020-10-19 | 2020-12-25 | 南通赛洋电子有限公司 | 一种基于噪声分类优化imcra算法的语音增强方法 |
CN113808607B (zh) * | 2021-03-05 | 2024-07-16 | 北京沃东天骏信息技术有限公司 | 基于神经网络的语音增强方法、装置及电子设备 |
CN113008364B (zh) * | 2021-04-05 | 2022-07-15 | 中国人民解放军63966部队 | 装甲车辆隐蔽性噪声测量方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06506330A (ja) * | 1991-03-05 | 1994-07-14 | ピクチュアテル コーポレイション | 可変ビットレート音声エンコーダ |
JPH08506434A (ja) * | 1993-11-30 | 1996-07-09 | エイ・ティ・アンド・ティ・コーポレーション | 通信システムにおける伝送ノイズ低減 |
JP2008116952A (ja) * | 2006-10-31 | 2008-05-22 | Harman Becker Automotive Systems Gmbh | 音声信号のモデルベース強化 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6263307B1 (en) * | 1995-04-19 | 2001-07-17 | Texas Instruments Incorporated | Adaptive weiner filtering using line spectral frequencies |
US7177805B1 (en) * | 1999-02-01 | 2007-02-13 | Texas Instruments Incorporated | Simplified noise suppression circuit |
US7454332B2 (en) * | 2004-06-15 | 2008-11-18 | Microsoft Corporation | Gain constrained noise suppression |
-
2009
- 2009-03-06 EP EP09154541.8A patent/EP2226794B1/en active Active
-
2010
- 2010-01-22 JP JP2010012611A patent/JP5439200B2/ja active Active
- 2010-03-05 US US12/718,473 patent/US8422697B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06506330A (ja) * | 1991-03-05 | 1994-07-14 | ピクチュアテル コーポレイション | 可変ビットレート音声エンコーダ |
JPH08506434A (ja) * | 1993-11-30 | 1996-07-09 | エイ・ティ・アンド・ティ・コーポレーション | 通信システムにおける伝送ノイズ低減 |
JP2008116952A (ja) * | 2006-10-31 | 2008-05-22 | Harman Becker Automotive Systems Gmbh | 音声信号のモデルベース強化 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015129760A1 (ja) * | 2014-02-28 | 2015-09-03 | 日本電信電話株式会社 | 信号処理装置、方法及びプログラム |
JPWO2015129760A1 (ja) * | 2014-02-28 | 2017-03-30 | 日本電信電話株式会社 | 信号処理装置、方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP5439200B2 (ja) | 2014-03-12 |
EP2226794B1 (en) | 2017-11-08 |
US20100226501A1 (en) | 2010-09-09 |
EP2226794A1 (en) | 2010-09-08 |
US8422697B2 (en) | 2013-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5439200B2 (ja) | バックグラウンドノイズ推定 | |
JP5184896B2 (ja) | アクティブノイズ制御システム | |
KR100860805B1 (ko) | 음성 강화 시스템 | |
EP2056296B1 (en) | Dynamic noise reduction | |
EP1580882B1 (en) | Audio enhancement system and method | |
US8744844B2 (en) | System and method for adaptive intelligent noise suppression | |
US8170221B2 (en) | Audio enhancement system and method | |
EP2283484B1 (en) | System and method for dynamic sound delivery | |
US10262673B2 (en) | Soft-talk audio capture for mobile devices | |
KR20070068379A (ko) | 오디오 신호의 인식된 라우드니스 및/또는 인식된 스펙트럼밸런스 계산 및 조정 | |
JP2007011330A (ja) | スピーチ信号の適合する強化のためのシステム | |
EP2234105A1 (en) | Background noise estimation | |
JP2020028016A (ja) | 残響抑制装置及び補聴器 | |
JP4522509B2 (ja) | オーディオ装置 | |
EP3531719B1 (en) | Dereverberation device and hearing aid | |
KR102403996B1 (ko) | 보청기의 채널영역 방식, 채널영역 방식을 이용한 보청기의 피팅방법, 그리고 이를 통해 피팅된 디지털 보청기 | |
EP4258263A1 (en) | Apparatus and method for noise suppression | |
CN117425812A (zh) | 一种音频信号处理方法、装置、存储介质和车辆 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130110 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130614 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130619 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130731 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131008 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131111 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20131205 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20131216 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5439200 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |