基于Matlab的FIR滤波器设计毕业设计

1、基于 matlab 的 fir 滤波器设计 摘 要 介绍了利用 matlab 信号处理工具箱进 fir 滤波器设计的三种方法：窗函数法、频 率采样法和最优化设计，给出了详细的设计步骤，并将设计的滤波器应用到一个混和 正弦波信号，以验证滤波器的性能。 关键词：matlab，数字滤波器，有限冲激响应，窗函数，仿真 the fir filter design based on matlab author：liu junyue tutor：huang liqun abstract introduced into the signal processing toolbox using matlab fi

2、r filter design in three ways: window function method, frequency sampling method and the most optimal design, detailed design steps are given, and the design of the filter is applied to a sine wave mixing signal, to verify the filter performance. key words: matlab, digital filter, finite impulse res

3、ponse, the window function, simulation 目目 录录 1 绪论.1 1.1 数字滤波器的研究背景和意义.1 1.2 数字滤波器的应用及现状.1 1.3 数字滤波器的设计方法.1 2 数字滤波器的结构.1 2.1 直接型.1 2.2 级联型.1 2.3 频率取样型.1 3 fir 滤波器.1 3.1 fir 滤波器的基本概念.1 3.2 fir 滤波器的特点.1 3.3 fir 滤波器的种类 .1 3.3.1 数字集成电路 fir 滤波器 .1 3.3.2 dsp 芯片 fir 滤波器 .1 3.3.3 可编程 fir 滤波器.1 4 fir 数字滤波器设计.

4、1 4.1 窗函数法.1 4.1.1 矩形窗.1 4.1.2 海宁(hanning)窗与汉明(hamming)窗.1 4.1.3 布来克曼窗(blackman).1 4.1.4 凯泽窗(kaiser).1 4.1.5 窗函数法设计步骤及举例.1 4.2 频率采样法.1 4.2.1 基本思想.1 4.2.2 约束条件.1 4.2.3 设计误差.1 4.3 最优化设计.1 4.3.1 等波纹切比雪夫逼近准则.1 4.3.2 仿真函数.1 5 结果分析.1 5.1 窗函数法仿真结果.1 5.2 频率采样法.1 5.3 最优化设计.1 结 论.1 致 谢.1 参考资料.1 附 录.1 附录 a .1

5、1 绪论 1.1 数字滤波器的研究背景和意义 当今，数字信号处理技术正飞速发展，它不但自成一门学科，更是以不同形式影响 和渗透到其他学科；它与国民经济息息相关，与国防建设紧密相连；它影响或改变着 我们的生产、生活方式，因此受到人们普遍的关注。 数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网 络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、 通信信号、导航信号等等。上述这些信号大部分是模拟信号，也有小部分是数字信号。 模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。大多 数情况下一维模拟信号的自变量是时间，经过时间上的离散

6、化（采样）和幅度上的离 散化（量化） ，这类模拟信号便成为一维数字信号。因此，数字信号实际上是数字序列 表示的信号，语音信号经采样和量化后，得到的数字信号是一个一维离散时间序列； 而图像信号的经采样和量化后，得到的是数字信号是一个二维离散空间序列。 数字滤波技术是数字信号分析、处理技术的重要分支。无论是信号的获取、传输， 还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传输时至 关重要的。在所有的电子系统中，使用最多技术最复杂的要算数字滤波器了。数字滤 波器的优劣直接决定产品的优劣。 1.2 数字滤波器的应用及现状 数字滤波器精确度高、使用灵活、可靠性高，具有模拟设备所没有的

7、许多优点， 已广泛应用于各个学科技术领域，例如数字电视、语音、通信、雷达、声纳、遥感、 图像、生物医学以及许多工程应用领域。随着信息时代数字时代的到来，数字滤波技 术已经成为一门极其重要的学科和技术领域。以往的滤波器大多采用模拟电路技术， 但是，模拟电路技术存在很多难以解决的问题，例如，模拟电路元件对温度的敏感性， 等等。而采用数字技术则避免很多类似的难题，当然数字滤波器在其他方面也有很多 突出的优点，所以采用数字滤波器对信号进行处理是目前的发展方向。 1.3 数字滤波器的设计方法 数字滤波器按照单位取样响应的时域特性可以分为无限脉冲响应（iir）系统( )h n 和有限脉冲响应（fir）系统

8、。fir 数字滤波器的优点在于它可以做成具有严格线性相 位,而同时可以具有任意的幅度特性;它的传递函数没有极点;这保证了设计出的fir 数 字滤波器一定是平稳的。 所谓数字滤波器设计，简单地说，就是要找到一组能满足特定滤波要求的系数向 量a和b。而滤波器设计完成后还需要进一步考虑如何将其实现，即选择什么样的滤波 器结构来完成滤波运算。fir数字滤波器的设计方法很多，其中较为常用的是窗函数设 计法、频率采样设计法和最优化设计法。本文讨论利用窗函数法、频率采样法和等波 纹切比雪夫逼近法（调用remez函数）来分别实现各种fir滤波器的设计。 窗函数法设计的基本思想是把给定的频率响应通过 idtft

9、（inverse discrete time fourier transform），求得脉冲响应，然后利用加窗函数对它进行截断和平滑，以实现 一个物理可实现且具有线性相位的 fir 数字滤波器的设计目的。其核心是从给定的频 率特性，通过加窗确定有限长单位取样响应；频率采样法设计的基本思想是把给( )h n 出的理想频率响应进行取样，通过 idft 从频谱样点直接求得有限脉冲响应；最优化 设计方法是指采用最优化准则来设计的方法，在 fir df 的最优化设计中 ,最优化准则 有均方误差最小化准则和等波纹切比雪夫逼近(也称最大误差最小化)准则两种。 2 数字滤波器的结构 数字滤波器可以用差分方程、

10、单位取样响应、以及系统函数等表示。对于研究系统 的实现方法，即它的运算结构来说，用流图表示最为直接。一个给定的输入输出关系， 可以用多种不同的数字网络来实现。在不考虑量化影响时，这些不同的实现方法是等 效的；但在考虑量化影响时，这些不同的实现方法性能上就有差异。因此，运算结构 是很重要的，同一系统函数 h(z)，运算结构的不同，将会影响系统的精度、误差、稳 定性、经济性以及运算速度等许多重要性能。 无限长单位脉冲响应滤波器与有限长单位脉冲响应滤波器在结构上有各自不同的特 点，必须分别讨论。 有限长单位脉冲响应数字滤波器的结构 有限长单位脉冲响应滤波器的系统函数为 (2.1) 1 0 ( )(

11、) n n n h zh n z 其差分方程为 (2.2) 1 0 ( )( ) () n k y nh k x nk 其基本结构形式有下述几种。 2.1 直接型 由公式（2.1）可得出图 2.1 所示的直接型结构。我们也称这种结构为抽头延迟线 滤波器或横向型结构。由于上述公式就是信号的卷积型式，故还称为卷积型结构。 图 2.1 fir 滤波器直接型结构 2.2 级联型 将(2.1)式的系统函数分解成若干一阶和二阶多项式的乘积( )h z 1 z ( )x n ( )y n (0)h 1 z 1 z (1)h(2)h(2)h n (1)h n (2.3) 21 12 012 01 ( )( )

12、() n n n kkk nk h zh n zzz 其中表示取的整数部分。若 n 为偶数，则 n-1 为奇数，故系数中有一个为零， 2 n 2 n 2k 图 2.2 fir 滤波器的级联型结构（n 为奇数） 这是因为，这时有奇数个根，其中复数根成共轭对，必为偶数，必然有奇数个实根。 图 2.2 示出了 n 为奇数时，fir 滤波器的级联结构，其中每一个二阶因子用直接型结构 实现。 直接型结构和级联型结构在雷达信号处理中作为相关器和对消器等获得了广泛的 应用。 2.3 频率取样型 系统函数在单位圆上作 n 等分取样的取样值就是的离散傅里叶变换值( )h z( )h n h(k)，用表示的内插公

13、式为( )h k( )h z (2.4) 1 1 0 1( ) ( )(1) 1 n n k k n h k h zz nwz 由上式可见，fir 系统可用一个子 fir 系统和一个子 iir 系统1 n z 实现，子 fir 系统是一个由 n 节延迟单元组成的梳状滤波器，如图 1 1 0 ( ) 1 n k k n h k wz 1 n z 2.3 所示。在单位圆上有 n 个等分的零点1 n z (2.5)10 n z (2.6) 2 0,1,1 n jk k knze 梳状滤波器的频率响应 1 z ( )x n( )y n 01 11 21 1 z 1 z 02 12 22 1 z 1 z

14、 2 0 n 2 1 n 2 2 n 1 z (2.7)()11 j jnjn z e h eze 其幅度特性为 (2.8)()2 sin 2 j n h e 其幅频响应如图 2.3 所示，由于其频谱形状像梳子，因此称其为梳状滤波器。 ( )y n n z ( )x n () j c h e 0 2 n 4 n 图 2.3 梳状滤波器结构及频率响应幅度 子 iir 系统是 n 个型的分式和的形式，每一个一阶网络在单 1 ( ) 1 k n h k wz 1 ( ) 1 k n h k wz 位圆上有一极点 (2.9) 2 k nk n jk zwe 因此网络对频率为的响应为，是一个谐振频率为的

15、无耗谐振器。并联 2 k n 2 k n 谐振器的极点正好各自抵消一个梳状滤波器的零点，从而使在频率点处的响 2 k n 应就是。因此控制滤波器的响应很直接，这正是频率取样型结构的特点。( )h k r im jz re z 11 图 2.4 抽样点改在 r passrad=0.2*pi; w1=boxcar(81); w2=triang(81); w3=hanning(81); w4=hamming(81); n=1:1:81; hd=sin(passrad*(n-41)./(pi*(n-41); hd(41)=passrad/pi; h1=hd.*rot90(w1); h2=hd.*rot

16、90(w2); h3=hd.*rot90(w3); h4=hd.*rot90(w4); mag1,rad=freqz(h1); mag2,rad=freqz(h2); mag3,rad=freqz(h3); mag4,rad=freqz(h4); subplot(2,2,1); plot(rad,20*log10(abs(mag1); grid on; subplot(2,2,2); plot(rad,20*log10(abs(mag2); grid on; subplot(2,2,3); plot(rad,20*log10(abs(mag3); grid on; subplot(2,2,4)

17、; plot(rad,20*log10(abs(mag4); grid on; 图 4.3 窗函数设计的低通滤波器的幅度频率特性 窗函数主要用来减少序列因截断而产生的gibbs效应。但当这个窗函数为矩形时， 得到的fir滤波器幅频响应会有明显的gibbs效应，并且任意增加窗函数的长度（即fir 滤波器的抽头数）gibbs效应也不能得到改善。为了克服这种现象,窗函数应该使设计的 滤波器： （1） 频率特性的主瓣宽度应尽量窄,且尽可能将能量集中在主瓣内； （2） 窗函数频率特性的旁瓣趋于 的过程中，其能量迅速减小为零。 4.2 频率采样法 频率采样法是从频域出发,根据频域采样定理，对给定的理想滤波

18、器的频率响应 加以等间隔的抽样 ，得到：() jw h e( ) d h k k=0,1,n-1 (4.15) (2 ) ( )() k n jw ddw hkhe 再利用可求得 fir 滤波器的系统函数及频率响应。( ) d hk( )h z() jw h e 而在各采样点间的频率响应则是其的加权内插函数延伸叠加的结果。但对于一个无 限长的序列，用频率采样法必然有一定的逼近误差,误差的大小取决于理想频响曲线的 形状, 理想频响特性变换越平缓, 则内插函数值越接近理想值,误差越小。为了提高 逼近的质量，可以通过在频率相应的过渡带内插入比较连续的采样点，扩展过渡带使 其比较连续，从而使得通带和阻

19、带之间变换比较缓慢，以达到减少逼近误差的目的。 选取 w0,2内 n 个采样点的约束条件为： (4.16) ( )() ()() h kh n k mn m 01kn (1) 增大阻带衰减三种方法： 1) 加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。 2) 过渡带的优化设计,利用线性最优化的方法确定过渡带采样点的值，得到要求 的滤波器的最佳逼近（而不是盲目地设定一个过渡带值） 。 3)增大 n。如果要进一步增加阻带衰减，但又不增加过渡带宽，可增加采样点数 n。代价是滤波器阶数增加，运算量增加。 直接从频域进行设计，物理概念清楚，直观方便；适合于窄带滤波器设计，这时 频率响应只有少数几个非零值

20、，但是截止频率难以控制。 典型应用：用一串窄带滤波器组成多卜勒雷达接收机，覆盖不同的频段，多卜勒 频偏可反映被测目标的运动速度。 4.2.1 基本思想 使所设计的 fir 数字滤波器的频率特性在某些离散频率点上的值准确地等于所需 滤波器在这些频率点处的值，在其他频率处的特性则有较好的逼近。 4.2.2 约束条件约束条件 为了设计线性相位的 fir 滤波器，采样值 h(k)要满足一定的约束条件。前已指出, 具有线性相位的 fir 滤波器，其单位采样响应 h(n)是实序列，且满足 h(n)=+/-h(n-1-n) (4.17) 由此得到的幅频和相频特性，就是对 h(k)的约束。 4.2.3 设计误

21、差设计误差 fir 设计步骤： 给定指标 k,hk h(k) h(ej)或 h(z)关心的是，由上述设计过程得到的 h(ej) 与 h（k）的逼近程度，以及 h(ej)与 h(k)的关系。 令 ，则 单位圆上的频响为： 这是一个内插公式，式中为内插函 数。 令则 所以，在每个采样点上，频响 h(ej) 严格地与理想特性 h(k)一致，在采样点之 间，频响由各采样点的内插函数延伸迭加而形成，因而有一定的逼近误差，误差大小 与理想频率响应的曲线形状有关，理想特性平滑，则误差小；反之，误差大，在理想 频率响应的不连续点会产生肩峰和波纹。n 增加，则采样点变密，内插误差减小。 图图 4.4 频率采样的

22、响应频率采样的响应 4.3 最优化设计 最优化设计方法是指采用最优化准则来设计的方法。在 fir df 的最优化设计中 , 最优化准则有均方误差最小化准则和等波纹切比雪夫逼近(也称最大误差最小化)准则两 种。实际设计中 ,只有采用窗函数法中的矩形窗 才能满足前一种最优化准则 ,但由于 吉布斯 (gibbs )效应的存在，使其根本不能满足设计的要求。为了满足设计的要求 ,可 以采用其它的窗函数来消除吉布斯效应 ,但此时的设计已经不能满足该最优化准则了。 因此 ,要完成 fir df 的最优化设计 ,只能采用后一种优化准则来实现。 前面介绍了 fir 数字滤波器的两种逼近设计方法，即窗口法（时域逼

23、近法）和频 率采样法（频域逼近法） ，用这两种方法设计出的滤波器的频率特性都是在不同意义上 对给定理想频率特性 hd(ej)的逼近。 说到逼近，就有一个逼近得好坏的问题，对“好”“坏”的恒量标准不同，也会得出不 同的结论，我们前面讲过的窗口法和频率采样法都是先给出逼近方法，所需变量，然 后再讨论其逼近特性，如果反过来要求在某种准则下设计滤波器各参数，以获取最优 的结果，这就引出了最优化设计的概念，最优化设计一般需要大量的计算，所以一般 需要依靠计算机进行辅助设计。 最优化设计的前提是最优准则的确定，在 fir 滤波器最优化设计中，常用的准则 有 最小均方误差准则 最大误差最小化准则。 1) 均

24、方误差最小化准则 若以 e(ej)表示逼近误差，则 (4.18) 那么均方误差为 (4.19) 均方误差最小准则就是选择一组时域采样值，以使均方误差，这一方法注重的是在整 个- 频率区间内总误差的全局最小，但不能保证局部频率点的性能，有些频率点可 能会有较大的误差，对于窗口法 fir 滤波器设计，因采用有限项的 h(n)逼近理想的 hd(n)，所以其逼近误差为： (4.20) )()( jj d j ehehee）（ deedeheh jjj d 2 2 2 2 1 2 1 n d nhnh 2 2 )()( 如果采用矩形窗 (4.21) 则有 (4.22) 2) 最大误差最小化准则（也叫最佳

25、一致逼近准则） (4.23) 其中 f 是根据要求预先给定的一个频率取值范围，可以是通带，也可以是阻带。 最佳一致逼近即选择 n 个频率采样值（或时域 h(n) 值） ，在给定频带范围内使频响的 最大逼近误差达到最小,也叫等波纹逼近。 优点：可保证局部频率点的性能也是最优的，误差分布均匀，相同指标下，可用 最少的阶数达到最佳化。 例如，我们提到的频率采样最优化设计，它是从已知的采样点数 n、预定的一组 频率取样和已知的一组可变的频率取样（即过渡带取样）出发，利用迭代法（或解析 法）得到具有最小的阻带最大逼近误差（即最大的阻带最小衰减）的 fir 滤波器。但 它只是通过改变过渡带的一个或几个采样

26、值来调整滤波器特性。如果所有频率采样值 （或 fir 时域序列 h(m)）都可调整，显然，滤波器的性能可得到进一步提高。可以证 明，这是一个最小均方误差。 所以，矩形窗窗口设计法是一个最小均方误差 fir 设计，根据前面的讨论，我们 知道其优点是过渡带较窄，缺点是局部点误差大,或者说误差分布不均匀。 4.3.1 等波纹切比雪夫逼近准则 在滤波器的设计中 ,通常情况下通带和阻带的误差要求是不一样的。等波纹切比雪 夫逼近准则就是通过对通带和阻带使用不同的加权函数 ,实现在不同频段(通常指的是 通带和阻带)的加权误差最大值相同 ,从而实现其最大误差在满足性能指标的条件下达 到最小值。 尽管窗函数法与

27、频率采样法在 fir 数滤波器的设计中有着广泛的应用, 但两者不 是最优化的设计 。通常线性相位滤波在不同的频带内逼近的最大容许误差要求不同。 等波纹切比雪夫逼近准则就是通过通带和阻带使用不同的加权函数,实现在不同频段(通 常指的是通带和阻带) 的加权误差最大值相同,从而实现其最大误差在满足性能指标的 条件下达到最小值，即使得 和之间的最大绝对误差最小。() jw d he() jw h e 其它0 1)( )( nnonh nh d 1 222 | )()(| )()(| nnn dd nhnhnhnh max|()| min j e ef 等波纹切比雪夫逼近是采用加权逼近误差，它可以表示为

28、： () jw e e (4.24) ()()()() jwjwjwjw d e ew eheh e 其中，为逼近误差加权函数在误差要求高的频段上，可以取较大的加权值，() jw w e 否则，应当取较小的加权值。 尽管按照 fir 数字滤波器单位取样响应 h(n)的对称性和 n 的奇、偶性，fir 数 字滤波器可以分为 4 种类型，但滤波器的频率响应可以写成统一的形式： (4.25) (1) 22 () ()( ) j nw jwjk h eeeh w 其中，k0 ,1 , h ()为幅度函数，且是一个纯实数，表达式也可以写成统一的 形式： (4.26)()( ) ( ) jw d heqp

29、 其中，为 的固定函数，为 m 个余弦函数的线性组合。( )q( )p 4.3.2 仿真函数 利用数字信号处理工具箱中的 remezord 和 remez 函数可以实现 firdf 的最优化 设计。在此先介绍这两个函数: (1) n ,fo ,ao ,weights =remezordf ,a ,dev 功能:利用 remezord 函数可以通过估算得到滤波器的近似阶数 n ,归一化频率带边 界 fo ,频带内幅值 ao 及各个频带内的加权系数 weights。输入参数 f 为频带边缘频率 ,a 为各个频带所期望的幅度值 ,dev 是各个频带允许的最大波动。 (2) h =remez(n ,f

30、o ,ao ,weights, ftype) 功能:利用 remez 函数可以得到最优化设计的 fir df 的系数 ,输入参数 n 是( )h n 滤波器的阶数 ,fo ,ao ,weights 参数含义说明同 (1)。ftype 是所设计的滤波器类型 ,它 除了可以设计普通的滤波器外 ,它还可以设计数字希尔钞特变换器以及数字微分器。实 际设计中 ,由于 remezord 函数可跑高估或低估滤波器的阶数 n ,因此在得到滤波器的 系数后 ,必须检查其阻带最小衰减是否满足设计要求。如果此时的技术指标不能满足设 计要求 ,则必须提高滤波器的阶数到 n +1 ,n +2 等。故等波纹切比雪夫逼近法

31、设计 fir 数字滤波器的步骤是： 给出所需的频率响应，加权函数和滤波器的单位取样响应() jw d he() jw w e 的长度 n。 ( )h n 由中给定的参数来形成所需的、和的表达式。 ( )w( ) d h( )p 根据 remez 算法,求解逼近问题。 利用傅立叶逆变换计算出单位取样响应 。( )h n 5 结果分析 5.1 窗函数法仿真结果 采用特殊的窗函数如 hamming 窗，可以减小 gibbs 效应，但同时也会使滤波器的 过度带变宽。波动幅度取决于窗函数幅度频谱旁瓣的相对幅度, 而波纹的多少取决于 窗函数旁瓣的多少，如图 5.1 所示。以上两点是就是窗函数直接截断hd(

32、n)引起的截断 效应在频域的反映, 截断效应直接影响滤波器的性能, 因为通带内的波动会影响滤波 器痛带中的平稳性, 阻带内的波动则影响阻带最小衰减,因此,减少截断效应也是 fir 数字滤波器设计的关键之一。 图5.1 窗函数设计的fir低通滤波器频率响应 5.2 频率采样法 图5.2为在间断点处增加一个过渡点后的情况。从图5.3中可以看出滤波器的带外衰 减指标有了明显的改善!但这同时增加了滤波器的过渡带宽。所以，在带外衰减和过渡 带宽这两个指标之间需要有一个折衷。因频率取样点都局限在2/n的整数倍点上，所 以在指定通带和阻带截止频率时，这种方法受到限制，比较死板。充分加大n，可以接 近任何给定

33、的频率，但计算量和复杂性增加。频率采样法偏离设计指标明显，阻带衰 减最小，只有适当选取过渡带样点值，才会取得较好的衰耗特性。 图 5.2 fir 的单位取样响应 图 5.3 fir 的低通衰减幅频特性 5.3 最优化设计 在设计中 ,如果该滤波器的特性不满足要求 ,那么 ,原有参数必须作适当调整。这 在程序中很容易实现 ,只需对参数进行重新设定 ,就可以得到新条件下滤波器的特性。 采用最优化设计方法时大大减小了滤波器的阶数,从而减小了滤波器的体积,并最终降 低了滤波器的成本。这样使得设计出来的滤波器更为简单经济。因而在实际的滤波器 设计中,这种最优化方法是完全可行的。在实际应用中 ,如果需要对

34、某一信号源进行特 定的滤波 ,并要检验滤波效果 ,应用传统方法实施起来比较繁琐。在matlab环境下,可先 用软件模拟产生信号源 ,再设计滤波器对其进行滤波 。 图 5.4 滤波器输出的幅频及相频响应特性 同样是设计一个 fir 低通数字滤波器，综合分析可以看出： (1) 窗函数法在阶数较低时，阻带特性不满足设计要求，只有当滤波器阶数较高时， 使用海明窗和凯塞窗基本可以达到阻带衰耗要求； (2) 频率采样法偏离设计指标最明显，阻带衰减最小，而且设计比采用窗函数法复 杂。只有适当选取过渡带样点值，才会取得较好的衰耗特性； (3) 利用等波纹切比雪夫逼近法则的设计可以获得最佳的频率特性和衰耗特性，

35、具 有通带和阻带平坦，过渡带窄等优点。 综上所述，fir 滤波器很容易实现具有严格线性相位的系统, 使信号经过处理后不 产生相位失真,舍入误差小,而且稳定，因此越来越受到广泛的重视。matlab 软件的 诞生, 使数字信号处理系统的分析与设计得简单,它已经成为电子工程师必备的一个工 具软件。 结 论 本文通过一个设计实例，介绍了利用 matlab 实现 fir 滤波器设计与滤波的三 种方法，从仿真结果可以看出它们均可以达到技术指标要求，而且方法简单、快捷， 大大减轻了工作量。滤波器的设计工作完成后，可以借助于 matlab 的 export 操作导 出所设计滤波器的系统函数 h(z)。由于 m

36、atlab 具有强大的接口功能，仿真后的结果可 以很方便的移植到 dsp、cpld 或 fpga 等器件中。在实际应用中,只需按要求修改滤波 器参数,并对程序作较少的改动,即可实现不同截止频率的 fir 滤波器,实用性较强 参考资料 1 董长虹等. matlab信号处理与应用.北京：国防工业出版社，2005. 2 美 m.h.海因斯 著，张建华等 译.数字信号处理.北京：科学出版社，2002. 3 张葛祥，李娜. matlab仿真技术与应用.北京：清华大学出版社，2003. 4 楼顺天，李博菡. 基于matlab的系统分析与设计.西安：西安电子科技大学出版社， 1998. 5 马昌风. 最优化

37、设计法及matlab设计.北京：科学出版社，2010 6 张磊. matlab实用教程.北京：人民邮电出版社，2008 7 陈龙等. 数字信号处理的matlab实现.北京：科学出版社，2010 8 甘本祓,吴万春. 现代微波滤波器的结构与设计.北京 科学出版社，2008 9 美 恒里y-f拉姆 著，模拟和数字滤波器设计与实现.北京：人民邮电出版社， 1985 附 录 附录 a 窗函数实现程序： passrad=0.4*pi; w1=boxcar(61); w2=hamming(61) n=1:1:61; hd=sin(passrad*(n-31)./(pi*(n-31); hd(31)=pas

38、srad/pi; h1=hd.*rot90(w1); h2=hd.*rot90(w2); mag1,rad=freqz(h1); mag2,rad=freqz(h2); subplot(2,2,1); plot(rad,20*log10(abs(mag1); grid on; title(designed by rectangular window); subplot(2,2,2); plot(rad,20*log10(abs(mag2); grid on; title(designed by hamming window); h1,w1=freqz(h1,1,100,2); subplot(

39、2,2,3); plot(w1,unwrap(angle(h1); grid on; h2,w2=freqz(h2,1,100,2); subplot(2,2,4); plot(w2,unwrap(angle(h2); grid on; 频率采样法实现程序：频率采样法实现程序： f=0.2 0.3; a=1 0; dev=0.04 0.02; %给出滤波器的参数 n f0 a0 w=remezord(f,a,dev); n=n; alpha=(n-1)/2; k=0:n-1; wp=0.2*pi; ws=0.3*pi; %计算理想低通滤波器的截止频率 wc=(wp+ws)/2; m=fix(w

40、c*n/(2*pi)+1); %在两边过渡带取值为 0.5 的采样点 t = 0.5; hrs=ones(1,m),t,zeros(1,n-2*m-1),t,ones(1,m-1); k1 = 0:floor(alpha); k2 = floor(alpha+1):n-1; phai=-alpha*(2*pi)/n*k1,alpha*(2*pi)/n*(n -k2); h =hrs.*exp(j*phai); %计算单位冲激响应 h =ifft(h,n); h1,w1 = freqz(h,1,256,1); hr = abs( h1); h1 = 20* log10(hr); %画出 fir

41、df 的单位取样响应 figure(1); k=0:n-1; stem(k,h,k.) axis(0,n-1,1.1*min(real(h),1.1*max(real(h); xlabel(n); ylabel(h(n); grid on; %画出 fir df 的低通衰减幅频特性 figure(2); plot(w1,h1); xlabel(normalized) frequency(rad/sample); ylabel(magnitude(db); grid on; 最优化设计程序实现： fs =2000; %采样频率 rp =3; %通带波纹 rs =40; %阻带波纹 f = 500

42、 600; %截止频率 a = 1 0; %期望的幅度 dev= (10(rp/20) - 1)/ (10(rp/20) + 1) 10( - rs/20) ; n ,fo ,ao ,w = remezord(f ,a ,dev,fs) ; b = remez(n ,fo ,ao ,w) ; %调用最优设计法中 remez 函数 reqz(b ,1 ,1024 ,fs);附录：外文翻译 in wang zuoliangs translation practices, he translated many poems, especially the poems written by robert

43、 burns. his translation of burns “a red, red rose” brought him fame as a verse translator. at the same time, he published about ten papers on the translation of poems. some argue that poems cannot be translated. frost stresses that poetry might get lost in translation. according to wang, verse trans

44、lation is possible and necessary, for “the poet-translator brings over some exciting work from another culture and in doing so is also writing his own best work, thereby adding something to his culture. in this transmission and exchange, a richer, more colorful world emerges. ”(wang, 1991:112). then

45、 how can we translate poems? according to wangs understanding, the translation of poems is related to three aspects: a poems meaning, poetic art and language. （1）a poems meaning “socio-cultural differences are formidable enough, but the matter is made much more complex when one realizes that meaning does not consist in the meaning of words only, but also in syntactical structures, speech rhythms, levels of style.” (wang, 1991:93). （2）poetic art according to wang, “blys point