在ＦＦＴ变换中出现虚假频率是怎么回事？
我做的是转子故障信号特征分析方面的课题，对原始信号进行高频和低频系数分解后，对高频和低频系数进行ＦＦＴ变换，但发现个别频谱线所对应频率过大，即产生虚假频率，这是怎么回事呢？下面是程序:
load ('C:\Program Files\MATLAB71\work\wd5.mat');
f=1280;
s=data(1:1024);
figure(1);
plot(s);
xlabel('时间t/s');
ylabel('幅值Amm');
title('原始信号');
=wavedec(s,5,'db3');
d5=wrcoef('d',c,l,'db3',5);
d4=wrcoef('d',c,l,'db3',4);
d3=wrcoef('d',c,l,'db3',3);
d2=wrcoef('d',c,l,'db3',2);
d1=wrcoef('d',c,l,'db3',1);
figure(2);
subplot(5,1,1);plot(d5);ylabel('d5');
title('高频系数d');
subplot(5,1,2);plot(d4);ylabel('d4');
subplot(5,1,3);plot(d3);ylabel('d3');
subplot(5,1,4);plot(d2);ylabel('d2');
subplot(5,1,5);plot(d1);ylabel('d1');
xlabel('时间t/s');
a5=wrcoef('a',c,l,'db3',5);
a4=wrcoef('a',c,l,'db3',4);
a3=wrcoef('a',c,l,'db3',3);
a2=wrcoef('a',c,l,'db3',2);
a1=wrcoef('a',c,l,'db3',1);
figure(3);
subplot(5,1,1);plot(a5);ylabel('a5');
title('低频系数a');
subplot(5,1,2);plot(a4);ylabel('a4');
subplot(5,1,3);plot(a3);ylabel('a3');
subplot(5,1,4);plot(a2);ylabel('a2');
subplot(5,1,5);plot(a1);ylabel('a1');
xlabel('时间t/s');
figure(4);
s=s-mean(s);
Y=fft(s,1024);
pyy=Y.*conj(Y)/1024;
plot((1:512)*f/1024,pyy(1:512));
xlabel('频率Hz');ylabel('幅值Amm');title('s的FFT频谱');
figure(5);
a4=a4-mean(a4);
Y=fft(a4,1024);
pyy=Y.*conj(Y)/1024;
subplot(3,1,1);plot((1:512)*f/1024,pyy(1:512));ylabel('幅值Amm');title('a4的FFT频谱');
Y=fft(d4,1024);
pyy=Y.*conj(Y)/1024;
subplot(3,1,2);plot((1:512)*f/1024,pyy(1:512));ylabel('幅值Amm');title('d4的FFT频谱');
Y=fft(d1,1024);
pyy=Y.*conj(Y)/1024;
subplot(3,1,3);plot((1:512)*f/1024,pyy(1:512));ylabel('幅值Amm');title('d1的FFT频谱');
xlabel('频率Hz');
下面是FFT变换后的频谱图:
[ 本帖最后由 lhy 于 2007-6-14 12:53 编辑 ]
回复 #1 lhy 的帖子
您能不能把图贴出来，最好带上您的主要程序，要不然这很难回答，一般情况下分解的时候容易出现虚假频率。 有些高频成分从时域看起来比较小,但是到了频谱可能很大. 局部放大了，你得看看它的幅值到底变大了没．有可能信号低频部分能量较大，高频小但是经过滤波把高低频分开了，实际上还是没改变只是视觉效果，数值没变．[ 本帖最后由 zhlong 于 2007-6-14 19:59 编辑 ] 大家能再说的清楚点么,谢谢 就是相同大小的加速度信号在低频时表现为位移大、振幅大，而在高频时表现为位移小，振幅较小。但是对它们做频谱分析，其在频谱里的值是一致的。
因此你分析得到和信号可能是真实存在的，不是虚假的。 楼主的图应和s的谱图相比较(figure4）。 这个可能跟滤波器的选择有关吧. 我现在也遇到这样的困惑，请教过测试老师，他说ＦＦＴ对等时间采样信号有可能出现虚假频率，如一个方波信号，ＦＦＴ后有１，３，５，７倍频率。但如何处理，他说需要采用其他分析方法？？
我还不知道，楼主解决了，在这里说一声。 最后一个图的y轴范围为10^4，较前两个小3个数量级。用相同的 y 轴范围，这个高频就不明显了。
页:
[1]