|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
从电厂bently上采集来的数据,文件名0-12.txt,晚上零点到中午12点的数据,跑了下面的程序:
权限不够,只能传200KB的,我将数据文件分拆一下:
0-2.txt
(160.26 KB, 下载次数: 98)
3-5.txt
(159.98 KB, 下载次数: 54)
6-8.txt
(159.93 KB, 下载次数: 39)
9-11.txt
(160 KB, 下载次数: 51)
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([0,1000,0,300]);
xlabel('频率(Hz)');
ylabel('幅值');
title('幅频谱图');
grid;
频域图和时域图如下,
时域图跟bently采集的趋势图一致,不明白为什么自己FFT变换的频域图,会成上面的样子,求高手解决!多谢!
频域应该表现在50HZ,100HZ和200HZ振幅大的,结果频域图完全不跟想象的一样,想知道哪里出问题了,求高手!
|
|