基于矩形窗设计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数字滤波器
郑州航空工业管理学院《电子信息系统仿真》课程设计14级专业班级题目基于矩形窗设计FIR数字滤波器姓名学号二О一六年十一月二十五日第一章 FIR 滤波器的设计原理及方法FIR 滤波器通常采用窗函数方法来设计。
窗设计的基本思想是，首先选择一个适当的理想选频滤波器（它总是具有一个非因果，无限持续时间脉冲响应），然后街区（加窗）它的脉冲响应得到线性相位和因果FIR 滤波器。
我们用Hd （e^jw ）表示理想的选频滤波器，它在通带上具有单位增益和线性相位，在阻带上具有零响应。
一个带宽wc<pi 的低通滤波器由下式给定：⎩⎨⎧=<-⋅<==<=wc w jaw e w wc jw e Hd ||),(^1||,0)^(π为了从hd(n)得到一个FIR 滤波器，必须同时在两边截取hd （n ）。
而要得到一个因果的线性相位滤波器，它的h(n)长度为N ，必须有：⎩⎨⎧-<=<==其他,010),()(N n n hd n h 2/)1(-=N α这种操作叫做加窗，h(n)可以看做是hd(n)与窗函数w （n ）的乘积： h(n)=hd(n)w(n) 其中⎩⎨⎧-≤≤=，其他对称，关于01n 0)(N n w α根据w （n ）的不同定义，可以得到不同的窗结构。
在频域中，因果FIR 滤波器响应H(e^jw)由Hd(e^jw)和窗响应W(e^jw)的周期卷积得到，即⎰--=⊗=ππλλωπd j e W jw e W jw e Hd jw e H )(^(2/1)^()^()^(常用的窗函数有矩形窗、巴特利特（BARTLETT）窗、汉宁（HANNING）窗、海明（HAMMING）窗、布莱克曼（BLACKMAN）窗、凯泽（KAISER）窗等。
FIR滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。
窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列，主要设计步骤为：(1) 通过傅里叶逆变换获得理想滤波器的单位脉冲响应hd(n)。
信号实验三 用窗函数设计FIR滤波器
实验三用窗函数设计FIR滤波器一、实验目的1.熟悉FIR滤波器设计的基本方法。
2.掌握用窗函数设计FIR数字滤波器的原理及方法，熟悉相应的计算机高级语言编程。
3.熟悉线性相位FIR滤波器的幅频特性和相位特性。
4.了解各种不同窗函数对滤波器性能的响应。
二、实验内容1、以下是几种常用的窗函数：1.矩形窗：2.Hanning（汉宁）窗：3.Hamming（海明）窗：4.Blackman窗：5.Kaiser窗：窗函数法设计线性相位FIR滤波器可以按如下步骤进行：1．确定数字滤波器的性能要求。
确定各临界频率{错误！未找到引用源。
}和滤波器单位脉冲响应长度N。
2.根据性能要求和N值，合理地选择单位脉冲响应h(n)有奇偶对称性，从而确定理想频率响应H d(Ω)的幅频特性和相位特性。
3.用傅里叶反变换公式，求得理想单位脉冲响应h d(n)。
4.选择适当的窗函数W(n)根据式h(n)=h d(n)w(n)求得所设计的FIR滤波器单位脉冲响应。
5.用傅理叶变换求得其频率响应H(Ω)，分析它的幅频特性，若不满足要求，可适当改变窗函数形式或长度N，重复上述过程，直至得到满意的结果。
2. 生成四种窗函数：矩形窗、三角窗、汉宁窗、海明窗，并观察其频率响应。
n=30;%矩形窗及其频响window1=rectwin(n);[h1,w1]=freqz(window1,1);subplot(4,2,1);stem(window1);title('矩形窗');subplot(4,2,2);plot(w1/pi,20*log(abs(h1))/abs(h1(1)));title('矩形窗频响');%三角窗及其频响window2=triang(n);[h2,w2]=freqz(window2,1);subplot(4,2,3);stem(window2);title('三角窗');subplot(4,2,4);plot(w2/pi,20*log(abs(h2))/abs(h2(1)));title('三角窗频响');%汉宁窗及其频响window3=hann(n);[h3,w3]=freqz(window3,1);subplot(4,2,5);stem(window3);title('汉宁窗');subplot(4,2,6);plot(w3/pi,20*log(abs(h3))/abs(h3(1)));title('汉宁窗频响');%汉明窗频响window4=hamming(n);[h4,w4]=freqz(window4,1);subplot(4,2,7);stem(window4);title('汉明窗');subplot(4,2,8);plot(w4/pi,20*log(abs(h4))/abs(h4(1)));title('汉明窗频响');3、（1）用Hanning窗设计一个线性相位带通滤波器，其长度N=15，上下边带宽截止频率分别为ω=0.3π错误！未找到引用源。
用窗口法设计FIR数字滤波器
信号、系统与信号办理实验Ⅱ实验报告实验名称：用窗口法设计FIR数字滤波器一、实验目的认识一个实质滤波器的设计过程，加深掌握用窗口法设计FIR数字滤波器的原理和窗函数对数字滤波器性能的影响二、实验内容与要求1．编写用改良的升余弦窗设计FIR线性相应低通数字滤波器的程序,已知wc=，N=21。
这里wc为低通滤波器的截止频次，N为滤波器的长度,滤波器的阶数为N-1。
2．调试运转程序,要求在屏幕上显示出单位脉冲响应h(n)的数值。
画出其幅度响应|H(ejw)|及20log10|H(ejw)|的曲线。
3.画出窗函数w(n)及其频谱|W(ejw)|和20log10|W(ejw)|的曲线。
三、实验程序与结果clear all;N=21;wc=*pi;b=fir1(N-1,wc/pi,hamming(N))freqz(b,1,512);h(n)的数值即b：b=Columns1through11Columns12through21 clear all;N=21;wc=*pi;r=(N-1)/2;n=0:N-1;hdn=sin(wc*(n-r))/pi./(n-r);if rem(N,2)~=0hdn(r+1)=wc/pi;endwn=hamming(N);h=hdn.*wn';H=fft(h,512);W=fft(wn,512);w=2*[0:511]/512;figure;subplot(3,1,1);plot(wn);subplot(3,1,2);plot(w,abs(W));ylabel('|W(eiw)|') subplot(3,1,3);plot(w,20*log10(abs(W))); ylabel('20lg|W(eiw)|'); figure;subplot(2,1,1);plot(w,abs(H));ylabel('|H(eiw)|') subplot(2,1,2);plot(w,20*log10(abs(H))); ylabel('20lg|H(eiw)|');四、仿真结果剖析图一是FIR低通滤波器的幅频和相频特征图图二是滤波器设计过以后的幅度特征图三是窗函数的幅度特征FIR滤波器的设计过程在于运用海明窗口进行设计，依据定义h(n)=hd(n)*w(n)出发来求出w(n)的值。
用窗口法设计FIR数字滤波器
实验8用窗口法设计FIR 数字滤波器8.1实验目的了解一个实际滤波器设计过程，加深掌握用窗口法设计FIR 数字滤波器的原理和窗函数对数字滤波器性能的影响。
8.2基本原理设所希望得到的滤波器的理想频率响应为H d （e jw ）。
那FI 滤波器的设计就在于寻找一个频率响应为()()∑=--=N n jwn jwH 0e n h e 去逼近H d （e jw ）。
在这种逼近中最直接的一种方法是从单位冲激响应序列h(n)着手，使h(n)逼近理想的单位冲激响应h d (n)。
我们知道h d (n)可以从理想频率响应H d （e jw ）通过傅里叶反变换得到，即： ()()∑∞∞=-=n jwn jwH e n h e dd ()⎰-=ππωπd H jwn jw de )(e 21n h d但是一般来说，这样得到的单位冲激响应h d (n)往往都是无限长序列，而且是非因果的。
我们以一个截止频率为ωC 的线性相位理想低通滤波器为例来说明。
设低边滤波器的时延为ɑ，即()πωωωωα≤≤≤⎩⎨⎧-cc jw jw H 0e ed (8-2) 则 ()()[]()απαωωπωωα--==⎰--n n d c jwn jw sin e e 21n h c c d 这是一个以ɑ为中心的偶对称的无限长非因果序列。
这样一个无限长的序列怎样用一个有限长主序列去近似呢?最简单的办法就是直接截取它的一段来代替它。
例如n=0到n=N —1的一段截取来作为h(n)，但是为了保证所得到的是线性相位滤波器，必须满足h(n)的对称性，所以时延ɑ应该取h(n)长度的一半，即ɑ=(N —1)/2。
()()n100h n h d 其他-≤≤⎩⎨⎧N n n (8-3) 这种直接截取的办法可以形象地想象为，h(n)好比是通过一个窗口所看到的一段h d (n)。
h(n)为h d (n)和一个"窗口函数"的乘积。
在这里，窗口函数就是矩形序列R N (n)，即h(n)=h d (n)·R N (n) (8-4)但是一般来说，窗口函数并不一定是矩形序列，可以在矩形以内对h d (n)作一定的加权处理，因此，一般可以表示为h(n)=h d (n)·ω(n) (8-5)这里ω(n)就是窗口函数。
fir滤波器窗函数设计法
fir滤波器窗函数设计法
FIR (Finite Impulse Response)滤波器的窗函数设计法是一种经典的数字滤波器设计方法。
它通过选择一个合适的窗函数来对滤波器的频率响应进行加权，从而实现对信号的滤波。
窗函数设计法的基本步骤如下：
1. 确定滤波器的设计规格：包括截止频率、通带和阻带的幅频响应要求等。
2. 根据设计规格，计算出滤波器的理想频率响应：可以使用理想滤波器的频率响应作为目标。
3. 选择一个合适的窗函数：常用的窗函数有矩形窗、汉宁窗、汉明窗等。
选择窗函数的关键是考虑到主瓣宽度和副瓣衰减的平衡。
4. 将选择的窗函数应用于理想频率响应上，得到加权后的频率响应。
5. 对加权后的频率响应进行反变换，得到滤波器的时域响应。
6. 根据需要，对时域响应进行截断或零增益处理，以满足设计规格。
7. 最后，根据计算得到的滤波器系数，可以通过巴特沃斯频率抽样公式将其转换为巴特沃斯直接型或传输函数型，以便在数字系统中实现滤波。
需要注意的是，FIR滤波器的窗函数设计法是一种近似方法，设计的滤波器无法完全符合理想要求。
设计过程中需要权衡主瓣宽度和副瓣衰减等因素，以及选择合适的截断或零增益处理方式，以获得满
足实际需求的滤波器性能。
实验六 用窗函数法设计FIR数字滤波器
实验六 用窗函数法设计FIR 数字滤波器一，实验目的1.掌握窗函数法设计FIR 数字滤波器的原理及具体方法；2.深入理解吉布斯现象，理解不同窗函数的特点。
二，实验内容例6.1 利用fir1函数和矩形窗设计一个N=51，截止频率ωc =0.5π的低通滤波器，画出幅频特性。
MATLAB 程序： clearN=51;wc=0.5;h=fir1(50,wc,boxcar(N)) [H,W]=freqz(h,1) plot(W/pi,abs(H));title(‘矩形窗振幅特性/dB ’);xlabel(‘相对频率’);ylabel(‘H(w)’)00.10.20.30.40.50.60.70.80.910.20.40.60.811.21.4矩形窗振幅特性/dB相对频率H (w )例6.2 利用fir1函数和布莱克曼窗设计一个N=51，截止频率为1p ω=0.3π，2p ω=0.4π的带通滤波器。
MATLAB 程序：clearN=51;wc=[0.3,0.4];h=fir1(50,wc,'bandpass',blackman(N)) [H,W]=freqz(h,1) plot(W/pi,abs(H));title('布莱克曼窗带通振幅特性/dB'); xlabel('相对频率');ylabel('H(w)')0.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.70.80.91布莱克曼窗带通振幅特性/dB相对频率H (w )6.3.1窗函数法设计低通数字滤波器，{πωπωπωω4.0||0,||4.0,02/)(≤≤≤<-=N j e j d eH（1）N=26，分别利用矩形窗，汉宁窗和布莱克曼窗设计该滤波器，且滤波器具有线性相位。
绘出脉冲响应()n h 及滤波器的频率响应； （2）增加N ，观察过渡带和最大肩峰值的变化。
窗函数设计fir滤波器的原理
窗函数设计fir滤波器的原理
FIR滤波器是数字信号处理中常用的一种滤波器，其特点是具有线性相位和稳定性。
在FIR滤波器中，窗函数是一种常用的设计方法，它可以用来控制滤波器的频率响应和滤波器的截止频率等参数。
窗函数的设计原理是基于信号的截断和补零，即将信号限制在一个有限的时间段内，并在信号的两端补零，使得信号在有限时间内变得平滑。
这样做的目的是为了避免信号在频域上出现不必要的波动，从而影响滤波器的性能。
在FIR滤波器中，窗函数的设计可以分为两个步骤：首先选择一个合适的窗函数，然后将该窗函数与理想滤波器的频率响应相乘，得到实际滤波器的频率响应。
常用的窗函数包括矩形窗、汉明窗、汉宁窗、布莱克曼窗等。
其中，矩形窗是最简单的窗函数，其频率响应为常数，但其截止频率较高，不适合用于滤波器的设计。
汉明窗和汉宁窗的频率响应较为平滑，但其截止频率较低，适合用于低通滤波器的设计。
布莱克曼窗的频率响应最为平滑，但其计算较为复杂，适合用于高精度的滤波器设计。
在选择窗函数后，需要将其与理想滤波器的频率响应相乘，得到实际
滤波器的频率响应。
理想滤波器的频率响应可以通过傅里叶变换得到，其截止频率和通带宽度可以根据滤波器的设计要求进行调整。
将窗函
数与理想滤波器的频率响应相乘后，得到的实际滤波器的频率响应可
以通过傅里叶反变换得到时域响应，从而得到滤波器的系数。
总之，窗函数是一种常用的FIR滤波器设计方法，其原理是通过信号
的截断和补零来控制滤波器的频率响应和截止频率等参数。
在设计过
程中，需要选择合适的窗函数，并将其与理想滤波器的频率响应相乘，得到实际滤波器的频率响应，从而得到滤波器的系数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
郑州航空工业管理学院《电子信息系统仿真》课程设计14级专业班级题目基于矩形窗设计FIR数字滤波器姓名学号二О一六年十一月二十五日第一章 FIR 滤波器的设计原理及方法FIR 滤波器通常采用窗函数方法来设计。
窗设计的基本思想是，首先选择一个适当的理想选频滤波器（它总是具有一个非因果，无限持续时间脉冲响应），然后街区（加窗）它的脉冲响应得到线性相位和因果FIR 滤波器。
我们用Hd （e^jw ）表示理想的选频滤波器，它在通带上具有单位增益和线性相位，在阻带上具有零响应。
一个带宽wc<pi 的低通滤波器由下式给定：⎩⎨⎧=<-⋅<==<=wc w jaw e w wc jw e Hd ||),(^1||,0)^(π为了从hd(n)得到一个FIR 滤波器，必须同时在两边截取hd （n ）。
而要得到一个因果的线性相位滤波器，它的h(n)长度为N ，必须有：⎩⎨⎧-<=<==其他,010),()(N n n hd n h 2/)1(-=N α这种操作叫做加窗，h(n)可以看做是hd(n)与窗函数w （n ）的乘积：h(n)=hd(n)w(n) 其中⎩⎨⎧-≤≤=，其他对称，关于01n 0)(N n w α根据w （n ）的不同定义，可以得到不同的窗结构。
在频域中，因果FIR 滤波器响应H(e^jw)由Hd(e^jw)和窗响应W(e^jw)的周期卷积得到，即⎰--=⊗=ππλλωπd j e W jw e W jw e Hd jw e H )(^(2/1)^()^()^(常用的窗函数有矩形窗、巴特利特（BARTLETT ）窗、汉宁（HANNING ）窗、海明（HAMMING ）窗、布莱克曼（BLACKMAN ）窗、凯泽（KAISER ）窗等。
FIR 滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。
窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列，主要设计步骤为：(1) 通过傅里叶逆变换获得理想滤波器的单位脉冲响应hd(n)。
(2) 由性能指标确定窗函数W(n)和窗口长度N 。
(3) 求得实际滤波器的单位脉冲响应h(n)， h(n)即为所设计FIR 滤波器系数向量b(n)。
第二章窗函数：。
1矩形窗矩形窗属于时间变量的零次幂窗。
矩形窗使用最多，习惯上不加窗就是使信号通过了矩形窗。
这种窗的优点是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了高频干扰和泄漏，甚至出现负谱现象。
2 三角窗三角窗亦称费杰（Fejer）窗，是幂窗的一次方形式。
与矩形窗比较，主瓣宽约等于矩形窗的两倍，但旁瓣小，而且无负旁瓣。
3 汉宁（Hanning）窗汉宁窗又称升余弦窗，汉宁窗可以看作是3个矩形时间窗的频谱之和，或者说是 3个 sine（t）型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了π/T，从而使旁瓣互相抵消，消去高频干扰和漏能。
可以看出，汉宁窗主瓣加宽并降低，旁瓣则显著减小，从减小泄漏观点出发，汉宁窗优于矩形窗．但汉宁窗主瓣加宽，相当于分析带宽加宽，频率分辨力下降。
4 海明（Hamming）窗海明窗也是余弦窗的一种，又称改进的升余弦窗。
海明窗与汉宁窗都是余弦窗，只是加权系数不同。
海明窗加权的系数能使旁瓣达到更小。
分析表明，海明窗的第一旁瓣衰减为一42dB．海明窗的频谱也是由3个矩形时窗的频谱合成，但其旁瓣衰减速度为20dB／（10oct），这比汉宁窗衰减速度慢。
海明窗与汉宁窗都是很有用的窗函数。
5 高斯窗高斯窗是一种指数窗。
高斯窗谱无负的旁瓣，第一旁瓣衰减达一55dB。
高斯富谱的主瓣较宽，故而频率分辨力低．高斯窗函数常被用来截断一些非周期信号，如指数衰减信号等。
对于窗函数的选择，应考虑被分析信号的性质与处理要求。
如果仅要求精确读出主瓣频率，而不考虑幅值精度，则可选用主瓣宽度比较窄而便于分辨的矩形窗，例如测量物体的自振频率等；如果分析窄带信号，且有较强的干扰噪声，则应选用旁瓣幅度小的窗函数，如汉宁窗、三角窗等；对于随时间按指数衰减的函数，可采用指数窗来提高信噪比。
不同的窗函数对信号频谱的影响是不一样的，这主要是因为不同的窗函数，产生泄漏的大小不一样，频率分辨能力也不一样。
信号的截断产生了能量泄漏，而用FFT算法计算频谱又产生了栅栏效应，从原理上讲这两种误差都是不能消除的，但是我们可以通过选择不同的窗函数对它们的影响进行抑制。
(矩形窗主瓣窄，旁瓣大，频率识别精度最高，幅值识别精度最低；布莱克曼窗主瓣宽，旁瓣小，频率识别精度最低，但幅值识别精度最高)第三章课程程序1. %程序功能：产生一个包含低频、中频、高频的连续信号进行取样并画出其时域图。
Fs=200;t=(1:100)/Fss1=sin(2*pi*t*5)s2=sin(2*pi*t*15)s3=sin(2*pi*t*30);s=s1+s2+s3;plot(t,s); % 画出信号的时域波形xlabel('Time(seconds)')ylabel('Time waveform')title('原始信号的时域波形')2. %程序功能：画出信号的频谱图。
S=fft(s,1024) % 对s进行快速傅立叶变换w=(0:255)/256*(Fs/2)plot(w,abs(S(1:256))) % 画出信号的幅度图xlabel('Frequency (Hz)')ylabel('幅度')title('幅度谱')axis([0 35 0 60]);gridplot(w,angle(S(1:256))) % 画出信号的相位图xlabel('Frequency (Hz)')ylabel('相位')title('相位谱')3.% 程序功能：设计低通滤波器并画出其频谱图：fb=10; fc=13 % 设置滤波器截止频率fs=100; wb=2*pi*fb/fsws=2*pi*fc/fs; wc=0.5*(wb+ws)tr_width=ws-wb % 过渡带宽M=ceil(1.8*pi/tr_width)n=0:M;r=M/2;hn=fir1(M,wc/pi,boxcar(M+1)); % 矩形窗加窗法figure(1);freqz(hn,1,100) %矩形窗FIR的频响曲线figure(2);subplot(1,2,1);stem(n,hn,'k.') %显示矩形窗FIR的系数title('矩形窗FIR')4.% 程序功能：画出滤波后信号的时域波形和频谱sf=filter(hn,[1],s) % sf为滤滤波后的信号plot(t,sf) % 画出滤波后信号的时域波形xlabel('Time(seconds)')ylabel('Time waveform')axis([0 1 -1 1])title('滤波后信号的时域波形')SF=fft(sf,512) % 对sf进行快速傅里叶变换w=(0:255)/256*(Fs/2);plot(w,abs(SF(1:256))) % 画出滤波后信号的幅度图xlabel('Frequency (Hz)')ylabel('幅度谱')title('滤波后信号的幅度谱'); gridplot(w,angle(SF(1:256))) % 画出滤波后信号的相位图xlabel('Frequency (Hz)')ylabel('相位谱')title('滤波后信号的相位谱')4.% 程序功能：对滤波前后信号进行比较plot(w,abs([S(1:256)' SF(1:256)'])) % 将滤波前后信号的幅度谱画在一起xlabel('Frequency (Hz)')ylabel('Mag.of Fourier transform')grid;legend({'before','after'}) %对两个曲线进行区分命名title('滤波前后信号对比')第四章仿真结果分析4.1产生原始信号并分析频谱产生原始信号，用plot函数画出其时域波形如图4.1所示：图4.1原始信号时域波形图用fft函数进行傅里叶变换后，如图4.2所示：图4.24.2将不同特性的数字滤波器用矩形窗表现出来在低通滤波器中，低频，中频，高频是信号含有的三个频率分量。
步骤如下: (1)取通带截止频率为fb=40Hz ，阻带截止频率为fs=46Hz ; (2)取样频率为fs=200Hz ，由wb=2*pi*fb/fs 、ws=2*pi*fc/fs ； (3)过渡带宽：s w —p w =1.8M，得出M ； (4)取一理想低通脉冲响应hd=ideal_lp(wc,M)，其中wc=0.5*(wb+ws)；(5)由于矩形窗用函数w_box=(boxcar(M))'，于是对hd 加矩形窗w_box ，得到一长度为M 的因果且线性相位FIR 滤波器h=hd.*w_box ；（6）用函数[db,mag,pha,grd,w]=freqz_m(h,[1])得出其频谱特性。
（7）画出理想脉冲响应hd 、矩形窗w_boxcar 、实际脉冲响应h 和滤波器的幅度响应(dB)，如图4.3所示。
图4.3滤波器的相关图4.3滤波后的信号处理用函数sf=filter(h,[1],s)分别得出滤波后的信号，其中h对应不同滤波器而不同，通过带通滤波器后的的信号时域波形和频谱图如图4.4(a)、4.4(b)、4.4(c)所示。
图4.4(a)滤波后信号的时域波形图图4.4(b)滤波后信号的幅度谱共享知识分享快乐卑微如蝼蚁、坚强似大象图4.4(c)滤波后信号的相位谱4.4 结果分析函数plot(w,abs([S(1:256)' SF(1:256)']))和legend({'before','after'})后，可在一张图中清晰的看出滤波前后的不同：第5章课程总结通过本次实验我掌握了MATLAB工具软件在工程设计中的使用并掌握窗函数法设计线性相位型FIR滤波器的方法。
同时也知道了不同窗函数对所设计数字滤波器性能的影响，可以清楚的对滤波输出结果作出分析，此次试验也让我认识到了我自己在设计滤波器这方面存在的不足，我会努力学习，在设计实验这方面多努力，这次实验对我今后的学习有很大帮助，也让我找到了今后学习的方向及目标。