ML

図6.58:MLの接続例 ノードの入出力とプロパティ

入力

INPUT FRAMES :Matrix<complex<float> >型．マルチチャネル複素スペクトル．行がチャネル，つまり，

各マイクロホンから入力された波形の複素スペクトルに対応し，列が周波数ビンに対応する．

INPUT SOURCES :Vector<ObjectRef>型．音源定位結果等が格納されたSource型オブジェクトのVector配 列である．典型的には，SourceTrackerノード，SourceIntervalExtenderノードと繋げ，その出力を用いる．

INPUT NOISE CM :Vector<ObjectRef>型．INPUT SOURCESと同じSource型オブジェクトのVector配 列である．雑音方向の情報のオプション入力である．

出力

OUTPUT :Map<int, ObjectRef>型．分離音の音源IDと，分離音の1チャネル複素スペクトル (Vector<complex<float> >型)のペア．

パラメータ

LENGTH :int型．分析フレーム長[samples]．前段階における値（AudioStreamFromMic，MultiFFTノード など）と一致している必要がある．デフォルト値は512[samples]．

ADVANCE :int型．フレームのシフト長[samples]．前段階における値（AudioStreamFromMic，MultiFFTノー ドなど）と一致している必要がある．デフォルト値は160[samples]．

SAMPLING RATE :int型．入力波形のサンプリング周波数[Hz]．デフォルト値は16000[Hz]．

DECOMPOSITION ALGORITHM : 型．音源分離で用いる演算アルゴリズムの選択．GEVDは一般化固有 値分解を，GSVDは一般化特異値分解を表す．GEVDはGSVDに比べて雑音抑制性能が良好だが計算時 間がかかる．使用目的や計算機環境に応じて演算アルゴリズムの使い分けができる．

ALPHA :float型．フィルタ更新係数．デフォルト値は0.99．

ENABLE DEBUG :bool型．デフォルトはfalse．trueが与えられると,分離状況が標準出力に出力される．

表6.48: MSNRで利用するパラメータ

パラメータ名 型 デフォルト値 単位 説明

LENGTH int 512 [pt] 分析フレーム長

ADVANCE int 160 [pt] フレームのシフト長

SAMPLING RATE int 16000 [Hz] サンプリング周波数

DECOMPOSITION ALGORITHM string GEVD 演算アルゴリズム

ALPHA float 0.99 フィルタ更新係数．

ENABLE DEBUG bool false デバッグ出力の可否

ノードの詳細

技術的な詳細: 基本的に詳細は下記の参考文献を参照されたい．

音源分離概要: 音源分離問題で用いる記号を表6.49にまとめる．演算はフレーム毎に周波数領域において行わ れるため，各記号は周波数領域での，一般には複素数の値を表す．音源分離はK個の周波数ビン(1≤k≤K) それぞれに対して演算が行われるが，本節ではそれを略記する．N,M, f をそれぞれ，音源数，マイク数，フ レームインデックスとする．

表6.49:変数の定義

変数 説明

S(f)=[

S₁(f), . . . ,S_N(f)]T

f フレーム目の音源の複素スペクトル X(f)=[

X1(f), . . . ,XM(f)]T

マイクロホン観測複素スペクトルのベクトル．INPUT FRAMES入力に対応．

N(f)=[

N₁(f), . . . ,N_M(f)]T

加法性雑音

H=[H₁, . . . ,H_N]∈C^M×N 1≤n≤N番目の音源から1≤m≤M番目のマイクまでの伝達関数行列   K(f)∈C^M×M 既知雑音相関行列

W(f)=[W₁, . . . ,W_M]∈C^N×M 分離行列 Y(f)=[

Y1(f), . . . ,YN(f)]T

分離音複素スペクトル 音のモデルは以下の一般的な線形モデルを扱う．

X(f) = HS(f)+N(f) (6.63)

分離の目的は，

Y(f) = W(f)X(f) (6.64)

として，Y(f)がS(f)に近づくように，W(f)を推定することである．

最尤法に基づく分離行列WMLは次式であらわされる．

WML(f) = K˜⁻¹(f)H

H^HK˜⁻¹(f)H (6.65)

ここで，

K(˜ f) = K(f)+||K(f)||FαI (6.66) であり，ここで||K(f)||Fは既知雑音相関行列K(f)のフロベニウスノルム，αはパラメータREG FACTOR，

Iは単位行列である.

トラブルシューティング: 基本的にはGHDSSノードのトラブルシューティングと同じ．

参考文献

[1 ] F. Asano: ’Array signal processingfor acoustics —Localization, tracking and separation of sound sources—, The Acoustical Society of Japan, 2011.