b = fir1(n,wn,window) window:指定窗函数 b = fir1(n,wn,'ftype',window) 注意：设计高通和带阻FIR滤波器时，阶数n只能取偶数以保证h(n)长度为n+1为奇数。但当用户将n设置为奇数时，fir1函数会自动对n加1。 例7.2.1 例7.2.2 例7.2.3 凯塞窗例1 凯塞窗例3
2应用MATLAB进行频率采样法设计FIRDF 例7.3.1 例7.3.2
13
实验七 综合实验-数字滤波器设计
一、实验目的
1.掌握FIR数字滤波器的综合设计 2.掌握IIR数字滤波器的综合设计
要求：分别设计FIR数字滤波器和IIR数字滤波器。
一个含有5Hz、20Hz和30Hz的混和正弦波信号，设采样频率为fs=150Hz，阻带衰减至少为50dB，绘制滤波前和滤波后的频谱图进行比较，要求如下： （1）只保留5Hz的正弦信号，即滤除20Hz和30Hz正弦信号； （2）只保留30Hz的正弦信号，即滤除5Hz和20Hz正弦信号； （3）只保留20Hz的正弦信号，即滤除5Hz和30Hz正弦信号； （4）保留5Hz和30Hz的正弦信号，即滤除20Hz正弦信号。
14
实验八 时域离散系统的实现
一、实验目的
1.掌握FIR数字滤波器的结构 2.掌握IIR数字滤波器的结构 3.掌握数字滤波器的格型结构 二、基本内容
1 FIR数字滤波器的结构 MATLAB函数：
直接型到级联型的转换：[sos,g] = tf2sos(b,1) 级联型到直接型的转换：[b,a]=sos2tf(sos,g)
B(z)b1?b2z?1?...bn?1z?n H(z)??A(z)a1?a2?...?am?1z?m
b0k?b1kz?1?b2kz?2 H(z)?g?Hk(z)?g??1?2?a2kzk?1k?11?a1kzLL例1：FIR数字滤波器的系统函数为：H(z)?1?16（1）直接型结构； （2）线性相位结构； （3）级联型结构 2 IIR数字滤波器的结构 MATLAB函数：
15
1?4z?z?8，画出： 16直接型到级联型的转换：[sos,g] = tf2sos(b,a) 级联型到直接型的转换：[b,a]=sos2tf(sos,g)
B(z)b1?b2z?1?...bn?1z?n H(z)???mA(z)a1?a2?...?am?1z
b0k?b1kz?1?b2kz?2 H(z)?g?Hk(z)?g??1?21?az?azk?1k?11k2kLL直接型到并联型的转换：[C,B,A]=dir2par(b,a)
523?z?1?z?233例2：已知某三阶数字滤波器的系统函数为：H(z)?，试画出其直1?11?11?2(1?z)(1?z?z)322接型、级联型和并联型结构。 3.数字滤波器的格型结构 (1) 全零点格型结构
直接型到FIR格型的转换函数：k= tf2latc(b) 例8.4.1
FIR格型到直接型的转换函数：b= latc2tf(k) 例8.4.2 （2）全极点格型结构 k = tf2latc(1,a)
例4：全极点IIR滤波器的系统函数为：H(z)?1，求其格型结构。
13?15?21?31?z?z?z2483（3）零极点格型梯形结构 [k,v]=tf2latc(b,a)
16
1?2z?1?2z?2?z?3例5：零极点IIR滤波器的系统函数为：H(z)?，求其格型梯形结构。
13?15?21?31?z?z?z2483
17
