MATLAB程序处理的结果和其他软件处理的结果不一样

luyong85

以下是我的matlab程序，实测的振动加速度数据见test1，请教一下程序哪里出现了问题，导致matlab处理的结果与其他软件处理的结果相差这么大。请指教！（频谱相差很大）
close all
clc
clear
format long
fs=500;        %%%%%%%%%%%注意输入采样频率
freq_cutoff=40;%%%%%%%%%%%注意输入采样频率
delta=1/fs;%定义周期
%一下为设计滤波器，以对数据进行40Hz低通滤波
[b,a]=butter(6,freq_cutoff/fs,'low'); %1/500交频butter巴特沃斯
load F:\matlab\matlab-temp\test1.txt;  %导入数据
res=test1;
[n0,m0] = size(res);
t0=res(:,1);
x1=res(:,2);
x2=res(:,5);
z1=res(:,4);
z2=res(:,7);
s1=size(x2);
N=s1(1,1);
L=length(x1);
time_input=[0:1/fs:(N-1)*(1/fs)]';%型成时间序列，中间的数据应与。。
data1=filter(b,a,x1);%对数据进行filter数字滤波
temp=[time_input,data1];
subplot(3,1,1)
plot(t0,x1,'black');
subplot(3,1,2);
plot(time_input,data1,'black');
NFFT = 2^nextpow2(L);%采样点数，为了满足2的整数倍，不足则在后面补零
Y1=fft(data1.*blackman(L),NFFT);%取值以符合采样定量的要求
Y=2*sqrt(Y1.*conj(Y1))/NFFT;
f0=fs/2*linspace(0,1,NFFT/2+1);
subplot(3,1,3)
plot(f0,abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

其他软件处理的40Hz低通滤波频谱图

                               
登录/注册后可看大图

其他软件处理的40Hz低通滤波时域图

                               
登录/注册后可看大图

matlab处理的40Hz一下信号时频图

                               
登录/注册后可看大图