求大虾 FFT变换 600MW汽轮机真实故障数据
从电厂bently上采集来的数据,文件名0-12.txt,晚上零点到中午12点的数据,跑了下面的程序:权限不够,只能传200KB的,我将数据文件分拆一下:
m=load('0-12.txt');
N=length(m);
fs=3200;%设定采样频率
n=0:N-1;
for i=1:N;
q(i,1)=m(i,2);
end;
d=q';
s0=d;
figure(1);
plot(1:N,s0(1:N));
title('原始信号');
grid;
%进行FFT变换并做频谱图
y=fft(s0);%进行fft变换
mag=abs(y);
f=(0:length(y)-1)*fs/length(y);
figure(2);
plot(f,mag);%做频谱图
%axis();
xlabel('频率(Hz)');
ylabel('幅值');
title('幅频谱图');
grid;
频域图和时域图如下,
时域图跟bently采集的趋势图一致,不明白为什么自己FFT变换的频域图,会成上面的样子,求高手解决!多谢!
频域应该表现在50HZ,100HZ和200HZ振幅大的,结果频域图完全不跟想象的一样,想知道哪里出问题了,求高手!
不好意思,时域图:
C:\Users\nanyi\Desktop\时域.jpg 怎么没人回帖……求助!! 会FFT的 进来教教我啊…… 你的原始数据有问题,请把第一列去掉,在做傅立叶变换 这是倒频分析的吧? 不好意思呀,我用你的软件分析也C:\Documents and Settings\Administrator\桌面是fft分析 本帖最后由 impulse 于 2011-4-29 09:32 编辑
当然不一样,它的频谱并不是取你看到的8周期(1024点)振动波形数据进行计算的,它的频谱数据在采集阶段就完成计算的,样本点也比1024多。但你这个频谱图也太离谱了吧。 不知道你用的是408还是208,如果真感兴趣,用408采集Raw Waveform,然后取出来进行信号分析,如果算法没问题,可以得到和它差不多的值。 回复 4 # 独孤一1984 的帖子
你那哪是时域图呀,趋势图吧,dt为12秒? 回复 10 # impulse 的帖子
我这里还有一套在bently实验台上的碰磨数据,采集来的是毫伏值,没有对应到振幅(据说是大致的线性对应),但是出来的趋势图(1024个点)确实是削波曲线,还是频域出问题了,出来的频域图里的振动大征兆是以32Hz为基频的,但是数据采集的人说当时的转速是3000r, 怎么会出现这种错误呢? 谢谢楼主了 谢谢楼主,先把数据下回去分析! 学习学习下 学习学习下
页:
[1]
2