JPH10161694A - 帯域分割型雑音低減方法 - Google Patents
帯域分割型雑音低減方法Info
- Publication number
- JPH10161694A JPH10161694A JP8317786A JP31778696A JPH10161694A JP H10161694 A JPH10161694 A JP H10161694A JP 8317786 A JP8317786 A JP 8317786A JP 31778696 A JP31778696 A JP 31778696A JP H10161694 A JPH10161694 A JP H10161694A
- Authority
- JP
- Japan
- Prior art keywords
- noise
- signal
- band
- input signal
- bands
- 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
Landscapes
- Noise Elimination (AREA)
Abstract
(57)【要約】
【課題】 帯域分割数を少なくしても品質劣化を伴わな
い。 【解決手段】 入力信号を帯域分割し、その各帯域ごと
の信号Ｘk (n) についてパワーＰX,k (n) を求め（２
４）、これより雑音平均バワーＰavN,k (n) を推定し
（５１）、ＰX,k (n) ，ＰavN,k (n) によりＳＮＲ
k (n) を推定し（２７）、ＳＮＲk (n) によりゲインフ
ァクターＧ（ＳＮＲk (n) ）を求め（３０）、ＳＮＲk
(n) により音声と雑音との比率に応じた加算率αを、Ｓ
／Ｎが悪い時大となるように決定し、Ｘk (n) にＧ（Ｓ
ＮＲk (n) ）を掛けて、雑音低減した信号Ｐ k (n) ′を
得（２８）、αX,k (n) ＋（１−α）Ｙk (n) ′でαに
応じてＹk (n) ′にＸk (n) を加算したＹk (n) を得
（５３）、その全帯域を合成して時間領域に変換する。
い。 【解決手段】 入力信号を帯域分割し、その各帯域ごと
の信号Ｘk (n) についてパワーＰX,k (n) を求め（２
４）、これより雑音平均バワーＰavN,k (n) を推定し
（５１）、ＰX,k (n) ，ＰavN,k (n) によりＳＮＲ
k (n) を推定し（２７）、ＳＮＲk (n) によりゲインフ
ァクターＧ（ＳＮＲk (n) ）を求め（３０）、ＳＮＲk
(n) により音声と雑音との比率に応じた加算率αを、Ｓ
／Ｎが悪い時大となるように決定し、Ｘk (n) にＧ（Ｓ
ＮＲk (n) ）を掛けて、雑音低減した信号Ｐ k (n) ′を
得（２８）、αX,k (n) ＋（１−α）Ｙk (n) ′でαに
応じてＹk (n) ′にＸk (n) を加算したＹk (n) を得
（５３）、その全帯域を合成して時間領域に変換する。
Description
【０００１】
【発明の属する技術分野】この発明は、音声会議装置・
ＴＶ会議装置等の音声／音響装置等において、目的とな
る音声信号と不要な雑音等の信号が混在する入力信号か
ら、雑音を低減した音声信号を出力する雑音低減方法に
関する。
ＴＶ会議装置等の音声／音響装置等において、目的とな
る音声信号と不要な雑音等の信号が混在する入力信号か
ら、雑音を低減した音声信号を出力する雑音低減方法に
関する。
【０００２】
【従来の技術】音声会議・ＴＶ会議等の拡声通話系で
は、マイクロホンで受音し、相手側に送出される送話信
号に目的となる音声以外の周囲雑音等が混入すると、音
声の明瞭性が損なわれ通話品質が著しく劣化する。この
為、送話信号に含まれる目的音声以外の周囲雑音を低減
する事が強く求められている。
は、マイクロホンで受音し、相手側に送出される送話信
号に目的となる音声以外の周囲雑音等が混入すると、音
声の明瞭性が損なわれ通話品質が著しく劣化する。この
為、送話信号に含まれる目的音声以外の周囲雑音を低減
する事が強く求められている。
【０００３】雑音低減方法とは、目的となる音声信号と
不要な周囲雑音等の信号が混在する入力信号から、雑音
を低減した信号を出力する技術である。図２Ａは収音シ
ステムを示すもので、これを用いて従来の雑音低減方法
を説明する。この明細書においては、信号の時間表現は
離散時間を表わす整数値ｎを用いて、例えばＸ（ｎ）と
表わす。今発声者１２が発声した目的とする音声信号１
３をＳ（ｎ）、空調などの不要な周囲雑音１４をＮ
（ｎ）、これら音声信号１３と雑音１４とがマイクロホ
ン１１で受音されて雑音低減装置１６へ入力される入力
信号１５をＸ（ｎ）、雑音低減装置１６の出力信号１７
をＹ（ｎ）とする。雑音低減装置１６への入力信号Ｘ
（ｎ）には、目的となる音声信号Ｓ（ｎ）以外に周囲雑
音Ｎ（ｎ）が混入している。即ち Ｘ（ｎ）＝Ｓ（ｎ）＋Ｎ（ｎ） （１） と表わされる。この時、入力信号Ｘ（ｎ）中の雑音Ｎ
（ｎ）を低減し、目的となる音声信号Ｓ（ｎ）に近い信
号を出力信号Ｙ（ｎ）として取り出す装置を雑音低減装
置と呼ぶ。
不要な周囲雑音等の信号が混在する入力信号から、雑音
を低減した信号を出力する技術である。図２Ａは収音シ
ステムを示すもので、これを用いて従来の雑音低減方法
を説明する。この明細書においては、信号の時間表現は
離散時間を表わす整数値ｎを用いて、例えばＸ（ｎ）と
表わす。今発声者１２が発声した目的とする音声信号１
３をＳ（ｎ）、空調などの不要な周囲雑音１４をＮ
（ｎ）、これら音声信号１３と雑音１４とがマイクロホ
ン１１で受音されて雑音低減装置１６へ入力される入力
信号１５をＸ（ｎ）、雑音低減装置１６の出力信号１７
をＹ（ｎ）とする。雑音低減装置１６への入力信号Ｘ
（ｎ）には、目的となる音声信号Ｓ（ｎ）以外に周囲雑
音Ｎ（ｎ）が混入している。即ち Ｘ（ｎ）＝Ｓ（ｎ）＋Ｎ（ｎ） （１） と表わされる。この時、入力信号Ｘ（ｎ）中の雑音Ｎ
（ｎ）を低減し、目的となる音声信号Ｓ（ｎ）に近い信
号を出力信号Ｙ（ｎ）として取り出す装置を雑音低減装
置と呼ぶ。
【０００４】図３は、Spectrum Subtraction（S.F.Bol
l, IEEE Trans. on ASSP, vol.27, no.2, pp.113-120,
Apr(1979)). Wiener Filer(J.S.Lim. & A.V.Oppenheim,
in Proc. IEEE, vol.67, no.12, pp.1586-1604, Dec(1
979)). Maximum Likelihood Envelop（R.J.McAulay &
M.L.Malpass, IEEE Trans. on ASSP, vol.28, no.2, p
p.137-145, Apr(1980)). minimum mean squared error
method(MMSE)(Y.Ephraim& D.Malah, IEEE Trans. on AS
SP, vol.32, no.6, pp.1109-1121, Dec(1984)).等の短
時間スペクトラル振幅（ＳＴＳＡ）評価（Short Time S
pectral Amplitude(STSA) Estimation）を基礎とした雑
音低減方式で、従来使用されている方法の機能構成を示
すものである。これを用いて従来の雑音低減方法を説明
する。図２Ａと同一の要素には共通の記号を用いた。
l, IEEE Trans. on ASSP, vol.27, no.2, pp.113-120,
Apr(1979)). Wiener Filer(J.S.Lim. & A.V.Oppenheim,
in Proc. IEEE, vol.67, no.12, pp.1586-1604, Dec(1
979)). Maximum Likelihood Envelop（R.J.McAulay &
M.L.Malpass, IEEE Trans. on ASSP, vol.28, no.2, p
p.137-145, Apr(1980)). minimum mean squared error
method(MMSE)(Y.Ephraim& D.Malah, IEEE Trans. on AS
SP, vol.32, no.6, pp.1109-1121, Dec(1984)).等の短
時間スペクトラル振幅（ＳＴＳＡ）評価（Short Time S
pectral Amplitude(STSA) Estimation）を基礎とした雑
音低減方式で、従来使用されている方法の機能構成を示
すものである。これを用いて従来の雑音低減方法を説明
する。図２Ａと同一の要素には共通の記号を用いた。
【０００５】まず、マイクロホン１１により受音され
た、目的信号と不要な雑音とが混入する入力信号１５を
Ａ／Ｄ変換部２１においてデジタル化し、周波数帯域分
割部２２と雑音判別部２３に転送する。周波数帯域分割
部２２では、転送された信号が複数の周波数帯域に分割
される。周波数帯域への分割は、例えば離散的フーリエ
変換等を用いて行う。ここで、帯域分割された信号は一
般に複素数であるが、分割方法によっては実数となる場
合もある。ここでは一般的に、複素数を仮定して議論す
るが、実数の場合も同じ議論が可能である。周波数帯域
に分割されたｋ番目の周波数帯域の信号を Ｘk (n) ＝Ｘk,r (n) ＋ｊＸk,i (n) （２） （Ｘk,r ，Ｘk,i はそれぞれＸk (n) の実数部分と虚数
部分）とすると、Ｘk (n) は、入力信号パワー計算部２
４、入力信号位相計算部２５、雑音パワー計算部２６に
転送される。入力信号パワー計算部２４では、各帯域ご
との入力信号のパワーレベル ＰX,k (n) ＝Ｘk,r (n)2＋Ｘk,i (n)2 （３） が、入力信号位相計算部２５では各帯域ごとの位相 Φk (n) ＝tan -1［Ｘk,i (n) ／Ｘk,r (n) ］ （４） それぞれが計算される。その後ＰX,k (n) はＳ／Ｎ比推
定部２７およびゲインファクター挿入部２８に転送さ
れ、Φk (n) は時間領域変換部２９に転送される。一
方、雑音判別部２３ではＡ／Ｄ変換部２１から転送され
たきたＸ（ｎ）に対して、まずパワーレベル ＰX (n) ＝Σ｛Ｘ(n-k) ｝2 （５） Σはｋ＝０からＬ−１までが計算される。ここで、Ｌは
積分時間を表わす。次に例えば予め決められたしきい値
Ｐthに対し、 ＰX (n) ＜Ｐth （６） の判定が行われ、この条件式を満たした場合には、雑音
であると判別する。雑音パワー計算部２６では、雑音判
別部２３において入力信号Ｘ（ｎ）が雑音であると判定
された時のみ、雑音の各帯域ごとのパワーレベルを ＰN,k (n) ＝Ｘk,r (n)2＋Ｘk,i (n)2 （７） として計算し、その時間平均ＰavN,k (n) をＳ／Ｎ比推
定部２７に転送する。時間平均は、例えば ＰavN,k (n) ＝（１／Ａ）Σm γm ＰN,k （ｎ−ｍ） （８） と計算される。ここでγm は例えば、 γm ＝（γ）m （９） と表わされるような指数重みづけの係数で（γ＜１）、
Ａは （１／Ａ）Σm γm ＝１ （10） となる正規化の為の定数である。
た、目的信号と不要な雑音とが混入する入力信号１５を
Ａ／Ｄ変換部２１においてデジタル化し、周波数帯域分
割部２２と雑音判別部２３に転送する。周波数帯域分割
部２２では、転送された信号が複数の周波数帯域に分割
される。周波数帯域への分割は、例えば離散的フーリエ
変換等を用いて行う。ここで、帯域分割された信号は一
般に複素数であるが、分割方法によっては実数となる場
合もある。ここでは一般的に、複素数を仮定して議論す
るが、実数の場合も同じ議論が可能である。周波数帯域
に分割されたｋ番目の周波数帯域の信号を Ｘk (n) ＝Ｘk,r (n) ＋ｊＸk,i (n) （２） （Ｘk,r ，Ｘk,i はそれぞれＸk (n) の実数部分と虚数
部分）とすると、Ｘk (n) は、入力信号パワー計算部２
４、入力信号位相計算部２５、雑音パワー計算部２６に
転送される。入力信号パワー計算部２４では、各帯域ご
との入力信号のパワーレベル ＰX,k (n) ＝Ｘk,r (n)2＋Ｘk,i (n)2 （３） が、入力信号位相計算部２５では各帯域ごとの位相 Φk (n) ＝tan -1［Ｘk,i (n) ／Ｘk,r (n) ］ （４） それぞれが計算される。その後ＰX,k (n) はＳ／Ｎ比推
定部２７およびゲインファクター挿入部２８に転送さ
れ、Φk (n) は時間領域変換部２９に転送される。一
方、雑音判別部２３ではＡ／Ｄ変換部２１から転送され
たきたＸ（ｎ）に対して、まずパワーレベル ＰX (n) ＝Σ｛Ｘ(n-k) ｝2 （５） Σはｋ＝０からＬ−１までが計算される。ここで、Ｌは
積分時間を表わす。次に例えば予め決められたしきい値
Ｐthに対し、 ＰX (n) ＜Ｐth （６） の判定が行われ、この条件式を満たした場合には、雑音
であると判別する。雑音パワー計算部２６では、雑音判
別部２３において入力信号Ｘ（ｎ）が雑音であると判定
された時のみ、雑音の各帯域ごとのパワーレベルを ＰN,k (n) ＝Ｘk,r (n)2＋Ｘk,i (n)2 （７） として計算し、その時間平均ＰavN,k (n) をＳ／Ｎ比推
定部２７に転送する。時間平均は、例えば ＰavN,k (n) ＝（１／Ａ）Σm γm ＰN,k （ｎ−ｍ） （８） と計算される。ここでγm は例えば、 γm ＝（γ）m （９） と表わされるような指数重みづけの係数で（γ＜１）、
Ａは （１／Ａ）Σm γm ＝１ （10） となる正規化の為の定数である。
【０００６】Ｓ／Ｎ比推定部２７では、各帯域ごとに入
力信号パワー計算部２４で計算されたＰX,k (n) 及び雑
音パワー計算部２６で計算された雑音パワーＰav
N,k (n) を用いて目的音声信号対雑音信号の比率である
Ｓ／Ｎ比が推定される。Ｓ／Ｎ比推定部２７について
は、処理の流れ図を図４Ａに記した。詳しい処理の説明
はこの図を用いて行う。
力信号パワー計算部２４で計算されたＰX,k (n) 及び雑
音パワー計算部２６で計算された雑音パワーＰav
N,k (n) を用いて目的音声信号対雑音信号の比率である
Ｓ／Ｎ比が推定される。Ｓ／Ｎ比推定部２７について
は、処理の流れ図を図４Ａに記した。詳しい処理の説明
はこの図を用いて行う。
【０００７】まず、ステップ３１において、 ＳＮＲk (n) ′＝ＰX,k (n) ／ＰN,k (n) （11） で定義されるＳＮＲk (n) ′を計算する。次にステップ
３２において、決定したＳＮＲk (n) ′に対して、（１
３）式で表わされる一時刻前の推定値（雑音低減された
パワー）ＰY,k （ｎ−１）を用いて平均化してＳＮＲk
(n) とする。即ち ＳＮＲk (n) ＝（１−β）Ｐ［ＳＮＲk (n) ′−１］ ＋β［ＰY,k （ｎ−１）／ＰN,k （ｎ−１）］ （12） とする。Ｐ［＊］は＊が正なら＊を、＊が負なら０をと
る。このＳＮＲk (n) 、必要に応じてＳＮＲk (n) ′
は、ゲインファクター計算部３０に転送される。
３２において、決定したＳＮＲk (n) ′に対して、（１
３）式で表わされる一時刻前の推定値（雑音低減された
パワー）ＰY,k （ｎ−１）を用いて平均化してＳＮＲk
(n) とする。即ち ＳＮＲk (n) ＝（１−β）Ｐ［ＳＮＲk (n) ′−１］ ＋β［ＰY,k （ｎ−１）／ＰN,k （ｎ−１）］ （12） とする。Ｐ［＊］は＊が正なら＊を、＊が負なら０をと
る。このＳＮＲk (n) 、必要に応じてＳＮＲk (n) ′
は、ゲインファクター計算部３０に転送される。
【０００８】ゲインファクター計算部３０では、Ｓ／Ｎ
比推定部２７より転送されてきたＳＮＲk (n) 、場合に
よると、これとＳＮＲk (n) ′を用いて、各雑音低減方
式で定義されている各周波数におけるゲインファクター
Ｇ（ＳＮＲk (n) ）を計算する。このゲインファクター
Ｇ（ＳＮＲk (n) ）は、ゲインファクター挿入部２８に
転送される。図５に各手法によるゲインファクターを表
わす。つまり図５中の上３つの手法ではＳＮＲk (n) の
みを用いてゲインファクターを求めるが、最も下のＭＭ
ＳＥ法による時は、ＳＮＲk (n) の他にＳＮＲk (n) ′
を用いる。
比推定部２７より転送されてきたＳＮＲk (n) 、場合に
よると、これとＳＮＲk (n) ′を用いて、各雑音低減方
式で定義されている各周波数におけるゲインファクター
Ｇ（ＳＮＲk (n) ）を計算する。このゲインファクター
Ｇ（ＳＮＲk (n) ）は、ゲインファクター挿入部２８に
転送される。図５に各手法によるゲインファクターを表
わす。つまり図５中の上３つの手法ではＳＮＲk (n) の
みを用いてゲインファクターを求めるが、最も下のＭＭ
ＳＥ法による時は、ＳＮＲk (n) の他にＳＮＲk (n) ′
を用いる。
【０００９】ゲインファクター挿入部２８では、各帯域
ごとに、ゲインファクター計算部３０において計算され
たゲインファクターを用いて雑音低減を行う。即ち入力
信号パワー計算部２４より転送されてきた帯域信号Ｐ
X,k (n) に対して、 ＰY,k (n) ＝Ｇ（ＳＮＲk (n) ）×ＰX,k (n) （13） を行い雑音を低減した帯域出力ＰY,k (n) を出力する。
ＰY,k (n) は時間領域変換部２９に転送され、入力信号
位相計算部２５から送られてきたΦk (n) を用いて、 Ｙk (n) ＝Ｙk,r (n) ＋ｊＹk,i (n) 但し、Ｙk,r (n) ＝ＰY,k (n) cos ［Φk (n) ］ Ｙk,i (n) ＝ＰY,k (n) sin ［Φk (n) ］ （14） に変換され、全帯域信号に合成され、更に例えば逆離散
的フーリエ変換により時間領域信号に変換される。この
結果を、Ｄ／Ａ変換部３４でアナログ信号にして雑音を
低減した信号１７、Ｙ（ｎ）を出力する。なお、Spectr
al subtractionでは（１３）式の演算で入力信号パワー
ＰX,k (n) から雑音パワーＰN,k (n) を引算した結果と
なる。
ごとに、ゲインファクター計算部３０において計算され
たゲインファクターを用いて雑音低減を行う。即ち入力
信号パワー計算部２４より転送されてきた帯域信号Ｐ
X,k (n) に対して、 ＰY,k (n) ＝Ｇ（ＳＮＲk (n) ）×ＰX,k (n) （13） を行い雑音を低減した帯域出力ＰY,k (n) を出力する。
ＰY,k (n) は時間領域変換部２９に転送され、入力信号
位相計算部２５から送られてきたΦk (n) を用いて、 Ｙk (n) ＝Ｙk,r (n) ＋ｊＹk,i (n) 但し、Ｙk,r (n) ＝ＰY,k (n) cos ［Φk (n) ］ Ｙk,i (n) ＝ＰY,k (n) sin ［Φk (n) ］ （14） に変換され、全帯域信号に合成され、更に例えば逆離散
的フーリエ変換により時間領域信号に変換される。この
結果を、Ｄ／Ａ変換部３４でアナログ信号にして雑音を
低減した信号１７、Ｙ（ｎ）を出力する。なお、Spectr
al subtractionでは（１３）式の演算で入力信号パワー
ＰX,k (n) から雑音パワーＰN,k (n) を引算した結果と
なる。
【００１０】これらの従来の方式では、２５６〜１０２
４の周波数帯域に分割して処理が行われる事が多いが、
これを会議システムにそのまま用いると大きな遅延が生
じる。会議を行う際に信号に遅延が生じると、通話性能
が劣化するという問題がある。一方、遅延を小さくする
為に周波数分割の数を減らすと、以下の（１５）式で定
義されているＳ／Ｎ比、 Ｓ／Ｎ＝１０ log10（ＰS ／ＰN ） ただし、ＰS ：目的信号の平均パワー ＰN ：雑音信号の平均パワー （15） が１５ｄＢ程度以上ある場合には、雑音低減によって音
声の品質が向上するが、１０ｄＢ程度以下の場合には、
雑音は低減されるが、それに伴い音声信号に歪みが生じ
たり、消し残された雑音が時間的に変化するのが原因で
聴感的に悪くなり品質が劣化することが判明した。これ
は、従来の雑音低減方式では、少数の帯域に分割された
信号で処理を行うと、非線形処理に起因する歪が大きく
なる為である。
４の周波数帯域に分割して処理が行われる事が多いが、
これを会議システムにそのまま用いると大きな遅延が生
じる。会議を行う際に信号に遅延が生じると、通話性能
が劣化するという問題がある。一方、遅延を小さくする
為に周波数分割の数を減らすと、以下の（１５）式で定
義されているＳ／Ｎ比、 Ｓ／Ｎ＝１０ log10（ＰS ／ＰN ） ただし、ＰS ：目的信号の平均パワー ＰN ：雑音信号の平均パワー （15） が１５ｄＢ程度以上ある場合には、雑音低減によって音
声の品質が向上するが、１０ｄＢ程度以下の場合には、
雑音は低減されるが、それに伴い音声信号に歪みが生じ
たり、消し残された雑音が時間的に変化するのが原因で
聴感的に悪くなり品質が劣化することが判明した。これ
は、従来の雑音低減方式では、少数の帯域に分割された
信号で処理を行うと、非線形処理に起因する歪が大きく
なる為である。
【００１１】即ち、従来法では、遅延を小さいままで、
効果的に雑音低減する事が出来ないという問題がある。
この為、音声会議装置・ＴＶ会議装置等、受聴を目的と
し音質が重要であり、またリアルタイム性が要求される
収音においては、この方法をそのまま適用する事はでき
ない。
効果的に雑音低減する事が出来ないという問題がある。
この為、音声会議装置・ＴＶ会議装置等、受聴を目的と
し音質が重要であり、またリアルタイム性が要求される
収音においては、この方法をそのまま適用する事はでき
ない。
【００１２】
【発明が解決しようとする課題】周囲雑音が混入した入
力信号の雑音を低減する方法として雑音の振幅分を音声
信号から減算する従来の方法では、周波数帯域分割数が
多いと遅延が大きくなり、周波数帯域分割数を減らすと
音声信号に歪が生じたり、引き残された雑音が時間的に
変化する為に聴感上好ましくない音をたてるという問題
があった。この発明の目的は、処理遅延が少なく、聴感
上の音質の劣化が少ない雑音低減処理方法を提供する事
である。
力信号の雑音を低減する方法として雑音の振幅分を音声
信号から減算する従来の方法では、周波数帯域分割数が
多いと遅延が大きくなり、周波数帯域分割数を減らすと
音声信号に歪が生じたり、引き残された雑音が時間的に
変化する為に聴感上好ましくない音をたてるという問題
があった。この発明の目的は、処理遅延が少なく、聴感
上の音質の劣化が少ない雑音低減処理方法を提供する事
である。
【００１３】
【課題を解決するための手段】この発明は、目的となる
音声信号と周囲雑音などの混在したマイクロホンでの受
音信号を複数の帯域に分割し、各々の帯域別の信号に対
し雑音パワーを推定し、推定された雑音パワーと実際に
入力されてきた入力信号パワーとを比較して音声信号と
雑音信号の比率を推定すると共に音声信号の振幅を雑音
低減された信号として推定し、推定された音声信号と雑
音信号の比率に基づいて計算された帯域ごとの雑音抑圧
の為のゲインファクターを、対応する帯域の入力信号に
掛け合わせる事によって得る。雑音低減された信号の推
定はその際生じる歪を少量の各帯域別の信号を、音声信
号と雑音信号の比率に応じて雑音が多い程、多く加算す
る事によってマスキングし、その結果歪が少ない効果的
な雑音低減を可能にする。
音声信号と周囲雑音などの混在したマイクロホンでの受
音信号を複数の帯域に分割し、各々の帯域別の信号に対
し雑音パワーを推定し、推定された雑音パワーと実際に
入力されてきた入力信号パワーとを比較して音声信号と
雑音信号の比率を推定すると共に音声信号の振幅を雑音
低減された信号として推定し、推定された音声信号と雑
音信号の比率に基づいて計算された帯域ごとの雑音抑圧
の為のゲインファクターを、対応する帯域の入力信号に
掛け合わせる事によって得る。雑音低減された信号の推
定はその際生じる歪を少量の各帯域別の信号を、音声信
号と雑音信号の比率に応じて雑音が多い程、多く加算す
る事によってマスキングし、その結果歪が少ない効果的
な雑音低減を可能にする。
【００１４】この発明は、次の様な特徴を持つ。まず、
雑音低減によって生じた歪を、少量の各帯域別の信号を
加算しマスキングする事によって音質の劣化を防ぐ。ま
た、各帯域別の信号の加算率を音声信号と雑音信号の比
率によって変化させる事によって、雑音が少ない、つま
り、雑音低減による処理歪が小さい場合には、各帯域別
信号の加算量を小さくし、雑音低減効果を最大限にする
事が可能である。
雑音低減によって生じた歪を、少量の各帯域別の信号を
加算しマスキングする事によって音質の劣化を防ぐ。ま
た、各帯域別の信号の加算率を音声信号と雑音信号の比
率によって変化させる事によって、雑音が少ない、つま
り、雑音低減による処理歪が小さい場合には、各帯域別
信号の加算量を小さくし、雑音低減効果を最大限にする
事が可能である。
【００１５】
【発明の実施の形態】図１にこの発明の実施例を適用し
た雑音低減装置の機能構成を示し、図３と対応するもの
については共通の記号を用いた。雑音パワー推定には、
図３に示したものでもよいがここでは、特願平８−６８
５４８や佐々木、羽田“損失制御を用いた帯域分割型雑
音低減方式について”春季日本音響学会予稿５１５〜５
１６（１９９６）で提案された方式を用いて推定を行う
事にして、以下の説明を行う。この方法によると、音声
と雑音とが混在している状態とみなす区間でも時間的に
音声がない区間があれば、雑音パワーを検出することが
できる。
た雑音低減装置の機能構成を示し、図３と対応するもの
については共通の記号を用いた。雑音パワー推定には、
図３に示したものでもよいがここでは、特願平８−６８
５４８や佐々木、羽田“損失制御を用いた帯域分割型雑
音低減方式について”春季日本音響学会予稿５１５〜５
１６（１９９６）で提案された方式を用いて推定を行う
事にして、以下の説明を行う。この方法によると、音声
と雑音とが混在している状態とみなす区間でも時間的に
音声がない区間があれば、雑音パワーを検出することが
できる。
【００１６】まず、マイクロホン１１により受音され
た、目的信号と不要な雑音等の混入する入力信号１５を
Ａ／Ｄ変換部２１においてデジタル化し、周波数帯域分
割部２２に転送する。周波数帯域分割部２２では、転送
された信号が周波数帯域に分割される。分割された各帯
域信号は、入力信号パワー計算部２４、ゲインファクタ
ー挿入部２８に転送される。以降、入力信号のｋ番目の
帯域信号をＸk (n) として、Ｘk (n) に対する処理の流
れを説明する。
た、目的信号と不要な雑音等の混入する入力信号１５を
Ａ／Ｄ変換部２１においてデジタル化し、周波数帯域分
割部２２に転送する。周波数帯域分割部２２では、転送
された信号が周波数帯域に分割される。分割された各帯
域信号は、入力信号パワー計算部２４、ゲインファクタ
ー挿入部２８に転送される。以降、入力信号のｋ番目の
帯域信号をＸk (n) として、Ｘk (n) に対する処理の流
れを説明する。
【００１７】入力信号パワー計算部２４では、転送され
てきたＸk (n) のパワーレベルを前記（７）式で計算
し、Ｓ／Ｎ比推定部２７、雑音パワー推定部５１に転送
される。雑音パワー推定部５１では転送されてきたＰ
X,k (n) を用いて雑音パワーＰav N,k (n) の推定が行わ
れる。雑音パワー推定部５１の処理の流れ図を図２Ａに
示す。
てきたＸk (n) のパワーレベルを前記（７）式で計算
し、Ｓ／Ｎ比推定部２７、雑音パワー推定部５１に転送
される。雑音パワー推定部５１では転送されてきたＰ
X,k (n) を用いて雑音パワーＰav N,k (n) の推定が行わ
れる。雑音パワー推定部５１の処理の流れ図を図２Ａに
示す。
【００１８】まずステップ６１において、転送されてき
たＰX,k (n) の時間平均パワーレベルＰX,k (n) を前記
（８）式、（９）式、（１０）式を用いて計算する。こ
の場合の平均時間としては、例えば５〜６ｍｓｅｃをと
る。次に、ステップ６２である一定時間におけるＰav
X,k (n) のレベル分布のヒストグラムをとる。つまりＰ
avX,k (n) が属するパワー区間の数を１加算、即ち、 ｈk (int PavX,k (n))＝ｈk (int PavX,k (n) ）＋１ （16） を行う。ｉｎｔ（＊）は小数点以下を切り捨て整数化す
る事を示す。更にステップ６３でヒストグラムｈk (i)
のピーク区間が検出され記憶される。即ち前後の値に対
して ｈk (i′）≦ｈk (i) （17） となるｉを求める。このｉの中で最も小さなｉを雑音の
パワーＰavN,x (n) とする。つまりピーク値が複数得ら
れた時、最小の値のピーク値を雑音パワーＰavN, k (n)
とし、これはＳ／Ｎ比推定部２７に転送される。
たＰX,k (n) の時間平均パワーレベルＰX,k (n) を前記
（８）式、（９）式、（１０）式を用いて計算する。こ
の場合の平均時間としては、例えば５〜６ｍｓｅｃをと
る。次に、ステップ６２である一定時間におけるＰav
X,k (n) のレベル分布のヒストグラムをとる。つまりＰ
avX,k (n) が属するパワー区間の数を１加算、即ち、 ｈk (int PavX,k (n))＝ｈk (int PavX,k (n) ）＋１ （16） を行う。ｉｎｔ（＊）は小数点以下を切り捨て整数化す
る事を示す。更にステップ６３でヒストグラムｈk (i)
のピーク区間が検出され記憶される。即ち前後の値に対
して ｈk (i′）≦ｈk (i) （17） となるｉを求める。このｉの中で最も小さなｉを雑音の
パワーＰavN,x (n) とする。つまりピーク値が複数得ら
れた時、最小の値のピーク値を雑音パワーＰavN, k (n)
とし、これはＳ／Ｎ比推定部２７に転送される。
【００１９】Ｓ／Ｎ比推定部２７では、入力信号パワー
計算部２４で計算されたＰX,k (n)及び雑音パワー推定
部５１で推定された雑音パワーＰavN,k (n) を用いて、
図４Ａで示した方法でＳＮＲk (n）′及びＳＮＲk (n）
が推定される。Ｓ／Ｎ比推定部２７で推定されたＳＮＲ
k (n）はゲインファクター計算部３０及び入力信号加算
率決定部５２に転送される。ゲインファクター計算部３
０で用いる計算法によってはＳＮＲk (n）′も入力され
る。
計算部２４で計算されたＰX,k (n)及び雑音パワー推定
部５１で推定された雑音パワーＰavN,k (n) を用いて、
図４Ａで示した方法でＳＮＲk (n）′及びＳＮＲk (n）
が推定される。Ｓ／Ｎ比推定部２７で推定されたＳＮＲ
k (n）はゲインファクター計算部３０及び入力信号加算
率決定部５２に転送される。ゲインファクター計算部３
０で用いる計算法によってはＳＮＲk (n）′も入力され
る。
【００２０】ゲンイファクター計算部３０では、Ｓ／Ｎ
比推定部２０１より転送されてきたＳＮＲk (n）、必要
に応じてＳＮＲk (n）′を用いて、ゲインファクターＧ
（ＳＮＲk (n））が決定される。ここで、ゲインファク
ターの具体的な計算は、従来技術の項で図５に示した方
法などが用いられる。ゲインファクター計算部３０で推
定されたゲインファクターＧ（ＳＮＲk (n））は、ゲイ
ンファクター挿入部２８に転送される。なお、Ｓ／Ｎ比
推定部２７の方法およびゲインファクター計算部３０の
方法は、この明細書に記載したもの以外の方法でもよ
い。
比推定部２０１より転送されてきたＳＮＲk (n）、必要
に応じてＳＮＲk (n）′を用いて、ゲインファクターＧ
（ＳＮＲk (n））が決定される。ここで、ゲインファク
ターの具体的な計算は、従来技術の項で図５に示した方
法などが用いられる。ゲインファクター計算部３０で推
定されたゲインファクターＧ（ＳＮＲk (n））は、ゲイ
ンファクター挿入部２８に転送される。なお、Ｓ／Ｎ比
推定部２７の方法およびゲインファクター計算部３０の
方法は、この明細書に記載したもの以外の方法でもよ
い。
【００２１】ゲインファクター挿入部２８では、ゲイン
ファクター計算部３０において計算されたゲインファク
ターＧ（ＳＮＲk (n））を用いて雑音低減を行う。即ち
周波数帯域分割部２２より転送されてきた帯域信号Ｘk
(n) に対して（１８）式の演算を行う。 Ｙ′k (n) ＝Ｇ（ＳＮＲk (n））×Ｘk (n) （18） この雑音を低減した信号Ｙ′k (n) を入力信号加算部５
３に転送する。
ファクター計算部３０において計算されたゲインファク
ターＧ（ＳＮＲk (n））を用いて雑音低減を行う。即ち
周波数帯域分割部２２より転送されてきた帯域信号Ｘk
(n) に対して（１８）式の演算を行う。 Ｙ′k (n) ＝Ｇ（ＳＮＲk (n））×Ｘk (n) （18） この雑音を低減した信号Ｙ′k (n) を入力信号加算部５
３に転送する。
【００２２】入力信号加算率決定部５２では、転送され
てきたＳＮＲk (n）を用いて、Ｓ／Ｎ比に基づいた入力
信号の加算率α′を決定する。α′は、ＳＮＲk (n）が
大きい時は小さい値をとる事が望ましい。例えば、 １５≦１０ log10［ＳＮＲk (n）］ ⇒ α′＝０ ５＜１０ log10［ＳＮＲk (n）］＜１５ ⇒ α′＝０．１ １０ log10［ＳＮＲk (n）］＜５ ⇒ α′＝０．３（19） の様に決定する。更に、この例では周波数帯域分割部２
２において帯域分割数Ｍに応じて決定される帯域分割数
ファクターβをα′に掛け合わせて入力信号加算率αと
した場合である。例えば α＝β×α′ １０２４≦Ｍ ⇒ β＝０．０１ Ｍ＜１０２４ ⇒ β＝１ （20） とする。つまり帯域分割数Ｍが大きい場合は加算率αを
小さくする。入力信号加算率αは、入力信号加算部５３
に転送される。
てきたＳＮＲk (n）を用いて、Ｓ／Ｎ比に基づいた入力
信号の加算率α′を決定する。α′は、ＳＮＲk (n）が
大きい時は小さい値をとる事が望ましい。例えば、 １５≦１０ log10［ＳＮＲk (n）］ ⇒ α′＝０ ５＜１０ log10［ＳＮＲk (n）］＜１５ ⇒ α′＝０．１ １０ log10［ＳＮＲk (n）］＜５ ⇒ α′＝０．３（19） の様に決定する。更に、この例では周波数帯域分割部２
２において帯域分割数Ｍに応じて決定される帯域分割数
ファクターβをα′に掛け合わせて入力信号加算率αと
した場合である。例えば α＝β×α′ １０２４≦Ｍ ⇒ β＝０．０１ Ｍ＜１０２４ ⇒ β＝１ （20） とする。つまり帯域分割数Ｍが大きい場合は加算率αを
小さくする。入力信号加算率αは、入力信号加算部５３
に転送される。
【００２３】入力信号加算部５３では、周波数帯域分割
部２２、ゲインファクター挿入部２８、入力信号加算率
決定部５２から転送されてきた帯域信号Ｘk (n) 、雑音
を低減した信号Ｙ′k (n) 及び入力信号加算率αを用い
て帯域出力信号Ｙk (n) を出力する。即ち、 Ｙk (n) ＝α×Ｘk (n) ＋（１−α）×Ｙ′k (n) （21） とする。Ｙk (n) は時間領域変換部２９に転送され、全
体域が合成された後時間領域に変換される。この結果
を、Ｄ／Ａ変換部３４でアナログ信号にして雑音を低減
した信号１７、Ｙ(n) を出力する。
部２２、ゲインファクター挿入部２８、入力信号加算率
決定部５２から転送されてきた帯域信号Ｘk (n) 、雑音
を低減した信号Ｙ′k (n) 及び入力信号加算率αを用い
て帯域出力信号Ｙk (n) を出力する。即ち、 Ｙk (n) ＝α×Ｘk (n) ＋（１−α）×Ｙ′k (n) （21） とする。Ｙk (n) は時間領域変換部２９に転送され、全
体域が合成された後時間領域に変換される。この結果
を、Ｄ／Ａ変換部３４でアナログ信号にして雑音を低減
した信号１７、Ｙ(n) を出力する。
【００２４】この本発明の有効性を、実際のＴＶ会議シ
ステムの状況を想定したオピニオン評価法を用いて評価
した。評価実験は、体積８７ｍ3 、残響時間３００ｍｓ
ｅｃ、暗騒音レベルは４６ｄＢ（Ａ）の可変残響室で行
った。帯域分割数Ｍは６４とした（これはサブバンド方
式エコーキャンセラで好ましいとされている分割数であ
る）。評価は３０名の被験者数で行った。作成した刺激
音を、スピーカから流し、それを拡声会議における受話
音声として、被験者に品質評価を行ってもらった。８ｄ
Ｂ、６ｄＢ、１ｄＢの３種類のＳ／Ｎ比の雑音付加音声
を図５に示した４種類の雑音低減方式で処理し、雑音低
減処理の前後で評価がどの様に変るかを調べた。
ステムの状況を想定したオピニオン評価法を用いて評価
した。評価実験は、体積８７ｍ3 、残響時間３００ｍｓ
ｅｃ、暗騒音レベルは４６ｄＢ（Ａ）の可変残響室で行
った。帯域分割数Ｍは６４とした（これはサブバンド方
式エコーキャンセラで好ましいとされている分割数であ
る）。評価は３０名の被験者数で行った。作成した刺激
音を、スピーカから流し、それを拡声会議における受話
音声として、被験者に品質評価を行ってもらった。８ｄ
Ｂ、６ｄＢ、１ｄＢの３種類のＳ／Ｎ比の雑音付加音声
を図５に示した４種類の雑音低減方式で処理し、雑音低
減処理の前後で評価がどの様に変るかを調べた。
【００２５】図２Ｂ（ａ）は、各雑音低減方式におけ
る、ＭＯＳ値の変化を示したものである。雑音低減処理
前と処理後を比較すると、Ｓ／Ｎ比の良い時は、ほぼ、
処理によってＭＯＳ値が雑音付加音声（原声）のそれよ
り上がるのに対し、Ｓ／Ｎ比が悪い時は、下がってしま
う。これは、帯域分割数が少ない場合、雑音低減による
品質改善よりも、処理による音声の歪や、消し残り雑音
による劣化が激しく、ＭＯＳ値を下げてしまう為だと考
えられる。この事から、Ｓ／Ｎ比が悪い場合には、従来
法では雑音低減による音質向上がみられないばかりか、
却って低下している事がわかる。
る、ＭＯＳ値の変化を示したものである。雑音低減処理
前と処理後を比較すると、Ｓ／Ｎ比の良い時は、ほぼ、
処理によってＭＯＳ値が雑音付加音声（原声）のそれよ
り上がるのに対し、Ｓ／Ｎ比が悪い時は、下がってしま
う。これは、帯域分割数が少ない場合、雑音低減による
品質改善よりも、処理による音声の歪や、消し残り雑音
による劣化が激しく、ＭＯＳ値を下げてしまう為だと考
えられる。この事から、Ｓ／Ｎ比が悪い場合には、従来
法では雑音低減による音質向上がみられないばかりか、
却って低下している事がわかる。
【００２６】図２Ｂ（ｂ）は、この発明により雑音低減
音に対し、原音を足し合わせた結果である。ここでは、
Ｓ／Ｎ比によらずα＝０．３とした。Ｓ／Ｎ比が悪い場
合に原音声を付加すると、付加前に比べＭＯＳ値がほぼ
０．５〜１上がり、その結果、雑音低減処理による音質
向上が見られる。また、Ｓ／Ｎ比に応じて加算率αを変
化させ、Ｓ／Ｎ比がいい場合には、付加率αを小さくす
ることで、更に効果的な音質向上が期待出来る事がわか
る。
音に対し、原音を足し合わせた結果である。ここでは、
Ｓ／Ｎ比によらずα＝０．３とした。Ｓ／Ｎ比が悪い場
合に原音声を付加すると、付加前に比べＭＯＳ値がほぼ
０．５〜１上がり、その結果、雑音低減処理による音質
向上が見られる。また、Ｓ／Ｎ比に応じて加算率αを変
化させ、Ｓ／Ｎ比がいい場合には、付加率αを小さくす
ることで、更に効果的な音質向上が期待出来る事がわか
る。
【００２７】ゲインファクター挿入を図３に示したよう
に、つまり（１３）式により行い、これより（１４）式
によりＹk (n) ′を得てもよい。先の実施例から理解さ
れるように、各帯域ごとの入力信号パワーと、対応帯域
ごとの雑音パワーとから、その帯域の目的音声信号に対
する雑音の比率を推定すると共に、その帯域の前記目的
音声信号の振幅を雑音低減された信号として推定すれば
よい。
に、つまり（１３）式により行い、これより（１４）式
によりＹk (n) ′を得てもよい。先の実施例から理解さ
れるように、各帯域ごとの入力信号パワーと、対応帯域
ごとの雑音パワーとから、その帯域の目的音声信号に対
する雑音の比率を推定すると共に、その帯域の前記目的
音声信号の振幅を雑音低減された信号として推定すれば
よい。
【００２８】
【発明の効果】以上述べたようにこの発明は、目的音声
と雑音が混在した信号を周波数帯域に分割し、帯域毎の
目的信号に対する雑音信号の比率を推定し、この推定結
果に基づいて帯域毎に雑音を低減するが、周波数帯域分
割数が少ない場合に生じる目的音声の歪を低減する為
に、雑音と音声の比率に応じて雑音低減音声に原信号を
その雑音が多い程多く加算することによって、歪の少な
い雑音低減音声が得られることができる。その際に帯域
分割数が少ない程、加算する比率を大とする。
と雑音が混在した信号を周波数帯域に分割し、帯域毎の
目的信号に対する雑音信号の比率を推定し、この推定結
果に基づいて帯域毎に雑音を低減するが、周波数帯域分
割数が少ない場合に生じる目的音声の歪を低減する為
に、雑音と音声の比率に応じて雑音低減音声に原信号を
その雑音が多い程多く加算することによって、歪の少な
い雑音低減音声が得られることができる。その際に帯域
分割数が少ない程、加算する比率を大とする。
【００２９】この発明方法による受聴を目的とした雑音
低減により、遅延が少なく、聞き易い目的信号を得る事
が可能になる。その結果、音声会議・ＴＶ会議等の拡声
通話系において、マイクロホンで受音し、相手側に送出
される送話信号に、目的となる音声以外の周囲雑音が混
入した場合でも、この発明方法による雑音の低減により
音声の明瞭性を保つ事が可能になり、通信品質が向上す
る。
低減により、遅延が少なく、聞き易い目的信号を得る事
が可能になる。その結果、音声会議・ＴＶ会議等の拡声
通話系において、マイクロホンで受音し、相手側に送出
される送話信号に、目的となる音声以外の周囲雑音が混
入した場合でも、この発明方法による雑音の低減により
音声の明瞭性を保つ事が可能になり、通信品質が向上す
る。
【図１】この発明の実施例を適用した雑音低減装置の機
能構成を示すブロック図。
能構成を示すブロック図。
【図２】Ａは図１中の雑音パワー推定部５１における処
理流れ図、Ｂ（ａ）は従来法による雑音低減効果を主観
評価した結果を示す図、（ｂ）はこの発明による雑音低
減効果を主観評価した結果を示す図である。
理流れ図、Ｂ（ａ）は従来法による雑音低減効果を主観
評価した結果を示す図、（ｂ）はこの発明による雑音低
減効果を主観評価した結果を示す図である。
【図３】従来の雑音低減方法を適用した装置の機能構成
を示すブロック図。
を示すブロック図。
【図４】Ａはゲインファクター計算部における処理の流
れ図、Ｂは雑音低減装置の原理説明図である。
れ図、Ｂは雑音低減装置の原理説明図である。
【図５】代表的なゲインファクターの計算式を示す図で
ある。
ある。
Claims (3)
- 【請求項１】 入力信号から雑音信号を除去した音声信
号を出力する雑音低減処理方法において、 前記入力信号を複数の周波数帯域に分割する帯域分割過
程と、 前記各帯域に分割された入力信号のパワーを計算する入
力信号パワー計算過程と、 前記各帯域に分割された入力信号中の雑音パワーの求め
る雑音パワー計算過程と、 前記各帯域毎の入力信号パワーと、前記各帯域毎の雑音
パワーとから前記各帯域毎の目的音声信号に対する雑音
の比率を推定すると共に、前記各帯域毎の前記目的音声
信号の振幅を雑音低域された信号として推定する推定過
程と、 前記推定された前記目的音声信号対雑音比率に基づい
て、前記帯域毎の前記雑音低減された信号と前記各帯域
に分割された入力信号との加算割合を決定する過程と、 前記決定された前記加算割合に基づいて、前記帯域毎の
前記雑音低減された信号と前記各帯域に分割された入力
信号を加算し、帯域出力信号を出力する過程と、 前記帯域出力信号を時間領域に変換し全帯域信号に合成
する過程と、 を有することを特徴とする帯域分割型雑音低減方法。 - 【請求項２】 請求項１記載の雑音低減方法において、 前記推定過程は前記各帯域毎の入力信号パワーと前記各
帯域毎の雑音パワーとを用いて、前記各帯域毎の前記目
的音声信号対雑音の比率を推定する過程と、 前記目的音声信号対雑音比率に基づいて前記各帯域毎の
ゲインファクターを決定するゲインファクター計算過程
と、 前記各帯域に分割された入力信号に対して前記各帯域毎
に決定された前記ゲインファクターを挿入して、前記各
帯域毎の前記雑音低減された信号を得るゲインファクタ
ー挿入過程とよりなることを特徴とする帯域分割型雑音
低減方法。 - 【請求項３】 請求項１又は２記載の雑音低減方法にお
いて、 前記帯域分割過程で分割する周波数帯域の数に応じて、
前記帯域毎の前記雑音低減された信号と前記各帯域に分
割された入力信号の加算割合を決定することを特徴とす
る帯域分割型雑音低減方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31778696A JP3454402B2 (ja) | 1996-11-28 | 1996-11-28 | 帯域分割型雑音低減方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31778696A JP3454402B2 (ja) | 1996-11-28 | 1996-11-28 | 帯域分割型雑音低減方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10161694A true JPH10161694A (ja) | 1998-06-19 |
| JP3454402B2 JP3454402B2 (ja) | 2003-10-06 |
Family
ID=18092036
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP31778696A Expired - Lifetime JP3454402B2 (ja) | 1996-11-28 | 1996-11-28 | 帯域分割型雑音低減方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3454402B2 (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001344000A (ja) * | 2000-05-31 | 2001-12-14 | Toshiba Corp | ノイズキャンセラとこのノイズキャンセラを備えた通信装置、並びにノイズキャンセル処理プログラムを記憶した記憶媒体 |
| WO2002101729A1 (fr) * | 2001-06-06 | 2002-12-19 | Mitsubishi Denki Kabushiki Kaisha | Attenuateur de bruit |
| JP2005257748A (ja) * | 2004-03-09 | 2005-09-22 | Nippon Telegr & Teleph Corp <Ntt> | 収音方法、収音装置、収音プログラム |
| JP2005321821A (ja) * | 2005-06-17 | 2005-11-17 | Nec Corp | ノイズ除去方法及び装置 |
| JP2008034928A (ja) * | 2006-07-26 | 2008-02-14 | Kenwood Corp | 音声入出力装置及びプログラム |
| JP2008306630A (ja) * | 2007-06-11 | 2008-12-18 | Fujitsu Ltd | 音量制御装置および方法 |
| WO2009047858A1 (ja) * | 2007-10-12 | 2009-04-16 | Fujitsu Limited | エコー抑圧システム、エコー抑圧方法、エコー抑圧プログラム、エコー抑圧装置、音出力装置、オーディオシステム、ナビゲーションシステム及び移動体 |
| JP2010102201A (ja) * | 2008-10-24 | 2010-05-06 | Yamaha Corp | 雑音抑圧装置及び雑音抑圧方法 |
| JP2013506878A (ja) * | 2009-10-01 | 2013-02-28 | クゥアルコム・インコーポレイテッド | オーディオ信号の雑音抑制 |
| JP2014085542A (ja) * | 2012-10-24 | 2014-05-12 | Fujitsu Ltd | オーディオ符号化装置、オーディオ符号化方法、オーディオ符号化用コンピュータプログラム |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4162604B2 (ja) | 2004-01-08 | 2008-10-08 | 株式会社東芝 | 雑音抑圧装置及び雑音抑圧方法 |
| JP4223441B2 (ja) * | 2004-06-03 | 2009-02-12 | 日本電信電話株式会社 | 雑音低減方法、この方法を実施する装置、プログラムおよびその記録媒体 |
-
1996
- 1996-11-28 JP JP31778696A patent/JP3454402B2/ja not_active Expired - Lifetime
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4580508B2 (ja) * | 2000-05-31 | 2010-11-17 | 株式会社東芝 | 信号処理装置及び通信装置 |
| JP2001344000A (ja) * | 2000-05-31 | 2001-12-14 | Toshiba Corp | ノイズキャンセラとこのノイズキャンセラを備えた通信装置、並びにノイズキャンセル処理プログラムを記憶した記憶媒体 |
| WO2002101729A1 (fr) * | 2001-06-06 | 2002-12-19 | Mitsubishi Denki Kabushiki Kaisha | Attenuateur de bruit |
| CN1308914C (zh) * | 2001-06-06 | 2007-04-04 | 三菱电机株式会社 | 噪声抑制装置和方法 |
| US7302065B2 (en) | 2001-06-06 | 2007-11-27 | Mitsubishi Denki Kabushiki Kaisha | Noise suppressor |
| JP4518817B2 (ja) * | 2004-03-09 | 2010-08-04 | 日本電信電話株式会社 | 収音方法、収音装置、収音プログラム |
| JP2005257748A (ja) * | 2004-03-09 | 2005-09-22 | Nippon Telegr & Teleph Corp <Ntt> | 収音方法、収音装置、収音プログラム |
| JP2005321821A (ja) * | 2005-06-17 | 2005-11-17 | Nec Corp | ノイズ除去方法及び装置 |
| JP2008034928A (ja) * | 2006-07-26 | 2008-02-14 | Kenwood Corp | 音声入出力装置及びプログラム |
| JP4580409B2 (ja) * | 2007-06-11 | 2010-11-10 | 富士通株式会社 | 音量制御装置および方法 |
| JP2008306630A (ja) * | 2007-06-11 | 2008-12-18 | Fujitsu Ltd | 音量制御装置および方法 |
| US8218777B2 (en) * | 2007-06-11 | 2012-07-10 | Fujitsu Limited | Multipoint communication apparatus |
| WO2009047871A1 (ja) * | 2007-10-12 | 2009-04-16 | Fujitsu Limited | エコー抑圧システム、エコー抑圧方法、エコー抑圧プログラム、エコー抑圧装置、音出力装置、オーディオシステム、ナビゲーションシステム及び移動体 |
| WO2009047858A1 (ja) * | 2007-10-12 | 2009-04-16 | Fujitsu Limited | エコー抑圧システム、エコー抑圧方法、エコー抑圧プログラム、エコー抑圧装置、音出力装置、オーディオシステム、ナビゲーションシステム及び移動体 |
| US8340963B2 (en) | 2007-10-12 | 2012-12-25 | Fujitsu Limited | Echo suppressing system, echo suppressing method, recording medium, echo suppressor, sound output device, audio system, navigation system and mobile object |
| JP2010102201A (ja) * | 2008-10-24 | 2010-05-06 | Yamaha Corp | 雑音抑圧装置及び雑音抑圧方法 |
| JP2013506878A (ja) * | 2009-10-01 | 2013-02-28 | クゥアルコム・インコーポレイテッド | オーディオ信号の雑音抑制 |
| US8571231B2 (en) | 2009-10-01 | 2013-10-29 | Qualcomm Incorporated | Suppressing noise in an audio signal |
| JP2014085542A (ja) * | 2012-10-24 | 2014-05-12 | Fujitsu Ltd | オーディオ符号化装置、オーディオ符号化方法、オーディオ符号化用コンピュータプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3454402B2 (ja) | 2003-10-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8010355B2 (en) | Low complexity noise reduction method | |
| US7454010B1 (en) | Noise reduction and comfort noise gain control using bark band weiner filter and linear attenuation | |
| US8521530B1 (en) | System and method for enhancing a monaural audio signal | |
| US8180064B1 (en) | System and method for providing voice equalization | |
| US9502048B2 (en) | Adaptively reducing noise to limit speech distortion | |
| AU771444B2 (en) | Noise reduction apparatus and method | |
| US9438992B2 (en) | Multi-microphone robust noise suppression | |
| US7555075B2 (en) | Adjustable noise suppression system | |
| US7649988B2 (en) | Comfort noise generator using modified Doblinger noise estimate | |
| JP3099870B2 (ja) | 音響エコーキャンセラ | |
| JP3459363B2 (ja) | 雑音低減処理方法、その装置及びプログラム記憶媒体 | |
| JPH09503590A (ja) | 会話の品質向上のための背景雑音の低減 | |
| JP3454403B2 (ja) | 帯域分割型雑音低減方法及び装置 | |
| US20180308503A1 (en) | Real-time single-channel speech enhancement in noisy and time-varying environments | |
| JP3454402B2 (ja) | 帯域分割型雑音低減方法 | |
| CN107045874B (zh) | 一种基于相关性的非线性语音增强方法 | |
| US20030033139A1 (en) | Method and circuit arrangement for reducing noise during voice communication in communications systems | |
| JP3309895B2 (ja) | 雑音低減方法 | |
| Yamashita et al. | Spectral subtraction iterated with weighting factors | |
| CN111933169B (zh) | 一种二次利用语音存在概率的语音降噪方法 | |
| Sasaki et al. | Noise reduction for subband acoustic echo canceller | |
| Sabnis et al. | Simulation of Speech Enhancement Algorithms and its Performance Comparison | |
| Wang et al. | Time-Frequency Thresholding: A new algorithm in wavelet package speech enhancement | |
| PORUBA | Subtractive-type algorithm utilizing the human ear masking characteristics | |
| Shruthi et al. | Speech intelligibility prediction and near end listening enhancement for mobile appliciation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080725 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080725 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090725 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090725 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100725 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100725 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110725 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120725 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130725 Year of fee payment: 10 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |