サンプルプログラム

プログラム例。
データを読んで窓をかけた後、フーリエ変換し、パワースペクトルを求める。

clear
cd '/exp/11/'
fget('fft_win.sci');

L=100000;
hanning = 0.5*cos(2*%pi*[0:L-1]/(L-1));
//pause
data=read("Seitai5.dat", 1000000, 1);
for i=1:10
data1=data(1+(i-1)*L:L*i);
data1=fft_win(hanning, data1);
powerspec1=power_sp(data1);
s(i)=sum(powerspec1);
c(i)=cent_sp(powerspec1)/s(i);
end

i=1:length(c);
subplot(211)
xtitle("", 'time(s)', 'power spectrum')
plot(i,c,'+')
//pause
//clf()
subplot(212)
xtitle("", 'time(s)', 'frequency')
plot(i,s,'+')
//plot(c,s)

//data2=data(L+1:2*L);
//data2=fft_win(hanning, data2);
//powerspec2=power_sp(data2);
//c2=cent_sp(powerspec2);
//s2=sum(powerspec2);

//xaxes=[0:length(powerspec1)-1];
//subplot(2,1,1)
//plot(xaxes, powerspec1)
//subplot(2,1,2)
//plot(xaxes, powerspec2)

//xtitle('', 'weight(kg)', 'EMG(V)');
//plot(x, value1, '+', x, value2, 'o', x, value3, '*')
//legend(['s';'m';'k'], 1)

//plot2d(t,y,1,rect=[xmin ymin xmax ymax])