数字信号实验报告殷光.docVIP

测试信号处理实验报告 班级：研097 姓名：殷光 学号：2009003018 实验一用脉冲响应不变法设计IIR滤波器 实验目的： 1.1掌握脉冲响应不变法设计IIR数字滤波器的原理及具体设计方法，熟悉实现这一设 计的计算机编程。 1.2观察用脉冲响应不变法设计的滤波器的时域特性和频域特性，比较所设计的数字滤 波器和相应模拟滤波器的频域特性，观察脉冲响应不变法设计屮产生的频率混淆现象。 1.3观察参数变化对滤波器性能的影响，通过了解脉冲响应不变法的优缺点，熟悉该方 法的应用范围。 实验原理与方法： 脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。它特别适用于那些 对滤波器的时域特性有一定要求的场合。具体地说：它可以保证所设计的IIR滤波器的脉冲 响应和相应模拟滤波器的冲激响应在采样点上完全一致。 一个模拟滤波器的传递函数可以川有理分式表示为： CIO + Cl\S +... + /?() +bs + ... + bNsN 通过反拉氏变换，我们就可以得到它的冲激响应h⑺，脉冲响应不变法就是要保持脉冲响 应不变，即 h⑻=ha(nT) 对上式中的冲激响应序列h(n)作z变换，就可以得到数字滤波器的传递函数H(z)。 一般来说，的分母多项式阶次总是大于分子多项式的阶次。数字滤波器的传递函 数H(z)经过合并简化，成为一般形式的有理分式传递函数 Co+AZ +“.+ cLz dQ + dz 1 +...4-dNz N 在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中，我们已经知道 H(z)\_=eXT = ^XHa(s^j^K) 按照z=/7’的关系，每一个S平面上宽度为2;r/T的水平条带将重迭映射到z平面上。因此 脉冲响应不变法将s平而映射到z平而，不是一个简单的一一对应的关系。对于髙采样频率 (T小)的情况，数字滤波器在频域可能有极高的增益。为此我们采用 N/=1 N /=1 N h(n) = Th(l(nT) = 在脉冲响应不变法设计屮，模拟频率与数字频率之间的转换关系是线性的 (^y = Q*Do同时，它可以保特脉冲响应不变，h(n)=/^nT)。因此，这一方法往往用于低 通吋域数字滤波器设计及相应的模拟系统数字仿真设计。 实验内容及步骤： 3.1复习脉冲响应不变法设计I1R滤波器和模拟滤波器没计的有关内界，认真阅读本实验 原理与方法的A容，熟悉脉冲响应不变法的设计步骤。 3.2编制一个用脉冲响应不变法设计IIR滤波器的程序。 3.3上机实验内容：完成以下两组数字滤波器的设计 采样频率f s=4Hz，周期T=0.25s，设计一个叫阶Butterworth低通滤波器，其3dB截止 频率 fe=lHz。 采样频率为1Hz,设计一个数字低通滤波器，要求其通带临界频率fp=0.2Hz，通带内 衰减小于ldB，阻带临界频率fs=0.3Hz,阻带袞减大于25dB。 用Butterworth逼近设计该滤波器。 用Chebyshev逼近设计该滤波器。 实验中绘制的幅频特性曲线是在0?31之间对幅频特性的等距采样，采样点数可选为 50。 实验主要程序(matlab): clear all;close all;clc; T=0.25;Fs=4;wc=2*pi:Omegac=wc/T;N=4;t=[0:0.01:3J;n=[0:l:3/T]; [b，a]=butter(N，Omegac，’s’); [c,d]=impinvar(b,a,T); rh,w]=freqs(b,a); fha，x，t|=impulse(b，a); mag = abs(h);mag = 20*logl0(mag); figure(l);impz(ha); title(’h(n)’);ylabel(’幅度’);grid on; figure(2) plot(t，ha);title(单位脉冲响应 ha⑴’);grid on; figure(3) subplot(l，2,l);plot(w，mag);title(?幅度响应 Ha(j/Omega));grid on; [h，w]=freqz(b，a);mag = abs(h);mag = 20*logl0(mag); subplot(l，2,2);plot(w，mag);title(’幅度响应 H(eAjw),);grid on; %用 Butterworth 设计 Fs=l ;Omegap=0.4;Omegar=0.6;Ap=l;Ar=25;T=l/Fs; [n,Wn] = buttord(Omegap,Omegar,Ap,Ar) [cs，ds]= butter(n,Wn); [b,a]=impinvar(cs,ds,T); [h,w]=freqz(b,a); mag = abs(h);db = 20*