自编两种滤波器结果比较函数----基于MATLAB

    该函数主要作用是，对同一段信号，用两种不同的滤波器进行滤波处理，然后比较两者在频域的结果。通过比较，可以选择合适的滤波器。这个函数中是对butter和fir两种进行了比较，读者要是想对比其他滤波器，或者同时比较三种及以上滤波器，只需要对代码稍加修改，就可以满足需求，这里不加赘述。
   还有一点，读者要是想输出不同滤波器的结果，只需将下述函数的第一行修改为：
   function  [filter1，filter2] = myfilter(Signal, Fs)   
   即可

function  myfilter(Signal, Fs)

%Signal----原始信号
%filter1---使用butter滤波器的滤波结果
%filter2---使用fir1滤波器的滤波结果
 
[b,a] = butter(5,[0.5*2/Fs 40*2/Fs],'bandpass');  %建立butter带通滤波器
filter1 = filtfilt(b,a,Signal);  %利用已建立滤波器进行滤波
NFFT1 = length(filter1);  %进行FFT的信号长度  
frequency1 = Fs*(0:NFFT1-1)/NFFT1;    %频率
filter1_FFT = abs(fft(filter1,NFFT1));   %FFT结果

[b1,a1] = fir1(20,[0.5*2/Fs 45*2/Fs],'bandpass');  %建立butter带通滤波器
filter2 = filtfilt(b1,a1,Signal);    %利用已建立滤波器进行滤波
NFFT2 = length(filter2);  %进行FFT的信号长度
frequency2 = Fs*(0:NFFT2-1)/NFFT2;   %频率
filter2_FFT = abs(fft