中国科学技术大学电子工程与信息科学系多媒体通信实验室(Copyright2000) axis(0 M-1 0 1.1);ylabel(w n); subplot(2,2,3):stem(h):tite(实际冲击响应)； axis([0 M-1 -0.1 0.3]);ylabel('h[n]); subplot(2,2,4)plot(wpi,db),title('衰减幅度)； axis([0 1-100 101);ylabel('Decibles'); 附录C4R滤波器的实现 1、freqs函数：模拟滤波器的频率响应 [例系统传递函数为H)=02s+03s+1 的模拟滤波器，在MATLAB中可以用以下 s2+0.4s+1 程序来实现： a=10.41;b=[0.20.31 %设置分子分母的系数 w=logspace(-l,l)月 %产生从10-1到10'之间地0个等间距点，即50个颜率点 freqs(b,a,w) %根据输入的参数绘制幅度谱和相位谱 2、freqz函数：数字波器的频率响应 [例]系统传递函数为H()= 0.2+0.3z+z2 1+0.4z+22 的模拟滤波器，在MATLAB中可以用以下 程序来实现： a=[10.41]b=[0.20.31; %根据输入的参数绘制幅度谱和相位谱，得到0到π之间128个点处的顿率响应 freqz(b,a,128) 3、ButterWorth模拟和数字滤波器 (1)butterd函数：Butter Worth滤波器阶数的选择。 调用格式：[n,Wn]=butterd(Wp,Ws,Rp,Rs),在给定滤波器性能的情况下（通带临界频率 Wp、阻带临界频率Ws、通带内最大衰减Rp和阻带内最小衰减Rs),计算Butter Worth滤波 器的阶数n和截止频率Wn。相同参数条件下的模拟滤波器则调用格式为： [n,Wn]=butterd(Wp,Ws,Rp,Rs,'s') (2)butter函数：ButterWorth滤波器设计。 调用格式：[b,a=butter(n,Wn),根据阶数n和截止频率Wn计算ButterWorth滤波器分 子分母系数(b为分子系数的矢量形式，a为分母系数的矢量形式)。相同参数条件下的模拟 滤波器则调用格式为：[b,a=butter(n,Wn,'s) [例]采样频率为1Hz,通带临界频率=0.2Hz,通带内衰减小于1dB(ap=1):阻带 如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系network@ustc.cdu.cn cxh@ustc.cdu.cn中国科学技术大学电子工程与信息科学系 多媒体通信实验室 （Copyright 2000） axis([0 M-1 0 1.1]);ylabel('w[n]'); subplot(2,2,3);stem(h);title('实际冲击响应'); axis([0 M-1 -0.1 0.3]);ylabel('h[n]'); subplot(2,2,4);plot(w/pi,db); title('衰减幅度'); axis([0 1 -100 10]);ylabel('Decibles'); 附录 C4 IIR 滤波器的实现 1、freqs 函数：模拟滤波器的频率响应 [例] 系统传递函数为 0.4 1 0.2 0.3 1 ( ) 2 2 + + + + = s s s s H s 的模拟滤波器，在 MATLAB 中可以用以下 程序来实现： a=[1 0.4 1]; b=[0.2 0.3 1]; %设置分子分母的系数 w=logspace(-1,1); %产生从 10−1 到101 之间地 0 个等间距点，即 50 个频率点 freqs(b,a,w) %根据输入的参数绘制幅度谱和相位谱 2、freqz 函数：数字滤波器的频率响应 [例] 系统传递函数为 2 2 1 0.4 0.2 0.3 ( ) − − + + + + = z z z z H z 的模拟滤波器，在 MATLAB 中可以用以下 程序来实现： a=[1 0.4 1];b=[0.2 0.3 1]; %根据输入的参数绘制幅度谱和相位谱，得到 0 到π之间 128 个点处的频率响应 freqz(b,a,128) 3、ButterWorth 模拟和数字滤波器 （1）butterd 函数：ButterWorth 滤波器阶数的选择。 调用格式：[n,Wn]=butterd(Wp,Ws,Rp,Rs)，在给定滤波器性能的情况下（通带临界频率 Wp、阻带临界频率 Ws、通带内最大衰减 Rp 和阻带内最小衰减 Rs），计算 ButterWorth 滤波 器的阶数 n 和截止频 率 Wn 。 相 同参数条 件下的模拟 滤波器则 调用格式为 ： [n,Wn]=butterd(Wp,Ws,Rp,Rs,’s’) （2）butter 函数：ButterWorth 滤波器设计。 调用格式：[b,a]=butter(n,Wn)，根据阶数 n 和截止频率 Wn 计算 ButterWorth 滤波器分 子分母系数（b 为分子系数的矢量形式，a 为分母系数的矢量形式）。相同参数条件下的模拟 滤波器则调用格式为：[b,a]=butter(n,Wn,’s’) [例] 采样频率为 1Hz，通带临界频率 fp =0.2Hz，通带内衰减小于 1dB（αp=1）；阻带 如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ustc.edu.cn cxh@ustc.edu.cn