第七章 有限脉冲响应数字滤波器的设计.ppt

1、2 3IIR数字滤波器：数字滤波器：可以利用模拟滤波器设计可以利用模拟滤波器设计但相位非线性，需要加相位校正网络但相位非线性，需要加相位校正网络FIR数字滤波器：数字滤波器： 可以严格线性相位，又可任意幅度特性可以严格线性相位，又可任意幅度特性因果稳定系统因果稳定系统但阶次比但阶次比IIR滤波器要高得多滤波器要高得多4 窗函数法窗函数法频率采样法频率采样法等纹波逼近法等纹波逼近法计算机优化设计方法计算机优化设计方法5 FIR滤波器的单位冲激响应：滤波器的单位冲激响应：( )01h nnN10( )( )NnnH zh n z系统函数：系统函数：在在 z 平面有平面有N 1 个零点个零点在在 z

2、 = 0 处是处是N 1 阶极点阶极点 6h(n)为实序列时，其频率响应：为实序列时，其频率响应：即群延时即群延时 是常数是常数( )dd 0( ) 第二类线性相位：第二类线性相位：( ) 第一类线性相位：第一类线性相位：()()jjH ee 10()( )Njj nnH eh n e()( )jHe 线性相位是指线性相位是指 是是 的线性函数。的线性函数。 H(wH(w):):幅度特性，幅度特性，(w):(w):相位特性。相位特性。7( ) ( )(1)01h nh NnnN 即：nh(n)为实序列；为实序列； nh(n) 关于关于 (N 1) /2 偶对称偶对称N为偶数时为偶数时N为奇数时

3、为奇数时80( ) ( )(1)01h nh NnnN 即：nh(n)为实序列；为实序列； nh(n) 关于关于 (N 1) /2 奇对称奇对称N为偶数时为偶数时N为奇数时为奇数时91100( )( )(1)NNnnnnH zh n zh Nn z 1(1)0( )NNmmh m z (1)1()NzH z 1mNn 令系统函数：系统函数：( )(1)01h nh NnnN 由1(1)0( )NNmmzh m z 10(1)11 ( )( )()2NH zH zzH z得11(1)001( )( )2NNnNnnnh n zzh n z1(1)01( )2NnNnnh nzzz11122120

4、( )2NNnnNNnzzzh n (1)1 ()NH zzH z 由1111221cos 221sin 2jNNnnz eNnzzNjn 11122120( )2NNnnNNnzzH zzh n12112011201( )cos2()( )1( )sin2jNNjnjz eNNjnNeh nnH eH zNjeh nn 11122120( )2NNnnNNnzzH zzh n13n频率响应：频率响应：( )(1)h nh Nn 11201()( )( )cos2jNNjjz enNH eH zeh nn12N1( )2N 101( )( )cos2NnNHh nn为第一类线为第一类线性相位性

5、相位14n频率响应：频率响应：( )(1)h nh Nn 11201()( )( )sin2jNNjjz enNH eH zjeh nn12N112201( )sin2NNjjnNeh nn0/21( )22N 101( )( )sin2NnNHh nn为第二类线为第二类线性相位性相位15n 当当h(n)为实序列为实序列n 且当且当h(n)奇或偶对称奇或偶对称且当且当h(nh(n) )偶对称时为第一类线性相位，偶对称时为第一类线性相位，当当h(nh(n) )奇对称时为第二类线性相位奇对称时为第二类线性相位时，其相位特性为线性相位时，其相位特性为线性相位16A、h(n)偶对称，偶对称，N为奇数为

6、奇数11cos(1)cos22NNNnn 11cos22NNn对呈偶对称101( )( )cos2NnNHh nn幅度函数：幅度函数：1cos2Nn17-32011( )2 ( )cos22NnNNHhh nn121112cos()22NmNNhhmm12Nnm令120( )( )cos()NnHa nn1(0)2Nah其中：其中：11,.,2Nn1( )22Na nhn18120( )( )cos()NnHa nn1(0)2Nah11,.,2Nn其中：其中：1( )22Na nhn19( )0, , 2 H对呈偶对称cos()0, 2 n对，呈偶对称120( )( )cos()NnHa nn

7、h(n)偶对称，偶对称，N为奇数为奇数可以实现各可以实现各种类型的滤种类型的滤波器波器20B、h(n)偶对称，偶对称，N为偶数为偶数12012 ( )cos2NnNh nn101( )( )cos2NnNHh nn幅度函数：幅度函数：21/211( )( )cos2NnHb nn( )22Nb nhn1,.,2Nn 其中：其中：22( )H对呈奇对称n不能设计成高通、带阻滤波器不能设计成高通、带阻滤波器 /211( )( )cos2NnHb nnh(n)偶对称，偶对称，N为偶数为偶数( )0, 2 H对呈偶对称23C、h(n)奇对称，奇对称，N为奇数为奇数11sin(1)sin22NNNnn

8、11sin22NNn对呈奇对称101( )( )sin2NnNHh nn幅度函数：幅度函数：1sin2Nn 24-3201( )2 ( )sin2NnNHh nn12112sin()2NmNhmm12Nnm令121( )( )sin()NnHc nn11,.,2Nn1( )22Nc nhn其中：1( )02Nh nNh奇对称且 为奇数25121( )( )sin()NnHc nn1( )22Nc nhn11,.,2Nn26( )0, 2H故对， 呈奇对称121( )( )sin()NnHc nnn sin()0, 2 n因对，呈奇对称h(n)奇对称，奇对称，N为奇数为奇数n只能实现带通滤波器只

9、能实现带通滤波器27/211( )( )sin2NnHd nn( )22Nd nhn1,.,2Nn D、h(n)奇对称，奇对称，N为偶数为偶数28( )H对呈偶对称/211( )( )sin2NnHd nnh(n)奇对称，奇对称，N为偶数为偶数( )0, 2 H对呈奇对称n不能实现低通和带阻滤波器不能实现低通和带阻滤波器293031(1)(1)要设计一个线性相位的低通要设计一个线性相位的低通DFDF 从幅度特性考虑，只能用第一种或第二种，从幅度特性考虑，只能用第一种或第二种， 即即h(nh(n) )为偶对称为偶对称(2)(2)要设计一个线性相位的高通要设计一个线性相位的高通DFDF 从幅度特性

10、考虑，可用第一种（从幅度特性考虑，可用第一种（h(nh(n) )偶，偶， N N为奇）或第四种（为奇）或第四种（h(nh(n) )为奇，为奇，N N为偶）为偶）(3)(3)要设计一个线性相位的带阻要设计一个线性相位的带阻DFDF 从幅度特性考虑只能用第一种（从幅度特性考虑只能用第一种（h(nh(n) ) 为偶，为偶，N N为奇）为奇）(4)(4)要设计一个线性相位的带通要设计一个线性相位的带通DFDF 从幅度特性考虑，可以四种中的任一种设计从幅度特性考虑，可以四种中的任一种设计32nh(n)为奇对称时，有为奇对称时，有900相移，适用于相移，适用于900移相器（也称正交变换网络），而移相器（也

11、称正交变换网络），而选频滤波器采用选频滤波器采用h(n)为偶对称时。为偶对称时。33( )0iH z*, 1/iizz即 也是零点(1)1( )()NH zzH z 得：得：由由1）若）若 z = zi 是是H(z)的零点，则的零点，则 z = zi-1 也是零点也是零点2）h(n)为实数，则零点共轭成对为实数，则零点共轭成对线性相位滤波器的零点是互为倒数的共轭对线性相位滤波器的零点是互为倒数的共轭对 即共轭成对且镜像成对即共轭成对且镜像成对(1)1()( )0NiiiH zzH z 3435 设设N为偶数，有为偶数，有1112002( )( )( )( )NNNnnnNnnnH zh n z

12、h n zh n z1122(1)0012(1)0( )( )(1)( )NNnN mnmNnN nnH zh n zh Nmzh nzz 令m=N-n-1,则有( )(1)h nh Nn 361112(1)201( )( )()2NNnN nnNH zh n zzhz 12(1)0( )( )NnN nnH zh n zz N为偶数：为偶数：N为奇数：为奇数：中间项h(N-1)/2单独列出仅需要仅需要N/2个乘法器个乘法器仅需要（仅需要（N N1 1）/2/2个乘法器个乘法器37第一类线性相位网络结第一类线性相位网络结构构x(n)y(n)z1z1z1z1z1z1z1h(0)h(1)h(2)h

13、(N/21)x(n)y(n)z1z1z1z1z1z1h(0)h(1)h(2)h(N1)/2)N 偶数N 奇数12(1)0( )( )NnN nnH zh n zz 38第二类线性相位网络结构第二类线性相位网络结构x(n)y(n)z1z1z1z1z1z1z1h(0)h(1)h(2)h(N/21)x(n)y(n)z1z1z1z1z1z1h(0)h(1)h(2)h(N1)/2)N 偶数N 奇数11111111112(1)0( )( )NnN nnH zh n zz 394010()( )Njj nddnHeh n e1( )2jj nddh nHeed( )( )dH zZT h n41线性相位理想

14、低通滤波器的频率响应：线性相位理想低通滤波器的频率响应：()0,jjccdcceHe 1sin()( )2()ccjj nccdcnh needn其理想单位抽样响应：其理想单位抽样响应：中心点为中心点为 的偶对称无限长非因果序列的偶对称无限长非因果序列42 h hd d(n(n) )为偶对称，无限长非因果序列为偶对称，无限长非因果序列43n将将n可以将可以将看作是看作是与矩形窗与矩形窗RN(n)相乘相乘( )01( )0dh nnNh n其它12N( )( )( )dNh nh n Rn444510()( )Njj nddnHeh n e1( )2jj nddh nHeed( ) ( )H z

15、ZT h n( )( )( )dNh nh n Rn窗函数法窗函数法4647n时域加窗处理时域加窗处理n频域卷积频域卷积( )( )( )dNh nh n Rn11()()()()22jjjjjdNdNH eHeReHeRed 48n时域表示时域表示n频域表示频域表示n幅度幅度n相位相位101( )0NnNRn其它1120sin(2)()sin(2)NNjjj nNnNReeesin(2)( )sin(2)NNR12N49nFIR滤波器的幅度特性滤波器的幅度特性11()()()( )()22jjjjdNdNH eHeReeHRd1( )( )()2dNHHRd5051(0)( )NHR近似于的

16、全部积分面积()0.5(0)cHH2cHN为最大值，正肩峰2cHN为最小值，负肩峰( )H随，绕零值波动( )(0)HH随，绕波动0n cn 2cNn 2cNn 2cNn 2cNn 1( )( )()2dNHHRd幅度函数：幅度函数：52n在理想特性不连续点在理想特性不连续点c c附近形成过渡带。过滤带的附近形成过渡带。过滤带的宽度近似等于宽度近似等于R RN N() )主瓣宽度，主瓣宽度，=4/N=4/N ，过渡，过渡带越小越好。带越小越好。 n通带内增加了波动，最大的峰值在通带内增加了波动，最大的峰值在c c- - 2/N2/N 处。处。阻带内产生了余振，最大的负峰在阻带内产生了余振，最大

17、的负峰在c c+ +2/N2/N处。处。n通带与阻带中波动的情况与窗函数的幅度谱有关。通带与阻带中波动的情况与窗函数的幅度谱有关。 R RN N() )波动愈快（加大时），通带与阻带内波动愈波动愈快（加大时），通带与阻带内波动愈快，快， R RN N() )旁瓣的大小直接影响波动的大小。旁瓣的大小直接影响波动的大小。n这些影响是对这些影响是对h hd d(n(n) )加矩形窗引起的，称之为加矩形窗引起的，称之为吉布斯吉布斯效应效应。 53时域截断，频域产生吉布斯现象时域截断，频域产生吉布斯现象频域截断，时域产生吉布斯现象频域截断，时域产生吉布斯现象信号与系统信号与系统541 1、在、在c c附

18、近形成过渡带，过渡带越小越好。附近形成过渡带，过渡带越小越好。 =4/N=4/N2 2、在通带和阻带内出现波动，使阻带衰减减、在通带和阻带内出现波动，使阻带衰减减 弱，甚至不满足阻带衰减要求。弱，甚至不满足阻带衰减要求。因此：减小过渡带，加大阻带衰减。因此：减小过渡带，加大阻带衰减。55n增加矩形窗口的宽度增加矩形窗口的宽度N，可以减少过渡带的宽度，可以减少过渡带的宽度， =4/N 。n但加大但加大N不能减少吉布斯效应的影响。不能减少吉布斯效应的影响。 N加大时，主、旁瓣幅度同时加高，相对值不变，总加大时，主、旁瓣幅度同时加高，相对值不变，总是最大肩峰比是最大肩峰比H(0)高高8.95%,最大

19、负峰比最大负峰比0值超过值超过8.95%。 不能改变主瓣和旁瓣幅度相对值，即不能不能改变主瓣和旁瓣幅度相对值，即不能加大阻带衰减。加大阻带衰减。 n选择合适的窗函数形状，使其谱函数的主瓣包含更选择合适的窗函数形状，使其谱函数的主瓣包含更多的能量，相应旁瓣幅度就变小了；旁瓣的减少可多的能量，相应旁瓣幅度就变小了；旁瓣的减少可使通带与阻带波动减少，从而加大阻带的衰减。但使通带与阻带波动减少，从而加大阻带的衰减。但这样总是以加宽过渡带为代价的。这样总是以加宽过渡带为代价的。sin/2sin/2( )NxNxRN（）56n矩形矩形 (Rectangle)窗窗n三角（三角（Bartlett）窗）窗 n汉

20、宁（汉宁（Hanning）窗）窗n汉明（汉明（Hamming）窗）窗 n布莱克曼（布莱克曼（Blackman）窗）窗n凯泽（凯泽（Kaiser）窗）窗574N主瓣宽度最窄：主瓣宽度最窄：旁瓣峰值幅度：旁瓣峰值幅度：13dB阻带最小衰减：阻带最小衰减：21dB 1120()( )NNjjj nRRnWew n eWe窗谱：窗谱：sin2( )sin2RNW幅度函数：幅度函数： 101( )0nNw n其它58旁瓣峰值幅度：旁瓣峰值幅度：20log第一旁瓣峰值主瓣峰值阻带最小衰减：阻带最小衰减：20log第一旁瓣面积主瓣面积5921012( )212112nNnNw nnNnNN8N主瓣宽度为：主

21、瓣宽度为：旁瓣峰值幅度旁瓣峰值幅度25dB阻带最小衰减阻带最小衰减25dB12()NjjW eWe窗谱：窗谱：2sin24( )sin2NWN幅度函数：幅度函数： 1N 6012( )1cos( )21Nnw nRnN8N主瓣宽度宽：主瓣宽度宽：旁瓣峰值幅度旁瓣峰值幅度31dB阻带最小衰减阻带最小衰减44dB22( )0.5( )0.25RRRWWWWNN幅度函数：幅度函数： 1N 612( )0.540.46cos( )1Nnw nRnN8N主瓣宽度宽：主瓣宽度宽：旁瓣峰值幅度：旁瓣峰值幅度：41dB阻带最小衰减：阻带最小衰减：53dB22( )0.54( )0.23RRRWWWWNN幅度函

22、数：幅度函数： 1N 6224( )0.420.5cos0.08cos( )11Nnnw nRnNN12N主瓣宽度最宽：主瓣宽度最宽：旁瓣峰值幅度：旁瓣峰值幅度：57dB阻带最小衰减：阻带最小衰减：74dB22( )0.42( )0.25RRRWWWWNN440.04RRWWNN幅度函数：幅度函数： 1N 6364窗型变化，旁瓣减小，阻带窗型变化，旁瓣减小，阻带衰减增大，主瓣宽度增加衰减增大，主瓣宽度增加6566 ：第一类变形零阶：第一类变形零阶 贝塞尔函数贝塞尔函数0( )I 2002111( )( )nINw nI01nN改变 可同时调整主瓣 宽度和旁瓣幅度 旁瓣幅度 但主瓣宽度670.4

23、0.112(8.7)500.5842(21)0.07886(21)2150021ssssssdBdBdBdB82.285stsptMBB和和M M的计算公式：的计算公式：6869阻带最小衰减只由窗形状决定阻带最小衰减只由窗形状决定过渡带宽则与窗形状和窗宽过渡带宽则与窗形状和窗宽N都有关都有关70对于窗函数法，加大对于窗函数法，加大N N只能减少过渡只能减少过渡带的宽度，不能加大阻带的衰减，要带的宽度，不能加大阻带的衰减，要加大阻带的衰减只能通过选择合适的加大阻带的衰减只能通过选择合适的窗函数来实现。窗函数来实现。基本原则：在保证阻带衰减的情况下，基本原则：在保证阻带衰减的情况下， 尽量选择主瓣

24、窄的窗函数。尽量选择主瓣窄的窗函数。 711.给出希望逼近的理想滤波器的频率响应函数给出希望逼近的理想滤波器的频率响应函数2.计算计算 如果如果不能用简单函数表示，可以用求和代替积不能用简单函数表示，可以用求和代替积分。分。3. 根据允许的过渡带宽及阻带衰减，选定窗函数和根据允许的过渡带宽及阻带衰减，选定窗函数和N值，值，见书表见书表7.2.1, 7.2.21( )()2jj nddh nHeed22101( )()MjkjknMMddkh nHeeM724. 将将与窗函数相乘得与窗函数相乘得FIR数字滤波器的数字滤波器的冲激响应冲激响应h(n)。5. 将将h(n)作作z变换，得到变换，得到H

25、(z)。( ) ( )H zZT h n( )( )( )dNh nh n Wn73/2/0.2ppspsf /2/0.4ststsstsf 50sdB解解：1）求数字频率）求数字频率设计一个线性相位设计一个线性相位FIR低通滤波器，低通滤波器，给定抽样频率为给定抽样频率为421.5 10 (/sec)srad 321.5 10 (/sec)prad通带截止频率为通带截止频率为3t23 10 (/sec)srad 阻带起始频率为阻带起始频率为阻带衰减不小于阻带衰减不小于-50dB，幅度特性如图所示，幅度特性如图所示742）求）求hd(n)()0,jjccdcceHe ccsf()11( )22

26、ccjj njndh needed1sin()()ccnnnn12N1/220.3psts 752( )0.540.46cos( )1Nnw nRnN20.2stptsB 6.6330.2tANB1162N4）确定）确定N 值值50sdB3）选择窗函数：由）选择窗函数：由 确定哈明窗（确定哈明窗（-53dB）6.6tBN哈明窗带宽：765）确定）确定FIR滤波器的滤波器的h(n)( )( ) ( )dh nh n w n33sin 0.3160.540.46cos( )1616nnRnn6）将将h(n)作作z变换，得变换，得到到H(z)。77对模拟信号进行低通滤波处理，要求通带对模拟信号进行低

27、通滤波处理，要求通带0 0f1.5kHzf1.5kHz内衰减小内衰减小于于1dB1dB，阻带，阻带2.5kHzf2.5kHzf上衰减大于上衰减大于40dB40dB。希望对模拟信号采。希望对模拟信号采样后用线性相位样后用线性相位FIRFIR数字滤波器实现上述滤波，采样频率数字滤波器实现上述滤波，采样频率Fs=10kHzFs=10kHz。用窗函数法设计满足要求的。用窗函数法设计满足要求的FIRFIR数字低通滤波器，求数字低通滤波器，求出出h(nh(n),),并画出损耗函数曲线。为了降低运算量，希望滤波器的并画出损耗函数曲线。为了降低运算量，希望滤波器的阶数尽量低。阶数尽量低。解解： 1）求数字频率

28、）求数字频率2150020.310000ppsfF2250020.510000sssfF0.42psc0.2tspB40sdB782）计算阶数）计算阶数M，窗长。，窗长。为了降低阶数选择凯塞窗为了降低阶数选择凯塞窗0.40.5842(21)0.07886(21)3.3953ss840822.28872.2852.285 0.2stMB取取M=23，N=M+1=243）得到）得到h(n)sin 0.4 ()( )( ) ( )( )()dnh nh n w nw nn79n函数函数fir1的调用格式为的调用格式为n b= fir1(n, Wc, ftype, Window)n其中，其中，n为滤波

29、器阶数，为滤波器阶数，Wc为截止频率为截止频率nftype决定滤波器类型，决定滤波器类型，ftype= high，设，设计高通计高通FIR滤波器，滤波器，ftype= stop，设计带，设计带通、带阻通、带阻FIR滤波器。滤波器。nWindow指定窗函数类型，默认为指定窗函数类型，默认为Hamming窗；可选窗；可选Hanning、Hamming、Blackman、bartlett和和boxcar窗，每种窗窗，每种窗都可以由都可以由Matlab的相应函数生成。的相应函数生成。80nwn=boxcar(N)nwn=bartlett(N)nwn=hanning(N)nwn=hamming(N)nw

30、n=blackman(N)nwn=kaiser(N，alph）81fp=1500;fs=2500;rs=40;Fs=10000；wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Bt=ws-wp; alph=0.5842*(rs-21)0.4+0.07886*(rs-21); M=ceil(rs-8)/2.285/Bt); wc=(wp+ws)/2/pi; hn=fir1(M,wc,kaiser(M+1,alph); 820102030-0.6nh(n)00.51-80-60-40-200/20lg|Hg()|830.350.20.15tlplsB用窗函数法设计一个线

31、性相位用窗函数法设计一个线性相位FIRFIR带阻滤波器。要求通带下截带阻滤波器。要求通带下截止频率止频率w wlplp=0.2=0.2，阻带下截止频率，阻带下截止频率w wlsls=0.35 =0.35 ，阻带上截止，阻带上截止频率频率w wusus=0.65 =0.65 ，通带上截止频率，通带上截止频率w wupup=0.8=0.8，通带最大衰减，通带最大衰减a ap p=1dB=1dB，阻带最小衰减，阻带最小衰减a as s=60dB=60dB。解解： 1）选窗型）选窗型2）确定窗长）确定窗长 Bt=12/N/N因为因为as=60dB，所以选择，所以选择blackman窗窗12800.15

32、tANB3）wc=(wlp+wls)/2/,(wup+wus)/2/84wlp=0.2*pi;wls=0.35*pi;wus=0.65*pi;wup=0.8*pi; B=wls-wlp; N=ceil(12*pi/B); wc=(wls+wlp)/2/pi,(wus+wup)/2/pi; hn=fir1(N,wc,stop,blackman(N+1); 85020406080-0.6nh(n)00.51-80-60-40-200/20lg|Hg()|86nFIRFIR线性相位滤波器的基本设计思路；线性相位滤波器的基本设计思路；n窗函数法设计线性相位数字滤波器的改进思路；窗函

33、数法设计线性相位数字滤波器的改进思路；n时域加窗，频域产生吉布斯效应；时域加窗，频域产生吉布斯效应；n减少吉布斯效应的措施；减少吉布斯效应的措施；窗型的选择、窗长的计算窗型的选择、窗长的计算n窗函数法设计窗函数法设计FIRFIR数字滤波器的步骤。数字滤波器的步骤。nMATLABMATLAB实现实现87n对理想滤波器的频率特性对理想滤波器的频率特性在在0,20,2 范围内等范围内等间隔地取样间隔地取样N N个点，得：个点，得：n 对对H Hd d(k(k) )作作N N点逆点逆DFTDFT，即，即IDFT,IDFT,得：得：n DFDF的系统函数为：的系统函数为：2101( )( )0,1,1N

34、jknNdkh nHk enNN2d( )()0,1,1jkNH kHekN( ) ( )H zZT h n适合直接型结构适合直接型结构88n对理想滤波器的频率特性对理想滤波器的频率特性在在0,20,2 范围范围内等间隔地取样内等间隔地取样N N个点，得：个点，得：n根据插值公式实现根据插值公式实现2d( )()0,1,1jkNH kHekN12011( )( )1NNjkkNzH kH zNez适合频率采样型结适合频率采样型结构构891）h(n)偶对称，偶对称，N为奇数为奇数 12NjjH eHe( )(2)HH 2kjkjNkH kH eH ekN kHH幅度偶对称：幅度偶对称：1 212

35、kNNkkNN 相位函数：相位函数：2kN90 设用理想低通作为希望设计的滤波器，截止频率为设用理想低通作为希望设计的滤波器，截止频率为c，采样点数采样点数N，N=奇数时：奇数时：( )()1,0,1,2,( )0,1,2,11( ),0,1,2,1ccccH kH NkkkH kkkkNkNkkkNN 912）h(n)偶对称，偶对称，N为偶数为偶数12()( )NjjH eHe( )(2)HH 2kjkjNkH kH eH ekN kHH 幅度奇对称：幅度奇对称：11kkN 相位函数：相位函数：2kN92 设用理想低通作为希望设计的滤波器，截止频率设用理想低通作为希望设计的滤波器，截止频率为

36、为c，采样点数，采样点数N，N=偶数时：偶数时：( )1,0,1,2,( )0,1,2,1()1,0,1,2,1( ),0,1,2,1gcgcccgcHkkkHkkkkNkHNkkkNkkkNN 93n从时域分析误差从时域分析误差 1( )2jj nddh nHeed( )()( )dNrh nh nrN Rn由频率域采样定理知道，在频域由频率域采样定理知道，在频域0202之间等间之间等间隔采样隔采样N N点，利用点，利用IDFTIDFT得到的得到的h(nh(n) )应是应是h hd d(n(n) )以以N N为周期，周期性延拓乘以为周期，周期性延拓乘以R RN N(n(n) )，即，即由：由

37、：若若H Hd d(e(ejwjw) )中有间断点，则中有间断点，则h hd d(n(n) )无限长，无限长，h(nh(n) )有时有时域混叠失真，当域混叠失真，当N N增大时，误差减小。增大时，误差减小。2101( )( )0,1,1NjknNdkh nHk enNN94n从频域分析误差从频域分析误差 12011( )( )1NNjkkNzH kH zNez102()( )NjkH eH kkN n 当当w=2k/N时，时，(w)=1,H(e(w)=1,H(ejwkjwk)=H(k)=H(k),),在采样点在采样点 处，完全恢复，误差为处，完全恢复，误差为0 0。n 当当w w 2k/N时，

38、时，H(ejw)插值实现，有误差。插值实现，有误差。 Hd(ejw)特性愈平滑的区域，误差愈小；特性愈平滑的区域，误差愈小； Hd(ejw)特性曲线间断点处，误差最大。特性曲线间断点处，误差最大。在间断点处形成过在间断点处形成过渡带，通带和阻带渡带，通带和阻带内产生震荡波纹，内产生震荡波纹，且间断点附近震荡且间断点附近震荡幅度最大，使阻带幅度最大，使阻带衰减减小甚至不满衰减减小甚至不满足指标要求。足指标要求。见书图见书图P219P21995n增加增加N N，使抽样点变密，可减小过渡带宽，使抽样点变密，可减小过渡带宽度，但间断点附近误差依然最大，不能加度，但间断点附近误差依然最大，不能加大阻带衰

39、减。大阻带衰减。n在频响间断点处区间内内插一个或几个过在频响间断点处区间内内插一个或几个过渡采样点渡采样点, ,可加大阻带衰减可加大阻带衰减, ,但导致过渡带但导致过渡带变宽。变宽。96增加过渡带抽样点，可加大阻带衰减增加过渡带抽样点，可加大阻带衰减不加过渡抽样点：不加过渡抽样点：s20dB 加一点：加一点：s44 54dB 加两点：加两点：s65 75dB 加三点：加三点：s85 95dB 102NjkH eH kkNn过渡带增加采样点阻带衰减加大过渡带增加采样点阻带衰减加大 97110Int824( )10Int191622ccNNkH kNNk 利用频率抽样法设计一个频率特性为矩形的理想

40、低通滤利用频率抽样法设计一个频率特性为矩形的理想低通滤波器，截止频率为波器，截止频率为0.50.5，抽样点数为，抽样点数为N N=33=33，要求滤波器，要求滤波器具有第一类线性相位。具有第一类线性相位。解：解：100cjdHe理想低通频率特性其它98得线性相位得线性相位FIR滤波器的频率响应：滤波器的频率响应：816133sin 33sin 33sin233233233sin33sin33sin2233233jjkkkH eekk2 /33过渡带宽：过渡带宽：阻带衰减：阻带衰减：-20dB99n增加一点过渡带抽样点增加一点过渡带抽样点令令H(9)=0.54 /33过渡带宽：过渡带宽：阻带衰减

41、：阻带衰减：-40dB100n增加两点过渡带抽样点增加两点过渡带抽样点令令H(17)=0.5886 H(18)=0.10656 /65过渡带宽：过渡带宽：阻带衰减：阻带衰减：-60dB增加抽样点数为增加抽样点数为N N=65=65101clear all;N=33;wc=0.5*pi;k=0:N-1;phase=(-pi*k*(N-1)/N; M= floor(wc/(2*pi/N);HK=ones(1,M+1),zeros(1,N-2*M-1),ones(1,M);HK1=HK.*exp(j*phase);hn=ifft(HK1,N);freqz(hn,1,512); 102103HK=on

42、es(1,M+1),zeros(1,N-2HK=ones(1,M+1),zeros(1,N-2* *M-1),ones(1,M);M-1),ones(1,M);改为：改为：HK=ones(1,M+1),0.5,zeros(1,N-2HK=ones(1,M+1),0.5,zeros(1,N-2* *M- M- 3,0.5,ones(1,M); 3,0.5,ones(1,M);104（1）根据阻带的最小衰减）根据阻带的最小衰减as选择过渡带采样点的个数选择过渡带采样点的个数m。（2）确定过渡带宽度）确定过渡带宽度Bt，估算频域采样点数，估算频域采样点数N(即窗长）。即窗长）。2(1)tNmB（3）

43、构造一个希望逼近的频率响应函数）构造一个希望逼近的频率响应函数Hd(ejw)。（4）进行频率采样得到）进行频率采样得到H(k)。 (5) 对对H(k)进行进行N点的点的IDFT，得到，得到h(n)。（6）检验设计结果，不过阻带衰减未达到指标要求，则要）检验设计结果，不过阻带衰减未达到指标要求，则要改改 变过渡带采样值，重复做，直到满足要求。变过渡带采样值，重复做，直到满足要求。105用频率采样法设计第一类线性相位低通用频率采样法设计第一类线性相位低通FIR数字滤波器，要数字滤波器，要求通带截止频率求通带截止频率wp=/3/3，阻带最小衰减大于，阻带最小衰减大于40dB40dB，过渡带，过渡带宽

44、度宽度B Bt t/16/16。解：解： s s=40dB=40dB时，查表时，查表.1，m=1m=1， Bt/16由由2(1)64tNmB取取N=65106T=input(T= )； %过渡带的采样值过渡带的采样值datB=pi/16;wc=pi/3; m=1;N=(m+1)*2*pi/datB+1; N=N+mod(N+1,2); Np=fix(wc/(2*pi/N);Ns=N-2*Np-1; Ak=ones(1,Np+1),zeros(1,Ns),ones(1,Np); Ak(Np+2)=T;Ak(N-Np)=T; thetak=-pi*(N-1)*(0:N-1)/N;

45、Hk=Ak.*exp(j*thetak);hn=real(ifft(Hk); Hw=fft(hn,1024); wk=2*pi*0:1023/1024;Hgw=Hw.*exp(j*wk*(N-1)/2); 10700.511.5200.51/Hg(k)0204060-nh(n)00.51-60-40-200/20lg|Hg()|108b=fir2(n,f,m,window)例：例： f=0 0.6 0.6f=0 0.6 0.6 1; 1; m=1 1 0 0 ; m=1 1 0 0 ; window=hanning(31); window=hanning(31); b=fi

46、r2(30,f,m,window); b=fir2(30,f,m,window); h,w h,w=freqz(b,1);figure(1);=freqz(b,1);figure(1); plot(f,m,w/pi,abs(h plot(f,m,w/pi,abs(h);); figure(2); figure(2); plot(w/pi,20 plot(w/pi,20* *log(abs(hlog(abs(h); ); legend(Ideal,fir2 Designed); legend(Ideal,fir2 Designed);1090100.20.4