MATLAB程序处理的结果和其他软件处理的结果不一样
本帖最后由 luyong85 于 2013-1-2 22:27 编辑以下是我的matlab程序，实测的振动加速度数据见test1，请教一下程序哪里出现了问题，导致matlab处理的结果与其他软件处理的结果相差这么大。请指教！（频谱相差很大）
close all
clc
clear
format long
fs=500; %%%%%%%%%%%注意输入采样频率
freq_cutoff=40;%%%%%%%%%%%注意输入采样频率
delta=1/fs;%定义周期
%一下为设计滤波器，以对数据进行40Hz低通滤波
=butter(6,freq_cutoff/fs,'low'); %1/500交频butter巴特沃斯
load F:\matlab\matlab-temp\test1.txt;%导入数据
res=test1;
= 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=';%型成时间序列，中间的数据应与。。
data1=filter(b,a,x1);%对数据进行filter数字滤波
temp=;
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低通滤波频谱图http://www.ilovematlab.cn/data/attachment/forum/201301/02/092018yrpep01p116qxpms.jpg其他软件处理的40Hz低通滤波时域图http://www.ilovematlab.cn/data/attachment/forum/201301/02/092019s1yqsivobglgygsb.jpg matlab处理的40Hz一下信号时频图http://www.ilovematlab.cn/data/attachment/forum/201301/02/092019oqtl5qlptoh57oit.jpg
页:
[1]