STFT分析图说明

16443

我对一组振动信号做STFT变换后，得到下边这个时频分布图，对图的意义不太了解，请指点一下。

zhlong

估计是你参数设置的问题，主要是窗宽和重叠。
另外你的STFT程序能否贴上来看看。

16443

clear
clc
A=textread('result.txt'); %read data
AA=A(:,2);
sig=AA(2:length(AA),1);
%For a faster computation, N should be a power of two
N=2048; % freqency number
M=length(sig);
if (N>M)
    error('sampling number N must smaller signal length M')
end
figure(1)
plot((1:length(sig)),sig)
title('Time-Amplitude');
xlabel('Time T');
ylabel('Amplitude A');
h = window(@gausswin,17);  %gausswin
[tfr,t,f]=tfrstft(sig,1:N,N,h); %time freqency
figure(2)
contour(t,f,abs(tfr))
title('time-freqency 17-dian gausswin');
xlabel('Time T');
ylabel('Freqency f');

songzy41

1，有可能楼主取的窗函数：
h = window(@gausswin,17);
窗长17太短，加长试试看；
2，从图中看能量都集中在0频附近，说明有较大的直流分量，建议楼主在做STFT之前先消除直流分量。

zhangnan3509

我也认为有相同的问题，17个跟2048个比起来是有点太小了

zhangnan3509

t = 0:.0001:.0255;
x1 = sin(2*pi*80*t);
x2= 0.8*exp(-(t-0.01875).^2*1000000).*sin(2*pi*1600*t);
%subplot 211; plot(t,x1,'r',t,x2,'b');
axis([0 0.0255 -1 1]); grid;
x=x1+x2;
%subplot 212; plot(t,x);
axis([0 0.0255 -1 1]); grid;
%figure
h=window('hamming',255);
[tfr,tt,f]=tfrstft(x',1:256,256,h);
imagesc(t,f(1:64)*20000,abs(tfr(1:64,:)));
axis('xy');
xlabel('Time'); ylabel('Frequency');
（你也可以看看songzy41 先生的例子）

16443

好的，我看看这个例子，顺便取大窗函数点数。
另外：用contour函数画出的图contour(t,f,abs(tfr))，图中的颜色值代表什么？

zhangnan3509

我是这样理解的，相同的颜色的图线可以在取相同的值时得到的不同的结果。类似于等势线。详细的请看 help里关于Contour的解释

[ 本帖最后由 zhangnan3509 于 2007-7-1 10:54 编辑 ]

16443

我把数据文件传上来，你帮我看看，为啥这个离散数据的信号时频处理后，频率基本集中在0附近。
这个信号是通过软件计算，然后按照等时间间隔t＝2E－4提取得到的。

[ 本帖最后由 16443 于 2007-7-1 12:48 编辑 ]

zhlong

　　楼主的信号本身频率就很低，接近0HZ。这个我们可以从图1楼主的信号时域图形来看，信号有6250个点，却只有大概10个周期。假定采样间隔为1s，那么一个周期为625s，信号频率为1/625=0.0016Hz。接近0。

　　上面的推导可以通过下面这个仿真信号来验证。
　　t=1:6250;
　　x=sin(2*pi*0.0016*t);
　　plot(x)

16443

谢谢，你的解答使我明白了为啥时频分布主要集中在0Hz附近。看来我的信号在采集中存在问题，需要重新来过 :@L

[ 本帖最后由 zhlong 于 2007-7-1 15:26 编辑 ]

zhlong

你做故障诊断的吗？你这是仿真信号？还是实验采集的？

16443

我这是个仿真信号 ，按时间等步长计算得到

zhlong

做的是什么故障的仿真呢？

16443

不是，做了一个简单悬臂梁（不是单一的梁）的动态响应。
在matlab中怎样间隔n行（比如n＝10），提取某一列的值组成一个新的列 ？