matlab FIR IIR 滤波器

主要是FIR IIR 滤波器的matlab实现

function y = zach_fir_bandpass_filter（x fs1 fp1 fp2 fs2 ap as Ft）
%Band Pass Filter
% ap = 1; as = 40;
% fp1 = 900;%通带截止频率；
% fp2 = 1300;%通带截止频率；
% fs1 = 600;%阻带截止频率
% fs2 = 1500;%阻带截止频率
% Ft = 4000;%采样率

%归一化角频率边界频率数字频率
wp1 = 2*pi*fp1/Ft;
wp2 = 2*pi*fp2/Ft;
ws1 = 2*pi*fs1/Ft;
ws2 = 2*pi*fs2/Ft;
delta_w = wp1 - ws1;%过渡带宽度
%求出M
[N Wn Bta filtype] = kaiserord（[ws1 wp1 wp2 ws2]/pi    [0 1 0]    [10^（-as/20）  1-10^（-ap/20）  10^（-as/20）]）
%设计滤波器
b = fir1（N Wn filtype hann（N+1） ‘noscale‘）;%用于与输入信号卷积
[hw] = freqz（b 1 512）;
subplot（311）;plot（ w*Ft/（2*pi）  20*log10（abs（h）） ）;
xlabel（‘Hz‘）; ylabel（‘dB‘）;

L = 512;
X = fft（x512）;%x为行向量
subplot（312）;plot（ Ft*（-L/2:L/2-1）/L abs（fftshift（X）））;
%subplot（312）;plot（（-L/2:L/2-1）*（2*pi/L） abs（fftshift（X）））;
title（‘原信号‘）;
xlabel（‘Hz‘）; ylabel（‘Mag‘）;

y = filter（b 1 x）;
subplot（313）;plot（  Ft*（-L/2:L/2-1）/L abs（fftshift（fft（y512））））;
title（‘高通滤波后的信号‘）;
xlabel（‘Hz‘）; ylabel（‘Mag‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5252  2017-05-04 07:38  tools\APPLE_chendanling2_女_20_1.txt

     文件       1041  2017-05-30 18:32  tools\zach_fir_bandpass_filter.asv

     文件       1072  2017-05-30 18:34  tools\zach_fir_bandpass_filter.m

     文件       1004  2017-05-30 18:34  tools\zach_fir_highpass_filter.m

     文件        992  2017-05-30 18:17  tools\zach_fir_lowpass_filter.asv

     文件       1001  2017-05-30 18:36  tools\zach_fir_lowpass_filter.m

     文件       1704  2017-05-30 18:39  tools\zach_iir_bandpass_filter.m

     文件       1308  2017-05-30 18:37  tools\zach_iir_highpass_filter.asv

     文件       1328  2017-05-30 18:39  tools\zach_iir_highpass_filter.m

     文件       1085  2017-05-30 17:17  tools\zach_iir_lowpass_filter.asv

     文件       1301  2017-05-30 18:39  tools\zach_iir_lowpass_filter.m

     目录          0  2017-05-30 18:37  tools

----------- ---------  ---------- -----  ----

                17088                    12