(57)【要約】【課題】所定のゲイン関数を利用して、ミュージカルノイズが低減された音声強調を行う。【解決手段】音声処理装置は、目的音および雑音が混入している入力音声の目的音を強調して音声周波数成分を取得する目的音強調部１０２と、入力音声の目的音を抑圧して雑音周波数成分を取得する目的音抑圧部１０４と、音声周波数成分および雑音周波数成分に応じた所定のゲイン関数を用いて音声周波数成分に乗算するゲイン値を算出するゲイン算出部１０６と、ゲイン算出部１０６により算出されたゲイン値を音声周波数成分に乗算するゲイン乗算部１０８と、を備え、ゲイン算出部１０６は、音声周波数成分と雑音周波数成分とのエネルギー比が所定値以下の場合にゲイン値およびゲイン関数の傾きが所定値より小さくなるゲイン関数を用いてゲイン値を算出する。【選択図】図３
2011191669.tif 000002
【特許請求の範囲】
【請求項１】
目的音および雑音が混入している入力音声の前記目的音を強調して音声周波数成分を取得する目的音強調部と、
前記入力音声の前記目的音を抑圧して雑音周波数成分を取得する目的音抑圧部と、
前記音声周波数成分および前記雑音周波数成分に応じた所定のゲイン関数を用いて前記音声周波数成分に乗算するゲイン値を算出するゲイン算出部と、
前記ゲイン算出部により算出されたゲイン値を前記音声周波数成分に乗算するゲイン乗算部と、
を備え、
前記ゲイン算出部は、前記音声周波数成分と前記雑音周波数成分とのエネルギー比が所定値以下の場合に前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなる前記ゲイン関数を用いて前記ゲイン値を算出する、音声処理装置。
【請求項２】
前記音声周波数成分には目的音成分と雑音成分が含まれており、前記ゲイン乗算部は、前記音声周波数成分に前記ゲイン値を乗算して前記音声周波数成分に含まれている前記雑音成分を抑圧する、請求項１に記載の音声処理装置。
【請求項３】
前記ゲイン算出部は、前記目的音抑圧部により取得された雑音周波数成分に雑音のみが含まれていると推定して、前記ゲイン値を算出する、請求項１に記載の音声処理装置。
【請求項４】
前記ゲイン関数は、前記音声周波数成分と前記雑音周波数成分とのエネルギー比において、雑音の比率が集中している雑音集中範囲の前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなるゲインカーブを有する関数である、請求項１に記載の音声処理装置。
【請求項５】
前記ゲイン関数は、前記雑音集中範囲以外で最も急斜である前記ゲイン関数の傾きよりも小さい傾きのゲインカーブを有する関数である、請求項４に記載の音声処理装置。
【請求項６】
前記入力音声に含まれる前記目的音が存在する区間を検出する目的音区間検出部を備え、
前記ゲイン算出部は、前記目的音区間検出部による検出結果に応じて、前記目的音強調部により取得された前記音声周波数成分のパワースペクトルおよび前記目的音抑圧部により取得された前記雑音周波数成分のパワースペクトルを平均化する、請求項１に記載の音声処理装置。
【請求項７】
前記ゲイン算出部は、前記目的音区間検出部による検出の結果、目的音が存在する区間であることが検出された場合に第１の平滑化係数を選択し、前記目的音が存在する区間であることが検出された場合に第２の平滑化係数を選択して、前記音声周波数成分および前記雑音周波数成分のパワースペクトルを平均化する、請求項６に記載の音声処理装置。
【請求項８】
前記ゲイン算出部は、平均化された前記音声周波数成分のパワースペクトルおよび前記雑音周波数成分のパワースペクトルを用いて算出されたゲイン値を平均化する、請求項６に記載の音声処理装置。
【請求項９】
前記目的音抑圧部により取得された雑音周波数成分の大きさを、前記目的音強調部により取得された音声周波数成分に含まれる雑音成分の大きさに対応させるように前記雑音周波数成分を補正する雑音補正部を備え、
前記ゲイン算出部は、前記雑音補正部により補正された前記雑音周波数成分に応じたゲイン値を算出する、請求項１に記載の音声処理装置。
【請求項１０】
前記雑音補正部は、ユーザ操作に応じて前記雑音周波数成分を補正する、請求項９に記載の音声処理装置。
【請求項１１】
前記雑音補正部は、検出された雑音の状態に応じて前記雑音周波数成分を補正する、請求項９に記載の音声処理装置。
【請求項１２】
目的音および雑音が混入している入力音声の前記目的音を強調して音声周波数成分を取得するステップと、
前記入力音声の前記目的音を抑圧して雑音周波数成分を取得するステップと、
前記音声周波数成分と前記雑音周波数成分とのエネルギー比が所定値以下の場合に前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなるゲイン関数を用いて前記音声周波数成分に乗算するゲイン値を算出するステップと、
前記ゲイン算出部により算出されたゲイン値を前記音声周波数成分に乗算するステップと、
を含む、音声処理方法。
【請求項１３】
コンピュータを、
目的音および雑音が混入している入力音声の前記目的音を強調して音声周波数成分を取得する目的音強調部と、
前記入力音声の前記目的音を抑圧して雑音周波数成分を取得する目的音抑圧部と、
前記音声周波数成分および前記雑音周波数成分に応じた所定のゲイン関数を用いて前記音声周波数成分に乗算するゲイン値を算出するゲイン算出部と、
前記ゲイン算出部により算出されたゲイン値を前記音声周波数成分に乗算するゲイン乗算部と、
を備え、
前記ゲイン算出部は、前記音声周波数成分と前記雑音周波数成分とのエネルギー比が所定値以下の場合に前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなる前記ゲイン関数を用いて前記ゲイン値を算出する、音声処理装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音声処理装置、音声処理方法およびプログラムに関する。
【背景技術】
【０００２】
従来から、雑音が混入している入力音声に対して、雑音を抑圧して目的音声を強調することが行われている（例えば、特許文献１〜３）。上記特許文献では、目的音声を強調した音声周波数成分には目的音声と雑音が含まれており、目的音声を抑圧した雑音周波数成分には雑音のみが含まれていると推定して、音声周波数成分のパワースペクトルから雑音周波数成分のパワースペクトルを減算することにより、入力音声から雑音音声を除去している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第３６７７１４３号公報
【特許文献２】特許第４１６３２９４号公報
【特許文献３】特許公開２００９−４９９９８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、上記特許文献では、処理後の音声信号にミュージカルノイズといわれる特有の歪みが生じたり、音声周波数成分に含まれる雑音と雑音周波数成分に含まれる雑音とが等しくない場合があったりするため、適切な雑音除去を行うことができないという問題があった。
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、所定のゲイン関数を利用して、ミュージカルノイズが低減された音声強調を行うことが可能な、新規かつ改良された音声処理装置、音声処理方法およびプログラムを提供することにある。
【課題を解決するための手段】
【０００５】
上記課題を解決するために、本発明のある観点によれば、目的音および雑音が混入している入力音声の前記目的音を強調して音声周波数成分を取得する目的音強調部と、前記入力音声の前記目的音を抑圧して雑音周波数成分を取得する目的音抑圧部と、前記音声周波数成分および前記雑音周波数成分に応じた所定のゲイン関数を用いて前記音声周波数成分に乗算するゲイン値を算出するゲイン算出部と、前記ゲイン算出部により算出されたゲイン値を前記音声周波数成分に乗算するゲイン乗算部と、を備え、前記ゲイン算出部は、前記音声周波数成分と前記雑音周波数成分とのエネルギー比が所定値以下の場合に前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなる前記ゲイン関数を用いて前記ゲイン値を算出する、音声処理装置が提供される。
【０００６】
また、前記音声周波数成分には目的音成分と雑音成分が含まれており、前記ゲイン乗算部は、前記音声周波数成分に前記ゲイン値を乗算して前記音声周波数成分に含まれている前記雑音成分を抑圧してもよい。
【０００７】
また、前記ゲイン算出部は、前記目的音抑圧部により取得された雑音周波数成分に雑音のみが含まれていると推定して、前記ゲイン値を算出してもよい。
【０００８】
また、前記ゲイン関数は、前記音声周波数成分と前記雑音周波数成分とのエネルギー比において、雑音の比率が集中している雑音集中範囲の前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなるゲインカーブを有する関数でもよい。
【０００９】
また、前記ゲイン関数は、前記雑音集中範囲以外で最も急斜である前記ゲイン関数の傾きよりも小さい傾きのゲインカーブを有する関数でもよい。
【００１０】
また、前記入力音声に含まれる前記目的音が存在する区間を検出する目的音区間検出部を備え、前記ゲイン算出部は、前記目的音区間検出部による検出結果に応じて、前記目的音強調部により取得された前記音声周波数成分のパワースペクトルおよび前記目的音抑圧部により取得された前記雑音周波数成分のパワースペクトルを平均化してもよい。
【００１１】
また、前記ゲイン算出部は、前記目的音区間検出部による検出の結果、目的音が存在する区間であることが検出された場合に第１の平滑化係数を選択し、前記目的音が存在する区間であることが検出された場合に第２の平滑化係数を選択して、前記音声周波数成分および前記雑音周波数成分のパワースペクトルを平均化してもよい。
【００１２】
また、前記ゲイン算出部は、平均化された前記音声周波数成分のパワースペクトルおよび前記雑音周波数成分のパワースペクトルを用いて算出されたゲイン値を平均化してもよい。
【００１３】
また、前記目的音抑圧部により取得された雑音周波数成分の大きさを、前記目的音強調部により取得された音声周波数成分に含まれる雑音成分の大きさに対応させるように前記雑音周波数成分を補正する雑音補正部を備え、前記ゲイン算出部は、前記雑音補正部により補正された前記雑音周波数成分に応じたゲイン値を算出してもよい。
【００１４】
また、前記雑音補正部は、ユーザ操作に応じて前記雑音周波数成分を補正してもよい。
【００１５】
また、前記雑音補正部は、検出された雑音の状態に応じて前記雑音周波数成分を補正してもよい。
【００１６】
また、上記課題を解決するために、本発明の別の観点によれば、目的音および雑音が混入している入力音声の前記目的音を強調して音声周波数成分を取得するステップと、前記入力音声の前記目的音を抑圧して雑音周波数成分を取得するステップと、前記音声周波数成分と前記雑音周波数成分とのエネルギー比が所定値以下の場合に前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなるゲイン関数を用いて前記音声周波数成分に乗算するゲイン値を算出するステップと、前記ゲイン算出部により算出されたゲイン値を前記音声周波数成分に乗算するステップと、を含む、音声処理方法が提供される。
【００１７】
また、上記課題を解決するために、本発明の別の観点によれば、コンピュータを、目的音および雑音が混入している入力音声の前記目的音を強調して音声周波数成分を取得する目的音強調部と、前記入力音声の前記目的音を抑圧して雑音周波数成分を取得する目的音抑圧部と、前記音声周波数成分および前記雑音周波数成分に応じた所定のゲイン関数を用いて前記音声周波数成分に乗算するゲイン値を算出するゲイン算出部と、前記ゲイン算出部により算出されたゲイン値を前記音声周波数成分に乗算するゲイン乗算部と、を備え、前記ゲイン算出部は、前記音声周波数成分と前記雑音周波数成分とのエネルギー比が所定値以下の場合に前記ゲイン値および前記ゲイン関数の傾きが所定値より小さくなる前記ゲイン関数を用いて前記ゲイン値を算出する、音声処理装置として機能させるためのプログラムが提供される。
【発明の効果】
【００１８】
以上説明したように本発明によれば、所定のゲイン関数を利用して、ミュージカルノイズが低減された音声強調を行うことができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の実施形態の概要を説明する説明図である。
【図２】本発明の実施形態の概要を説明する説明図である。
【図３】本発明の第１の実施形態にかかる音声処理装置の機能構成を示すブロック図である。
【図４】同実施形態にかかるゲイン算出部の機能構成を示すブロック図である。
【図５】同実施形態にかかるゲイン算出部による平均化処理を示すフローチャートである。
【図６】同実施形態にかかる目的音区間検出部の機能構成を示すブロック図である。
【図７】同実施形態にかかる目的音の検出処理について説明する説明図である。
【図８】同実施形態にかかる目的音の検出処理について説明する説明図である。
【図９】同実施形態にかかる目的音区間の検出処理を示すフローチャートである。
【図１０】同実施形態にかかる目的音の検出処理について説明する説明図である。
【図１１】同実施形態にかかる白色化について説明する説明図である。
【図１２】同実施形態にかかる雑音補正部の機能構成を示すブロック図である。
【図１３】同実施形態にかかる雑音補正の処理を示すフローチャートである。
【図１４】同実施形態にかかる雑音補正部の機能構成を示すブロック図である。
【図１５】同実施形態にかかる雑音補正の処理を示すフローチャートである。
【図１６】同実施形態にかかる音声処理装置の機能構成を示すブロック図である。
【図１７】同実施形態にかかる定式化による出力信号の差を説明する説明図である。
【図１８】本発明の第２の実施形態にかかる機能構成を示すブロック図である。
【図１９】同実施形態にかかる目的音強調前後の雑音スペクトルを説明する説明図である。
【図２０】同実施形態にかかる目的音強調前後の目的音スペクトルを説明する説明図である。
【図２１】従来の技術を説明する説明図である。
【図２２】従来の技術を説明する説明図である。
【発明を実施するための形態】
【００２０】
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００２１】
また、以下に示す順序に従って、当該「発明を実施するための形態」を説明する。
１．本実施形態の目的
２．第１実施形態
３．第２実施形態
【００２２】
＜１．本実施形態の目的＞
まず、本実施形態の目的について説明する。従来から、雑音が混入している入力音声に対して、雑音を抑圧して目的音声を強調することが行われている（例えば、上記特許文献１〜３）。特許文献１では、複数個のマイクを利用して、目的音声を強調した信号（以降、音声周波数成分と称する。）と、目的音声を抑圧した信号（以降、雑音周波数成分と称する。）が取得される。
【００２３】
そして、音声周波数成分には目的音声と雑音が含まれ、雑音周波数成分には雑音のみが含まれていると推定して、両者を利用してスペクトルサブトラクションが行われる。特許文献１におけるスペクトルサブトラクション処理においては、処理後の音声信号にミュージカルノイズといわれる特有の歪みが生じてしまうという問題があった。また、音声周波数成分に含まれる雑音と雑音周波数成分に含まれる雑音は等しいと仮定して処理しているが、実際には等しくない場合があるという問題があった。
【００２４】
ここで、一般的なスペクトルサブトラクションの処理について説明する。一般に、スペクトルサブトラクションでは、信号に含まれる雑音成分を推定し、パワースペクトル上で引き算が行われる。以下では、音声周波数成分Ｘに含まれる目的音成分をＳ、雑音成分をＮ、雑音周波数成分をＮ′とする。処理後周波数成分Ｙのパワースペクトルは以下の式により得られる。
【００２５】
【数１】
2011191669.tif 000003
一般には、入力信号の位相を利用して復元するので、以下のように引き算であってもＸにある値（以下、ゲイン値）を乗じることにより雑音成分を抑圧することができる。
【００２６】
【数２】
2011191669.tif 000004
【００２７】
Ｗs(h)をＸとＮ′の比ｈの関数とみなすと、その外形は図２１に示した外形となる。ｈ＜１の範囲はフロアリングといわれ、一般には、Ｗs(h)＝０．０５など適当な小さい値に置き換えられる。図２１に示したように、Ｗs(h)の外形は、ｈが小さいところで非常に大きな傾きを持っている。したがって、ｈが、ｈの小さい範囲（例えば１＜ｈ＜２）で少し振動すると、その結果得られるゲイン値が大きく振動することとなる。これにより、周波数成分に対して、時間−周波数ごとに変号の大きな値が乗じられることになり、いわゆるミュージカルノイズが生じると考えられる。
【００２８】
ｈが小さい値をとる場合とは、音声周波数成分Ｘにおいて、Ｓが非常に小さい場合もしくは、Ｓ＝０となる非音声区間であり、この区間での音質の劣化が著しくなる。また、Ｎ＝Ｎ′と仮定しているが、この仮定が正しくない場合に、特に非音声区間でゲイン値が大きく振動し、音質劣化の要因となる。
【００２９】
また、上記した特許文献３では、音声周波数成分（Ｘ＝Ｓ＋Ｎ）と雑音周波数成分Ｎ′に対して、出力の適応において音声周波数成分に含まれる雑音成分Ｎと雑音周波数成分Ｎ′の大きさをそろえている。しかし、ポストフィルタリング手段でＭＡＰ最適化などを行っているものの、Ｗｉｅｎｅｒ  Ｆｉｌｔｅｒに基づいた手法になっており、出力の適応の効果を十分に活かすことができない。
【００３０】
Ｗｉｅｎｅｒ  Ｆｉｌｔｅｒは、目的音成分Ｓと雑音成分Ｎに対して、以下で与えられる値を音声周波数成分に乗じることにより雑音の抑圧をおこなう。
【００３１】
【数３】
2011191669.tif 000005
実際にはＳとＮは観測できないため、観測可能な音声周波数成分Ｘと雑音周波数成分Ｎ′を利用し、以下のように求める。
【００３２】
【数４】
2011191669.tif 000006
【００３３】
これを、前述のスペクトルサブストラクションと同様にｈの関数と考えると、その外形は図２２に示した外形となる。図２１のスペクトルサブトラクションと同様に、ｈの値が小さい範囲において、Ｗ(h)の傾きが大きくなっている。出力の適応により、非音声区間では、ｈ自体の散らばりは小さくなり（１の付近に集まる）、従来と比べると、乗じるゲイン値の変動を小さく抑えることが可能となっている。しかし、傾き自体が大きいところにｈの値が集中するのは望ましくない。
【００３４】
そこで、上記のような事情を一着眼点として、本実施形態にかかる音声処理装置が創作されるに至った。本実施形態にかかる音声処理装置によれば、所定のゲイン関数を利用して、ミュージカルノイズが低減された音声強調を行うことが可能となる。
【００３５】
＜２．第１実施形態＞
次に、第１実施形態について説明する。図１および図２を参照して、第１実施形態の概要について説明する。第１実施形態では、雑音抑圧に利用するゲイン関数Ｇ(r)が以下の特徴を有する。
（１）ｒが小さい値の範囲Ｒ１（例えばｒ＜２）では、なるべく小さな値かつ、小さな傾きを有する。
（２）ｒが中程度の範囲Ｒ２（例えば２＜ｒ６）では、大きな正の傾きを有する。
（３）ｒが十分大きい範囲Ｒ３（例えばｒ≧６）では、傾きは小さくなり、1に収束する。
（４）Ｇ(r)は変曲点に対して非対称。
【００３６】
図１のグラフ３００は、上記（１）〜（４）の条件を満たす関数Ｇ(r)の外形を示している。図２は、実際に観測されたデータにおいて、雑音のみが存在する区間でのｈの値の分布をグラフ化したものである。ヒストグラム３０１に示したように、実際に観測されたデータにおいて、雑音のみが存在する区間でのｈの値のほとんど（８０％）は、０〜２に集中している。したがって、上記（１）の条件におけるｒが小さい範囲とは、雑音のみが存在する区間において、雑音の比率（ｈ）のヒストグラムを算出したときに、８０％のデータが含まれる範囲とすることができる。以下では、ｒ＜２の範囲Ｒ１において、なるべく小さな値かつ、小さな傾きを有するゲイン関数Ｇ(r)を用いて雑音抑圧をおこなっている。
【００３７】
また、本実施形態では、目的音区間か否かを検出して、時間方向のパワースペクトルの平均化を行う。例えば、目的音が存在しない区間で大きく平均化することにより、時間方向の分散を小さくする。これにより、上記したゲイン関数によりｒが小さい範囲Ｒ１において変動が少ない値を出力し、かつ、時間方向にも変動の少ない値を得ることが可能となり、更に、ミュージカルノイズを低減することができる。
【００３８】
また、本実施形態では、音声周波数成分に含まれる雑音成分Ｎと、雑音周波数成分Ｎ′の比がＧ(r)のＲ１の範囲に収まるように、周波数特性の補正を行う。これにより、さらに、ゲイン値の算出において、ｈを小さくし、さらに分散を小さくすることが可能となり、大きな雑音抑圧および大幅なミュージカルノイズの低減を実現することができる。
【００３９】
次に、図３を参照して、音声処理装置１００の機能構成について説明する。図３は、音声処理装置１００の機能構成を示すブロック図である。音声処理装置１００は、目的音強調部１０２と、目的音抑圧部１０４と、ゲイン算出部１０６と、ゲイン乗算部１０８と、目的音区間検出部１１０と、雑音補正部１１２などを備える。
【００４０】
目的音強調部１０２は、目的音および雑音が混入している入力音声の目的音を強調して、音声周波数成分Ｙempを取得する機能を有する。本実施形態では、複数本のマイクロホンから音声Ｘｉが入力されるが、かかる例に限定されず、一本のマイクロホンから音声Ｘｉが入力されてもよい。目的音強調部により取得された音声周波数成分Ｙempは、ゲイン算出部１０６、ゲイン乗算部１０８、目的音区間検出部１１０に提供される。
【００４１】
目的音抑圧部１０４は、目的音および雑音が混入している入力音声の目的音を抑圧して、雑音周波数成分Ｙsupを取得する機能を有する。目的音抑圧部１０４により目的音を抑圧して、雑音成分が推定される。目的音抑圧部１０４により取得された雑音周波数成分Ｙsupは、ゲイン算出部１０６、目的音区間検出部１１０、雑音補正部１１２に提供される。
【００４２】
ゲイン算出部１０６は、目的音強調部１０２により取得された音声周波数成分および目的音抑圧部１０４により取得された雑音周波数成分に応じた所定のゲイン関数を用いて音声周波数成分に乗算するゲイン値を算出する機能を有する。所定のゲイン関数とは、図１に示したように、音声周波数成分と雑音周波数成分とのエネルギー比が所定値以下の場合にゲイン値およびゲイン関数の傾きが所定値より小さくなるゲイン関数である。
【００４３】
ゲイン乗算部１０８は、ゲイン算出部１０６により算出されたゲイン値を目的音強調部１０２により取得された音声周波数成分に乗算する機能を有する。音声周波数成分に図１に示したゲイン関数を用いたゲイン値が乗算されることにより、ミュージカルノイズを低減して、雑音を抑圧することが可能となる。
【００４４】
目的音区間検出部１１０は、入力音声に含まれる目的音が存在する区間を検出する機能を有する。目的音区間検出部１１０は、目的音強調部１０２により提供される周波数スペクトルＹempと、目的音抑圧部１０４から得られる周波数スペクトルＹsupから振幅スペクトルを算出して、それぞれ入力音声Ｘｉとの相関を求めることにより目的音の区間を検出する。目的音区間検出部１１０による目的音の検出処理については後で詳細に説明する。
【００４５】
ゲイン算出部１０６は、目的音区間検出部１１０による検出結果に応じて、目的音強調部１０２により取得された音声周波数成分のパワースペクトルおよび目的音抑圧部１０４により取得されたパワースペクトルを平均化する。ここで、図４を参照して、目的音区間検出部１１０による検出結果に応じたゲイン算出部１０６の機能について説明する。
【００４６】
図４に示したように、ゲイン算出部１０６は、演算手段１２２と、第１平均化手段１２４と、第１保持手段１２６と、ゲイン算出手段１２８と、第２平均化手段１３０、第２保持手段１３２などを有する。演算手段１２２は、目的音強調部１０２により取得された周波数スペクトルＹempおよび目的音抑圧部１０４により取得された周波数スペクトルＹsupに対して、パワースペクトルを算出する機能を有する。
【００４７】
そして、第１平均化手段１２４は、目的音区間検出部１１０により検出された目的音区間を示す制御信号に応じて、パワースペクトルの平均化を行う。第１平均化手段１２４では、例えば、一次の減衰を利用して、目的音区間検出部１１０の検出結果に応じて、パワースペクトルの平均化を行う。目的音が存在する区間では、以下の数式によりパワースペクトルの平均化を行う。
【００４８】
【数５】
2011191669.tif 000007
【００４９】
また、目的音が存在しない区間では、以下の数式によりパワースペクトルの平均化を行う。
【００５０】
【数６】
2011191669.tif 000008
【００５１】
上記では、ｒ１＜ｒ２で、例えば、ｒ１＝０．３、ｒ２＝０．９などの値を利用する。また、ｒ３は、例えば、ｒ２と同程度の値を利用することが望ましい。また、目的音の存在に応じて、ｒ１、ｒ２を切り替えるのではなく、連続的に変化させてもよい。ｒ１およびｒ２を連続的に変化させる方法については後で詳細に説明する。また、上記では１次の減衰を利用した平滑化を行っているが、かかる例に限定されない。例えば、Ｎフレームを平均して、そのＮをｒ同様に制御してもよい。すなわち、目的音が存在するときは、過去３フレームの平均値を利用し、目的音が存在しないときは過去７フレームの平均値を利用するなどの制御を行う。
【００５２】
上記では、目的音が存在しない区間で、なるべく、ＰｘおよびＰｎを大きく平均化することにより、時間方向の分散を小さくすることができる。本実施形態にかかるゲイン関数では、図１に示したように、ｒが小さい範囲（Ｒ１）において、変動が少ない値を出力することができる。つまり、ゲイン関数Ｇ（ｒ）を利用することにより、ｒが小さい範囲においてミュージカルノイズを発生しにくくしているが、パワースペクトルの平均化により時間方向にも変動が少ない値を得ることが可能となる。これにより、ミュージカルノイズを更に低減することが可能となる。一方、目的音が存在する区間で大きな平均化を行うとエコー感の原因となるため、目的音の有無に応じて平滑化係数ｒの制御を行う。
【００５３】
ゲイン算出手段１２８は、ｈ＝Ｐｘ／Ｐｎに応じて、図１に示した外形を持つ値を算出する。このとき、あらかじめ保持したテーブルの値を利用してもよいし、図１の外形を持つ以下の関数を利用してもよい。
【００５４】
【数７】
2011191669.tif 000009
例えば、ｂ＝０．８、ｃ＝０．４とする。
【００５５】
第２平均化手段１３０は、第１平均化手段１２４と同様の平均化処理をゲイン値に対して行う。平均化の係数は、ｒ１、ｒ２、ｒ３と同じ値でもよいし、異なる値であってもよい。次に、図５を参照して、ゲイン算出部１０６による平均化処理について説明する。図５は、ゲイン算出部１０６による平均化処理を示すフローチャートである。
【００５６】
図５に示したように、まず、目的音強調部１０２および目的音抑圧部１０４から周波数スペクトル（Ｙemp、Ｙsup）を取得する（Ｓ１０２）。そして、パワースペクトル（Ｙ２emp、Ｙ２sup）を算出する（Ｓ１０４）。そして、第１保持手段１２６から、過去の平均化されたパワースペクトル（Ｐｘ、Ｐｎ）を取得する（Ｓ１０６）。そして、目的音が存在する区間であるか否かを判定する（Ｓ１０８）。
【００５７】
ステップＳ１０８において、目的音が存在する区間であると判定された場合には、平滑化係数にｒ＝ｒ１を選択する（Ｓ１１０）。ステップＳ１０８において、目的音が存在しない区間であると判定された場合には、平滑化係数にｒ＝ｒ２を選択する。そして、以下の数式によりパワースペクトルの平均化を行う（Ｓ１１４）。
【００５８】
【数８】
2011191669.tif 000010
【００５９】
そして、Ｐｘ、Ｐｎを利用して、ゲイン値ｇを算出する（Ｓ１１６）。そして、第２保持手段１３２から、過去のゲイン値Ｇを取得する（Ｓ１１８）。ステップＳ１１８において取得したゲイン値Ｇを以下の数式により平均化する。
【００６０】
【数９】
2011191669.tif 000011
【００６１】
ステップＳ１２０において平均化されたゲイン値Ｇをゲイン乗算部１０８へ送る（Ｓ１２２）。そして、第１保持手段１２６にＰｘおよびＰｎを保持し（Ｓ１２４）、第２保持手段にゲイン値Ｇを保持する（Ｓ１２６）。上記処理は、すべての周波数域に対して実行される。また、上記処理では、パワースペクトルの平均化とゲインの平均化において、同じ平均化係数を用いているが、かかる例に限定されず、それぞれ異なる平均化係数を用いてもよい。
【００６２】
次に、図６を参照して、目的音区間検出部１１０による目的音の検出処理について説明する。図６に示したように、目的音区間検出部１１０は、演算手段１３２と、相関算出手段１３４と、比較手段１３６と、判定手段１３８などを有する。
【００６３】
演算手段１３２には、目的音強調部１０２から提供される周波数スペクトルＹempと、目的音抑圧部１０４から提供される周波数スペクトルＹsupと、入力信号のうち一つの周波数スペクトルＸiが入力される。周波数スペクトルＸiの選択については、どのマイクロホンを選択してもよいが、目的音が入力される位置が予めわかっている場合には、目的音に最も近い位置のマイクロホンを利用することが望ましい。これにより、最も大きな音で目的音を入力することができる。
【００６４】
演算手段１３２は、入力された各周波数スペクトルに対して、振幅スペクトルもしくは、パワースペクトルを算出する。そして、相関算出手段１３４は、ＹempとＸiの振幅スペクトルの相関Ｃ１と、ＹsupとＸiの振幅スペクトルの相関Ｃ２を求める。比較手段１３６は、相関算出手段１３４により算出された相関Ｃ１と相関Ｃ２とを比較する。判定手段１３８は、比較手段１３６による比較結果に応じて、目的音が存在するか否かを判定する。
【００６５】
判定手段１３８は、振幅スペクトルの相関から以下の手法により目的音が存在するか否かを判定する。まず、演算手段１３２に入力される信号に含まれる成分を以下に示す。
目的音強調部１０２から得られる周波数スペクトルＹemp：目的音声＋抑圧された雑音成分
目的音抑圧部１０４から得られる周波数スペクトルＹsup：雑音成分
入力信号のうち一つの周波数スペクトルＸi：目的音声＋抑圧された雑音成分
【００６６】
振幅スペクトルの相関は、二つのスペクトルが似ているときに大きな値をとる。図７のグラフ３１０に示したように、目的音が存在する区間では、Ｘiの形状は、ＹsupよりもＹempに似ているスペクトルとなることがわかる。また、図７のグラフ３１２に示したように、目的音が存在しない区間では、雑音のみとなる。このため、Ｘiの形状は、ＹsupとＹempでは同程度となり、明確な差のないスペクトルとなることがわかる。
【００６７】
よって、ＸiとＹempの相関値Ｃ１は、ＸiとＹsupの相関値Ｃ２に比べて、目的音が存在する区間では大きくなる。また、目的音が存在しない区間では、Ｃ１とＣ２は同程度の値となる。図８のグラフ３１４に示したように、相関値Ｃ１から相関値Ｃ２を減算した値は、実際の目的音の存在区間と同程度の値となっていることがわかる。このように、振幅スペクトルの相関を比較することにより、目的音が存在する区間と目的音が存在しない区間とを区別することが可能となる。
【００６８】
次に、図９を参照して、目的音区間検出部１１０による目的音区間の検出処理について説明する。図９は、目的音区間検出部１１０による目的音区間の検出処理を示すフローチャートである。図９に示したように、まず、目的音強調部１０２から周波数スペクトルＹemp、目的音抑圧部１０４から周波数スペクトルＹsup、マイクロホンの入力から周波数スペクトルＸiを取得する（Ｓ１３２）。
【００６９】
ステップＳ１３２において取得した周波数スペクトルから振幅スペクトルを算出する（Ｓ１３４）。そして、ＸiとＹempの振幅スペクトルの相関Ｃ１、ＸiとＹsupの振幅スペクトルの相関Ｃ２を算出する（Ｓ１３６）。そして、相関Ｃ１から相関Ｃ２を減算した値（Ｃ１−Ｃ２）がＸiの閾値Ｔｈより大きいかを判定する（Ｓ１３８）。
【００７０】
ステップＳ１３８において、ＴｈよりＣ１−Ｃ２が大きいと判定された場合には、目的音が存在すると判断する（Ｓ１４０）。ステップＳ１３８において、ＴｈよりＣ１−Ｃ２が小さいと判定された場合には、目的音が存在しないと判断する（Ｓ１４２）。以上、目的音区間検出部１１０による目的音区間の検出処理について説明した。
【００７１】
次に、目的音区間検出部１１０が、数式により目的音区間を算出する場合について説明する。まず、各振幅スペクトルを以下のように定義する。
【００７２】
【数１０】
2011191669.tif 000012
【００７３】
Ａxiの平均値を用いて、以下の白色化を行う。
【００７４】
【数１１】
2011191669.tif 000013
【００７５】
そして、ＡＷxiとの相関を取る。ここで、ｐ(k)は周波数ごとの重みである。
【００７６】
【数１２】
2011191669.tif 000014
【００７７】
上記した重みｐ(k)は、例えば、図１０の関数３１６で示される。音声は主として低域に強いエネルギーが集中し、雑音は広い帯域に渡ってエネルギーが存在する。このため、主として音声の強い帯域のみを利用することで精度を上げることが可能となる。例えば、Ｎ＝５１２（ＦＦＴサイズ）に対して、Ｎｏ＝４０、Ｌ＝３などを利用することができる。
【００７８】
ここで、図１１を参照して、上記した白色化について説明する。図１１のグラフ３１８に示したように、振幅スペクトルは正の値しかもたない。このため、相関値も正の値しかもたず、値のレンジが小さくなってしまう。実際には０．６〜１．０程度のレンジとなる。そこで、基準となる直流成分を減算することにより、正・負両方の値をとるようにする操作を行っている。この操作を本実施形態では白色化と呼んでいる。このように、白色化することにより、相関値についても、−１〜１のレンジの値をもつことが可能となる。これにより、目的音検出の精度を上げることが可能となる。
【００７９】
また、上記で平滑化係数ｒ１およびｒ２は連続的に変化させてもよいとしたが、以下ではｒ１およびｒ２を連続的に切り替える場合について説明する。以下では、目的音区間検出部１１０により算出されるＣ１、Ｃ２および閾値Ｔｈを利用する。これらの値を利用して、以下の数式により１以下の値を算出する。例えば、β＝１または２とする。ｍｉｎは二つのｔの値のうち小さいほうを選択する関数である。
【００８０】
【数１３】
2011191669.tif 000015
【００８１】
上記数式において、ｖは目的音が存在するときに１に近い値をとる。このことを利用して、平滑化係数を連続的に以下のように求めることができる。目的音が存在するときには、ｒ≒ｒ１で、それ以外ではｒ≒ｒ２と制御される。
【００８２】
【数１４】
2011191669.tif 000016
【００８３】
図３に戻り、音声処理装置１００の機能構成の説明を続ける。雑音補正部１１２は、目的音抑圧部１０４により取得された雑音周波数成分の大きさを、目的音強調部１０２により取得された音声周波数成分に含まれる雑音成分の大きさに対応させるように雑音周波数成分を補正する機能を有する。これにより、ゲイン算出部１０６によるゲイン値の算出において、ｈを小さくし、さらに分散を小さくすることが可能となり、大きな雑音抑圧および大幅なミュージカルノイズの低減を実現することができる。
【００８４】
まず、雑音補正部１１２による雑音補正の考え方について説明する。以下の処理は各周波数成分に同様に施されるが、説明を容易にするため、周波数インデックスは省略して記載する。
目的となる音源のスペクトルをＳとし、目的音源からマイクロホンまでの伝達特性をＡとし、各マイクロホンに観測される雑音成分をＮとする。このとき、マイクロホンに観測される信号Ｘは、以下のように記載することができる。Ｍはマイクロホン数である。
【００８５】
【数１５】
2011191669.tif 000017
【００８６】
目的音強調部１０２および目的音抑圧部１０４は、それぞれＸに対してある重みをかけて足す処理をおこなっているため、各部の出力信号は以下のように与えられる。Ｘに対してかけられる重みの作り方により、目的音を小さくしたり大きくしたりすることができる。
【００８７】
【数１６】
2011191669.tif 000018
【００８８】
したがって、Ｗemp、Ｗsupが一致しない限り、目的音強調部１０２の出力に含まれる雑音成分と、目的音抑圧部１０４の出力は異なる。具体的には、パワースペクトル上で雑音抑圧をおこなうため、各周波数ごとに雑音の大きさのレベルが一致しないこととなる。そこで、Ｗemp、Ｗsupを補正することにより、ゲイン値算出におけるｈの値を１に近づけることが可能となる。すなわち、ゲイン値において小さい値かつ傾きの小さいところに値を集中することができる。ｈは以下の数式により表される。
【００８９】
【数１７】
2011191669.tif 000019
【００９０】
例えば、
【００９１】
【数１８】
2011191669.tif 000020
の場合は、補正を行うことにより、ｈは１より大きい値から１に近づく。よって、雑音抑圧量を向上することができる。また、
【００９２】
【数１９】
2011191669.tif 000021
の場合は、補正を行うことにより、ｈは１より小さい値から１に近づく。よって、音声の劣化を低減することができる。
【００９３】
ｈが１付近の小さい値に集中すると、ゲイン関数の最小値を小さくすることができる。これにより、雑音抑圧量の向上に寄与することが可能となる。Ｗemp、Ｗsupは既知の値であるため、雑音スペクトルＮの共分散Ｒnがわかれば、以下の数式により雑音補正を行うことができる。
【００９４】
【数２０】
2011191669.tif 000022
【００９５】
次に、図１２を参照して、雑音補正部１１２による雑音補正処理について説明する。図１２に示したように、雑音補正部１１２は、演算手段１４０と保持手段１４２などを有する。演算手段１３０には、目的音抑圧部１０４により取得された周波数スペクトルＹsupが入力される。そして、保持手段１４２を参照し補正係数を算出して、入力された周波数スペクトルＹsupに乗じて雑音スペクトルＹcompを算出する。算出されたＹcompは、ゲイン算出部１０６に提供される。保持手段１４２には、雑音の共分散、目的音強調部１０２および目的音抑圧部１０４で用いられる係数が保持されている。
【００９６】
次に、図１３を参照して、雑音補正部１１２による雑音補正の処理について説明する。図１３は、雑音補正部１１２による雑音補正の処理を示すフローチャートである。図１３に示したように、まず、目的音抑圧部１０４から周波数スペクトルＹsupを取得する（Ｓ１４２）。そして、保持手段１４２から共分散、目的音強調の係数、目的音抑圧の係数を取得する（Ｓ１４４）。そして、周波数毎に補正係数Ｇcompを算出する（Ｓ１４６）。
【００９７】
そして、周波数毎に周波数スペクトルにステップＳ１４６において算出された補正係数Ｇcompを乗じる（Ｓ１４８）。
【００９８】
【数２１】
2011191669.tif 000023
【００９９】
そして、ゲイン算出部１０６にステップＳ１４８における算出結果Ｙcompを送る（Ｓ１５０）。雑音補正部１１２による上記処理は、すべての周波数域に対して繰り返し実行される。
【０１００】
上記した雑音の共分散Ｒnは、例えば、以下の数式により算出することができる（参照：Measurement of
Correlation Coefficients in Reverberant Sound Fields, Richard K. Cook et al THE
JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, VOLUME 26, NUMBER 6, NOVEMBER
1955）。
【０１０１】
【数２２】
2011191669.tif 000024
【０１０２】
直線状に並んだマイクに対して、拡散雑音場を仮定すると、
【０１０３】
【数２３】
2011191669.tif 000025
【０１０４】
直線状に並んだマイクに対して、全方位から互いに相関のない雑音が到来する場を仮定すると、
【０１０５】
【数２４】
2011191669.tif 000026
【０１０６】
また、雑音の共分散Ｒnは、数式により算出する以外にも、例えば、あらかじめ大量のデータを収録して、その平均値を求めることにより得ることができる。この場合、マイクロホンに観測されるのは雑音のみとなるため、以下の数式により雑音の共分散を得ることができる。
【０１０７】
【数２５】
2011191669.tif 000027
【０１０８】
また、目的音強調部１０２、上述した伝達特性Ａ、共分散Ｒnを用いて以下のような係数を構築することができる。一般的に、最尤ビームフォーミングと呼ばれる（参照：アダプティブアンテナ技術  菊間信良著  オーム社）。
【０１０９】
【数２６】
2011191669.tif 000028
【０１１０】
また、最尤ビームフォーミング手法に限定されず、遅延和ビームフォーミングと呼ばれる手法を用いてもよい。この場合、上記において、Ｒnが単位行列であることと同義となる。また、目的音抑圧部１０４では、上記したＡとＡ以外の伝達特性を利用して以下のような係数が構築される。以下は、目的音とは別の方向に対して１、目的音の方向の信号をゼロとするような係数となる。
【０１１１】
【数２７】
2011191669.tif 000029
【０１１２】
また、雑音補正部１１２は、制御部（図示せず）からの選択信号に基づいて、補正係数を変更するようにしてもよい。例えば、図１４に示したように、雑音補正部１１２は、演算手段１５０と、選択手段１５２と、複数の保持手段（第１保持手段１５４、第２保持手段１５６、第３保持手段１５８）を有してもよい。複数の保持手段には、それぞれ異なる補正係数が保持されている。選択手段１５２は、制御部から提供される選択信号に基づいて、第１保持手段１５４、第２保持手段１５６、第３保持手段１５８に保持されている補正係数のいずれかの補正係数を取得する。
【０１１３】
制御部は、例えば、ユーザ入力に応じて動作するか、雑音の状態に応じて動作して、雑音補正部の選択手段１５２に選択信号を提供する。そして、演算手段１５０は、選択手段１５２により選択された補正係数を用いて、入力された周波数スペクトルＹsupに当該補正係数を乗じて雑音スペクトルＹcompを算出する。
【０１１４】
次に、図１５を参照して、選択信号に基づいて補正係数を取得する場合の雑音補正処理について説明する。図１５に示したように、まず、目的音抑圧部１０４から周波数スペクトルＹsupを取得する（Ｓ１５２）。そして、制御部から選択信号を取得する（Ｓ１５４）。そして、取得した選択信号の値が現在の値と異なっているか否かを判定する（Ｓ１５６）。
【０１１５】
ステップＳ１５６において、取得した値が現在の値と異なっていると判定された場合には、取得した選択信号の値を利用して、選択信号の値に対応する保持手段からデータを取得する（Ｓ１５８）。そして、周波数毎に補正係数Ｇcompを算出する（Ｓ１６０）。そして、以下の数式により、周波数毎に周波数スペクトルに補正係数を乗じる（Ｓ１６２）。
【０１１６】
【数２８】
2011191669.tif 000030
【０１１７】
ステップＳ１５６において、取得した値が現在の値と同じであると判定された場合には、ステップＳ１６２の処理を実行する。そして、ゲイン算出部１０６にステップＳ１６２における算出結果Ｙcompを送る（Ｓ１６４）。雑音補正部１１２による上記処理は、すべての周波数域に対して繰り返し実行される。
【０１１８】
また、図１６に示したように、音声処理装置２００のように、雑音補正部２０２が目的音区間検出部１１０の検出結果を利用して雑音共分散の算出を行ってもよい。雑音補正部２０２は、目的音抑圧部１０４から出力された周波数スペクトルＹsupだけでなく、目的音強調部１０２から出力された周波数スペクトルＹempおよび目的音区間検出部１１０により検出された検出結果を利用して、雑音補正を行う。
【０１１９】
以上、第１実施形態について説明した。第１実施形態によれば、図１の特徴を有するゲイン関数Ｇ（ｒ）を利用して雑音を抑圧することができる。すなわち、音声の周波数成分と雑音の周波数成分のエネルギー比に応じたゲイン値を音声の周波数成分に乗算して適切に雑音を抑圧することができる。
【０１２０】
また、目的音区間か否かを検出し、スペクトル時間方向の平均化制御を行うことにより、時間方向の分散を小さくして、時間方向に変動が少ない値を得ることが可能となり、ミュージカルノイズの発生を更に低減することが可能となる。また、音声周波数成分に含まれる雑音成分Ｎと、雑音周波数成分Ｎ′の比がＧ(r)のＲ１の範囲に収まるように、周波数特性の補正を行う。これにより、さらに、ゲイン値の算出において、ｈを小さくしさらに分散を小さくすることが可能となり、大きな雑音抑圧および大幅なミュージカルノイズの低減を実現することができる。
【０１２１】
本実施形態にかかる音声処理装置１００または２００は、携帯電話やＢｌｕｅｔｏｏｔｈのヘッドセットや、コールセンターやＷｅｂ会議に用いられるヘッドセット、ＩＣレコーダやビデオ会議システム、ノートＰＣの本体に付加されたマイクを用いたＷｅｂ会議やボイスチャットに利用することができる。
【０１２２】
＜３．第２実施形態＞
次に、第２実施形態について説明する。第１実施形態では、ゲイン関数を利用して、大きな雑音抑圧を実現しつつ、ミュージカルノイズを低減する方法について説明した。以下では、複数マイクロホンを利用することにより、スペクトルサブストラクション（以降、ＳＳとも称する）を利用して、非常に簡易にミュージカルノイズを低減し、目的音声を強調する方法について説明する。ＳＳベースの場合、以下の数式が成立する。
【０１２３】
【数２９】
2011191669.tif 000031
ＳＳの定式化として、フロアリングの行い方によって２通りの記述が可能である。
【０１２４】
＜定式化１＞
【０１２５】
【数３０】
2011191669.tif 000032
【０１２６】
＜定式化２＞
【数３１】
2011191669.tif 000033
【０１２７】
定式化１では、Ｇが負にならない限りはフロアリングが生じないが、定式化２では、Ｇthより小さい場合はＧthという一定の利得を掛けることが差となる。定式化１では、Ｇは非常に小さい値までとることが可能となり、雑音自体の抑圧量が大きくなる。しかし、第１実施形態で説明したように、ＳＳは、ゲインという観点から見ると、時間−周波数的に不連続な値をとる可能性が高いため、ミュージカルノイズを発生させる。
【０１２８】
また、定式化２では、Ｇth（例えば０．１）より小さな値は乗じられないため、雑音自体の抑圧量は小さい。しかし、多くの時間−周波数において、一定のＧthが乗じられることにより、ミュージカルノイズ自体の発生を抑えることが可能となる。例えば、雑音を小さくする方法として、音量を下げることが考えられる。上記現象は、例えば、ラジオに雑音が乗っているときに音量を下げると雑音は小さくなり、変な歪みを持った音が出てこないことからもわかる。すなわち、違和感の少ない音声を提供するためには、雑音抑圧を大きくするよりも、雑音の変形を一定にすることが有効であることがわかる。
【０１２９】
ここで、図１７を参照して、上記した定式化によるＳＳの出力信号の差について説明する。図１７は、定式化によるＳＳの出力信号の差を説明する説明図である。図１７のグラフ４０１は、マイクロホンから出力された音声周波数Ｘである。グラフ４０２は、定式化１により、Ｇが乗じられた場合である。この場合、レベル自体を下げることができるが、周波数の形が崩れてしまう。また、グラフ４０３は、定式化２により、Ｇが乗じられた場合である。この場合、周波数の形は保持されたまま、レベルが下がる。
【０１３０】
以上から、音声の成分はなるべくＧthより大きな値が乗算され、雑音の成分はすべてＧthの値が乗算されるようにすればよいことがわかる。
【０１３１】
【数３２】
2011191669.tif 000034
【０１３２】
一般的には、αを２程度に設定し、大きめに雑音成分を減算することで上記処理を実現する。しかし、一般的に推定した雑音成分Ｎが正しくなければ意味をなさない。
【０１３３】
また、本実施形態の第２のポイントは、複数マイクロホンを用いた処理を利用することである。上記処理に適した雑音成分を効率的に見つけ、一定の値Ｇthを乗算できるようにしたものである。図１８を参照して、本実施形態にかかる音声処理装置３００の機能構成について説明する。図１８に示したように、音声処理装置３００は、目的音強調部１０２、目的音抑圧部１０４、目的音区間検出部１１０、雑音補正部３０２、ゲイン算出部３０４などを備える。以下では、第１実施形態と異なる機能について特に詳細に説明し、第１実施形態と同様の機能については詳細な説明は省略する。
【０１３４】
第１実施形態では、雑音補正部１１２によりＹsupとＹempのパワーが等しくなるように補正が行われていた。つまり、目的音強調後の雑音パワーを推定していた。しかし、本実施形態では、ＹsupとＸiのパワーが等しくなるような補正を行う。すなわち、目的音強調前の雑音のパワーを推定する。
【０１３５】
目的音強調前の雑音を推定するには、雑音補正部３０２で算出される値
【０１３６】
【数３３】
2011191669.tif 000035
を以下の数式のように変形する。
【０１３７】
【数３４】
2011191669.tif 000036
【０１３８】
これにより、目的音強調前のマイクロホンiに含まれる雑音成分を推定することが可能となる。実際に、目的音強調後の雑音スペクトルと推定された目的音強調前の雑音スペクトルを比較すると、図１９のグラフ４１０に示したようになる。グラフ４１０に示したように、目的音強調前の雑音は、目的音強調後の雑音より大きく、特に、低域で顕著に現れている。
【０１３９】
また、実際に、目的音強調後の目的音スペクトルとマイクに入力された目的音スペクトルを比較すると、図２０のグラフ４１２に示したようになる。グラフ４１２に示したように、目的音強調後の目的音スペクトルと、マイクに入力された目的音スペクトルとを比較すると、目的音強調後と目的音強調前とで目的音成分は大きく変化していないことがわかる。
【０１４０】
以上から、ＳＳにおける雑音成分Ｎとして、目的音強調前の推定雑音を利用すると、多くの時間−周波数において、Ｇは負の値となる（ここではα＝１とした。）。なぜならば、推定雑音（Ｎ）の方が実際に含まれる雑音成分（Ｘ）より大きいからである。目的音強調とは、雑音を抑圧することであるので、目的音強調前の方が雑音自体の大きさは目的音強調後よりも大きくなっている。これは、複数マイクロホンを利用した処理によって得られるものである。
【０１４１】
また、雑音成分には一定のゲインＧthが乗算される。一方、目的音については、多少劣化があるものの、Ｇthに比べられ１に近い値が乗算される。よって、ＳＳに基づくゲイン関数を利用したとしても、ミュージカルノイズの発生の少ない音声を得ることが可能となる。このように、マイクロホンアレイ処理の特徴を生かし、目的音強調前の雑音成分を推定し、この雑音成分を利用することによりスペクトルサブストラクションベースの手法であっても、簡易にミュージカルノイズを低減して音声強調を行うことができる。
【０１４２】
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。
【０１４３】
例えば、本明細書の音声処理装置１００、２００、３００の処理における各ステップは、必ずしもフローチャートとして記載された順序に沿って時系列に処理する必要はない。すなわち、音声処理装置１００、２００、３００の処理における各ステップは、異なる処理であっても並列的に実行されてもよい。
【０１４４】
また、音声処理装置１００、２００、３００に内蔵されるＣＰＵ、ＲＯＭおよびＲＡＭなどのハードウェアを、上述した音声処理装置１００、２００、３００の各構成と同等の機能を発揮させるためのコンピュータプログラムも作成可能である。また、該コンピュータプログラムを記憶させた記憶媒体も提供される。
【符号の説明】
【０１４５】
１００、２００、３００    音声処理装置
１０２    目的音強調部
１０４    目的音抑圧部
１０６    ゲイン算出部
１０８    ゲイン乗算部
１１０    目的音区間検出部
１１２    雑音補正部
【図１】
2011191669.tif 000037
【図２】
2011191669.tif 000038
【図３】
2011191669.tif 000039
【図４】
2011191669.tif 000040
【図５】
2011191669.tif 000041
【図６】
2011191669.tif 000042
【図７】
2011191669.tif 000043
【図８】
2011191669.tif 000044
【図９】
2011191669.tif 000045
【図１０】
2011191669.tif 000046
【図１１】
2011191669.tif 000047
【図１２】
2011191669.tif 000048
【図１３】
2011191669.tif 000049
【図１４】
2011191669.tif 000050
【図１５】
2011191669.tif 000051
【図１６】
2011191669.tif 000052
【図１７】
2011191669.tif 000053
【図１８】
2011191669.tif 000054
【図１９】
2011191669.tif 000055
【図２０】
2011191669.tif 000056
【図２１】
2011191669.tif 000057
【図２２】
2011191669.tif 000058