FIR滤波器窗函数设计
fir滤波器的主要设计方法 -回复
fir滤波器的主要设计方法-回复fir滤波器是一种基本的数字滤波器，主要用于数字信号处理中的滤波操作。
它的设计方法有很多种，包括频率采样法、窗函数法、最优权系数法等。
本文将一步一步回答"[fir滤波器的主要设计方法]"，让我们一起来了解一下吧。
一、频率采样法频率采样法是fir滤波器设计的最基本方法之一。
它的主要思想是在频域中对滤波器的频响特性进行采样，然后通过反变换得到滤波器的冲激响应。
这种方法的优点是设计简单，适用于各种滤波器的设计。
1. 确定滤波器的截止频率和通带、阻带的要求。
根据应用的具体需求，确定滤波器的频率范围和滤波特性。
2. 设计理想的滤波器频率响应。
根据频率范围和滤波特性的要求，设计所需的滤波器频率响应。
常见的有低通、高通、带通、带阻等类型。
3. 进行频率采样。
根据滤波器频率响应的要求，在频域中进行一系列均匀或者非均匀的采样点。
4. 反变换得到滤波器的冲激响应。
对采样得到的频率响应进行反傅里叶变换，得到滤波器的冲激响应。
5. 标准化处理。
对得到的冲激响应进行标准化处理，使得滤波器的增益等于1。
6. 实现滤波器。
根据得到的冲激响应，使用差分方程或者卷积的方法实现fir滤波器。
二、窗函数法窗函数法是一种常用的fir滤波器设计方法，它主要是通过在频域中将理想的滤波器乘以一个窗函数来实现滤波器的设计。
1. 确定滤波器的截止频率和通带、阻带的要求，根据具体应用的需求确定滤波器的频率范围和滤波特性。
2. 设计理想的滤波器频率响应。
根据频率范围和滤波特性要求，设计所需的滤波器频率响应。
3. 选择窗函数。
根据滤波器的频率响应和窗函数的性质，选择合适的窗函数。
4. 计算窗函数的系数。
根据选择的窗函数，计算窗函数的系数。
5. 实现滤波器。
将理想滤波器的频率响应与窗函数相乘，得到实际的滤波器频率响应。
然后使用反变换将频率响应转换为滤波器的冲激响应。
6. 标准化处理。
对得到的冲激响应进行标准化处理，使得滤波器的增益等于1。
窗函数法设计FIR滤波器
FIR 数字滤波器的设计方法IIR 数字滤波器最大缺点：不易做成线性相位，而现代图像、语声、数据通信对线性相位的要求是普遍的。
正是此原因，使得具有线性相位的FIR 数字滤波器得到大力发展和广泛应用。
1． 线性相位FIR 数字滤波器的特点FIR DF 的系统函数无分母，为∑∑-=--=-==11)()(N n n N i ii z n h zb z H ，系统频率响应可写成：∑-=-=10)()(N n jwn jwe n h e H ，令)(jw e H ＝)()(w j e w H Φ，H(w)称为幅度函数，)(w Φ称为相位函数。
这与模和幅角的表示法有所不同，H(w)为可正可负的实数，这是为了表达上的方便。
如某系统频率响应)(jw e H ＝wj we34sin -，如果采用模和幅角的表示法，w 4sin 的变号相当于在相位上加上)1(ππj e =-因，从而造成相位曲线的不连贯和表达不方便，而用)()(w j e w H Φ这种方式则连贯而方便。
线性相位的FIR 滤波器是指其相位函数)(w Φ满足线性方程：)(w Φ＝βα+-w （βα,是常数）根据群时延的定义，式中α表示系统群时延，β表示附加相移。
线性相位的FIR 系统都具有恒群时延特性，因为α为常数，但只有β＝0的FIR 系统采具有恒相时延特性。
问题：并非所有的FIR 系统都是线性相位的，只有当它满足一定条件时才具有线性相位。
那么应满足什么样的条件？从例题入手。
例题：令h(n)为FIR 数字滤波器的单位抽样相应。
N n n ≥<或0时h(n)=0，并假设h(n)为实数。
（a ） 这个滤波器的频率响应可表示为)()()(w j jwew H e H Φ=（这是按幅度函数和相位函数来表示的，不是用模和相角的形式），)(w H 为实数。
（N 要分奇偶来讨论） （1） 当h(n)满足条件)1()(n N h n h --=时，求)(w H 和)(w Φ（π≤≤w 0） （2） 当h(n)满足条件)1()(n N h n h ---=时，求)(w H 和)(w Φ（π≤≤w 0）（b ） 用)(k H 表示h(n)的N 点DFT（1） 若h(n)满足)1()(n N h n h ---=，证明H(0)=0； （2） 若N 为偶数，证明当)1()(n N h n h --=时，H(N/2)=0。
(整理)FIR数字滤波器的(海明)窗函数法设计.
FIR数字滤波器的（海明）窗函数法设计1．课程设计目的（1）熟悉并掌握MATLAB中有关声音（wave）录制、播放、存储和读取的函数。
（2）加深对FIR数字滤波器设计的理解，并用窗函数法进行FIR数字滤波器的设计。
（3）将设计出来的FIR数字滤波器利用MATLAB进行仿真。
（4）对一段音频文件进行加入噪声处理，对带有噪声的文件进行滤波处理。
2.设计方案论证2.1 Matlab语言概述MATLAB是一种以矩阵运算为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求。
随着版本的不断升级，内容不断扩充，功能更加强大，从而被广泛应用于仿真技术、自动控制和数字信号处理领域。
此高级语言可用于技术计算此开发环境可对代码、文件和数据进行管理交互式工具可以按迭代的方式探查、设计及求解问题数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以及数积分等二维和三维图形函数可用于可视化数据各种工具可用于构建自定义的图形用户界面各种函数可将基于MATLAB的算法与外部应用程序和语言（如C、C++、Fortran、Java、COM 以及Microsoft Excel）集成不支持大写输入，内核仅仅支持小写2.2声音处理语音是人类获取信息的重要来源和利用信息的重要手段。
语音信号处理是一门发展十分迅速，应用非常广泛的前沿交叉学科，同时又是一门跨学科的综合性应用研究领域和新兴技术。
声音是一种模拟信号，而计算机只能处理数字信息0和1。
因此，首先要把模拟的声音信号变成计算机能够识别和处理的数字信号，这个过程称为数字化，也叫“模数转换”。
在计算机对数字化后的声音信号处理完后，得到的依然是数字信号。
必须把数字声音信号转变成模拟声音信号，然后再图1 选择windows下的录音机”或是点击快捷按钮图5 加噪后语音信号和频谱图7 滤波器幅频特性与相频特性设计的滤波器是用单位采样响应h(n)表示的，可以利用带噪声语音图8滤波器系统函数。
实验7 窗函数法设计FIR数字滤波器
实验7窗函数法设计FIR数字滤波器一、实验目的掌握窗函数法设计F1R数字滤波器的原理和具体方法二、实验设备与环境计算机、Mat1ab软件环境三、实验基础理论1>基本原理窗函数设计法的基本思想为，首先选择一个适当的理想的滤波器Hd(,3),然后用窗函数截取它的单位脉冲响应％(九)，得到线性相位和因果的FIR滤波器，这种方法的重点是选择一个合适的窗函数和理想滤波器,使设计的滤波器的单位脉冲响应逼近理想滤波器的单位脉冲响应。
2、设计步骤(1)给定理想滤波器的频率响应Hd("3),在通带上具有单位增益和线性相位，在阻带上具有零响应。
一个带宽为g(3c<Tr)的低通滤波器由下式给定h(e j^=(eW∣ω∣≤ωc虱)一1Oωc<∣ω∣<π其中α为采样延迟，其作用是为了得到因果的系统。
(2)确定这个滤波器的单位脉冲响应为了得到一个h(n)长度为N的因果的线性相位FIR滤波器，我们令N-Ia=-2-(3)用窗函数截取hd(τι)得到所设计FIR数字滤波器h(n)h(n)=h d(n)w(n)3、窗函数的选择常用的窗函数有矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗等。
Mat1ab提供了一些函数用于产生窗函数，如下表所示：在设计过程中我们需要根据给定的滤波器技术指标,选择滤波器长度N 和窗函数3(n)°表7.2列出了常用的窗函数的一些特性，可供设计时参考。
其中幻是修正的零阶贝塞尔函数，参数B 控制最小阻带衰减，这种窗函数对于相同的N 可以提供不同的过渡带宽。
由于贝塞尔函数比较更杂，这种窗函数的设计方程很难推导,然而幸运的是，有一些经验设计方程可以直接使用。
已知给定的指标叫Msc,Rp 和4,滤波器长度N 和凯瑟窗参数B 可以按如下凯瑟窗设计方程给出过渡带宽：∆ω=ωst -ωp入一7.95 2.285∆ω_(0.1102(4-8.7) ,P=iθ.5842(4-21)04+0.07886(4-21), 四、实验内容1、设计一个数字低通FIR 滤波器，其技术指标如下ωp =0.2τr,RP=0.25dBωst =0.3τr,A s =50dB分别采用矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗设计该滤波器。
fir滤波器窗函数设计法
fir滤波器窗函数设计法
FIR (Finite Impulse Response)滤波器的窗函数设计法是一种经典的数字滤波器设计方法。
它通过选择一个合适的窗函数来对滤波器的频率响应进行加权，从而实现对信号的滤波。
窗函数设计法的基本步骤如下：
1. 确定滤波器的设计规格：包括截止频率、通带和阻带的幅频响应要求等。
2. 根据设计规格，计算出滤波器的理想频率响应：可以使用理想滤波器的频率响应作为目标。
3. 选择一个合适的窗函数：常用的窗函数有矩形窗、汉宁窗、汉明窗等。
选择窗函数的关键是考虑到主瓣宽度和副瓣衰减的平衡。
4. 将选择的窗函数应用于理想频率响应上，得到加权后的频率响应。
5. 对加权后的频率响应进行反变换，得到滤波器的时域响应。
6. 根据需要，对时域响应进行截断或零增益处理，以满足设计规格。
7. 最后，根据计算得到的滤波器系数，可以通过巴特沃斯频率抽样公式将其转换为巴特沃斯直接型或传输函数型，以便在数字系统中实现滤波。
需要注意的是，FIR滤波器的窗函数设计法是一种近似方法，设计的滤波器无法完全符合理想要求。
设计过程中需要权衡主瓣宽度和副瓣衰减等因素，以及选择合适的截断或零增益处理方式，以获得满
足实际需求的滤波器性能。
FIR滤波器设计分析
FIR滤波器设计分析FIR（Finite Impulse Response）滤波器是一类数字滤波器，其输出只取决于输入信号的有限数量的过去样本。
FIR滤波器的设计分析主要包括滤波器的设计目标、设计方法、设计参数选择、滤波器性能评估等方面。
首先，FIR滤波器的设计目标是根据特定的应用需求，设计一个能够满足给定要求的滤波器。
比如，在音频信号处理中，常见的设计目标包括降低噪声、增强语音清晰度等。
接下来，FIR滤波器的设计方法主要有窗函数法和频率采样法。
窗函数法是通过选择合适的窗函数来设计FIR滤波器，常见的窗函数有矩形窗、汉宁窗、汉明窗等。
频率采样法是通过在频域上选择一组等间隔的频率样点，然后通过频域设计方法将这些样点连接起来，得到FIR滤波器的频响。
设计参数选择是FIR滤波器设计的重要环节。
常见的设计参数包括滤波器阶数、截止频率、过渡带宽等。
滤波器阶数决定了滤波器的复杂度，一般情况下，滤波器阶数越高，滤波器的性能也会越好。
截止频率是指滤波器的频段边界，过渡带宽是指频域中通过频样点与阻带频样点之间的频带范围。
最后，FIR滤波器的性能评估主要包括幅频响应、相频响应、群延迟等指标。
幅频响应可以用来评估滤波器的频率特性，相频响应则描述了信号在滤波过程中的相对延迟。
群延迟是指信号通过滤波器时的延迟时间，对于实时信号处理应用非常重要。
总结起来，FIR滤波器设计分析主要涉及设计目标、设计方法、设计参数选择和滤波器性能评估四个方面。
通过合理选择设计方法和参数，并对滤波器的性能进行评估，可以设计出满足特定要求的FIR滤波器，从而实现信号处理、噪声降低等应用。
基于汉明窗函数的FIR低通滤波器的设计
基于汉明窗函数的FIR低通滤波器的设计1. 滤波器的概念和作用滤波是信号处理中重要的一个部分，可以去除噪声、修复信号、提取特定频率成分等。
滤波器则是一种用来实现滤波功能的电路或软件，其主要作用是选择性地影响输入信号的不同频率成分。
在数字信号处理中，滤波器通常是由一组数字滤波器系数组成的，其中FIR滤波器（Finite Impulse Response）是最简单和最容易实现的一种。
FIR滤波器的基本原理是将输入信号与滤波器系数的乘积相加，从而产生输出信号。
而滤波器系数的选择则决定了滤波器的特性。
在本文中，我们将介绍一种基于汉明窗函数的FIR低通滤波器的设计方法。
2. 汉明窗函数的基本原理汉明窗函数是一种广泛应用于数字信号处理中的窗函数，其形式为：$$w(n) = 0.54 - 0.46\cos\frac{2\pi n}{N-1}, \quad 0 \leq n \leq N-1 $$其中，$N$为窗口长度，决定了窗口的主瓣宽度和副瓣抑制程度。
汉明窗函数的特点是在窗口内具有平坦的频率响应和较高的副瓣抑制能力，而主瓣宽度较大。
3. FIR低通滤波器的设计过程在设计FIR低通滤波器时，我们的目标是将信号中高于一定频率的部分滤除，从而达到降噪的效果。
因此，我们需要选择适当的滤波器系数来实现这一目标。
具体设计过程如下：3.1 确定滤波器参数在设计FIR滤波器时，需要确定一些关键参数，包括：（1）采样频率：即采样器的采样率，一般情况下为信号最高频率的两倍以上。
（2）截止频率：即希望滤除的信号部分的最高频率。
（3）滤波器阶数：阶数越高，滤波器的频率响应越陡峭。
（4）窗口长度：即用于窗函数的样本点数。
3.2 选择窗函数在滤波器设计中，窗函数的选择对于滤波器性能有很大的影响。
通常，我们可以选择常用的窗函数，如矩形窗函数、汉明窗函数、升余弦窗函数等。
在本文中，我们将选择汉明窗函数作为窗口函数。
3.3 计算滤波器系数通过选择合适的窗口函数，我们可以得到对应的窗口系数，然后将其与所需滤波器类型（低通、高通、带通等）的理想频率响应进行卷积，即可得到所需的滤波器系数。
用窗函数法设计FIR数字低通滤波器要点
河北科技大学课程设计报告学生姓名: 学号:专业班级：课程名称:学年学期指导教师:20年月课程设计成绩评定表目录1. 窗函数设计低通滤波器1.1设计目的 (1)1.2设计原理推导与计算 (1)1．3设计内容与要求 (2)1.4设计源程序与运行结果 (3)１.5思考题……………………………………………………………………１01．6心得体会 (14)参考文献……………………………………………………………………… 1５１.窗函数设计低通滤波器1.1设计目的1. 熟悉设计线性相位数字滤波器的一般步骤。
2. 掌握用窗函数法设计FIR 数字滤波器的原理和方法。
3. 熟悉各种窗函数的作用以及各种窗函数对滤波器特性的影响。
4. 学会根据指标要求选择合适的窗函数。
1.2设计原理推导与计算如果所希望的滤波器的理想的频率响应函数为()ωj d e H ，则其对应的单位脉冲响应为()()ωπωωππd e e H n h j j d d ⎰-=21 （４.1）窗函数设计法的基本原理是设计设计低通FIR 数字滤波器时，一般以理想低通滤波特性为逼近函数()ωj e H ，即()⎪⎩⎪⎨⎧≤<≤=-πωωωωωαωc c j jd ，，e e H 0，其中21-=N α()()()[]()a n a n d e e d e eH n h c j j j j d d cc--===⎰⎰---πωωπωπωαωωωαωππωsin 2121用有限长单位脉冲响应序列()n h 逼近()n h d 。
由于()n h d 往往是无限长序列，而且是非因果的，所以用窗函数()n ω将()n h d 截断,并进行加权处理，得到:()()()n n h n h d ω=(4．2)()n h 就作为实际设计的FI Ｒ数字滤波器的单位脉冲响应序列，其频率响应函数()ωj e H 为()()nj N n j en h eH ωω∑-==1ﻩ ﻩ（4.3)式中,N 为所选窗函数()n ω的长度。
利用窗函数设计FIR数字滤波器
设 计 分析从 性 能表 可 以看 出 汉 宁 窗，海 明窗和 凯瑟 窗可 提供 的衰 减
大 于 ４ｄ ，但 是汉 宁 窗的 旁瓣 峰值 较 小 ，而 主瓣 宽度 和海 明窗一 样 ，可 以 ０Ｂ 使 滤波 器具 有较 少 的阶数 ，所 以选择 汉 宁窗 。用汉 宁窗 函数 设计低 通Ｆ Ｒ Ｉ滤
波 器 的过程 如下 ：
ＦＥ Ｚ Ｒ Ｑ Ｍ． Ｍ
Ｉ
ｙ ｎ ＝ ．（ ｘｎ ） （） ｈ ｉ （ —ｆ ）
ｉ０ ＝
（） １
ｆ ｎ ｔ ｏ ［ ｂ ａ ，ｐ ａ ｒ ，Ｗ ＝ ｒ ｑ ｍ（，ａ ｕ ｃ ｉ ｎ ｄ ，ｍ ｇ ｈ ，ｇ ｄ ］ ｆ ｅ ｚ ｂ ）
％ｄ 为相 对 振 幅 （Ｂ ：ｍｇ 绝对 振 幅 ：ｐ ａ ｂ ｄ） ａ 为 ｈ 为相 位 响应 ；ｇｄ 群 延 ｒ为 时 ；Ｗ 为频 率样 本 点向量 。
关键词： Ｆ Ｒ 字滤 波器设计 ；窗函数 ；实例 分析 Ｉ数 中图分类号 ：Ｔ ７ 文献标 识码：Ａ 文 章编号：１ ７ －７ ９ ２ １ ）０ １ ０ ４ ０ Ｎ ６ １ ５ ７（ ０ ０ ８ ０ ４ － １
数 字信 号处理 在通信 、语音 、 图像 ， 自动 控 制 、雷达 、军 事 、航 空航
ｄ ＝ ０ ｌ ｇ Ｏ （ａ ＋ ｐ ） ａ （ ａ ） ； ｂ ２ ＊ ｏ ｌ （ｍ ｇ ｅ ｓ ／ｍ ｘ ｍ ｇ ） ｐａ ａｇ ｅＨ ： ｈ ＝ ｎ ｌ （）
最后 利用 计算 机仿 真 实现 设计 的滤 波器 ， 再分析 滤 波结 果来 判断 设计
是否 满足 要求 。工程 实 际中常 用 的窗 函数有 ５ ， 即矩形 窗 、三角 窗、汉 宁 种
１ Ｉ数 字滤 波器 的介 绍 １Ｆ Ｒ
实验6FIR滤波器设计
实验6FIR滤波器设计FIR (Finite Impulse Response)滤波器是一种数字滤波器，其输出信号仅取决于振荡器的输入以前的有限个值。
FIR滤波器设计的目的是通过调整滤波器的系数以实现所需的频率响应。
在FIR滤波器设计中，首先确定滤波器的类型和频率响应的规格。
常见的滤波器类型有低通滤波器、高通滤波器、带通滤波器和带阻滤波器。
频率响应的规格由滤波器的截止频率、通带增益和阻带衰减等参数决定。
FIR滤波器的设计步骤如下：1.确定滤波器的类型和频率响应规格。
根据应用的需求，选择适当的滤波器类型和定义频率响应的参数。
2.确定滤波器的阶数。
阶数决定了滤波器的复杂度和性能。
一般而言，阶数越高，滤波器的性能越好，但计算复杂度也越高。
3.根据频率响应规格和系统设计的约束，选择一种滤波器设计方法。
常见的设计方法有窗函数法、频率采样法、最小均方误差法等。
4.设计滤波器的理想频率响应。
根据所选的设计方法，确定滤波器的理想频率响应。
这通常是一个分段线性函数，其中包括通带增益和阻带衰减。
5.将理想频率响应转换为时域的冲激响应。
这可以通过将理想频率响应进行反傅里叶变换来实现。
6.通过选择合适的窗函数，对冲激响应进行窗函数变换。
窗函数的选择是设计滤波器性能的重要因素。
7.通过窗函数变换得到滤波器的系数。
通过将窗函数变换应用于冲激响应，可以得到设计滤波器的系数。
这些系数确定了滤波器的时间响应和频率响应。
8.可选地，通过优化算法对滤波器的系数进行优化。
优化算法可以用来进一步改善滤波器的性能。
常用的优化算法包括加权最小二乘方法、梯度下降法等。
9.实现滤波器。
将设计好的滤波器系数应用于输入信号，得到滤波器输出。
可以使用编程语言或滤波器设计工具来实现滤波器。
10.验证滤波器的性能。
通过将滤波器应用于不同的输入信号，检验滤波器输出是否符合设计要求。
可以使用频谱分析工具和滤波器性能评估指标来评估滤波器的性能。
FIR滤波器设计是数字信号处理中重要的课题之一、设计一个性能良好的FIR滤波器需要对滤波器原理和设计方法有深入的了解，以及熟练的使用滤波器设计工具和编程工具。
用窗函数法设计FIR滤波器
1．用窗函数法设计一线性相位FIR低通滤波器，要求通带截止频率 ，
（1）选择一个合适的窗函数（如hamming窗），取单位冲击响应h（n）的长度N＝15，观察所设计滤波器的幅频特性，分析是否满足设计要求；
（2）取N＝45，重复上述设计，观察幅频和相频特性的变化，分析长度N变化的影响；
（3）保持N＝45不变，改变窗函数（如hamming窗变为blackman窗），观察并记录窗函数对滤波器幅频特性的影响。
xlabel('n');ylabel('h(n)');
title('hamming窗设计的h（n）'2);
hw=fft(hn,512);
w=2*[0:511]/512;
plot(w,20*log10(abs(hw)));
xlabel('w/pi');ylabel('Magnitude(dB)');
（4）由 ，得出单位脉冲响应 ;
（5）对 作离散时间傅立叶变换，得到 。
2．在MATLAB中，可以用b=fir1(N,Wn,’ftype’,taper)等函数辅助设计FIR数字滤波器。N代表滤波器阶数；Wn代表滤波器的截止频率(归一化频率)，当设计带通和带阻滤波器时，Wn为双元素相量；ftype代表滤波器类型，如’high’高通，’stop’带阻等；taper为窗函数，默认为海明窗，窗函数实现需要用窗函数blackman,hamming,hanningchebwin,kaiser产生。
用窗函数法设计FIR滤波器是在时域进行的,先用傅里叶变换求出理想滤波器单位抽样相应hd(n),然后加时间窗w(n)对其进行截断,以求得FIR 滤波器的单位抽样响应h(n)。
数字信号处理第五章 FIR滤波器、窗函数法
(5.23)
0
h(n)偶对称N为偶数时的频 率响应特性
不适用高通或带阻滤波器
2）H(w)以
为奇对称，以
偶对称
以
为奇对称，以
为偶对称
22
1. 成立条件
5.2.3 恒群延时单独成立
解线性微分方程
＝
的取值条件？
0
的取值条件？
有什么特性？
(5.24)
23
离散傅立叶变换： 由相位定义： 化简得：
(5.25)
n
28
相位函数：
(5.30) 线性相位
幅度函数：
(5.31)
其中
1）H(w)在
关于
处奇对称
奇对称
2）传输函数H(z)在
处有零点
0
h(n)奇对称N为奇数时的频 率响应特性
不适用于低通、高通、带阻滤波器
29
3. h(n)奇对称N为偶数（线性相位IV型）时的频率响应
所以
456 7
0123
n
(5.36)
数字信号处理 Digital Signal Processing
第五章 FIR滤波器设计和实现
Finite Impulse Response Filter Design and Implementation
2
FIR数字滤波器
• Finite Impulse Response (FIR)：有限冲激响应，非递归 滤波器，输出只与当前和有限个过去输入有关
偶对称
0 h(n)偶对称N为奇数时的 频率响应特性
20
3. h(n)偶对称N为偶数（线性相位II型）时的频率响应
分拆
0123456 7 n
变量代换及
FIR数字滤波器设计窗函数法与频率抽样法比较
ＦＩＲ数字滤波器设计窗函数法与频率抽样法比较（安阳工学院，河南安阳４５５０００）姚海燕摘要：用ＭＡＴＬＡＢ分析实例来比较窗函数设计法与频率取样法，并通过比较结果给出在ＦＩＲ数字滤波器的设计中选择设计方法的一个参考。
关键词：ＦＩＲ数字滤波器；窗函数；频率抽样；ＭＡＴＬＡＢ中图分类号：ＴＮ９１１文献标识码：Ａ文章编号：１６７３－２９２８（２００７）０６－００５１－０３收稿日期：２００７－０１－１７作者简介：姚海燕（１９８１－），女，安阳市人，安阳工学院教师。
研究方向：数字信号处理。
１ＦＩＲ数字滤波器数字滤波器是对频率有一定筛选作用的离散时间系统。
分为ＩＩＲ数字滤波器与ＦＩＲ数字滤波器两种。
其中ＩＩＲ数字滤波器可以借助比较成熟模拟滤波器设计的方法来设计，设计方法比较简单，又由于其零极点分布比较任意，所以阶次较低，但是由模拟到数字的变换导致其相位非线性，另外由于其冲激响应无限长，因而不能用ＦＦＴ来计算。
ＦＩＲ数字滤波器与ＩＩＲ的相比，由于极点固定在原点，使其阶次较高，但是ＦＩＲ数字滤波器能做到严格的线性相位，可以设计任意的频率特性，又因为其为有限长的冲激响应，可以用ＦＦＴ计算，再配合计算机辅助设计，其应用越来越广泛。
２窗函数设计法与频率抽样法的分析比较窗函数设计法与频率抽样法是ＦＩＲ数字滤波器设计的两种常用方法。
窗函数设计法比较简单，有闭合形式的公式可循［１］。
频率抽样法可以在频域直接设计，在抽样点处与理想滤波器严格相等，可以设计任意幅度响应的滤波器。
设计ＦＩＲ数字滤波器时，在选择设计方法时，设计出的滤波器的阶次也是我们作出选择时要着重考虑的。
那么就需要对这两种方法设计的ＦＩＲ数字滤波器进行比较。
下面通过一个实例来说明。
例：设计一个线性相位的ＦＩＲ数字低通滤波器，要求通带频率０．３π，衰减不超过５ｄＢ，阻带频率０．４π阻带衰减超过５０ｄＢ。
２．１采用窗函数设计法窗函数设计法是根据设计的性能要求，选择一个理想滤波器，然后用一个合适的窗函数，与理想滤波器在时域中的单位冲激响应相乘，即所谓加窗，得到一个有限长的冲激响应的数字系统，通过调整窗函数的参数来逼近理想滤波器的性能参数，从而达到设计要求。
数字信号处理实验——用窗函数设计FIR滤波器
实验四 用窗函数设计FIR 滤波器一、 实验目的1、熟悉FIR 滤波器设计的基本方法。
2、掌握用窗函数设计FIR 数字滤波器的原理及方法，熟悉相应的计算机高级语言编程。
3、熟悉线性相位FIR 滤波器的幅频特性和相位特性。
4、了解各种不同窗函数对滤波器性能的响应。
二、 实验原理和方法窗函数法设计的任务在于寻找一个可实现有限长单位脉冲响应的传递函数H(e jw )=∑-=10N n h(n)e -jwn 去逼近h d (n)=1/2π⎰π20H d (e jw )e jwn dw即h(n)=h d (n)w （n ） (一)几种常用的窗函数1、矩形窗 w(n)=R N (n)2、Hanning 窗 w(n)=0.5[1-cos(2πn ／N-1)]R N (n)3、Hamming 窗 w(n)=[0.54-0.46cos(2πn ／N-1)]R N (n)4、Blackman 窗 w(n)=[0.42-0.5 cos(2πn ／N-1)+0.08 cos(4πn ／N-1)] R N (n)5、Kaiser 窗 w(n)=I 0(β（1-[(2n ／(N-1))-1]2）½)／I 0(β）（二）窗函数法设计线性相位FIR 滤波器的步骤1、确定数字滤波器的性能要求。
确定各临界频率{w k }和滤波器单位脉冲响应长度N 。
2、根据性能要求和N 值，合理地选择单位脉冲响应h(n)有奇偶对称性，从而确定理想频率响应h d (e jw )的幅频特性和相位特性。
3、用傅里叶反变换公式求得理想单位脉冲响应h d (n)。
4、选择适当的窗函数W （n ），求得所设计的FIR 滤波器单位脉冲响应。
5、用傅里叶变换求得其频率响应H (e jw )，分析它的幅频特性，若不满足要求，可适当改变窗函数形式或长度N ，重复上述过程，直至得到满意的结果。
三、实验内容和步骤1、分别用矩形窗、Hanning 窗、Hamming 窗、Blackman 窗、Kaiser 窗（β=8.5）设计一个长度N=8的线性相位FIR 滤波器。
FIR滤波器设计与实现实验报告
FIR滤波器设计与实现实验报告目录一、实验概述 (2)1. 实验目的 (3)2. 实验原理 (3)3. 实验设备与工具 (4)4. 实验内容与步骤 (6)5. 实验数据与结果分析 (7)二、FIR滤波器设计 (8)1. 滤波器设计基本概念 (9)2. 系数求解方法 (10)频谱采样法 (11)最小均方误差法 (14)3. 常用FIR滤波器类型 (15)线性相位FIR滤波器 (16)非线性相位FIR滤波器 (18)4. 设计实例与比较 (19)三、FIR滤波器实现 (20)1. 硬件实现基础 (21)2. 软件实现方法 (22)3. 实现过程中的关键问题与解决方案 (23)4. 滤波器性能评估指标 (25)四、实验结果与分析 (26)1. 实验数据记录与处理 (27)2. 滤波器性能测试与分析 (29)通带波动 (30)虚部衰减 (31)相位失真 (32)3. 与其他设计方案的对比与讨论 (33)五、总结与展望 (34)1. 实验成果总结 (35)2. 存在问题与不足 (36)3. 未来发展方向与改进措施 (37)一、实验概述本次实验的主要目标是设计并实现一个有限脉冲响应（Finite Impulse Response，简称FIR）滤波器。
FIR滤波器是数字信号处理中常用的一种滤波器，具有线性相位响应和易于设计的优点。
本次实验旨在通过实践加深我们对FIR滤波器设计和实现过程的理解，提升我们的实践能力和问题解决能力。
在实验过程中，我们将首先理解FIR滤波器的基本原理和特性，包括其工作原理、设计方法和性能指标。
我们将选择合适的实验工具和环境，例如MATLAB或Python等编程环境，进行FIR滤波器的设计。
我们还将关注滤波器的实现过程，包括代码编写、性能测试和结果分析等步骤。
通过这次实验，我们期望能够深入理解FIR滤波器的设计和实现过程，并能够将理论知识应用到实践中，提高我们的工程实践能力。
本次实验报告将按照“设计原理设计方法实现过程实验结果与分析”的逻辑结构进行组织，让读者能够清晰地了解我们实验的全过程，以及我们从中获得的收获和启示。
Matlab中的滤波器设计方法详解
Matlab中的滤波器设计方法详解滤波器在信号处理中起着至关重要的作用，能够对信号进行去噪、去除干扰、增强所需频谱等操作。
Matlab作为一种强大的数学计算工具，提供了丰富的滤波器设计方法。
本文将详细介绍Matlab中常用的滤波器设计方法，包括FIR和IIR滤波器的设计原理和实现步骤。
一、FIR滤波器的设计方法FIR滤波器全称为有限脉冲响应滤波器，其特点是具有线性相位和稳定性。
在Matlab中，常用的FIR滤波器设计方法有窗函数法、最小二乘法和频率抽取法。
1. 窗函数法窗函数法是最简单直观的FIR滤波器设计方法。
它的基本思想是，在频域上用一个窗函数乘以理想滤波器的频率响应，再进行频域到时域的转换，得到滤波器的冲激响应。
常用的窗函数有矩形窗、汉宁窗、汉明窗等。
Matlab中，我们可以使用fir1函数进行窗函数法滤波器设计。
该函数的输入参数包括滤波器阶数、归一化截止频率和窗函数类型。
通过设计不同阶数和不同窗函数的滤波器，可以得到不同性能的滤波器。
2. 最小二乘法最小二乘法是一种优化方法，通过最小化滤波器的输出与目标响应之间的均方误差来设计滤波器。
在Matlab中，我们可以使用fir2函数进行最小二乘法滤波器设计。
该函数的输入参数包括滤波器阶数、频率向量和响应向量。
通过调整频率向量和响应向量，可以实现对滤波器的精确控制。
3. 频率抽取法频率抽取法是一种有效的FIR滤波器设计方法，能够实现对特定频带的信号进行滤波。
在Matlab中，我们可以使用firpm函数进行频率抽取法滤波器设计。
该函数的输入参数包括滤波器阶数、频率向量、增益向量和权重向量。
通过调整频率向量、增益向量和权重向量，可以实现对滤波器的灵活控制。
二、IIR滤波器的设计方法IIR滤波器全称为无限脉冲响应滤波器，其特点是具有非线性相位和多项式系数。
在Matlab中，常用的IIR滤波器设计方法有巴特沃斯滤波器、切比雪夫滤波器和椭圆滤波器。
1. 巴特沃斯滤波器巴特沃斯滤波器是一种最常用的IIR滤波器，其特点是具有最平滑的通带和最陡峭的阻带。
FIR数字滤波器两种设计方法的比较
Science &Technology Vision 科技视界0前言FIR 滤波器能严格做到线性相位和群时延为常数,并且不会因为滤波运算的舍入误差而产生极限环振荡现象。
但它与IIR 滤波器相比存在需要较高节数和滤波过程时间较长的缺点。
常见的FIR 滤波器设计方法有窗函数法,频率采样法,和切比雪夫等波纹优化设计法。
本文将对窗函数法和切比雪夫等波纹优化设计法进行比较。
[1]1设计方法比较1.1窗函数法窗函数法的设计是在时域中进行的。
设理想滤波器的单位脉冲响应为h d (n ),则有:h d(e jω)=∞n=-∞∑h d(n )e -jnω(1-1)h d (n )=12ππ-π∫H d (e jω)e jnωdω(1-2)因而,一旦H d (e jω)给定,就可求得h d (n ),但这样求得的h d (n )一般是无限长,且是非因果的。
为了得到一个因果的有限长的滤波器h (n ),要使用一个窗口函数ω(n )对h d (n )进行加窗处理,即:h (n )=h d (n )ω(n )(1-3)因此选择窗口函数的形状和长度就成了窗口函数法的关键。
表1列出了6种窗函数基本参数的比较。
表16种窗函数基本参数的比较利用窗函数设计FIR 滤波器的过程可设计如下:1)利用式(1-2),由给定的滤波器的幅频响应参数求出理想的单位脉冲响应h d (n )。
2)按允许的过渡带宽度Δω及其阻带衰减A s =-201g δs ,选择合适的窗函数ω(n ),并估计节数N,N =A Δω,其中A 由窗函数的类型决定。
3)确定延时值α=N -12(滤波器的中心)。
4)求h (n )=h d (n )ω(n )。
5)必要时验算FIR 滤波器的频率响应。
1.2等波纹设计法我们所希望设计的滤波器幅度响应为H d (ω),实际逼近的幅度响应为H g (ω)。
线性相位FIR DF 根据单位抽样响应h (n )的奇偶对称性以及h (n )的长度N 的奇偶性,总共可以分为四种类型。
数字信号处理---窗函数法设计FIR滤波器,窗函数选取原则
)
)]
W Hm ( ) 0.54W R ( ) 0.23[W R (
2 N 1
) W R (
2 N 1
)]
15 Logo
请加入标题 窗函数介绍
汉明窗
16 Logo
请加入标题 窗函数介绍
汉明窗
主瓣宽度 B m 8 / 23 1.092 -41dB
旁瓣峰值衰减
幅频响应窗函数介绍布莱克曼窗04205cos18频率响应窗函数介绍布莱克曼窗19窗函数介绍布莱克曼窗20主瓣宽度旁瓣峰值衰减57db阻带最小衰减74db精确过渡带12231638231502凯塞窗kaiserwindow窗函数介绍21其中22窗函数介绍23旁瓣峰值衰减57db阻带最小衰减80db精确过渡带1023136524各种窗函数性能比较窗类型主瓣宽度精确过旁瓣峰值衰减db阻带最小衰减db矩形窗18n1321三角窗61n252562n314466n4153布莱克曼窗12n11n577410n5780各种窗函数性能比较克曼窗阻带最小衰减较大凯塞窗的阻带最小衰减最大各窗函数的旁瓣峰值衰减规律与此关系相同
精确过渡带 Bt 6.6 / 23 0.9010 阻带最小衰减 53dB
17 Logo
布莱克曼窗
窗函数介绍
布莱克曼窗(Blankman Window) ，又
称二阶升余弦窗。
w Bl ( n ) [0.42 0.5 cos( 2 n N 1 ) 0.08 cos(
j
4 n N 1
各种窗函数性能比较
窗类型
矩形窗 三角窗
主瓣宽度
精确过
渡带 B t
1.8π/N 6.1π/N
旁瓣峰值
衰减/dB
-13 -25
用窗函数法设计 FIR 数字滤波器实验报告及代码展示
本科生实验报告数字信号处理 课 程 实 验 报 告实验名称 用窗函数法设计 FIR 数字滤波器 一、实验原理、目的与要求1. 实验原理如果所希望的滤波器的理想频率响应函数为，则其对应的单位脉冲响应为：用窗函数w(n)将)(d n h 截断，并进行加权处理，得到：h(n)就作为实际设计的FIR 数字滤波器的单位脉冲响应序列，其频率响应函数）（jw H e 为：式中，N 为所选窗函数w(n)的长度。
用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N 。
因此，在设计过程中，要根据对阻带最小衰减和过度带宽度的要求选择合适的窗函数类型和窗口长度N 。
选定窗函数了形和长度N 后，求出单位脉冲响应h(n)=hd(n)·w(n)，并可以求出）（jw H e 。
）（jw H e 是否满足要求，要进行验算。
一般在h(n)尾部加零使长度满足2的整数次幂，以便用FFT 计算）（jw H e 。
如果要观察细节，补零点数增多即可。
如果）（jw H e如果要求线性相位特性，则h(n)还必须满足：根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。
要根据所设计的滤波特性正确选择其中一类。
例如，要设计线性相位低通特性，可选择h(n)=h(N-1-n)一类，而不能选h(n)=-h(N-1-n)一类。
2. 实验目的（1）掌握用窗函数法设计 FIR 数字滤波器的原理和方法。
（2）熟悉线性相位 FIR 数字滤波器特性。
（3）了解各种窗函数对滤波特性的影响。
3. 实验要求（1）简述实验目的及原理。
（2）按照实验步骤及要求，比较各种情况下的滤波性能，说明窗口长度 N 和窗函数类型对滤波特性的影响。
（3）总结用窗函数法设计 FIR 滤波器的主要特点。
（4）简要回答思考题。
二、实验仪器设备（标注实验设备名称及设备号）Windows 计算机台号 22Matlab 软件三、实验内容步骤及结果分析1.用升余弦窗设计一线性相位低通 FIR 数字滤波器，截至频率wc = π/ 4 rad。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
课题名称：FIR滤波器窗函数设计FlR 滤波器窗函数设计引言：数字滤波器（DigitalFilter ）是指输入、输出都是离散时间信号，通过一定运算 关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。
在许 多数字信号处理系统中，如图像信号处理等，有限冲激响应（ FIR ）滤波器是最 常用的组件之一，它完成信号预调、频带选择和滤波等功能。
FIR 滤波器虽然在 截止频率的边沿陡峭性能上不及无限冲激响应 （IIR ）滤波器，但是却具有严格的 线性相位特性，稳定性好，能设计成多通带（或多阻带）滤波器组，所以能够在 数字信号处理领域得到广泛的应用。
数字滤波器的分类1） 根据系统响应函数的时间特性分为两类1. FIR (Finite ImPUISe Response 数字滤波器网络My[n]b k x[n k]k0特点：不存在反馈支路，其单位冲激响应为有限长2. IIR ( Infinite ImPUISe Response 数字滤波器网络特点：存在反馈支路，即信号流图中存在环路，其单位冲激响应为无限长（2） FIR 数字滤波器IIR 数字滤波器的区别1. 从性能上来说，IlR 滤波器传递函数包括零点和极点两组可调因素，对极点的惟一限制是在单位圆内。
因此可用较低的阶数获得高的选择性，所用的存储 单元少，计算量小，效率高。
但是这个高效率是以相位的非线性为代价的。
选择性越好，则相位非线性越严重。
FIR 滤波器传递函数的极点固定在原点， 是不能动的，它只能靠改变零点位置来改变它的性能。
所以要达到高的选择 性，必须用较高的阶数；对于同样的滤波器设计指标， FIR 滤波器所要求的 阶数可能比IIR 滤波器高5-10倍，但是FIR 滤波器可以得到严格的线性相位。
2. 从结构上看，IIR 滤波器必须采用递归结构，极点位置必须在单位圆内，否则 系统将不稳定。
相反，FIR 滤波器只要采用非递归结构，不论在理论上还是 在实际的有限精度运算中都不存在稳定性问题， 因此造成的频率特性误差也 较小。
此外FIR 滤波器可以采用快速傅里叶变换算法， 在相同阶数的条件下， 运算速度可以快得多。
3. 从设计工具看，IIR 滤波器可以借助于模拟滤波器的成果，因此一般都有有效 的圭寸闭形式的设计公式可供准确计算，计算工作量比较小，对计算工具的要 求不高。
hnb n, 0 n M 0, 其他 ny[n]b k x[n k]k0a k y[n k1k]从上面的简单比较可以看到IIR与FIR滤波器各有所长，所以在实际应用时应该从多方面考虑来加以选择。
从使用要求上来看，在对相位要求不敏感的场合，如语言通信等，选用IIR较为合适，这样可以充分发挥其经济高效的特点；对于图像信号处理，数据传输等以波形携带信息的系统，则对线性相位要求较高。
如果有条件，采用FIR滤波器较好。
另外，不论IIR和FIR阶数越高，信号延迟越大；同时在IIR滤波器中，阶数越高，系数的精度要求越高，否则很容易造成有限字长的误差使极点移到单位园外。
因此在阶数选择上是综合考虑的。
FIR数字滤波器的设计(一)FIR滤波器原理N 1y[n] h(k)x(n k)k 0FIR滤波器的系统输入输出差分方程为：所以FIR滤波器的系统函数为：H(Z)Y⑵X(Z)N 1h(n )z nn O由于FIR滤波器的单位脉冲响应h(n)是一个有限长序列，H (Z)是- 的(N-I)次多项式，它在Z平面上有(N-1)个零点，同时在原点有(N-I)阶重极点。
因此，H(Z) 永远稳定。
FIR滤波器设计的任务是选择有限长度的h( n),使传输函数覺肿笃满足一定的幅度特性和线性相位要求。
由于FIR⅛波器很容易实现严格的线性相位，所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。
(二) FlR滤波器的窗函数设计原理窗函数法的设计思想是按照所要求的理想滤波器频率响应IwP,设计个FIR⅛波器，使之频率响应来逼近%討逹O先由%廳严》的⅛-0傅里叶反变换导出理想滤波器的冲激响应序列沁阳,即:∖e^ωτιdω-π由于*加謝胛⅛是矩形频率特性，所以一，：一是一无限长的序列，且是非因果的，而要计的FIR 滤波器的冲激响应序列是有限长的，所以要用有限长的序列 h(n)来逼近无限长的序列 诂劇,最有效的方法是截断匸屈初，或者说用一个有限长度的窗 口函数w( n)序列来截取g 蜃胡，即：一。
按照复卷积公式，在时域中的乘积关系可表示成在频域中的周期性卷积关系， 即可得所设计的FIR 滤波器的频率响应：其中，箱?涉乜为截断窗函数的频率特性。
由此可见，实际的FIF 数字滤波器的频率响应盼笃逼近理想滤波器频率响应 碌贋卜F 的好坏，完全取决于窗函数的频 率特性…严》如果w(n)具有下列形式:0, n 0,n N w(n)1, OnNw( n)相当于一个矩形，我们称之为矩形窗。
即我们可采用矩形窗函数w( n)将无限 脉冲响应讪術截取一段脯爛来近似为乳述爛。
经过加矩形窗后所得的滤波器实际 频率响应能否很好地逼近理想频率响应呢？下图给出了理想滤波器加矩形窗后 的情况。
理想低通滤波器的频率响应'WJ 如图中左上角图，矩形窗的频率响 应•吩％为左下角图。
根据卷积定理，即得实际滤波器的频率响应 蹭加乍图形 为图中右图。
由图可看出，加矩形窗后使实际频率响应偏离理想频率响应， 主要影响有三 个方面：(1) 理想幅频特性陡直边缘处形成过渡带，过渡带宽取决于矩形窗函数频 率响应的主瓣宽度。
f H ci (評)Mz o I d)--T^)dθ(2)过渡带两侧形成肩峰和波纹，这是矩形窗函数频率响应的旁瓣引起的，旁瓣相对值越大，旁瓣越多，波纹越多。
(3)随窗函数宽度N的增大，矩形窗函数频率响应的主瓣宽度减小，但不改变旁瓣的相对值。
为了改善滤波器的性能，需使窗函数谱满足：主瓣尽可能窄，以使设计出来的滤波器有较陡的过渡带；第一副瓣面积相对主瓣面积尽可能小，即能量尽可能集中在主瓣，外泄少，使设计出来的滤波器的肩峰和余振小逼近于理想滤波器。
但是这两个条件是相互矛盾的，实际应用中，折衷处理，兼顾各项指标。
(三) MATLAB信号处理中提供的窗函数上边只考虑了矩形窗，如果我们使窗的主瓣宽度尽可能地窄，旁瓣尽可能地小，可以获得性能更好的滤波器，通过改变窗的形状来达到这个目的。
在数字信号处理的发展过程中形成了不同于矩形窗的很多窗函数，这些窗函数在主瓣和旁瓣特性方面各有特点，可满足不同的要求。
为此，用窗函数法设计器时，要根据给定的滤波器性能指标选择窗口宽度介绍几类类窗函数及其特性。
1. 矩形窗矩形窗函数的时域形式可以表示为:它的频域特性为:FIR数字滤波w(n)。
下面具体N和窗函数1, w(n) R N (n) O O n 其他它的频域特性为:.N 1 JWe j eSinSin 2. 汉宁窗函数汉宁窗函数的时域形式可以表示为:w(k) 0.5 1cos 2亠1,2, , N0.5W R0.25 W R 2πW R其中，w ()为矩形窗函数的幅度频率特性函数。
汉宁窗函数的最大旁瓣值比主瓣值低 31dB,但是主瓣宽度比矩形窗函数的主 瓣宽度增加了 1倍，为8π /N 。
3. 海明窗函数海明窗函数的时域形式可以表示为它的频域特性为：其中，W R ()为矩形窗函数的幅度频率特性函数。
海明窗函数的最大旁瓣值比主瓣值低 41dB,但它和汉宁窗函数的主瓣宽度是 一样大的。
4. 三角窗函数形式可以表示为:当n 为奇数时:当n 为偶数时:2k 1n, 1 k - n 22(n k 1) n-k n它的频域特性为:w(k)0.54 0.46COS 2 πk 1,2, ,NW( )0.54W R ( ) 0.23 W R2π N 1W R2π N 1三角窗是最简单的频谱函数 悴討'迪为非负的一种窗函数。
三角窗函数的时域w(k)2kn 11 k n 12 2(n k 1)n 1 Ik nn 1,2w(k)24Sin —2三角窗函数的主瓣宽度为 8π /N ,比矩形窗函数的主瓣宽度增加了一倍，但 是它的旁瓣宽度却小得多。
(四) 各窗函数的图形及幅频特性1.用MATLA 编程绘制各种窗函数的形状，窗函数的长度为21。
程序如下：clf;NWi n=21; n=0:NWi n-1; figure(1) for ii=1:4 SWitCh ii CaSe 1w=boxcar(Nw in); %矩形窗stext=' 矩形窗'； CaSe 2w=ha nnin g(Nw in); stext='汉宁窗'；CaSe 3% 汉宁窗w=hammi ng(Nwi n); stext='海明窗'；CaSe 4% 海明窗w=tria ng(Nwi n); stext='三角窗'；endPOSPlOt=['2,2,' i nt2str(ii)]; subplot(posplot); % 三角窗Stem (n ,w); % hold On绘出窗函数 plot (n ,w,'r'); %绘出包络线xlabel(' n');ylabel('w( n)');title(stext); hold off;grid on;end程序运行结果如下图:SinW R e j汉宁窗它⅛,f05n:严坪J…Φ≡ 0.50 5 10 15 20n°0 5 10 15n2.用MATLAB编程，clf;Nf=512;NWin=20; figure(1) for ii=1:4SWitCh iiCaSe 1w=boxcar(Nw in); stext=' 矩形窗'；CaSe 2w=ha nnin g(Nw in); stext=' 汉宁窗'；CaSe 3w=hammi ng(Nwi n); stext=' 海明窗'；CaSe 4w=tria ng(Nwi n); stext=' 三角窗'；end[y,f]=freqz(w,1,Nf); % mag=abs(y);POSPlOt=['2,2,' i nt2str(ii)];subplot(posplot);plot(f∕pi,20* Iog10(mag∕max(mag))); %采用512个频率点绘制各窗函数的幅频特性。
程序如下:%%窗函数复数频率特性的数据点数窗函数数据长度矩形窗汉宁窗海明窗三角窗求解窗函数的幅频特性%求得窗函数幅频特性绘制窗函数的幅频特性XIabe1（'归一化频率’）ylabel （' 振幅∕dB'）;title （stext ）;grid on; end程序运行结果如下图：窗函数主瓣频宽第旁瓣相对主瓣衰减（dB ）矩形窗 4τt ΛΓ-13 汉宁窗 8TFJV-31 海明窗 BπJV -41 三角窗8ττ JW F-25-50-100-1505J-15□l ---------------------- ; --------------------- -2000 05 1 0 0,51归一化頻率归一化频率由上图各窗函数的幅频特性可以看到各种窗函数都具有明显的主瓣和旁瓣。