sig = wavread('doremi.wav');
sig = sig(:,1);
fs = 44100;                # サンプリング周波数 （CD と同じ）
nf = 8196;		   # フーリエ変換するデータ数
f = (0:nf/2)/nf*44100;
for i = 1:4410:(size(sig)-nf/2)   # 0.1 sec おきに測定
	sp = fft(sig(i:i+nf-1).*hamming(nf));
#	sp = fft(sig(i:i+nf-1));
	[x, ix] = max(20*log10(abs(sp(1:nf/2+1))));
	peak = ix/nf*44100;
#	disp(i)
	disp(peak)
endfor
