(57)【要約】【課題】  ケプストラム解析により空間伝達関数の全域通過成分を推定することによって観測信号の残響を抑圧する残響抑圧方法とその装置を提供できるようにする。【解決手段】  観測信号の複素ケプストラムから観測信号の最小位相成分を求め、ショートパスリフタを用いて声道と声門の成分、並びに空間伝達関数の最小位相成分を抽出して、その平均値を逐次的に計算することによってガウス性白色雑音を除去すると同時に、ロングパスリフタを用いてピッチ成分を抽出して、ケフレンシ領域における観測信号から声道と声門の成分と空間伝達関数の最小位相成分の平均値とピッチ成分を減算することによって残響を抑圧させる。【選択図】図１
2015049502.tif 000002
【特許請求の範囲】
【請求項１】
音声信号と空間伝達関数の畳み込みによる残響が存在するＪ個のマイクロホンの観測信号からケプストラム解析を用いて残響を抑圧する残響抑圧方法であって、調音フィルタが全極形で近似でき、空間伝達関数が最小位相成分と全域通過成分の積で表現できることに基づき、ケフレンシ領域に変換された観測信号の最小位相成分を求め、ショートパスリフタを用いて声道と声門の成分、並びに空間伝達関数の最小位相成分を抽出して、その平均値を逐次的に計算することによってガウス性白色雑音を除去すると同時に、ロングパスリフタを用いてピッチ成分を抽出して、ケフレンシ領域における観測信号から声道と声門の成分と空間伝達関数の最小位相成分の平均値とピッチ成分を減算することによって残響を抑圧することを特徴とする残響抑圧方法。ここで、Ｊ≧１で与えられる。
【請求項２】
請求項１記載のケプストラム解析による残響抑圧方法において、ｉ番目のマイクロホン
2015049502.tif 000049
で求めることを特徴とする残響抑圧方法。ｍはフレーム時刻、２Ｎは短時間離散フーリエ変換の点数である。
【数１】
2015049502.tif 000050
【請求項３】
請求項１乃至請求項２のいずれか１項に記載のケプストラム解析による残響抑圧方法において、声道と声門の成分、並びに空間伝達関数の最小位相成分を式（２）のショートパ
2015049502.tif 000051
）のピッチ周期で、Ｎｄｉｖは０＜Ｎｄｉｖ＜Ｎ０（ｍ）＜Ｎを満足するように選択される。
【数２】
2015049502.tif 000052
【請求項４】
請求項１乃至請求項３のいずれか１項に記載のケプストラム解析による残響抑圧方法において、音声のピッチ成分を式（３）のロングパスリフタを用いて抽出することを特徴とする残響抑圧方法。
【数３】
2015049502.tif 000053
【請求項５】
請求項１乃至請求項４のいずれか１項に記載のケプストラム解析による残響抑圧方法において、声道と声門の成分、並びに空間伝達関数の最小位相成分からガウス性白色雑音を
2015049502.tif 000054
計算によって求めることを特徴とする残響抑圧方法。μは０＜μ≦１を満足するように選択される。
【数４】
2015049502.tif 000055
【請求項６】
請求項１乃至請求項５のいずれか１項に記載のケプストラム解析による残響抑圧方法に
2015049502.tif 000056
を特徴とする残響抑圧方法。
【数５】
2015049502.tif 000057
【請求項７】
請求項１乃至請求項６のいずれか１項に記載のケプストラム解析による残響抑圧方法に
2015049502.tif 000058
に変換した後、ｅｘｐを底とする指数関数によって残響抑圧信号Ｙ（ωｋ，ｍ）をＪ個求め、式（６）のベクトルから計算された相互相関行列Ｒ（ωｋ，ｍ）の固有値と固有ベクトルから得られた信号空間とこれに直交する雑音空間から音声信号の到来方向を推定することを特徴とする残響抑圧方法。
【数６】
2015049502.tif 000059
【請求項８】
請求項１乃至請求項７のいずれか１項に記載の残響抑圧方法を用いて残響抑圧を行うように構成されていることを特徴とする残響抑圧方法を用いた残響抑圧装置。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、インパルス系列、または白色雑音が調音フィルタを通過することにより生成される音声信号と空間伝達関数の畳み込みによる、残響が存在する観測雑音のみから空間伝達関数の全域通過成分を推定して残響を抑圧する残響抑圧方法に係わり、特に、ケプストラム解析を用いて観測信号の最小位相成分を使用することによって観測信号から残響を抑圧する残響抑圧方法及び残響抑圧装置に関する。
【背景技術】
【０００２】  複数のインパルス系列、または白色雑音が調音フィルタを通過することにより生成される音声信号と空間伝達関数の畳み込みにより残響が存在する信号が観測されるとき、観測信号のみから残響を抑圧する処理を残響抑圧という。【０００３】  ケプストラム解析に基づく残響抑圧方法が非特許論文１で提案されている。ケフレンシ領域において観測信号の最小位相成分を求めた後、そのフレーム平均によって空間伝達関数の最小位相成分を推定し、現フレーム時刻の観測信号から減算することによって空間伝達関数の全域通過成分を求め、観測信号の残響を抑圧する。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】  Ａ．Ｓｔｅｐｈｅｎｎｅ  ａｎｄ  Ｂ．Ｃｈａｍｐａｇｎｅ，  “Ｃｅｐｓｔｒａｌ  ｐｒｅｆｉｌｔｅｒｉｎｇ  ｆｏｒ  ｔｉｍｅ  ｄｅｌａｙ  ｅｓｔｉｍａｔｉｏｎ  ｉｎ  ｒｅｖｅｒｂｅｒａｎｔ  ｅｎｖｉｒｏｎｍｅｎｔｓ”，Ｐｒｏｃ．ｏｆ  ＩＣＡＳＳＰ９５，ｖｏｌ．５，ｐｐ．３０５５−３０５８，１９９５
【発明の概要】
【発明が解決しようとする課題】
【０００５】  従来の残響抑圧方法では、ケフレンシ領域において観測信号の最小位相成分を求めた後、有声音と無声音にそれぞれ対応するインパルス系列と白色雑音に無関係に、そのフレーム平均によって空間伝達関数の最小位相成分を推定する為、空間伝達関数の全域通過成分を求めるために、フレーム平均によって求めた空間伝達関数の最小位相成分を現在のフレーム時刻の観測信号から減算すると、インパルス系列と白色雑音の平均成分が空間伝達関数の全域通過成分に含まれ、残響抑圧効果が劣化する。【０００６】  本発明はこのような事情を鑑みてなされたものであり、ケフレンシ領域で空間伝達関数の全域通過成分を推定することによって能率的に観測信号の残響を抑圧することを目的とする。
【課題を解決するための手段】
【０００７】  このような目的に応えるために本発明に係る残響抑圧方法は、インパルス系列、または白色雑音が調音フィルタを通過することにより生成される音声信号と空間伝達関数の畳み込みによる残響が存在する観測信号のみからブラインドで空間伝達関数の全域通過成分をケプストラム解析により推定する方法であって、空間伝達関数が最小位相成分と全域通過成分の積で表現できることに基づき、ケフレンシ領域でフレーム毎に観測信号の最小位相成分を求めた後、ショートパスリフタを用いて声道と声門の成分、並びに空間伝達関数の最小位相成分を抽出して、その平均値を逐次的に計算することによってガウス性白色雑音を除去すると同時に、ロングパスリフタを用いてピッチ成分を抽出して、ケフレンシ領域における観測信号から声道と声門の成分と空間伝達関数の最小位相成分の平均値とピッチ成分を減算することによって残響を抑圧する。【０００８】  すなわち、本発明によれば空間伝達関数が最小位相成分と全域通過成分の積で表現できることに基づき、有声音と無声音に対応する音源信号の種類に無関係に空間伝達関数の全域通過成分を推定することによって観測信号の残響を抑圧させる。
【発明の効果】
【０００９】  本発明によれば、空間伝達関数が最小位相成分と全域通過成分の積で与えられることに基づき、空間伝達関数の全域通過成分を能率的に推定することによって、本発明に係る残響抑圧方法では、部屋の残響時間が長い場合であっても残響を十分に抑圧できるという効果がある。
【図面の簡単な説明】
【００１０】
【図１】  本発明に係る残響抑圧方法の実施の形態を示す図である。
【図２】  本発明に係る残響抑圧方法においてケプストラム解析の後、ショート・ロングパスリフタ、平均と減算を組み合わせて残響を抑圧する方法について説明するためのフローチャートである。
【図３】  本発明に係る残響抑圧方法の実施例１、並びに実施例２における話者（スピーカ）とマイクロホンの位置関係を表す平面図である。
【図４】  本発明に係る残響抑圧方法の実施例１における女性の音声信号の所望到来方向とその推定到来方向との差の絶対値と残響時間の関係を示す図である。
【図５】  本発明に係る残響抑圧方法の実施例１における男性の音声信号の所望到来方向とその推定到来方向との差の絶対値と残響時間の関係を示す図である。
【図６】  本発明に係る残響抑圧方法の実施例１における女性の音声信号の所望到来方向とその推定到来方向との差の絶対値の標準偏差と残響時間の関係を示す図である。
【図７】  本発明に係る残響抑圧方法の実施例１における男性の音声信号の所望到来方向とその推定到来方向との差の絶対値の標準偏差と残響時間の関係を示す図である。
【図８】  本発明に係る残響抑圧方法の実施例２における女性の音声信号の所望到来方向とその推定到来方向との差の絶対値と短時間フーリエ変換の点数との関係を示す図である。
【図９】  本発明に係る残響抑圧方法の実施例２における男性の音声信号の所望到来方向とその推定到来方向との差の絶対値の標準偏差と短時間フーリエ変換の点数との関係を示す図である。
【発明を実施するための最良の形態】
【００１１】  本発明に係る残響抑圧方法の実施の形態について図面を参照して説明する。【００１２】１．音源生成モデル
図１の話者１１では、声帯の振動、または声道の狭めに伴う乱流により発せられた音源信号が、調音器官により形成される調音フィルタを通過することにより音声信号として生成される。有声音として声帯の振動にインパルス系列、無声音として乱流に白色雑音をそれぞれ対応させ、調音フィルタの伝達関数を式（７）で表される全極形で仮定する。インパルス系列、または白色雑音を調音フィルタに入力すると、その出力信号が音声信号になるモデルが音声生成の信号モデルである。【００１３】
【数７】
2015049502.tif 000003
ここで、Ｇｖｔは調音フィルタの利得、２Ｌ（Ｌ＝１、２、…）は極の個数で、Ｌは十分に大きな値である。また、ｚ平面において極ｒｉは単位円の内側に存在し、｜ｒｉ｜＜１である。上付添字†は複素共役を示す。【００１４】  時刻ｔにおいてインパルス系列、または白色雑音ｕ（ｔ）は調音フィルタのインパルス応答ｈｖｔ（ｔ）と畳み込みされて、音声信号ｓ（ｔ）として出力される。音声信号ｓ（ｔ）は式（８）で表される。【００１５】
【数８】
2015049502.tif 000004
ここで、＊は畳み込み演算、インパルス応答ｈｖｔ（ｔ）は調音フィルタの伝達関数Ｈｖｔ（ｚ）を逆ｚ変換したものである。【００１６】２．畳み込みモデルと残響抑圧システム
図１に示すように、時刻ｔにおいて話者１１から発せられた音声信号ｓ（ｔ）が畳み込み演算されてｘｉ（ｔ）として観測される。ｎｉ（ｔ）はマイクロホン２ｉに加わる平均０、分散σ２のガウス性白色雑音で、音声信号と互いに統計的独立な信号である。時刻ｔにおいてＪ個のマイクロホン２１、２２、…、２Ｊで観測される観測信号ｘｉ（ｔ）は式（９）で表される。【００１７】
【数９】
2015049502.tif 000005
ここで、インパルス応答ｈｉ（ｔ）は話者１１からマイクロホン２ｉまでの空間伝達関数Ｈｉ（ｚ）を逆ｚ変換したものである。【００１８】  一般に、残響が存在する線形時不変な話者１１からマイクロホン２ｉまでの空間伝達関数Ｈｉ（ｚ）は非最小位相系で、式（１０）のように最小位相成分Ｈｉｍｉｎ（ｚ）と全域通過成分Ｈｉａｐ（ｚ）の積によって表される。【００１９】
【数１０】
2015049502.tif 000006
【００２０】  Ｈｉｍｉｎ（ｚ）とＨｉａｐ（ｚ）を式（１１）、（１２）に示す。ａｉｊ，ｂｉｊはｚ平面における零点または極を表す複素数である。Ｈｉｍｉｎ（ｚ）は最小位相成分であるので、｜ａｉｊ｜＜１のようにｚ平面において零点は単位円の内側に存在する。また、Ｈｉａｐ（ｚ）は安定で因果律を満足していると仮定すると、｜ｂｉｊ｜＜１より極は単位円の内側に存在する。【００２１】
【数１１】
2015049502.tif 000007
ここで、Ｇｉ＞０である。【００２２】  話者１１からマイクロホン２ｉまでの距離をｄｉ、音速をｃでそれぞれ表記すると、時間遅延τｉは式（１３）で表される。【００２３】
【数１２】
2015049502.tif 000008
【００２４】  次に、ケプストラム解析３ｉの手順について述べる。式（１４）で定義された窓関数ｗ（ｔ）を用いて式（１５）により観測信号ｘｉ（ｔ）を短時間フーリエ変換する。【００２５】
【数１３】
2015049502.tif 000009
ここで、ｍはフレーム時刻、ＴＳは隣接フレーム間のシフトサイズ、ωｋ＝πｋ／Ｎで、ｋ＝０、１、２、…、２Ｎ−１である。【００２６】  ２Ｎがｈｉ（ｔ）のインパルス応答長よりも十分に長ければ、時間周波数領域において観測信号は式（１６）のように近似することができる。【００２７】
【数１４】
2015049502.tif 000010
ここで、Ｓ（ωｋ，ｍ）とＮｉ（ωｋ，ｍ）はフレーム時刻ｍにおけるｓ（ｔ）とｎｉ（ｔ）の短時間フーリエ変換、Ｈｉ（ωｋ）はｈｉ（ｔ）の離散時間フーリエ変換をそれぞれ表す。【００２８】  式（１６）の両辺の複素自然対数をとると、式（１７）を得る。ただし、Ｖｉ（ωｋ，ｍ）は式（１８）によって与えられる。【００２９】
【数１５】
2015049502.tif 000011
【００３０】  逆短時間フーリエ変換によって式（１７）をケフレンシ領域に変換すると、式（１９）のように複素ケプストラムが求められる。【００３１】
【数１６】
2015049502.tif 000012
2015049502.tif 000013
ｇＳ（ωｋ，ｍ）、ｌｏｇＶｉ（ωｋ，ｍ）の逆短時間フーリエ変換、Ｈｉ（ｔ）はｌｏｇＨｉ（ωｋ）は逆離散フーリエ変換をそれぞれ表す。【００３２】  空間伝達関数Ｈｉ（ωｋ）は、式（２０）のように最小位相成分Ｈｉｍｉｎ（ωｋ）と全域通過成分Ｈｉａｐ（ωｋ）の積によって表される。【００３３】
【数１７】
2015049502.tif 000014
【００３４】  式（２０）の複素ケプストラムは式（２１）によって与えられる。【００３５】
【数１８】
2015049502.tif 000015
2015049502.tif 000016
に奇関数である。【００３６】
【数１９】
2015049502.tif 000017
【００３７】  式（１９）に式（２１）を代入すると、式（２４）を得る。【００３８】
【数２０】
2015049502.tif 000018
【００３９】  ケフレンシ領域における音声信号は、音声信号ｓ（ｔ）を短時間フーリエ変換した後、複素自然対数をとると、式（２５）によって求められる。【００４０】
【数２１】
2015049502.tif 000019
ここで、Ｈｖｔ（ωｋ，ｍ）とＳ（ωｋ，ｍ）はフレーム時刻ｍにおけるｈｖｔ（ｔ）とｕ（ｔ）の短時間フーリエ変換をそれぞれ表す。【００４１】  ｌｏｇＳ（ωｋ，ｍ）の逆短時間フーリエ変換を式（２６）に示す。【００４２】
【数２２】
2015049502.tif 000020
【００４３】  式（２４）に式（２６）を代入すると、式（２７）を得る。【００４４】
【数２３】
2015049502.tif 000021
2015049502.tif 000022
）の逆短時間フーリエ変換を表す。【００４５】  図１の３１、３２、…、３Ｊにおける観測信号をケプストラム解析した後、残響抑圧システム４０において本発明により残響抑圧方法４１、４２、…、４Ｊでマイクロホン毎に
2015049502.tif 000023
時間周波数領域に変換して、ｅｘｐを底とする指数関数７ｉを求める。フレーム時刻ｍにおけるＪ個の出力７ｉを重ねて式（２８）のように行列Ｙ（ωｋ，ｍ）を形成する。【００４６】
【数２４】
2015049502.tif 000024
【００４７】  Ｊ個の出力７ｉの相互相関行列Ｒ（ωｋ，ｍ）∈ＣＪ×Ｊを式（２９）によって計算する。【００４８】
【数２５】
2015049502.tif 000025
ここで、Ｅ［・］は期待値、Ｈは複素共役転置をそれぞれ表す。【００４９】  Ｒ（ωｋ，ｍ）の固有値と固有ベクトルにより、信号空間ＵＳ（ωｋ，ｍ）またはＨａｐ（ωｋ）Ｈａｐ（ωｋ）Ｈに直交する雑音空間Ｕｎ（ωｋ，ｍ）が求められる。ここで、Ｈａｐ（ωｋ）は式（３０）によって与えられる。【００５０】
【数２６】
2015049502.tif 000026
【００５１】  マイクロホン間隔１に比べｄｉが十分に大きいとすると、式（３１）を満足するθ（ωｋ）が音声信号の到来方向の推定値である。ここで、ステアリングベクトルａ（θ（ωｋ））は式（３２）で与えられる。式（３１）を満足するθ（ωｋ）は、ＭＵＳＩＣ法などを用いて推定することができる。【００５２】
【数２７】
2015049502.tif 000027
【００５３】３．残響抑圧方法
本発明に係る残響抑圧方法について、図１乃至図２を参照して詳細に説明する。図２は、図１の３１、３２、…、３Ｊにおける観測信号をケプストラム解析した後、残響抑圧システム４０において本発明によりマイクロホン毎に観測信号の残響を抑圧する手順をフローチャートで示したものである。【００５４】
2015049502.tif 000028
ｔ，ｍ）の最小位相成分を出力する。【００５５】
【数２８】
2015049502.tif 000029
【００５６】  有声音では連続フレーム上において声道と声門はゆっくりと変化する。これは、声道と声門はケプストラム領域では低い時間部分に現れることを意味する。一方、ピッチ成分は連続フレーム上において速く変化する。即ち、ケプストラム領域においてビッチ成分は高い時間部分に現れる。そこで、ステップＳ１０２において、式（３４）で与えられるショートパスリフタを用いてケプストラム領域において原点付近に現れる声道と声門の成分、
2015049502.tif 000030
【００５７】
【数２９】
2015049502.tif 000031
2015049502.tif 000032
足するように選択される。男性が発した音声信号なのか、女性が発した音声信号なのか特定できない場合、Ｎｄｉｖは女性の平均ピッチ周期よりも短い値に設定する。これは、一般に女性のピッチ周期が男性に比較して短いことに因る。【００５８】  他方、ピッチ成分を抽出するためには、ステップＳ１０３において式（３５）のロングパスリフタを用いる。【００５９】
【数３０】
2015049502.tif 000033
【００６０】
2015049502.tif 000034
ｏｒｔ（ｔ，ｍ）の平均値を式（３６）によって逐次的に計算する。すなわち、ステップＳ
2015049502.tif 000035
【００６１】
【数３１】
2015049502.tif 000036
2015049502.tif 000037
ある。【００６２】
2015049502.tif 000038
）によって与えられる。【００６３】
【数３２】
2015049502.tif 000039
【００６４】  無声音では白色雑音が調音フィルタに入力されるので、ステップＳ１０２においてショートパスリフタを用いると、ケプストラム領域において原点付近に現れる成分が式（３９）のように抽出される。【００６５】
【数３３】
2015049502.tif 000040
【００６６】  他方、無声音においてピッチ成分を抽出するためにはステップＳ１０３において式（４０）のロングパスリフタを用いる。【００６７】
【数３４】
2015049502.tif 000041
【００６８】
2015049502.tif 000042
なわち、ステップＳ１０４でフレーム時刻ｍが１と判定されると、ステップＳ１０６にお
2015049502.tif 000043
される。【００６９】
2015049502.tif 000044
2015049502.tif 000045
【００７０】
【数３５】
2015049502.tif 000046
【００７１】
2015049502.tif 000047
【実施例１】
【００７２】４．１  評価データ
図３のように、５．０６×３．４１×２．４４メートルの部屋に４センチメートル間隔で２個のマイクロホン２１、２２と、半径１メートルの円の円周上に話者（１個のスピーカ）１１をそれぞれ高さ１．５メートルに配置した。尚、図３は話者（１個のスピーカ）１１と２個のマイクロホン２１、２２の位置関係を示す平面図である。音声信号の所望到来方向は±６０度、±４５度、±３０度、±１５度とした。部屋の残響時間は１５０ミリ秒から８００ミリ秒に設定し、標本化周波数９６ｋＨｚ、量子化ビット数１６ビットで話者１１とマイクロホン２１、２２の間のインパルス応答ｈ１（ｔ）、ｈ２（ｔ）を人工的に発生させた。これを阻止域端周波数８ｋＨｚの低域通過フィルタに入力し、その出力信号を１／６倍にダウンサンプリングして標本化周波数１６ｋＨｚのインパルス応答を算出した。音声データは男性２名と女性２名から３．４５秒間のデータを作成した。実験条件は、２Ｎ＝４０９６点の短時間フーリエ変換、シフトサイズＴＳ＝６４、窓関数はα＝０．９９９８５の指数窓を用いた。ＳＮＲは２０ｄＢに設定した。本発明に係る残響抑圧方
2015049502.tif 000048
ｉｖ＝４ミリ秒、女性が発した音声信号ではＮｄｉｖ＝２ミリ秒に設定した。【００７３】４．２  評価指標
本発明に係る残響抑圧方法で観測信号から残響を抑圧した後、ＭＵＳＩＣ法により推定した音声信号の到来方向と音声信号の所望到来方向との差の絶対値と部屋の残響時間との関係、また、音声の所望到来方向と推定到来方向の差の絶対値の標準偏差と部屋の残響時間との関係によって、残響抑圧方法の性能を評価した。【００７４】４．３  評価対象
ケプストラム領域で空間伝達関数の平均最小位相成分を観測信号から減算することによって残響を抑圧する残響抑圧方法（非特許文献１）と、残響抑圧方法を使用しない場合を比較対象とする。【００７５】４．４  評価結果
音声信号の所望到来方向とその推定到来方向との差の絶対値と残響時間の関係を図４と図５、差の絶対値の標準偏差と残響時間の関係を図６と図７にそれぞれ示す。図４と図６は女性の音声信号、図５と図７は男性の音声信号の到来方向を推定した結果をそれぞれ示している。両図から音声信号の到来方向の推定に、本発明に係る残響抑圧方法が、従来の方法よりも効果があることが明らかとなった。
【実施例２】
【００７６】５．１  評価データ
実施例１と同じ方法で話者１１とマイクロホン２１、２２の間のインパルス応答ｈ１（ｔ）、ｈ２（ｔ）を人工的に発生させた。音声信号の所望到来方向、ＳＮＲ、シフトサイズは実施例１と同じ値に設定した。部屋の残響時間は２７０ミリ秒に設定した。音声データは男性２名と女性２名から、短時間フーリエ変換の点数２Ｎ＝１０２４では３．２６秒間、２Ｎ＝２０４８では３．３２秒間、２Ｎ＝４０９６では３．４５秒間、２Ｎ＝８１９２では３．７１秒間、２Ｎ＝１６３８４では４．２２秒間のデータを作成した。実験条件は、窓関数は短時間フーリエ変換の点数２Ｎ＝１０２４ではα＝０．９９４、２Ｎ＝２０４８ではα＝０．９９７、２Ｎ＝４０９６ではα＝０．９９８５、２Ｎ＝８１９２ではα＝０．９９９２、２Ｎ＝１６３８４ではα＝０．９９９５の指数窓をそれぞれ用いた。本発明に係る残響抑圧方法のμとＮｄｉｖも実施例１と同じ値に設定した。【００７７】５．２  評価対象
本発明に係る残響抑圧方法で観測信号から残響を抑圧した後、ＭＵＳＩＣ法により推定した音声信号の到来方向と音声信号の所望到来方向との差の絶対値と短時間フーリエ変換の点数との関係、また、音声の所望到来方向と推定到来方向の差の絶対値の標準偏差と短時間フーリエ変換の点数の関係によって、残響抑圧方法の性能を評価した。【００７８】５．３  評価対象
実施例１と同様に、ケプストラム領域で空間伝達関数の平均最小位相成分を観測信号から減算することによって残響を抑圧する残響抑圧方法（非特許文献１）と、残響抑圧方法を使用しない場合を比較対象とする。【００７９】５．４  評価結果
音声信号の所望到来方向とその推定到来方向との差の絶対値と短時間フーリエ変換の点数の関係を図８、差の絶対値の標準偏差と短時間フーリエ変換の点数の関係を図９にそれぞれ示す。両図は男性と女性の音声信号の推定到来方向の平均値をそれぞれ示している。両図から音声信号の到来方向の推定に、本発明に係る残響抑圧方法が、従来の方法よりも効果があることが明らかとなった。
【符号の説明】
【００８０】１１…話者、２１、２２、…、２Ｊ…マイクロホン、３１、３２、…、３Ｊ…ケプストラム解析、４０…残響抑圧システム、５１、５２、…、５Ｊ…残響抑圧方法、６１、６２、…、６Ｑ…短時間フーリエ変換、７１、７２、…、７Ｊ…ｅｘｐを底とする指数関数
【図１】
2015049502.tif 000060
【図２】
2015049502.tif 000061
【図３】
2015049502.tif 000062
【図４】
2015049502.tif 000063
【図５】
2015049502.tif 000064
【図６】
2015049502.tif 000065
【図７】
2015049502.tif 000066
【図８】
2015049502.tif 000067
【図９】
2015049502.tif 000068