|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我的采样数据为发动机振动的加速度,为何通过MATLAB后的加速度的时域图与FFT后的频谱图幅值相差很大?谢谢
data = load('C:\Documents and Settings\Administrator\桌面\盛鹏程毕业论文\2008.07.16\加速度1.txt');
dt = 0.0005;
x1 = data(:,1);
x2 = data(:,2);
x3 = data(:,3);
x4 = data(:,4);
N1 = length(x1);
N2 = length(x2);
N3 = length(x3);
N4 = length(x4);
t1 = [0:N1-1]*dt;
t2 = [0:N2-1]*dt;
t3 = [0:N3-1]*dt;
t4 = [0:N4-1]*dt;
figure(1)
subplot(4,1,1),plot(t1,x1);
xlabel('时间/s');ylabel('振幅');legend('前左加速度时域图');
xlim([0 0.5]);
ylim([-30 30]);
grid on;
subplot(4,1,2),plot(t2,x2);
xlabel('时间/s');ylabel('振幅');legend('前右加速度时域图');
grid on;
xlim([0 0.5]);
ylim([-30 30]);
subplot(4,1,3),plot(t3,x3);
xlabel('时间/s');ylabel('振幅');legend('后左加速度时域图');
grid on;
xlim([0 0.5]);
ylim([-30 30]);
subplot(4,1,4),plot(t4,x4);
xlabel('时间/s');ylabel('振幅');legend('后右加速度时域图');
grid on;
xlim([0 0.5]);
ylim([-30 30]);
figure(2)
y1 = fft(x1);
y2 = fft(x2);
y3 = fft(x3);
y4 = fft(x4);
mag1 =abs(y1)*2/N1;
f1 = [0:N1-1]/(N1*dt);
subplot(4,1,1),plot(f1(1:N1/2),mag1(1:N1/2));
xlabel('频率/Hz');ylabel('振幅');legend('前左加速度幅频图');
xlim([0 300]);
grid on;
mag2 =abs(y2)*2/N1;
f2 = [0:N2-1]/(N2*dt);
subplot(4,1,2),plot(f2(1:N2/2),mag2(1:N2/2));
xlabel('频率/Hz');ylabel('振幅');legend('前右加速度幅频图');
xlim([0 300]);
grid on;
mag3 =abs(y3)*2/N3;
f3 = [0:N3-1]/(N3*dt);
subplot(4,1,3),plot(f3(1:N3/2),mag3(1:N3/2));
xlabel('频率/Hz');ylabel('振幅');legend('后左加速度幅频图');
xlim([0 300]);
grid on;
mag4 =abs(y4)*2/N4;
f4 = [0:N4-1]/(N4*dt);
subplot(4,1,4),plot(f4(1:N4/2),mag4(1:N4/2));
xlabel('频率/Hz');ylabel('振幅');legend('后右加速度幅频图');
xlim([0 300]);
grid on; |
|