引用:
作者: lincolnlee0126
大侠们快快现身吧 小弟急等啊
|
fs=200;%设定采样频率,每秒钟采100个点数据
N=89; %采集数据样本数
n=0:N-1;
t=n/fs;
a=xlsread('NIUZHEN0105');
xx=a(:,1);
figure(1);
subplot(211);
plot(t,xx);%作时域波形
xlabel('时间(s)');
ylabel('微应变');
title('工况二第五道次上轴转矩时域波形');
grid;
%进行FFT变换并做频谱图
y=fft(xx,N);%进行fft变换
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
subplot(212);
plot(f,mag);%做频谱图
axis([3,100,0,6000]);
xlabel('频率(Hz)');
ylabel('幅值');
title('幅频谱图');
grid;
供参考