求大虾 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