STFT

王瀚

load('default.txt');
x=default(:,4);
figure(1);
fs=200;
n=length(x);
t=0:1/fs:(n-1)/fs;
plot(t,x);
hold on;
%nfft=2^nextpow2(n);
[tfr,t,f]=tfrstft(x,1:length(x),256);%语句中512也可以改为256、1024等，越大的话tfr行数（列数）越多，图形显示困难，建议512或1024，不能取信号一样长度，要不然一般计算机显示了时频图象，就会出现
%% 时频表示
figure(2);
subplot(211);imagesc(abs(tfr));
subplot(212);imagesc(angle(tfr));
xlabel('时间t');
ylabel('频率f');

问题：数据特征，5HZ基频 50HZ高频。但上列程序无法识别该特征。为什么？？

王瀚

数据在文件中，最后一列

华电机械

王瀚 发表于 2010-3-28 22:17

                               
登录/注册后可看大图

数据在文件中，最后一列

您能给我发一下STFT的工具包吗？现在急用，，在此感激！  邮箱：354773461@qq.com