プロッターの作成2
学会発表用に整理します。
stacksize(50000000);
loadmatfile(xgetfile('*.mat'));
getf('../CalcSSP.sce');
MEG=ave_0001(:,1:306)*CalcSSP([PCA_0001,PCA_0002,...PCA_000n]);
MEG(:,3:3:306)=MEG(:,3:3:306)/5;
getf('../FiltfiltIIR.sce');
MEG=FiltfiltIIR(MEG,ave_0001(:,$),4,'bp','ellip',[0,5,40],[0.1,0.1]);
// 0.5〜40Hzの周波数帯域通過フィルターを使いました。
loadmatfile('../vv_lout.mat');
smp=size(MEG,1);
T=linspace(0,6,smp);
X=vv_lout(:,2)*ones(1,smp)+ones(1,smp)+ones(306,1)*T;
Y=vv_lout(:,3)*ones(1,smp)+MEG'/100*5;
scf();plot2d(X',Y');
xset('use color',0);//白黒表示をお勧めします
hf=gcf();hf.visible='off';
hf.children.axes.visible='off';//Figure.Axes
hf.children.margins=[0,0,0,0];//Figure.Axes
//hf.children.children(1).children(2).foreground=2;//青色表示
X=X(:,[1,$]);
Y=vv_lout(:,3)*ones(1,2);
plot2d(X',Y');//横軸表示
//hf.children.children(1).children(:).foreground=1;
for t=1:smp;if ave_0001(t,$)>=0;break;end;end;
X=vv_lout(:,2)*ones(1,2)+ones(306,1)*[T(t),T(t)];
Y=vv_lout(:,3)*ones(1,2)+ones(306,1)*[-100,100]/100*5;
plot2d(X',Y');//0msecの縦軸表示 100fTとしている
//hf.children.children(1).children(:).foreground=1;
for t=1:smp;if ave_0001(t,$)>=500;break;end;end;
X=vv_lout(:,2)*ones(1,2)+ones(306,1)*[T(t),T(t)];
Y=vv_lout(:,3)*ones(1,2)+ones(306,1)*[-50,50]/100*5;
plot2d(X',Y');//500msecの縦軸表示 100fTとしている
//hf.children.children(1).children(:).foreground=1;
hf.visible='on';
理屈では
hf.children.children(1)...とすれば色を変えられるはずですが、
カラー表示にすると絶望的な遅さになります。
改めてMATLABの偉大さがわかりました。