数字信号出来实验报告--matlab滤波器设计.docVIP

广 西 工 学 院 实 验 报 告 用 纸 实验名称 IIR数字滤波器的设计 实验成绩 指导老师 陈艳 系（院） 计算机工程系 班级 学号 学生姓名 一、实验目的 1、掌握双线性变换法和脉冲响应不变法设计IIR数字滤波器的具体设计方法及原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。 2、观察双线性变换及脉冲响应不变法设计滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。 3、熟悉Btterworth滤波器、Chebyshev滤波器和椭圆滤波器的频率特性。 二、实验原理与方法 1、冲击响应不变法 用数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲击响应ha(t)，让h(n)正好等于ha(t)的采样值，即：h(n)=ha(nT)，其中T为采样间隔，如果以Ha(s)及H(z)分别表示ha(t)的拉氏变换及h(n)的z变换，则 2、双线性变换法 S平面与Z平面之间满足以下映射关系： S平面的虚轴单值地映射于Z平面的单位圆上，S平面的左半平面完全映射到Z平面的单位圆内，双线性变换不存在混叠问题。双线性变换时一种非线性变换，这种非线性引起的幅频特性畸变可通过预畸而得到校正。 IIR低通、高通、带通数字滤波器设计采用双线性原型变换公式： 变换类型 变换关系式 备注 低通 高通 带通 为通带的上下边带临界频率 数字低通滤波器的设计步骤如下： 1、确定数字滤波器的性能指标：通带临界频率fp、阻带临界频率fr；通带内最大衰减ap；阻带内最小衰减ar；采样周期T； 2、确定相应的数字角频率， 3、计算经过预畸的相应模拟低通原型的频率：； 共 张 第 张 广 西 工 学 院 实 验 报 告 用 纸 实验名称 实验成绩 指导老师 系（院） 班级 学号 学生姓名 4、根据和计算模拟低通原型滤波器的阶数N，并求得低通原型的传递函数ha(s)； 5、用上面的双线性变换公式带入ha(s)，求出所设计的传递函数H(z)； 6、分析滤波器的特性，检查其指标是否满足要求。 三、实验内容及步骤 （1）、fp=0.3KHz,ap=0.8dB,fr=0.2KHz,ar=20dB,T=1ms;设计一Chebyshev高通滤波器；观察其是否满足要求。 （2）、fp=0.2KHz,ap=1dB,fr=0.3KHz,ar=25dB,T=1ms;分别用脉冲响应不变法和双线性变换法设计一Butterworth数字低通滤波器，观察其幅频特性曲线，比较两种方法的优缺点。 （3）、利用双线性变换法设计满足下列指标的Butterworth型、Chebyshev型和椭圆型滤波器。 fp=1.2KHz,ap0.5≤dB,fr=2KHz,ar≥40dB,fs=8KHz （4）、利用双线性变换法设计一Butterworth型数字滤波器，已知fs=30KHz,其等效模拟滤波器指标为 ap＜3dB,2KHz＜f≤3KHz,ar≥5dB,fr≥6KHz,ar≥20dB,f≤1.5KHz 实验程序及结果： （1）、fp=300;ap=0.8;fr=200;ar=20;T=0.001; fs=1/T; [n,wn]=cheb1ord(fp*2/fs,fr*2/fs,ap,ar); [b,a]=cheby1(n,ap,wn,high); [h,w]=freqz(b,a,512,fs); plot(w,abs(h));axis([0 500 0 1]);grid on; ---------------------------------------------------------------------------------------------------------------------------------- （2）、fp=200;ap=1;fr=300;ar=25; T=0.001; fs=1/T; wp=fp*2*pi/fs; wr=fr*2*pi/fs; wp1=2*fs*tan(wp/2); wr1=2*fs*tan(wr/2); [n,wn]=buttord(wp1,wr1,ap,ar,s); [z,p,k]=buttap(n); [bap,aap]=zp2tf(z,p,k); [b,a]=lp2lp(bap,aap,wn); [bz1,az1]=impinvar(b,a,fs); [h1,w1]=freqz(bz1,az1); plot(w1*fs/(2*pi),abs(h1),-r);hold on; [bz2,az2]=bilinear(b,a,fs); [h2,w2]=freqz(bz2,az2); plot(w2*fs/(2*pi),abs(h2),:b);grid on;