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

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);
=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

把数据上传下

数据在文件中,最后一列

华电机械 发表于 2013-5-20 15:47

王瀚 发表于 2010-3-28 22:17 static/image/common/back.gif
数据在文件中,最后一列

您能给我发一下STFT的工具包吗?现在急用,,在此感激!邮箱:354773461@qq.com
页: [1]
查看完整版本: STFT