mscohereは welchのオーバーラップ付き平均ピリオドグラムを用いてcoherenceを推定します。
例えば、welchのオーバーラップ付き平均ピリオドグラムは、時系列信号に対して、
pwelch(x, hann(256), 128, 256, fs);
と定義した場合、
x(1:256) の信号に hann窓を掛けて、256点でFFT
x(129:384) の信号に hann 窓を掛けて、256点でFFT
x(257:512) の信号に hann 窓を掛けて、256点でFFT
・・・・・・・・・・・・・・・・
128サンプルオーバーラップさせながら区分ごとに窓関数付きfft を実行し、最後に平均化して、パワースペクトル密度 pxx が得られることになります。
オーバーラップの値を窓のサイズに近づければ近づけるほど、少しずつ窓をスライドしていくことになりますので平均化の回数が多くなります。
pwelch(x, hann(256), 128, 512, fs);
とすると、 256点の窓に対して 更に 256点分の零値をパディングして 512点でFFTを行うことになります。
※ 上記は pwelch の例となりますが、mscohere についても概念は同じことになります。
