高通数字滤波器的设计及DSP实现.doc

高通数字滤波器的设计及DSP实现 班级： 学号： 姓名： 现场答辩（20分） 书面报告（80分） 成绩 Matlab设计（20分） DSP程序（25分） 输入输出分析（20分） 格式规范（15分） 总分 等级 1、设计目的 1.熟悉FIR高通数字滤波器工作原理。 2.熟悉线性相位FIR高通数字滤波器的特性。 3.深入了解利用Matlab设计FIR高通数字滤波器的基本方法。 4.在设计的过程中掌握Matlab和CCS程序编译及软件设计的基本方法。 5.锻炼自己通过网络及各种资料解决实际问题的能力。 6.提高自己对于新知识的学习能力及进行实际操作的能力。 2、设计原理 滤波器是指用来对输入信号进行滤波的硬件和软件。所谓数字滤波器，是指输入、输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器一般可用两种方法实现：一种是根据描述的数字滤波器的数学模型或信号流图，用数字硬件装配成一台专门的设备，构成专用的信号处理机，这就是硬件实现方式；另一种方法就是直接利用通用计算机，将所需要的运算编成程序来让计算机来执行，这也就是软件实现方式。在硬件实现方式中，是以数字组件如延迟器、加法器和乘法器作为基本部件构成专用的数字信号处理系统。在软件实现方式中，它是借助于通用计算机机器语言、汇编语言或高级语言程序来做数字滤波器的运算过程。 滤波器性能最容易通过它的幅频响应的形状来描述?。滤波器在某个频率的幅度增益决定了滤波器对此频率输入的放大因子，增益可取任意值。增益高的频率范围，信号可以通过，称之为滤波器的通带，增益低的频率范围，滤波器对信号有衰减和阻塞作用，称这位滤波器的阻带。按通频带分类，滤波器可分为：低通滤波器、高通滤波器、带通滤波器、带阻滤波器等。本次设计的是高通滤波器。 数字信号处理就是在有限区间使用所观测到的信号序列进行各种各样的处理。截取信号区间两端的位置微小变化，将会导致结果产生很大变化，这种现象一般在信号处理中是必须避免的。为此，不仅应注意截取信号区间的长度，而且要注意截取区间两端不应造成急剧的变化。?截取持续信号中部分信号的工作，可以看作是通过一个窗函数采集所看到的信号序列。这种为截取信号所使用的窗口称为窗函数(window?function)。窗函数不仅适用于离散时间信号，也适用于连续信号。?实际上，设观测到的信号为x(t)，窗函数为w(t)，从窗口观测的信号为y(t)，则?y(t)=w(t)*x(t)。窗口函数的种类很多：矩形窗、汉宁窗、汉明窗、布莱克曼窗等。本次设计采用的是布莱克曼窗。 3、设计要求 1. A/D采样率 f s 为 100KHZ；采用窗函数设计法；数字角频率ω = ΩT = 2π f / fs （3-1），设计一个通带截止频率为 15KHZ 的高通滤波器。 2. 利用 matlab 进行相应滤波器设计，给出窗口选择、设计方法、设计过程及设计结果。 3. 利用上述设计所计算的滤波器系数，编写相应的 DSP 实现程序， 4. 输入信号为 1.5KHz、10KHz、20KHz 的叠加信号，即x(t)=sin(3000πt)+sin(20000πt)+ sin(40000πt)，观察其经过滤波器后的响应，记录输入、输出的时域、频域波形图，并对结果作相应分析。 5. 书写报告，要求整洁、规范，标题用四号宋体、正文用小四号宋体/Times New Roman，图表内文字、程序语句用五号宋体/Times New Roman；对所列图、表均给出标题说明，关键程序语句作注释。 4、主要设计内容  4.1MATLAB设计滤波器参数及实现 利用窗函数法设计FIR高通数字滤波器，绘制出高通数字滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进行处理，对比滤波前后的信号时域和频域图，验证滤波器的效果。 此次设计的FIR高通数字滤波器通带截止频率为fc=15KHZ ，设阻带边缘频率为f1=10KHZ，通带边缘频率为f2=20KHZ，根据公式（3-1）得，通带截止数字角频率为wc=2*pi*fc/fs；阻带边缘数字角频率ws=2*pi*f1/fs，通带边缘数字角频率wp=2*pi*f2/fs，过滤带宽wb=wp-ws，计算布莱克曼窗的长度N=11*pi/wb=55。 利用fir1函数设计线性相位RIRDF：b=fir1(n,Wn,’ftype’,window)n为滤波器的阶数；Wn为滤波器的截止频率，它是一个0到1的数。ftype为滤波器的类型，ftype=’high’时，设计的是高通滤波器；window为指定的窗函数，布莱克曼窗为blackman(n)。本次设计b=fir1(54,wc