Matlab实现IIR滤波器：低通、高通、带通、带阻设计

说明：（1）为了使滤波器阶数尽可能低，每个滤波器的边界频率选择原则是尽量使滤波器过渡带宽尽可能宽。
（2）与信号产生函数mstg相同，采样频率Fs=10kHz。
（3）为了滤波器阶数最低，选用椭圆滤波器。（之后，再依次实现巴特沃斯、切比雪夫1、切比雪夫2数字滤波器）
按照图2 所示的程序框图编写的实验程序为exp1.m。
2、实验程序清单
%实验1程序exp1.m
% IIR数字滤波器设计及软件实现
clear all;close all
Fs=10000;T=1/Fs;   %采样频率
%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st 
st=mstg;
%低通滤波器设计与实现=========================================
fp=280;fs=450;
wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;   %DF指标（低通滤波器的通、阻带边界频）
[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp
[B,A]=ellip(N,rp,rs,wp);      %调用ellip计算椭圆带通DF系统函数系数向量B和A
y1t=filter(B,A,st);     %滤波器软件实现
% 低通滤波器设计与实现绘图部分
figure(5);
subplot(2,1,1);
myplot(B,A);  %调用绘图函数myplot绘制损耗函数曲线
yt='y_1(t)';
subplot(2,1,2);
tplot(y1t,T,yt); %调用绘图函数tplot绘制滤波器输出波形
%带通滤波器设计与实现====================================================
fpl=440;fpu=560;fsl=275;fsu=900;
wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];rp=0.1;rs=60; 
[N,wp]=ellipord(wp,ws,rp,rs);    %调用ellipord计算椭圆DF阶数N和通带截止频率wp
[B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和A
y2t=filter(B,A,st);     %滤波器软件实现