Hilbert変換
後頭葉のα帯域の波形の大きさを包絡線envelopeを使って表示したいと思います。
包絡線を描くのに便利なツールとしてHilbert変換があります。
安静閉眼時のデータを読みます。
stacksize(50000000);
loadmatfile(xgetfile('*.mat'));
meg=diag(CalXRnge(1:306))*double(raw_data(1:306,1:6000));
getf('.../CalcSSP2.sce');
meg=CalcSSP2(meg);
time=((1:6000)-1)/CalXRnge($);
scf();plot2d(time,meg(235,:));//MEG2113
getf('.../FiltfiltIIR.sce');
x=FiltfiltIIR(meg(157,:)',time'*1000,4,'bp','ellip',[8,14],[1,1]);
scf();plot2d(time,x);
安静時の後頭部の10秒間の8-14Hzの波形です。Hilbert変換を使って包絡線envelopeを描きます。
以下の
hilbert.sceファイルを作成します。
行列には対応していません。
getf('.../hilbert.sce');
y=hilbert(x);
plot2d(time,abs(y));
h=gca();h.children(1).children.foreground=2;
これで包絡線envelope(青線)がかけました。