下载文档原格式
基于MATLAB的数字带通滤波器课程设计报告1西安文理学院机械电子工程系课程设计报告专业班级08级电子信息工程1班题目基于MATLAB的数字带通滤波器学号学生姓名指导教师2011 年12 月西安文理学院机械电子工程系课程设计任务书学生姓名_______专业班级________ 学号______指导教师______ 职称副教授教研室电子信息工程课程数字信号处理题目基于MATLAB 的数字带通滤波器设计任务与要求设计任务：要求设计一个IIR 带通滤波器，其中通带的中心频率为πω5.0=po ，通带的截止频率πω4.01=p ，πω6.02=p ，通带最大衰减dB p 3=α；阻带最小衰减dB s 15=α，阻带截止频率πω3.01=s ，πω7.02=s 。
设计要求：1. 根据设计任务要求给出实现方案及实现过程。
2. 给出所实现的滤波器幅频特性及相频特性曲线并加以分析。
3. 论文要求思路清晰，结构合理，语言流畅，书写格式符合要求。
开始日期2011.12.19 完成日期2011.12.302011年12月18 日一、设计任务设计一数字带通滤波器，用IIR 来实现，其主要技术指标：通带边缘频率：wp 1=0.4π，wp2=0.6π通带最大衰减：Ap=3dB阻带边缘频率：ws 1=0.3π，ws2=0.7π阻带最小衰减：As=15dB设计总体要求：用MATLAB 语言编程进行设计，给出IIR 数字滤波器的参数，给出幅度和相位响应曲线，对IIR 实现形式和特点等方面进行讨论。
二、设计方法IIR 数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR 滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。
比较常用的原型滤波器有巴特沃什滤波器（Butterworth ）、切比雪夫滤波器（Chebyshev ）、椭圆滤波器（Ellipse ）和贝塞尔滤波器（Bessel ）等。
他们有各自的特点，巴特沃什滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带和阻带里有波动，可以提高选择性；贝塞尔滤波器通带内有较好的线性相位特性；椭圆滤波器的选择性最好。
基于MATLAB信号处理工具箱的数字滤波器设计摘要：传统数字滤波器工作起来，需要计算量很大，对于设计要求来说，又比较复杂，需要调整滤波特性，有很大的难度，并且影响它本身的应用以及功能。
本文中，主要介绍了利用MATLAB信号处理工具箱的功能，设计常规数字滤波器，设计过程中使用一些常用的专用软件。
对于，使用MATLAB语言程序设计以及信号处理工具箱的FDATool工具，各种界面的详细步骤都有详细的介绍。
对于MATLAB设计滤波器，为了减轻庞大的工作量，也为了简便直观，我们对比设计要求，以及滤波器特性的参数，这样才能更有利于，优化滤波器本身。
在本文中，主要解读了仿真软件Simulink对设计出来的滤波器仿真，设计过程中利用工具MATLAB。
我们验证结果的准确性，是通过软件仿真实验，从而达到对滤波器的设计。
关键词：MATLAB；模拟滤波器；数字滤波器目录1 引言 (1)2 模拟滤波器的设计 (1)2.1 几种常用的模拟滤波器 (2)2.2 巴特沃斯低通滤波器设计 (3)3 数字滤波器的设计 (5)4 结束语 (10)参考文献 (10)1引言基于对正弦基波，以及整词谐波的分析，可以得知系统电压，电流信号中，成分比较复杂，因此，在电力系统中，最关键的部件还是滤波器。
MATLAB和Mathematica、Maple共同称为三大数学软件。
MATLAB在数学类科技应用软件中，在数值计算方面地位不可代替。
MATLAB可以进行矩阵运算，绘制函数和数据，实现算法，创建用户界面，连接其他编程语言的程序等功能，它的应用方面主要在于，工程计算；控制设计；信号处理与通讯；图像处理信号检测；金融建模设计与分析等领域。
MATLAB的基本数据单位是矩阵，它的指令表达式与工程中常用的形式十分相似，故用MATLAB来解算问题要比其他语言完成事情简捷，因此可以直接使用,我们可以把自己的程序导入进MATLAB函数库中，方便以后的使用。
因此我们设计滤波时，根据指标，直接使用CAD程序，就可以得到高通、带通和带阻滤波器系统函数。
目录1 引言 (1)2 设计任务 (3)2.1设计内容 (3)2.2设计要求 (3)3 语音信号的采集及时频分析 (4)3.1语音信号的采集 (4)3.2语音信号的时频分析 (4)4 基于MATLAB的数字滤波器的设计 (5)4.1数字滤波器的设计 (5)4.1.1数字滤波器的基本概念 (5)4.1.2 IIR滤波器设计思想 (5)4.2IIR数字滤波器设计 (5)4.2.1 IIR低通滤波器设计 (5)4.2.2 IIR带通滤波器设计 (7)4.2.3 IIR高通滤波器设计 (7)5 合成信号及其滤波 (9)5.1合成信号 (9)5.2合成信号滤波 (9)6 设计系统界面 (11)6.1系统界面设计工具—GUI概述 (11)6.2界面设计及使用说明 (11)7 心得体会 (13)参考文献 (14)1 引言数字滤波在通信、图像编码、语音编码、雷达等许多领域中有着十分广泛的应用。
目前，数字信号滤波器的设计在图像处理、数据压缩等方面的应用取得了令人瞩目的进展和成就。
近年来迅速发展起来。
MATLAB 是美国 MathWorks 公司推出的一套用于工程计算的可视化高性能语言与软件环境。
MATLAB 为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。
它以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。
MATLAB 推出的工具箱使各个领域的研究人员可以直观方便地进行科学研究、工程应用，其中的信号处理(signalproeessing)、图像处理(imageproeessing)、小波(wavelet)等工具箱为数字滤波研究的蓬勃发展提供了有力的工具。
IIR滤波器采用递归型结构，即结构上带有反馈环路。
IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。
同时，IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。
实验6 基于MATLAB 的数字滤波器设计实验目的：加深对数字滤波器的常用指标和设计过程的理解。
实验原理：低通滤波器的常用指标：PP j P for e G ωωδδω≤+≤≤-,1)(1πωωδω≤≤≤S S j for eG ,)(通带边缘频率：Pω，阻带边缘频率：Sω ，通带起伏：P δ，通带峰值起伏：])[1(log2010dB p pδα--=，阻带起伏：sδ，最小阻带衰减：])[(log 2010dB s S δα-=。
数字滤波器有IIR 和FIR 两种类型，它们的特点和设计方法不同。
在MATLAB 中，可以用[b ，a]=butter （N,Wn ）等函数辅助设计IIR 数字滤波器, 也可以用b=fir1(N,Wn,’ftype’) 等函数辅助设计FIR 数字滤波器。
实验内容： 利用MATLAB 编程设计一个数字带通滤波器，指标要求如下：通带边缘频率：10.4P ωπ=，20.6P ωπ=，通带峰值起伏：][1dB p≤α。
阻带边缘频率：πω3.01=S ，20.7S ωπ=，最小阻带衰减：][40dB S ≥α。
分别用IIR 和FIR 两种数字滤波器类型进行设计。
实验要求：给出IIR 数字滤波器参数和FIR 数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。
实验步骤：①Butterworth 滤波器的设计（IIR ）>>wp=[0.4*pi,0.6*pi]; >>wr=[0.3*pi,0.7*pi]; >>Ap=1; >>Ar=40;Sω-Pω-P SωPassband StopbandTransition bandFig 7.1 Typical magnitudespecification for a digital LPF>>[N,Wn]=buttord(wp/pi,wr/pi,Ap,Ar)N =7Wn =0.3854 0.6146>> [b,a]=butter(N,Wn,'bandpass')b =Columns 1 through 120.0002 0 -0.0014 0 0.0042 0 -0.0071 0 0.0071 0 -0.0042 0Columns 13 through 150.0014 0 -0.0002a =Columns 1 through 121.0000 0.0000 3.7738 0.0000 6.5614 0.0000 6.6518 0.0000 4.2030 0.0000 1.6437 0.0000Columns 13 through 150.3666 0.0000 0.0359>> [H,w]=freqz(b,a);>>mag=abs(H);>>plot(w/pi,mag);>>xlabel('角频率(\Omega)');>>ylabel('幅度|Ha(j\Omega)|');>> title('数字butterworth带通滤波器幅度响应|Ha(j\Omega)|');>>phase=angle(H);>>plot(w/pi,phase);>>xlabel('角频率(\Omega)');>>ylabel('相位');>> title('数字butterworth带通滤波器相位响应曲线');>>②FIR滤波器的设计：>> wp1 = 0.4*pi; wp2 = 0.6*pi;>> ws1 = 0.3*pi; ws2 = 0.7*pi;>>tr_width = min((wp1-ws1),(ws2-wp2)) tr_width =0.3142>> M = ceil(6.2*pi/tr_width) + 1M =63>> n=[0:1:M-1];>> wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2; >>wc=[wc1/pi,wc2/pi];>>window= hanning(M);>> [h1,w]=freqz(window,1);>>figure(1);>>subplot(2,1,1)>>stem(window);>>axis([0 60 0 1.2]);>>grid;>>xlabel('n');>> title('Hanning窗函数');>>subplot(2,1,2)>>plot(w/pi,20*log(abs(h1)/abs(h1(1)))); >>axis([0 1 -350 0]);>>grid;>>xlabel('w/\pi');>>ylabel('幅度(dB)');>> title('Hanning窗函数的频谱');>>hn = fir1(M-1,wc, hanning (M));>> [h2,w]=freqz(hn,1,512);>>figure(2);>>subplot(2,1,1)>>stem(n,hn);>>axis([0 60 -0.25 0.25]);>>grid;>>xlabel('n');>>ylabel('h(n)');>> title('Hanning窗函数的单位脉冲响应'); >>subplot(2,1,2)>>plot(w/pi,20*log(abs(h2)/abs(h2(1)))); >>grid;>>xlabel('w/\pi');>>ylabel('幅度(dB)');>>figure(3);>>phase=angle(h1);>>plot(phase);>>axis([1 pi -1 0]);>>xlabel('w/\pi');>>ylabel('线性相位');>> title('Hanning窗函数相位特性曲线'); >>实现形式及特点分析：1.在本例中，相同的技术指标下，IIR滤波器实现的阶数为N=7，而FIR滤波器的阶数N=63。
第一章绪论1.1引言随着信息时代与数字技术的发展，数字信号处理己逐渐发展成为当今极其重要的学科与技术领域之一。
数字信号处理在通信、语音、图像、自动控制雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。
在数字信号处理的基本方法中，通常会涉与到变换、滤波、频谱分析、调制解调和编码解码等处理。
其中，滤波是应用非常广泛的一个环节，数字滤波器的相关理论也一直都是人们研究的重点之一。
数字滤波器是数字信号处理的重要基础，在对信号的滤波、检测与参数的估计等信号应用中，数字滤波器是使用最为广泛的一种线性系统。
数字滤波器根据其单位冲击响应函数的时域特性可分为两类：无限冲击响应(IIR)数字滤波器和有限冲击响应(FIR)数字滤波器。
与IIR数字滤波器相比,FIR数字滤波器的实现是非递归的，稳定性好，精度高；更重要的是FIR数字滤波器在满足幅度响应要求的同时，可以获得严格的线性相位。
因此，它在高保真的信号处理中，如数字音频、图像处理、数据传输和生物医学等领域得到广泛应用。
1.2数字滤波器的研究背景与意义滤波在通信、图像编码、语音编码、雷达等许多领域中有着十分广泛的应用。
目前，数字信号滤波器的设计在图像处理、数据压缩等方面的应用取得了令人瞩目的进展和成就。
它是数字信号处理理论的一部分。
数字信号处理主要是研究用数字或符号的序列来表示信号波形，并用数字的方式去处理这些序列，以便估计信号的特征参量，或削弱信号中的多余分量和增强信号中的有用分量。
具体来说，凡是用数字方式对信号进行滤波、变换、调制、解调、均衡、增强、压缩、固定、识别、产生等加工处理，都可纳入数字信号处理领域。
数字信号处理学科的一项重大进展是关于数字滤波器设计方法的研究。
关于数字滤波器，50年代已有人讨论过数字滤波器，但直到60年代中期，才开始形成关于数字滤波器的一整套完整的正规理论。
在这一时期，提出了各种各样的数字滤波器结构，有的以运算误差最小为特点，有的则以运算速度高见长，而有的则二者兼而有之。
基于matlab的fir数字滤波器仿真与设计(开题报告)(共9篇)：开题仿真报告设计数字滤波器matlab fir滤波器函数fir1 fir滤波器阶数计算篇一：基于MATLAB的FIR数字滤波器的设计与仿真(开题报告) 毕业（设计）论文开题报告系别电子信息工程系专业电子信息工程班级学生姓名学号指导教师XXX 报告日期毕业（设计）论文开题报告表篇二：基于MATLAB的FIR数字滤波器仿真与设计(开题报告)1、课题来源2、研究目的和意义3、国内外研究现状和发展趋势及综述4、本课题的主要研究内容及拟采取的技术路线、试验方案5、研究基础6、预期达到的目标及进度安排7、参考文献篇三：基于MATLAB的FIR数字滤波器仿真与设计(开题报告)1、课题来源2、研究目的和意义3、国内外研究现状和发展趋势及综述4、本课题的主要研究内容及拟采取的技术路线、试验方案5、研究基础6、预期达到的目标及进度安排7、参考文献篇四：基于matlab的FIR滤波器设计开题报告一、研究背景数字滤波器和模拟滤波器有着相同的滤波概念，根据其频率响应特性可分为低通、高通、带通、带阻等类型，与模拟滤波器相比，数字滤波器除了具有数字信号处理的固有优点外，还有滤波精度高、稳定性好、灵活性强等优点。
数字滤波器可以过滤时间离散信号，通过对抽样数据进行数字处理来达到频域滤波的目的。
由于计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现。
FIR滤波器具有严格的线性相位，可以具有任意的幅度特性，且是可物理实现的因果系统，因此被广泛地应用在现代通信技术当中，如解调器中的位同步与位定时提取、自适应均衡去码间串扰以及话音的自适应编码等。
可见对基于MATLAB的FIR滤波器的研究是具有非常重要的现实意义的。
在高度信息化的今天，随着计算机技术，基于视觉的图像处理技术，模式识别等领域的不断研究与发展，图像识别技术得到更快的发展，并在检测的实时性和准确性方面有更新的突破，所以将数字滤波器应用在工业缺陷零件的检测上也成为当今一热点课题。
毕业论文(设计)学院：计算机科学学院专业: 自动化年级：2007题目：基于MA TLAB的数字滤波器的设计学生姓名：伍杰杰学号：07064098指导教师：程立职称：讲师2010年5 月11日中南民族大学本科论文（设计）原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。
除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。
本人完全意识到本声明的法律后果由本人承担。
作者签名：年月日目录摘要 (4)关键词 (4)Abstract (4)Keywords (4)1 数字滤波器的概述 (4)1.1 数字滤波器的基本结构 (5)1.1.1 IIR 滤波器的基本结构 (5)1.1.2 FIR 滤波器的基本结构 (6)1.2 数字滤波器的设计原理 (6)1.2.1 滤波器的性能指标 (7)1.2.2 IIR 数字滤波器的设计方法 (7)1.2.3 FIR 数字滤波器的设计方法 (8)2.MA TLAB中数字滤波器的设计 (10)2.1 IIR数字滤波器设计 (11)2.1.1. IIR滤波器阶数的选择 (11)2.1.2. IIR滤波器的设计 (11)2.2 FIR数字滤波器设计 (13)2.2.1常用窗函数 (13)2.2.2窗函数法设计FIR数字滤波器 (14)3. 基于MA TLAB GUI的数字滤波器设计思路及实现 (15)3.1 GUI界面设计概述 (15)3.2 “滤波器设计软件”设计所实现任务 (16)3.3基于MATLAB GUI的数字滤波器设计实现 (17)3.3.1 “滤波器设计软件”GUI界面设计 (17)3.3.1 “滤波器设计软件”回调函数编写 (17)3.3.2 AutoChoose.m程序的编写 (21)3.4 运行结果及显示 (25)5. 设计总结和心得 (26)参考文献 (26)基于MATLAB的数字滤波器的设计摘要：数字滤波器是数字信号处理的基础，用来对信号进行过滤、检测与参数估计等处理，在通信、图像、语音等许多领域有着十分广泛的应用，尤其在图像处理、数据压缩等方面取得了令人瞩目的进展与成就。
基于MATLAB的数字滤波器的设计1 绪论1.1 本课题的研究目的及意义数字滤波是数字信号处理的重要基础，在对信号的滤波、检测及参数的估计等信号应用中，数字滤波是使用最为广泛的一种技术。
数字滤波是语音和图象处理、模式识别、谱分析等应用中的一个基本的处理技术，它是数字信号处理中最重要的组成部分之一，几乎出现在所有的数字信号处理系统中[1]。
数字滤波的实现主要是依靠数字滤波器，数字滤波器是由数字乘法器、加法器和延时单元组成的一种装置，其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。
数字滤波器在语音信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用，例如电视、VCD、音响、医学图像处理等。
由于电子计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。
数字滤波器没有漂移，能够处理低频信号，频率响应特性可做成非常接近于理想的特性，且具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点，它可以满足滤波器对幅度和相位特性的严格要求，可以避免模拟滤波器所无法克服的电压漂移、温度漂移和噪声等问题[2]。
而且数字滤波器的种类繁多，有低通、高通、带通、带阻和全通等类型，它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的，便于选择和使用，这些优势决定了数字滤波器的应用越来越广泛，因此滤波器的研究具有极为重要的意义。
1.2 国内外发展现状关于数字滤波器，早在上世纪40年代末期就有人讨论过它的可能性问题，在50年代也有人讨论过数字滤波器的实现问题。
但直到60年代中期，才开始形成关于数字滤波器的一套完整的正规理论。
这一时期，提出了各种各样的数字滤波器结构，有的以运算误差最小为特点，有的则以运算速度快见长。
同时也出现了数字滤波器的各种逼近方法和实现方法，对递归和非递归两类滤波器作了全面的比，统一了数字滤波器的基本概念和理论，开始形成有限冲激响应（IIR）滤波器和无限冲激响应（FIR）滤波器的初步理论。
实验报告姓名：李鹏博 实验名称： IIR 数字滤波器设计 学号：2011300704 课程名称： 数字信号处理 班级：03041102 实验室名称： 航海西楼303 组号： 1 实验日期： 2014.06.20一、实验目的、要求掌握IIR 数字滤波器设计的冲激响应不变法和双线性变换法。
掌握IIR 数字滤波器的计算机编程实现方法，即软件实现。
二、实验原理为了从模拟滤波器设计IIR 数字滤波器，必须先设计一个满足技术指标的模拟滤波器，然后将其数字化，即从s 平面映射到z 平面，得到所需的数字滤波器。
虽然IIR 数字滤波器的设计本质上并不取决于连续时间滤波器的设计，但是因为在许多应用中，数字滤波器就是用来模仿模拟滤波器功能的，所以由模拟滤波器转化为数字滤波器是很自然的。
因此，由模拟滤波器设计数字滤波器的方法准确、简便，是目前最普遍采用的方法。
三、实验环境PC 机，Windows XP ，office 2003，Matlab 软件。
四、实验过程、数据记录、分析及结论实验过程1．编程设计滤波器，用冲激响应不变法设计IIR 数字滤波器。
2．编程设计滤波器，用双线性变换法设计IIR 数字滤波器。
3．求脉冲响应、频率响应以及零极点。
4．编程滤波，求滤波器输出，完成对不同频率的多个正弦信号的滤波。
实验步骤根据所给定的技术指标进行指标转换。
112c c f πΩ=，222c c f πΩ=，112s s f πΩ=，222s s f πΩ=，21p c c B Ω==Ω-Ω，221222s s s s s B Ω-ΩΩΩ=Ω，3,18p s αα=-=-。
根据指标设计Butterworth 模拟低通滤波器。
调用函数[n,wn]=buttord(wp,ws,rp,rs,’s ’)确定阶次。
调用函数[zl,pl,kl]=buttap(n)，求低通原型的模型。
调用函数[bl,al]=zp2tf(zl,pl,kl)实现模型转换。
实验7\8基于MATLAB勺数字滤波器设计实验目的：加深对数字滤波器的常用指标和设计过程的理解。
实验原理：低通滤波器的常用指标：
1 一6P 兰G(e^) ≤ 1 + 6P , for 国≤ ωP
G(J") ≤ 6s, for 国s ≤ ⑷≤ ∏
通带边缘频率：'P ,阻带边缘频率：'s，
通带起伏：J P，通带峰值起伏：
C(P= —20 IOg io (^-OP )【d B 】阻带起伏.冠S
PaSSband StOPband
Tran Siti on
band Fig 7.1 TyPiCaI magn itude SPeCifiCati On for a digital LPF
:S = -20 log ιo(r)[dB 】
O
数字滤波器有IIR和FlR两种类型，它们的特点和设计方法不同。
在MATLAB^,可以用［b , a］=butter ( N,Wr)等函数辅助设计IIR数字滤波器，也可以用b=fir1(N,Wn, 'type ')等函数辅助设计FIR数字滤波器。
实验内容：利用MATLAB编程设计一个数字带通滤波器，指标要求如下：
通带边缘频率：∙∙P1=0.45^,∙∙ P2=0∙65 二，通带峰值起伏：［dB】O
阻带边缘频率：'s1
0.3…,'s2 0.75…，最小阻带衰减：-S 4°[dB] O
分别用IIR和FlR两种数字滤波器类型进行设计。
实验要求：给出IIR数字滤波器参数和FIR数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。
实验内容：
IRR代码：
wp=[0.45*pi,0.65*pi];
ws=[0.3*pi,0.75*pi];
Ap=1;
A S=40;
[N,Wc]=buttord(wp∕pi,ws∕pi,Ap,As);
[b,a]=butter(N,Wc)%[b,a] = butter( n, Wn,'ftype')
最小阻带衰减:
I -W llrreqz(b-a=
SUbP-Of
(211
=
FnagHabS(H)-
P-
Of(WHLmag)-
x ωb e 一
(-
>5
(W))
y ωb e
一
(-
≡M 一 HamW)
口
≡e (--R 一 HamW)
一) grid On-
SUbPOf(212=
PhaSeHang-e(H=
P -
OfWPLPhaSe=
x ωb e 一
(-
>5
(W))
y ωb e -(
- W⅛) ≡e (-一R
盘
a ≡m s
B )
≡-x ⅛x -
幅度IHaOMI
≡
相拉
S
S 5 CU
S S S 0⅛
口 g >酉
tt 2 =R ⅛⅛⅛
J ≡B ⅛
O S F NJ £Q 4
S S U7 CJ
CD S >⅛⅛≡
F-R
V V
WPl H
O∙4*pr
WP2
H
0∙6*pr
V V WSl H 0∙3*pr
ws2 H
0∙7*pr
V
V =Γl w i d f h H
min((wprwsu(ws2lwp2))
>> tr_width =0.3142
>> M = ceil(6.2*pi/tr_width) + 1
>> M = 63
>> n=[0:1:M-1];
>> wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2; >> wc=[wc1/pi,wc2/pi];
>> window= hanning(M);
>> [h1,w]=freqz(window,1);
>> figure(1);
>> subplot(2,1,1)
>> stem(window);
>> axis([0 60 0 1.2]);
>> grid;
>> xlabel('n');
>> title('Hanning 窗函数');
>> subplot(2,1,2)
>> plot(w/pi,20*log(abs(h1)/abs(h1(1)))); >> axis([0 1 -350 0]);
>> grid;
>> xlabel('w/\pi');
>> ylabel('幅度(dB)');
>> title('Hanning 窗函数的频谱');
>> hn = fir1(M-1,wc, hanning (M));
>> [h2,w]=freqz(hn,1,512);
>> figure(2);
>> subplot(2,1,1)
>> stem(n,hn);
>> axis([0 60 -0.25 0.25]);
>> grid;
>> xlabel('n');
>> ylabel('h(n)');
>> title('Hanning 窗函数的单位脉冲响应'); >> subplot(2,1,2)
>> plot(w/pi,20*log(abs(h2)/abs(h2(1)))); >> grid;
>> xlabel('w/\pi');
>> ylabel(' 幅度(dB)');
>> figure(3);
>> phase=angle(h1);
>> plot(phase);
>> axis([1 pi -1 0]);
>> xlabel('w/\pi');
>> ylabel(' 线性相位');
>> title('Hanning 窗函数相位特性曲线')
回
FiIe Edit VieW InSert TOOlS DeSktOP WindOW HeIP
◎ A 聾紳⑥毘謠▼层□ Ξ∣ ■ 0
M FigUre 1 Hanning 窗函数
O
C S I I * i
Hanning®函数的频谱
-100
-200 -300 0
0.1 0.2 0.3 0.4
0.5 0.6 0.7 08 0.9 1
w/x
(8P)
置
Q FigUre 2 口 回 耳
i
Hanni叩窗函数相位特性曲线O I I I I I I I。
