(57)【要約】      （修正有）【課題】事前に観測、或いは予測した各マイクロホンの方向毎の感度等の先験情報を用いることなく、複数個のマイクロホンを用いて収音した複数方向の音源信号を個別に低雑音化して抽出する方向別収音装置及びプログラムを提供する。【解決手段】本発明の方向別収音装置１は、複数個のマイクロホン１０を用いて収音した複数方向の音源信号についてそれぞれ周波数領域の観測信号に変換するＤＦＴ処理部２０、各観測信号についてＰＳＤを算出するＰＳＤ算出部３０、感度比推定部４３により推定算出した所定の複素感度比及びＰＳＤ和推定部５０により推定算出したＰＳＤ和を用いてｍ番目のマイクロホンの観測信号から雑音成分を除去するための第ｍウィーナーフィルタを構成し、該観測信号に適用して雑音除去信号を生成する第ｍ雑音抑圧部４０及びその雑音除去信号に対しＩＤＦＴを施し時間領域信号に変換するＩＤＦＴ処理部６０を備える。【選択図】図２
2020010196.tif 000002
【特許請求の範囲】
【請求項１】
複数個のマイクロホンを用いて収音した複数方向の音源信号について個別に低雑音化して抽出する方向別収音装置であって、
前記複数個のマイクロホンを用いて収音した複数方向の音源信号についてそれぞれ離散フーリエ変換を施し周波数領域の観測信号に変換するＤＦＴ処理手段と、
それぞれの当該観測信号についてパワースペクトル密度を算出するＰＳＤ算出手段と、
推定算出した所定の複素感度比、及び推定算出した各マイクロホンの正面方向の音源信号のパワースペクトル密度の総和であるＰＳＤ和を用いて、ｍ番目のマイクロホンの観測信号から他の特定のマイクロホンの観測信号の成分を除去するためのクロススペクトル密度に基づくウィーナーフィルタを個別に算出して総乗し、当該ｍ番目のマイクロホンの観測信号から他の全てのマイクロホンの観測信号の成分を除去するための第ｍウィーナーフィルタを構成するウィーナーフィルタ構成手段と、
当該ｍ番目のマイクロホンの観測信号に対し前記第ｍウィーナーフィルタを施して該観測信号に関する雑音除去信号を生成するウィーナーフィルタ処理手段と、
前記雑音除去信号に対し逆離散フーリエ変換を施し、時間領域信号に変換して出力するＩＤＦＴ処理手段と、
を備えることを特徴とする方向別収音装置。
【請求項２】
前記ｍ番目のマイクロホンの観測信号と、該観測信号に関する前記雑音除去信号とを用いて単位時間毎に前記所定の複素感度比を推定算出し、当該個別に算出したウィーナーフィルタを単位時間毎に更新する感度比推定手段を更に備えることを特徴とする、請求項１に記載の方向別収音装置。
【請求項３】
前記感度比推定手段は、前記複数個（Ｍ個）のマイクロホンのうち周波数ωにおける当該ｍ番目のマイクロホンの正面方向の感度とｌ番目のマイクロホンの正面方向における該ｍ番目のマイクロホンの感度との複素感度比について当該推定算出した推定複素感度比ｇｌｍ(ω) （ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ）を用いて、単位時間τ毎の更新を考慮した重み係数αを、
【数１】
2020010196.tif 000031
に基づいて推定算出することを特徴とする、請求項２に記載の方向別収音装置。
【請求項４】
前記複数個のマイクロホンの全てに関する当該複素感度比を用いて重み係数を算出し、該重み係数による重み付き線形和により、各マイクロホンの正面方向の音源信号のＰＳＤの総和であるＰＳＤ和を推定算出するＰＳＤ和推定手段を更に備えることを特徴とする、請求項３に記載の方向別収音装置。
【請求項５】
前記ＰＳＤ和推定手段は、各マイクロホンで収音される信号のパワースペクトル密度ＰＸｉ(ω，τ) （ｉ＝１，２，…，Ｍ）に対し当該重み係数αｉを用いて
【数２】
2020010196.tif 000032
で表される重み付き線形和により、前記ＰＳＤ和を推定算出することを特徴とする、請求項４に記載の方向別収音装置。
【請求項６】
コンピューターを、請求項１から５のいずれか一項に記載の方向別収音装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数個のマイクロホンを用いて収音した複数方向の音源信号についてウィーナーフィルタを応用し個別に低雑音化して抽出する方向別収音装置及びプログラムに関する。
【背景技術】
【０００２】
従来、任意方向の音源信号を強調して抽出するために指向性収音技術とウィーナーフィルタを組み合わせた技法が知られている（例えば、非特許文献１、特許文献１参照）がある。これらの技法では、複数個の指向性マイクロホンよりなるマイクロホンアレイでビームフォーミングを行った後に、ウィーナーフィルタを用いて、目的方向外からの雑音を除去することで方向別収音を実現している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１７−１４３３２５号公報
【非特許文献】
【０００４】
【非特許文献１】Y.Hioka and K.Furuya, K.Kobayashi, K.Niwa, Y.Haneda,’ Underdetermined Sound Source Separation Using Power Spectrum Density Estimated by Combination of Directivity Gain,’ IEEE Transaction on audio, speech, and language processing, Vol.21, No.6, June 2013
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来技術における、複数個の指向性マイクロホンよりなるマイクロホンアレイと、ウィーナーフィルタを用いて、方向別収音を行う技法では、そのウィーナーフィルタのフィルタ設計に、マイクロホンアレイ周囲方向毎の感度が既知である必要があり、事前にそれらを計測しておかなければならなかった。
【０００６】
そのため、各マイクロホンの方向毎の感度が未知の場合には直ちに適応できないという問題がある。
【０００７】
従って、複数個のマイクロホンを用いて収音した複数方向の音源信号を個別に抽出する際に、事前に観測、或いは予測した各マイクロホンの方向毎の感度等の先験情報を用いることなく、当該複数方向の音源信号を個別に低雑音化して抽出する技法が望まれる。
【０００８】
本発明の目的は、上述の問題に鑑みて、事前に観測、或いは予測した各マイクロホンの方向毎の感度等の先験情報を用いることなく、複数個のマイクロホンを用いて収音した複数方向の音源信号を個別に低雑音化して抽出する方向別収音装置及びプログラムを提供することにある。
【課題を解決するための手段】
【０００９】
本発明に係る方向別収音装置においては、複数個のマイクロホンを用いて収音した複数方向の音源信号を個別に抽出する際に、各マイクロホンの正面外に存在する雑音源方向の感度と正面に存在する所望音源方向の感度比を推定してウィーナーフィルタを適用することにより所望方向外からの雑音を除去し、当該複数方向の音源信号を個別に低雑音化して抽出する。
【００１０】
即ち、本発明の方向別収音装置は、複数個のマイクロホンを用いて収音した複数方向の音源信号について個別に低雑音化して抽出する方向別収音装置であって、前記複数個のマイクロホンを用いて収音した複数方向の音源信号についてそれぞれ離散フーリエ変換を施し周波数領域の観測信号に変換するＤＦＴ処理手段と、それぞれの当該観測信号についてパワースペクトル密度を算出するＰＳＤ算出手段と、推定算出した所定の複素感度比、及び推定算出した各マイクロホンの正面方向の音源信号のパワースペクトル密度の総和であるＰＳＤ和を用いて、ｍ番目のマイクロホンの観測信号から他の特定のマイクロホンの観測信号の成分を除去するためのクロススペクトル密度に基づくウィーナーフィルタを個別に算出して総乗し、当該ｍ番目のマイクロホンの観測信号から他の全てのマイクロホンの観測信号の成分を除去するための第ｍウィーナーフィルタを構成するウィーナーフィルタ構成手段と、当該ｍ番目のマイクロホンの観測信号に対し前記第ｍウィーナーフィルタを施して該観測信号に関する雑音除去信号を生成するウィーナーフィルタ処理手段と、前記雑音除去信号に対し逆離散フーリエ変換を施し、時間領域信号に変換して出力するＩＤＦＴ処理手段と、を備えることを特徴とする。
【００１１】
また、本発明の方向別収音装置において、前記ｍ番目のマイクロホンの観測信号と、該観測信号に関する前記雑音除去信号とを用いて単位時間毎に前記所定の複素感度比を推定算出し、当該個別に算出したウィーナーフィルタを単位時間毎に更新する感度比推定手段を更に備えることを特徴とする。
【００１２】
また、本発明の方向別収音装置において、前記感度比推定手段は、前記複数個（Ｍ個）のマイクロホンのうち周波数ωにおける当該ｍ番目のマイクロホンの正面方向の感度とｌ番目のマイクロホンの正面方向における該ｍ番目のマイクロホンの感度との複素感度比について当該推定算出した推定複素感度比ｇｌｍ(ω) （ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ）を用いて、単位時間τ毎の更新を考慮した重み係数αを、所定の算出法（後述する［数２５］）に基づいて推定算出することを特徴とする。
【００１３】
また、本発明の方向別収音装置において、前記複数個のマイクロホンの全てに関する当該複素感度比を用いて重み係数を算出し、該重み係数による重み付き線形和により、各マイクロホンの正面方向の音源信号のＰＳＤの総和であるＰＳＤ和を推定算出するＰＳＤ和推定手段を更に備えることを特徴とする。
【００１４】
また、本発明の方向別収音装置において、前記ＰＳＤ和推定手段は、各マイクロホンで収音される信号のパワースペクトル密度ＰＸｉ(ω，τ) （ｉ＝１，２，…，Ｍ）に対し当該重み係数αｉを用いて所定の算出法（後述する［数２６］）で表される重み付き線形和により、前記ＰＳＤ和を推定算出することを特徴とする。
【００１５】
また、本発明のプログラムは、コンピューターを、本発明の方向別収音装置として機能させるためのプログラムとして構成する。
【発明の効果】
【００１６】
本発明によれば、事前に観測、或いは予測した各マイクロホンの方向毎の感度等の先験情報を用いることなく、目的方向の音源信号だけを明瞭に収音することができる。
【図面の簡単な説明】
【００１７】
【図１】本発明による一実施形態の方向別収音装置に係るウィーナーフィルタを適用する音響モデルの概略図である。
【図２】本発明による一実施形態の方向別収音装置の概略構成を示すブロック図である。
【図３】（ａ）,（ｂ）は、それぞれ本発明による一実施形態の方向別収音装置に係る雑音抑圧処理の有無を比較した収音例を示す図である。
【発明を実施するための形態】
【００１８】
まず、図１を参照して本発明による一実施形態の方向別収音装置１に係るウィーナーフィルタを適用する音響モデルについて説明し、続いて図２を参照して、ウィーナーフィルタを応用した本発明による一実施形態の方向別収音装置１について説明する。
【００１９】
〔ウィーナーフィルタを適用する音響モデル〕
図１は、本発明による一実施形態の方向別収音装置１に係るウィーナーフィルタ１００を適用する音響モデルの概略図である。
【００２０】
図１に示すように、Ｍ個の音源からの音源信号を個別に抽出するために各音源方向に向けたＭ個のマイクロホン（指向性マイクロホン）１０で構成されたマイクロホンアレイにより収音することを想定する。そして、ｍ番目のマイクロホン正面に位置する互いに無相関な音源信号をＳｍ(ω，τ)とし、ｍ番目のマイクロホン１０の観測信号Ｘｍ(ω，τ)を式（１）のようにモデル化する。
【００２１】
【数１】
2020010196.tif 000003
【００２２】
ｈｉｍ（ω）は、 周波数ωにおけるｍ番目のマイクロホン１０の正面方向の感度とｉ番目のマイクロホン１０の正面方向におけるｍ番目のマイクロホン１０の感度との比（以下、複素感度比とする）を表し、｜ｈｉｍ(ω)｜＜１であると想定する。即ち、ｍ番目のマイクロホン１０の角度ごとの感度を角度θの関数としてｈｍ(θ，ω)とすると、ｍ番目のマイクロホン１０の正面方向（θ＝θｍ）とｉ番目のマイクロホン１０の正面方向との開き角がθｉのとき、ｈｉｍ（ω） ＝ ｈｍ(θ＝θｍ ＋θｉ, ω)／ｈｍ（θ＝θｍ, ω）となる。ここで、ω、τはそれぞれ周波数のインデックス、時間フレームのインデックスを表す。
【００２３】
このとき、ウィーナーフィルタ１００において、観測信号Ｘｍ(ω，τ)を、式（２）のウィーナーフィルタｗｌ(ω，τ)で式（３）のようにフィルタリングすることで、観測信号Ｘｍ(ω，τ)から不所望な音源信号Ｓｌ(ω，τ)のみを除去することができる。
【００２４】
【数２】
2020010196.tif 000004
【００２５】
【数３】
2020010196.tif 000005
【００２６】
ここで、Ｐｍ(ω，τ)はｍ番目のマイクロホン１０で収音する音源信号Ｓｍ(ω，τ)のパワースペクトル密度（ＰＳＤ）である。Ｅ［＊］は期待値演算を表すが、実用においては過去数フレームの平均によって推定されるため、τフレームにおける推定値をＰｍ(ω，τ)としている。
【００２７】
また、０≦｜ｋｉ｜≦１（ｉ＝１，２,…,Ｍ）となる任意の係数ｋｉについて、ウィーナーフィルタｗｌ(ω，τ)を式（４）とする。
【００２８】
【数４】
2020010196.tif 000006
【００２９】
式（４）でウィーナーフィルタｗｌ(ω，τ)を表わすと、特定の音源信号Ｓｉ(ω，τ)（ｉ≠ｌ，ｍ）に影響を与えるものの、ｍ番目のマイクロホンの観測信号Ｘｍ(ω，τ)から、所望の音源信号Ｓｍ(ω，τ)を損なうことなく、不所望な音源信号Ｓｌ(ω，τ)を除去することができる。
【００３０】
したがって、ウィーナーフィルタ１００において、観測信号Ｘｍ(ω，τ)から所望の音源信号Ｓｍ(ω，τ)のみを復元するには、式（５）で表されるウィーナーフィルタＷ (ω，τ)によって観測信号Ｘｍ(ω，τ)を式（６）のようにフィルタリングすればよい。
【００３１】
【数５】
2020010196.tif 000007
【００３２】
【数６】
2020010196.tif 000008
【００３３】
以下、ウィーナーフィルタ１００における、式（５）のウィーナーフィルタＷ(ω，τ)の設計手法について詳細に述べる。
【００３４】
各マイクロホン１０の正面方向に位置する音源が互いに無相関であると仮定すると、観測信号Ｘｍ(ω，τ)のＰＳＤを示すＰＸｍ(ω，τ) は、式（７）のように推定できる。
【００３５】
【数７】
2020010196.tif 000009
【００３６】
Ｅ［＊］は期待値演算を表し、実際の処理の際には複数フレームにおける平均で代用する。
【００３７】
また、ｌ番目のマイクロホン１０の観測信号Ｘｌ(ω，τ)とｍ番目のマイクロホン１０の観測信号Ｘｍ(ω，τ)のクロススペクトル密度（ＣＳＤ：Cross Spectral Density）であるＣｌｍ(ω，τ)は、式（８）のように推定できる。
【００３８】
【数８】
2020010196.tif 000010
【００３９】
ここで、観測信号Ｘｍ(ω，τ)のＰＳＤを示すＰＸｍ(ω，τ)から、観測信号Ｘｌ(ω，τ)のＰＳＤを示すＰｌ(ω，τ)の項を除去するために式（９）の計算を行う。
【００４０】
【数９】
2020010196.tif 000011
【００４１】
ここで、Ｒｅ［＊］は複素数の実部を示す。式（９）から式（１０）が得られる。
【００４２】
【数１０】
2020010196.tif 000012
【００４３】
2020010196.tif 000013
【００４４】
【数１１】
2020010196.tif 000014
【００４５】
そして、ウィーナーフィルタ１００では、Ｍ個の観測信号Ｘｉ(ω，τ)に基づくＰＸｉ(ω，τ) （ｉ＝１，２,…,Ｍ）の重み付け線形和の信号として処理する。
【００４６】
【数１２】
2020010196.tif 000015
【００４７】
ここで、ＰＸｉ(ω，τ)に対する重みα＝［α１α２α３…αＭ］Ｔが式（１３）を満たすとする。ここで、［＊］Ｔはベクトル・行列の転置、［＊］−１は逆行列を表す。
【００４８】
【数１３】
2020010196.tif 000016
【００４９】
式（１３）を満たすとき、式（１４）が得られる。
【００５０】
【数１４】
2020010196.tif 000017
【００５１】
ここで、ｍ番目のマイクロホン１０とｉ番目のマイクロホン１０の正面方向の複素感度比ｈｉｍ（ω）は未知であるが、雑音対象のｌ番目のマイクロホン１０について後述するように推定複素感度比ｇｌｍ(ω)が得られるとすると、ｈｌｍ(ω) ≒ｇｌｍ(ω)として扱うことができ、式（１３）は式（１５）のようにｈｌｍ(ω)をｇｌｍ(ω)で代用することができる。
【００５２】
【数１５】
2020010196.tif 000018
【００５３】
従って、式（１１）と式（１４），（１５）から、ウィーナーフィルタｗｌ(ω，τ)を式（１６）とすることができる。
【００５４】
【数１６】
2020010196.tif 000019
【００５５】
2020010196.tif 000020
【００５６】
続いて、ｇｌｍ(ω)≒ｈｌｍ(ω)の推定方法について詳細を述べる。
【００５７】
まず、式（１７）のように、ウィーナーフィルタ１００を用いて観測信号Ｘｌ(ω，τ)から音源信号Ｓｌ(ω，τ)を復元できたとする。
【００５８】
【数１７】
2020010196.tif 000021
【００５９】
ｈｌｍ（ω）は、ｌ番目のマイクロホン１０で観測される音源信号Ｓｌ(ω，τ)とｍ番目のマイクロホン１０で観測される観測信号Ｘｍ(ω，τ)に含まれるＳｌ(ω，τ)成分との感度比であり、適応アルゴリズムを用いて逐次的に同定することが可能である。適応アルゴリズムは様々なアルゴリズムが考案されているが、ここでは例えばＬＭＳアルゴリズムを用いると、ｇｌｍ(ω，τ)は、式（１８），（１９）のようにフレーム毎に更新することができる。
【００６０】
【数１８】
2020010196.tif 000022
【００６１】
【数１９】
2020010196.tif 000023
【００６２】
μはステップサイズパラメータと呼ばれる、学習の速度を調整するパラメータであり、ユーザーが任意に決定することができる。
【００６３】
このように、複数個のマイクロホン１０を用いて収音した複数方向の音源信号を個別に抽出する際に、各マイクロホン１０の雑音方向の感度と雑音方向の感度比を推定してウィーナーフィルタ１００を適用することにより所望方向外からの雑音を除去し、当該複数方向の音源信号を個別に低雑音化して抽出することができる。
【００６４】
以上を踏まえ、図２を参照して、ウィーナーフィルタを応用した本発明による一実施形態の方向別収音装置１について説明する。
【００６５】
〔方向別収音装置〕
図２に示すように、本発明による一実施形態の方向別収音装置１は、Ｍ（Ｍ＞１）個のマイクロホン１０を用いて収音した複数方向の音源信号についてウィーナーフィルタを応用し個別に低雑音化して抽出する装置であり、Ｍ個のＤＦＴ処理部２０、Ｍ個のＰＳＤ算出部３０、Ｍ個の雑音抑圧部４０（４０‐１，４０‐２，…，４０‐Ｍ）、１個のＰＳＤ和推定部５０、及びＭ個のＩＤＦＴ処理部６０を備える。
【００６６】
以下、方向別収音装置１が備える各構成要素について、Ｍ個のマイクロホン１０から出力される出力信号を処理する順に説明する。
【００６７】
まず、Ｍ個のマイクロホン１０は、それぞれ図１に示したように収音を所望する音源に向けられており、各マイクロホン１０の正面方向に当該音源が位置している。
【００６８】
Ｍ個のＤＦＴ処理部２０は、それぞれＭ個のマイクロホン１０の出力信号を入力し、それぞれ当該出力信号に対し離散フーリエ変換（ＤＦＴ）を施し周波数領域の観測信号Ｘｉ(ω，τ) （ｉ＝１，２,…,Ｍ）に変換して、Ｍ個のＰＳＤ算出部３０にそれぞれ出力する。ここで、ω、τはそれぞれ周波数のインデックス、時間フレームのインデックスを表す。
【００６９】
Ｍ個のＰＳＤ算出部３０は、それぞれＭ個のＤＦＴ処理部２０から入力された観測信号Ｘｉ(ω，τ) （ｉ＝１，２，…，Ｍ）についてパワースペクトル密度（ＰＳＤ）を算出し、Ｍ個の雑音抑圧部４０の全てに出力する。
【００７０】
尚、各観測信号Ｘｉ(ω，τ)のＰＳＤであるＰＸｉ(ω，τ)は、上述した式（７）と同様に次式から算出する。
【００７１】
【数２０】
2020010196.tif 000024
【００７２】
ここで、Ｅ［＊］は期待値演算を表し、実際の処理の際には複数フレームにおける平均で代用する。
【００７３】
これにより、着目するｍ番目のマイクロホン１０の観測信号Ｘｍ(ω，τ)のＰＳＤを示すＰＸｍ(ω，τ)、当該観測信号Ｘｍ(ω，τ)から除去したいｌ番目のマイクロホン１０の観測信号Ｘｌ(ω，τ)のＰＳＤを示すＰｌ(ω，τ)が得られる。
【００７４】
Ｍ個の雑音抑圧部４０（４０‐１，４０‐２，…，４０‐Ｍ）は、それぞれ同様の構成要素を有し、代表して図２に示すように、ｍ番目のマイクロホン１０の観測信号Ｘｍ(ω，τ)を雑音抑圧の処理対象とする第ｍ雑音抑圧部４０‐ｍについて説明する。
【００７５】
（ウィーナーフィルタによる雑音抑圧処理）
第ｍ雑音抑圧部４０‐ｍは、ウィーナーフィルタ構成部４１、ウィーナーフィルタ処理部４２、及び感度比推定部４３を備える。また、ウィーナーフィルタ構成部４１は、ＣＳＤ算出部４１１、第ｍ‐ｌウィーナーフィルタ算出部４１２、及びフィルタ乗算部４１３を有する。
【００７６】
第ｍ雑音抑圧部４０‐ｍにおいて、まずＣＳＤ算出部４１１は、第ｍマイクロホン１０の観測信号Ｘｍ(ω，τ)と、他のマイクロホン１０である第ｌマイクロホン１０の観測信号Ｘｌ(ω，τ)とのクロススペクトル密度（ＣＳＤ）であるＣｍｌ(ω，τ)を、上述した式（８）に基づいて算出し、第ｍ‐ｌウィーナーフィルタ算出部４１２に出力する。
【００７７】
尚、式（８）に基づいたＣＳＤの算出において、Ｅ［＊］は期待値演算を表し、実際の処理の際には複数フレームにおける平均で代用する。
【００７８】
第ｍ‐ｌウィーナーフィルタ算出部４１２は、後述する感度比推定部４３によって推定算出した複素感度比ｇｌｍ(ω，τ) 、及び後述するＰＳＤ和推定部５０によって推定算出した各マイクロホン正面方向の音源信号のＰＳＤの総和であるＰＳＤ和｛ΣＭｉ＝１Ｐｉ(ω，τ)｝を用いて、第ｍマイクロホン１０の観測信号Ｘｍ(ω，τ)から他のマイクロホン１０である第ｌマイクロホン１０の観測信号Ｘｌ(ω，τ)の成分を除去するためのクロススペクトル密度（ＣＳＤ）に基づくウィーナーフィルタを算出する機能部である（ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ）。
【００７９】
より具体的に、第ｍ‐ｌウィーナーフィルタ算出部４１２は、当該観測信号Ｘｍ(ω，τ)のＰＳＤを示すＰＸｍ(ω，τ)、観測信号Ｘｌ(ω，τ)のＰＳＤを示すＰｌ(ω，τ)、及び、観測信号Ｘｌ(ω，τ)とマイクロホンｍの観測信号Ｘｍ(ω，τ)のクロススペクトル密度（ＣＳＤ）であるＣｍｌ(ω，τ)と、後述する推定複素感度比ｇｌｍ(ω，τ)及びＰＳＤ和｛ΣＭｉ＝１Ｐｉ(ω，τ)｝とを用いて、次式によって第ｍ‐ｌウィーナーフィルタｗｍ−ｌ(ω，τ)を算出し、フィルタ乗算部４１３に出力する。
【００８０】
【数２１】
2020010196.tif 000025
【００８１】
フィルタ乗算部４１３は、各第ｍ‐ｌウィーナーフィルタｗｍ−ｌ(ω，τ) （ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ）を入力して、次式のように、得られた第ｍ‐ｌウィーナーフィルタ（ｌ≠ｍ）の関数をＭ個総乗することで第ｍウィーナーフィルタＷｍ(ω，τ)を算出し、ウィーナーフィルタ処理部４２に出力する。
【００８２】
【数２２】
2020010196.tif 000026
【００８３】
このようにして、第ｍ雑音抑圧部４０‐ｍは、当該観測信号Ｘｍ(ω，τ)のＰＳＤを示すＰＸｍ(ω，τ)、観測信号Ｘｌ(ω，τ)のＰＳＤを示すＰｌ(ω，τ)、及び、観測信号Ｘｌ(ω，τ)とマイクロホンｍの観測信号Ｘｍ(ω，τ)のクロススペクトル密度（ＣＳＤ）であるＣｍｌ(ω，τ)と、後述する推定複素感度比ｇｌｍ(ω，τ)及びＰＳＤ和｛ΣＭｉ＝１Ｐｉ(ω，τ)｝とを用いて、第ｍウィーナーフィルタＷｍ(ω，τ)を算出し、ウィーナーフィルタ処理部４２に出力する。
【００８４】
ウィーナーフィルタ処理部４２は、次式のように、当該観測信号Ｘｍ(ω，τ)に対し第ｍウィーナーフィルタＷｍ(ω，τ)を施して当該観測信号Ｘｍ(ω，τ)に関する雑音除去信号を生成し、ＩＤＦＴ処理部６０に出力する。
【００８５】
【数２３】
2020010196.tif 000027
【００８６】
ＩＤＦＴ処理部６０は、ウィーナーフィルタ処理部４２から得られる雑音除去信号に対し逆離散フーリエ変換を施し、時間領域信号に変換して出力する。
【００８７】
（ウィーナーフィルタ算出に用いる感度比推定処理）
続いて第ｍ雑音抑圧部４０‐ｍにおける感度比推定部４３の処理について説明する。
【００８８】
感度比推定部４３は、ｍ番目のマイクロホンで観測される観測信号Ｘｍ(ω，τ)とウィーナーフィルタ処理部４２から得られる雑音除去信号を基に、ｌ番目のマイクロホンで観測される音源信号Ｓｌ(ω，τ)と、ｍ番目のマイクロホンで観測される観測信号Ｘｍ(ω，τ)に含まれる音源信号Ｓｌ(ω，τ)の成分との感度比である推定複素感度比ｇｌｍ(ω，τ)を単位時間毎に算出し、第ｍ‐ｌウィーナーフィルタ算出部４１２に出力する。これにより、各第ｍ‐ｌウィーナーフィルタｗｍ−ｌ(ω，τ) （ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ）は単位時間毎に更新されるようになる。
【００８９】
推定複素感度比ｇｌｍ(ω，τ)は、適応アルゴリズムを用いて逐次的に同定することが可能である。適応アルゴリズムは様々なアルゴリズムが考案されているが、ここでは例えばＬＭＳアルゴリズムを用いると、推定複素感度比ｇｌｍ(ω，τ)は、次式のようにフレーム毎に更新することができる。
【００９０】
【数２４】
2020010196.tif 000028
【００９１】
μ はステップサイズパラメータと呼ばれる、学習の速度を調整するパラメータであり、ユーザーが任意に決定することができる。
【００９２】
Ｍ個の雑音抑圧部４０（４０‐１，４０‐２，…，４０‐Ｍ）は、それぞれ上記の第ｍ雑音抑圧部４０‐ｍと同様の処理を並行して行うことで全てのマイクロホン１０に関する複素感度比ｇｌｍ(ω，τ)（ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ）を推定算出することができる。
【００９３】
（ウィーナーフィルタ算出に用いるＰＳＤ和推定処理）
ＰＳＤ和推定部５０は、重み乗算部５１及びＰＳＤ和算出部５２を有する。
【００９４】
重み乗算部５１は、Ｍ個の雑音抑圧部４０（４０‐１，４０‐２，…，４０‐Ｍ）から、全てのマイクロホン１０に関する推定複素感度比ｇｌｍ(ω，τ)を入力し、ｇｌｍ(ω，τ)（ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ；ｍ＝１，２，…，Ｍ）を用いて、上述した式（１５）に示したものに単位時間毎の更新を考慮し、次式のように重み係数αを推定算出する。
【００９５】
【数２５】
2020010196.tif 000029
【００９６】
ここで、重み乗算部５１は、周波数ωにおける当該ｍ番目のマイクロホン１０の正面方向の感度とｌ番目のマイクロホン１０正面方向におけるｍ番目のマイクロホン１０の感度との複素感度比について、Ｍ個の雑音抑圧部４０によって当該推定算出した推定複素感度比ｇｌｍ(ω，τ)に関する重み係数αｉをＰＳＤ和算出部５２に出力する。
【００９７】
ＰＳＤ和算出部５２は、重み係数αｉを用いて、各マイクロホン１０の正面方向の音源信号のＰＳＤの総和であるＰＳＤ和｛ΣＭｉ＝１Ｐｉ(ω，τ)｝について推定算出し、Ｍ個の雑音抑圧部４０（４０‐１，４０‐２，…，４０‐Ｍ）における第ｍ‐ｌウィーナーフィルタ算出部４１２に出力する。
【００９８】
ＰＳＤ和｛ΣＭｉ＝１Ｐｉ(ω，τ)｝は、次式のように各マイクロホン１０で収音される観測信号のＰＳＤであるＰＸｉ(ω，τ) （ｉ＝１，２，…，Ｍ）に対する重み係数αｉによる重み付き線形和で得ることができる。
【００９９】
【数２６】
2020010196.tif 000030
【０１００】
これにより、推定複素感度比ｇｌｍ(ω，τ)によって各第ｍ‐ｌウィーナーフィルタｗｍ−ｌ(ω，τ) （ｌ＝１，２，…，ｍ−１，ｍ＋１，…，Ｍ）について単位時間毎に更新した場合でも、第ｍウィーナーフィルタＷｍ(ω，τ)についてＰＳＤ和で正規化することができ、当該観測信号Ｘｍ(ω，τ)に対し第ｍウィーナーフィルタＷｍ(ω，τ)を施して得られる雑音除去信号を安定化させることができる。
【０１０１】
以上のように、本実施形態の方向別収音装置１では、複数個のマイクロホン１０を用いて収音した複数方向の音源信号を個別に抽出する際に、各マイクロホン１０の雑音方向の感度と雑音方向の感度比を推定してウィーナーフィルタＷｍを構成し、これを適用することにより所望方向外からの雑音を除去し、当該複数方向の音源信号を個別に低雑音化して抽出するようにした。
【０１０２】
これにより、本実施形態の方向別収音装置１によれば、事前に観測、或いは予測した各マイクロホン１０の方向毎の感度等の先験情報を用いることなく、適応的に目的音方向にマイクロホン１０を向けた場合でも、目的方向の音源信号だけを明瞭に収音することができるようなる。
【０１０３】
図３（ａ）,（ｂ）には、それぞれ本発明による一実施形態の方向別収音装置１に係る雑音抑圧処理の有無を比較した収音例を示している。より具体的には、図３（ａ）,（ｂ）にそれぞれ示す収音例１，２は周波数毎の雑音レベルを示すスペクトラム図であり、或る音源Ｍに向けられた第ｍマイクロホンの出力信号の特性（雑音抑圧処理無し）と、当該第ｍマイクロホンの出力信号に対し本実施形態の方向別収音装置１により雑音除去を行って抽出した出力信号の特性（雑音抑圧処理有り）とを比較して示している。図３（ａ）,（ｂ）から理解されるように、本発明による一実施形態の方向別収音装置１によれば全体として雑音除去が有効に作用していることが分かる。
【０１０４】
以上の実施形態における方向別収音装置１は、コンピューターにより構成することができ、方向別収音装置１の各処理部を機能させるためのプログラムを好適に用いることができる。具体的には、方向別収音装置１の各処理部を制御するための制御部をコンピューター内の中央演算処理装置（ＣＰＵ）で構成でき、且つ、各処理部を動作させるのに必要となるプログラムを適宜記憶する記憶部を少なくとも１つのメモリで構成させることができる。即ち、そのようなコンピューターに、ＣＰＵによって該プログラムを実行させることにより、方向別収音装置１の各処理部の有する機能を実現させることができる。更に、方向別収音装置１の各処理部の有する機能を実現させるためのプログラムを、前述の記憶部（メモリ）の所定の領域に格納させることができる。そのような記憶部は、装置内部のＲＡＭ又はＲＯＭなどで構成させることができ、或いは又、外部記憶装置（例えば、ハードディスク）で構成させることもできる。また、そのようなプログラムは、コンピューターで利用されるＯＳ上のソフトウェア（ＲＯＭ又は外部記憶装置に格納される）の一部で構成させることができる。更に、そのようなコンピューターに、方向別収音装置１の各処理部として機能させるためのプログラムは、コンピューター読取り可能な記録媒体に記録することができる。また、方向別収音装置１の各処理部をハードウェア又はソフトウェアの一部として構成させ、各々を組み合わせて実現させることもできる。
【０１０５】
以上、特定の実施形態の例を挙げて本発明を説明したが、本発明は前述した例に限定されるものではなく、その技術思想を逸脱しない範囲で種々変形可能である。例えば、上述した例では、特定の音源方向に向けたマイクロホン１０について指向性マイクロホンを想定して説明したが、任意の音源方向を収音する任意の指向性のマイクロホンに対して本発明に係る方向別収音装置１を適用できる。従って、本発明に係る方向別収音装置１は、上述した実施形態の例に限定されるものではなく、特許請求の範囲の記載によってのみ制限される。
【産業上の利用可能性】
【０１０６】
本発明によれば、事前に観測、或いは予測した各マイクロホンの方向毎の感度等の先験情報を用いることなく、適応的に目的音方向にマイクロホンを向けた場合でも、目的方向の音源信号だけを明瞭に収音することができるので、マイクロホンアレイを用いた音源信号を収音する用途に有用である。
【符号の説明】
【０１０７】
１  方向別収音装置
１０  マイクロホン（指向性マイクロホン）
２０  ＤＦＴ処理部
３０  ＰＳＤ算出部
４０，４０‐ｍ，４０‐Ｍ  雑音抑圧部
４１  ウィーナーフィルタ構成部
４２  ウィーナーフィルタ処理部
４３  感度比推定部
５０  ＰＳＤ和推定部
５１  重み乗算部
５２  ＰＳＤ和算出部
６０  ＩＤＦＴ処理部
１００  ウィーナーフィルタ
４１１  ＣＳＤ算出部
４１２  ウィーナーフィルタ算出部
４１３  フィルタ乗算部
【図１】
2020010196.tif 000033
【図２】
2020010196.tif 000034
【図３】
2020010196.tif 000035