scilabでお手軽に周波数解析と時間周波数解析

昨日のエントリーでscilabで周波数解析をする例を紹介したが、デモを見たところ既に用意されているようだ。

analyze //周波数解析
mapsound //時間周波数解析


//テストプログラム
Fs=22050;
Tl=1;
t=soundsec(Tl,Fs);

f1=1000;
s=(sin(2*%pi*t.*(f1+t*f1)));

scf
subplot(211)
analyze(s,100,5000,Fs);
subplot(212)
mapsound(s,0.05,100,5000,1,Fs);

かなり簡易的なもので、特にmapsoundの方は実用には耐えないと思われる。