JP5609157B2 - 係数設定装置および雑音抑圧装置 - Google Patents
係数設定装置および雑音抑圧装置 Download PDFInfo
- Publication number
- JP5609157B2 JP5609157B2 JP2010041950A JP2010041950A JP5609157B2 JP 5609157 B2 JP5609157 B2 JP 5609157B2 JP 2010041950 A JP2010041950 A JP 2010041950A JP 2010041950 A JP2010041950 A JP 2010041950A JP 5609157 B2 JP5609157 B2 JP 5609157B2
- Authority
- JP
- Japan
- Prior art keywords
- suppression
- power
- coefficient
- noise
- setting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2227/00—Details of public address [PA] systems covered by H04R27/00 but not provided for in any of its subgroups
- H04R2227/001—Adaptation of signal processing in PA systems in dependence of presence of noise
Description
|Y(f)|K=|X(f)|K−a|N(f)|K ……(A)
図1は、本発明の第1実施形態に係る雑音抑圧装置100のブロック図である。雑音抑圧装置100には信号供給装置12と放音装置14と入力装置16とが接続される。信号供給装置12は、音響信号x(t)を雑音抑圧装置100に供給する。音響信号x(t)は、以下の数式(1)で示すように、目的音成分(例えば音声や楽音等の音響)s(t)と雑音成分n(t)との混合音の波形を示す時間領域の信号である。
周囲の音響を収音して音響信号x(t)を生成する収音機器や、可搬型または内蔵型の記録媒体から音響信号x(t)を取得して雑音抑圧装置100に出力する再生装置や、通信網から音響信号x(t)を受信して雑音抑圧装置100に出力する通信装置が信号供給装置12として採用され得る。
数式(2)の記号jは虚数単位を意味し、記号θx(f,τ)は音響信号x(t)の位相角(位相スペクトル)を意味する。音響信号y(t)の振幅(振幅スペクトル)|Y(f,τ)|は、以下の数式(3A)および数式(3B)で定義されるように、音響信号x(t)(振幅|X(f,τ)|)から雑音成分n(t)(振幅|N(f,τ)|)を抑圧することで算定される。
抑圧処理前の確率密度関数P(x)の確率分布D1は、数式(3A)の累乗処理(K乗)で図2の部分(B)の確率分布D2に変化する。確率変数xに対する写像gを想定すると、変化後の確率分布D2を表す確率密度関数P(y)(y=g(x))は、以下の数式(8)で表現される。
数式(8)の記号|J|は、以下の数式(9)で定義されるヤコビアンを意味する。
累乗処理後の確率密度関数P(y)の確率分布D2は、数式(3A)および数式(3B)の減算処理で図2の部分(C)の確率分布D3に変化する。確率分布D3は、図2の部分(C)に矢印で示すように、雑音成分n(t)の期待値E[y]と抑圧係数βとの乗算値に応じた度合で確率分布D2を確率変数yの負側に平行移動し(数式(3A))、移動後に負数となる確率変数yの確率(度数)の総和を確率変数yのゼロの確率に累算した(数式(3B))形状となる。したがって、確率分布D3の確率密度関数Pss(y)は、以下の数式(15A)および数式(15B)で表現される。
数式(15A)および数式(15B)の確率密度関数Pss(y)は、数式(3A)の累乗根処理で、パワーに相当する確率変数で定義される確率密度関数Pss(x)に変換される。累乗根処理後の確率密度関数Pss(x)は、数式(15)および数式(15B)の変数yを累乗処理と同様の方法で変数x(x=|y(f,τ)|2)に置換した以下の数式(16A)および数式(16B)で表現される。
数式(18)の記号Γ(α,w)は、以下の数式(19)で定義される第2種不完全ガンマ関数を意味する。
数式(22)の記号sは信号成分(強調の目的となる音声成分)を意味し、記号nは雑音成分を意味する。また、添字inは雑音抑圧前を意味し、添字outは雑音抑圧後を意味する。すなわち、数式(22)の分母が雑音抑圧前のSN比に相当し、数式(22)の分子が雑音抑圧後のSN比に相当する。
次に、本発明の第2実施形態について説明する。第1実施形態では、音響信号x(t)(振幅|X(f,τ)|)から雑音成分n(t)(振幅|N(f,τ)|)を減算することで音響信号y(t)の振幅|Y(f,τ)|を算定したが、音響信号y(t)を生成するための演算は減算(スペクトル減算)に限定されない。第2実施形態では、音響信号x(t)の振幅|X(f,τ)|に対する所定の係数(ゲイン)の乗算で音響信号y(t)の振幅|Y(f,τ)|を算定する。なお、以下の各例示において作用や機能が第1実施形態と同等である要素については、以上と同じ参照符号を流用して各々の詳細な説明を適宜に省略する。
数式(25)から理解されるように、音響信号x(t)における雑音成分n(t)の振幅|N(f,τ)|が大きい周波数fほど係数値γ(f)は小さい数値に設定される。したがって、音響信号x(t)のうち雑音成分n(t)の振幅|N(f,τ)|が大きい周波数fほど振幅|X(f,τ)|を抑制した音響信号y(t)(第1実施形態と同様に雑音成分n(t)を抑圧した音響信号)が生成される。
以上の各形態は多様に変形され得る。具体的な変形の態様を以下に例示する。以下の例示から任意に選択された2以上の態様は適宜に併合され得る。
各変数の設定の方法は適宜に変更される。例えば、以上の各形態では冪指数Kを利用者が指示する構成を例示したが、指数設定部54が自動的に(すなわち利用者からの指示を必要とせずに)冪指数Kを設定する構成も採用され得る。例えば、指数設定部54は、演算処理装置22の演算性能に応じて冪指数K(例えば浮動小数点数等の演算性能の制約の範囲内で最小の冪指数K)を設定する。冪指数Kの設定の方法に関わらず、第1実施形態と同様に、指数設定部54が冪指数Kを0.1未満(更に好適には0.01未満)の正数に設定する構成が好適である。また、以上の各形態では、形状母数αと雑音抑圧率NRRの目標値N0とを可変に設定したが、形状母数αおよび目標値N0の少なくとも一方を所定値に固定した構成も採用され得る。したがって、母数設定部56や抑圧率設定部52は省略され得る。
以上の各形態では、係数設定部58が数式(24)の演算を実行することで抑圧係数βを算定したが、冪指数K(さらには形状母数αや雑音抑圧率NRR)に応じた抑圧係数βを特定する方法は適宜に変更される。例えば、冪指数Kと形状母数αと雑音抑圧率NRRの目標値N0との各数値の組合せに対して抑圧係数βを対応させたデータテーブルを記憶装置24に格納し、各変数の指示値(K,α,N0)に対応する抑圧係数βを、記憶装置24のデータテーブルから係数設定部58が検索して雑音抑圧部42に指示する構成も採用され得る。
第1実施形態の雑音抑圧(数式(3A))や第2実施形態の係数列Gの算定(数式(25))では、雑音成分n(t)の振幅|N(f,τ)|をK乗してから時間平均したが(Eτ[|N(f,τ)|K])、雑音成分n(t)の振幅|N(f,τ)|を時間平均してからK乗する構成({Eτ[|N(f,τ)|]}K)も採用され得る。すなわち、冪指数Kで累乗される雑音成分n(t)の振幅は、時間平均前の振幅|N(f,τ)|および時間平均後の振幅Eτ[|N(f,τ)|]の何れでもよい。なお、雑音成分n(t)の時間平均を省略した構成(例えば1個のフレームの振幅|N(f,τ)|のK乗を抑圧係数βに応じて振幅|X(f,τ)|から減算する構成)も採用され得る。
以上の各形態では、音響信号x(t)から雑音成分n(t)を減算した数値(|X(f,τ)|K−βEτ[|M(f,τ)|K])が負数となる場合に音響信号y(t)の振幅|Y(f,τ)|をゼロに設定(フロアリング処理)したが、フロアリング処理に適用される数値はゼロに限定されない。例えば、音響信号x(t)から雑音成分n(t)を減算した数値が負数となる周波数fの振幅|Y(f,τ)|を、振幅|X(f,τ)|や振幅|N(f,τ)|に応じた数値(例えば数値a1|X(f,τ)|や数値a2|N(f,τ)|(係数a1や係数a2は所定値に設定される))に設定する構成も採用され得る。
以上の各形態では、変数制御部44と雑音抑圧部42とを含む雑音抑圧装置100を例示したが、雑音抑圧に適用される抑圧係数βを設定する係数設定装置としても本発明は特定され得る。係数設定装置は、以上の各形態における変数制御部44を含んで構成される。係数設定装置が雑音抑圧部42と一体に構成される(すなわち以上の各形態の雑音抑圧装置100を構成する)か、係数設定装置と雑音抑圧部42(雑音抑圧装置)とが別体に構成されるかは本発明において不問である。
Claims (10)
- 雑音成分の各周波数の振幅のK乗(冪指数Kは正数)を音響信号の各周波数の振幅のK乗から抑圧する度合を示す抑圧係数を設定する係数設定装置であって、
冪指数Kを設定する指数設定手段と、
前記指数設定手段が設定した冪指数Kが小さいほど前記雑音成分の抑圧の度合が減少するように前記冪指数Kに応じて前記抑圧係数を可変に設定する係数設定手段と
を具備する係数設定装置。 - 雑音抑圧率の目標値を設定する抑圧率設定手段と、
音響信号の強度分布を近似する確率分布の形状母数を前記音響信号から算定する母数設定手段とを具備し、
前記係数設定手段は、前記指数設定手段が設定した冪指数Kと前記抑圧率設定手段が設定した雑音抑圧率の目標値と前記母数設定手段が算定した形状母数とに応じて前記抑圧係数を設定する
請求項1の係数設定装置。 - 前記係数設定手段は、前記母数設定手段が算定する形状母数が大きいほど前記雑音成分の抑圧の度合が減少するように前記形状母数に応じて前記抑圧係数を可変に設定する
請求項2の係数設定装置。 - 前記指数設定手段は、前記冪指数Kを0.1未満の数値に設定する
請求項1から請求項3の何れかの係数設定装置。 - 冪指数K(正数)を設定する指数設定手段と、
前記指数設定手段が設定した冪指数Kが小さいほど雑音成分の抑圧の度合が減少するように前記冪指数Kに応じて抑圧係数を可変に設定する係数設定手段と、
前記係数設定手段が設定した抑圧係数に応じた度合で雑音成分の各周波数の振幅のK乗を音響信号の各周波数の振幅のK乗から抑圧する処理を含む雑音抑圧で前記雑音成分の抑圧後の音響信号を生成する雑音抑圧手段と
を具備する雑音抑圧装置。 - 雑音抑圧率の目標値を設定する抑圧率設定手段と、
音響信号の強度分布を近似する確率分布の形状母数を前記音響信号から算定する母数設定手段とを具備し、
前記係数設定手段は、前記指数設定手段が設定した冪指数Kと前記抑圧率設定手段が設定した雑音抑圧率の目標値と前記母数設定手段が算定した形状母数とに応じて前記抑圧係数を設定する
請求項5の雑音抑圧装置。 - 雑音成分の各周波数の振幅のK乗を音響信号の各周波数の振幅のK乗から抑圧する処理を含む雑音抑圧で前記雑音成分の抑圧後の音響信号を生成する雑音抑圧手段を具備し、
前記冪指数Kは0.1未満の正数に設定される
雑音抑圧装置。 - 雑音成分の各周波数の振幅のK乗(冪指数Kは正数)を音響信号の各周波数の振幅のK乗から抑圧する度合を示す抑圧係数を設定するためのプログラムであって、
冪指数Kを設定する指数設定処理と、
前記指数設定処理で設定した冪指数Kが小さいほど前記雑音成分の抑圧の度合が減少するように前記冪指数Kに応じて前記抑圧係数を可変に設定する係数設定処理と
をコンピュータに実行させるプログラム。 - 冪指数K(正数)を設定する指数設定処理と、
前記指数設定処理で設定した冪指数Kが小さいほど雑音成分の抑圧の度合が減少するように前記冪指数Kに応じて抑圧係数を可変に設定する係数設定処理と、
前記係数設定処理で設定した抑圧係数に応じた度合で雑音成分の各周波数の振幅のK乗を音響信号の各周波数の振幅のK乗から抑圧する処理を含む雑音抑圧で前記雑音成分の抑圧後の音響信号を生成する雑音抑圧処理と
をコンピュータに実行させるプログラム。 - 雑音成分の各周波数の振幅のK乗を音響信号の各周波数の振幅のK乗から抑圧する処理を含む雑音抑圧で前記雑音成分の抑圧後の音響信号を生成する雑音抑圧処理をコンピュータに実行させるプログラムであって、
前記冪指数Kは0.1未満の正数に設定される
プログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010041950A JP5609157B2 (ja) | 2010-02-26 | 2010-02-26 | 係数設定装置および雑音抑圧装置 |
US12/932,473 US20110211711A1 (en) | 2010-02-26 | 2011-02-25 | Factor setting device and noise suppression apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010041950A JP5609157B2 (ja) | 2010-02-26 | 2010-02-26 | 係数設定装置および雑音抑圧装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011180219A JP2011180219A (ja) | 2011-09-15 |
JP5609157B2 true JP5609157B2 (ja) | 2014-10-22 |
Family
ID=44505267
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010041950A Expired - Fee Related JP5609157B2 (ja) | 2010-02-26 | 2010-02-26 | 係数設定装置および雑音抑圧装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20110211711A1 (ja) |
JP (1) | JP5609157B2 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5633673B2 (ja) * | 2010-05-31 | 2014-12-03 | ヤマハ株式会社 | 雑音抑圧装置およびプログラム |
US9318125B2 (en) * | 2013-01-15 | 2016-04-19 | Intel Deutschland Gmbh | Noise reduction devices and noise reduction methods |
JP6169849B2 (ja) * | 2013-01-15 | 2017-07-26 | 本田技研工業株式会社 | 音響処理装置 |
JP6559576B2 (ja) * | 2016-01-05 | 2019-08-14 | 株式会社東芝 | 雑音抑圧装置、雑音抑圧方法及びプログラム |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4628529A (en) * | 1985-07-01 | 1986-12-09 | Motorola, Inc. | Noise suppression system |
US4811404A (en) * | 1987-10-01 | 1989-03-07 | Motorola, Inc. | Noise suppression system |
IL84948A0 (en) * | 1987-12-25 | 1988-06-30 | D S P Group Israel Ltd | Noise reduction system |
JP3484757B2 (ja) * | 1994-05-13 | 2004-01-06 | ソニー株式会社 | 音声信号の雑音低減方法及び雑音区間検出方法 |
US5544250A (en) * | 1994-07-18 | 1996-08-06 | Motorola | Noise suppression system and method therefor |
US6122384A (en) * | 1997-09-02 | 2000-09-19 | Qualcomm Inc. | Noise suppression system and method |
US7392181B2 (en) * | 2004-03-05 | 2008-06-24 | Siemens Corporate Research, Inc. | System and method for nonlinear signal enhancement that bypasses a noisy phase of a signal |
WO2007089189A1 (en) * | 2006-01-31 | 2007-08-09 | Telefonaktiebolaget Lm Ericsson (Publ). | Non-intrusive signal quality assessment |
JP5152800B2 (ja) * | 2008-07-09 | 2013-02-27 | 国立大学法人 奈良先端科学技術大学院大学 | 雑音抑圧評価装置およびプログラム |
JP5152799B2 (ja) * | 2008-07-09 | 2013-02-27 | 国立大学法人 奈良先端科学技術大学院大学 | 雑音抑圧装置およびプログラム |
JP2010220087A (ja) * | 2009-03-18 | 2010-09-30 | Yamaha Corp | 音響処理装置およびプログラム |
-
2010
- 2010-02-26 JP JP2010041950A patent/JP5609157B2/ja not_active Expired - Fee Related
-
2011
- 2011-02-25 US US12/932,473 patent/US20110211711A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
US20110211711A1 (en) | 2011-09-01 |
JP2011180219A (ja) | 2011-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5641186B2 (ja) | 雑音抑圧装置およびプログラム | |
CA2732723C (en) | Apparatus and method for processing an audio signal for speech enhancement using a feature extraction | |
JP6169849B2 (ja) | 音響処理装置 | |
JP2012155339A (ja) | 音声状態モデルを使用したマルチセンサ音声高品質化 | |
JP6339896B2 (ja) | 雑音抑圧装置および雑音抑圧方法 | |
JP3960834B2 (ja) | 音声強調装置及び音声強調方法 | |
US9418677B2 (en) | Noise suppressing device, noise suppressing method, and a non-transitory computer-readable recording medium storing noise suppressing program | |
JP5609157B2 (ja) | 係数設定装置および雑音抑圧装置 | |
CN112712816A (zh) | 语音处理模型的训练方法和装置以及语音处理方法和装置 | |
JP5152799B2 (ja) | 雑音抑圧装置およびプログラム | |
JP5454330B2 (ja) | 音響処理装置 | |
JP5633673B2 (ja) | 雑音抑圧装置およびプログラム | |
JP5942388B2 (ja) | 雑音抑圧用係数設定装置、雑音抑圧装置および雑音抑圧用係数設定方法 | |
JP5152800B2 (ja) | 雑音抑圧評価装置およびプログラム | |
US20170323656A1 (en) | Signal processor | |
JP5728903B2 (ja) | 音響処理装置およびプログラム | |
JP4533126B2 (ja) | 近接音分離収音方法、近接音分離収音装置、近接音分離収音プログラム、記録媒体 | |
JP5772723B2 (ja) | 音響処理装置および分離マスク生成装置 | |
JP2015169901A (ja) | 音響処理装置 | |
JP2013250356A (ja) | 係数設定装置および雑音抑圧装置 | |
JP2013182161A (ja) | 音響処理装置およびプログラム | |
Fallis et al. | Efficiency of Spectral Subtraction Algorithms for an Urban Audio Acquisition System Using IoT Devices | |
JP2015169900A (ja) | 雑音抑圧装置 | |
JP2015004959A (ja) | 音響処理装置 | |
JP2014077916A (ja) | 音響処理装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130124 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20130128 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130418 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130925 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131118 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131203 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20140117 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140127 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20140120 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140218 |
|
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: 20140805 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140818 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5609157 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |