DSP课程设计FIR滤波器设计.docx

实用标准文档 DSP 课程设计报告 题目：FIR滤波器设计学院：电气信息学院专业：通信工程 姓名： 学号： 指导老师：曹玉英 文案大全 DSP DSP 课程设计报告 PAGE PAGE 10 一、设计目标 设计一个 FIR 高通滤波器，通带边界频率为 2khz，通带纹波小于 1dB，采样频率为 8khz，实现当多个频率的输入信号输入时只保留大于 2khz 的信号功能， 其中 FIR 滤波器的设计可以用 MATLAB 窗函数法进行。 二、算法研究 FIR 的原理和参数生成公式 FIR 数字滤波器是非递归性的线性时不变因果系统，这样的系统的差分方程可以表示为： 令输入信号 x(n)=δ(n),代入(2.1)式，有 =a δ(n)+a δ(n-1)+··+a (2.1) δ[n-(N-1)] (2.2) 0 1 n-1 这时的 y(n)即为冲激响应 h(n)。由(2.2)式很容易得到 h(0)=a ，h(1)=a , ··,h(N-1)=a ； 0 1 n-1 又由(2.2)式可知，当 n0 以及 nN-1 时，h(n)=0,即这个系统的冲激响应 是有限长度的，这样的滤波器就叫做有限冲激响应(FIR)滤波器。将 a =h(i)(i=0,1, ··，N-1)代入(2.2)式，得到 (2.3)i (2.3) 将(2.3)式的两边进行 z 变换后，可以得到 FIR 滤波器的系统函数： (2.4)? ? (2.4) 又由(2.4)式，有 ? ? ? ? 因此，FIR 滤波器的系统函数 H(z)的极点都位于 z=0 处，为 N-1 阶极点； 而 N-1 个零点由冲激响应 h(n)决定，一般来说，可以位于有限 z 平面的任何位置。 由于 FIR 数字滤波器的极点都集中在单位圆内的原点 z=0 处，与系数 h(n) 无关，因此 FIR 滤波器总是稳定的，这是 FIR 数字系统的一大优点。 利用 MATLAB 计算滤波系数 用来设计标准频率响应的基于窗函数的 FIR 滤波器，可实现加窗线性相位 FIR 数字滤波器的设计。语法：b=fir1(n,Wn) b=fir1(n,Wn,‘ftype’) b=fir1(n,Wn,Window) b=fir1(n,Wn,‘ftype’,window) n:滤波器的阶数；Wn：滤波器的截止频率； ftype：用来决定滤波器的类型， 当 ftype=high 时，可设计高通滤波器； 当 ftype=stop 时，可设计带阻滤波器； Window：用来指定滤波器采用的窗函数类型，其默认值为汉明窗。3.编写产生滤波器输入信号的程序 输入信号应该至少包含两种频率成分的正弦信号，一种信号频率小于2000hz，一种信号频率大于 2000hz。可以再MATLAB 中产生，也可编写DSP 程序产生。 三、开发平台 1.MATLAB MATLAB 是美国 MathWorks 公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB 和 Simulink 两大部分。 MATLAB 是matrixlaboratory 两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如 C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 2.CCS CCS(Code Composer Studio)是 TI 公司推出的针对 TMS320 系列 DSP 的集成开发环境。在CCS 下，开发者可对软件进行编辑，编译，调试，代码性能测试和项目管理等所有工作，并能将程序下载到目标 DSP 上运行调试。 在一个开放式的插件结构下，CCS 内部集成了以下软件工具： C5000 代码产生工具（包括 C5000 的编译器，汇编优化器，汇编器和连接器）； 软件模拟器（Simulator）； 3）实时基础软件 DSP/BIOSTM ； 主机与目标机之间的实时数据交换软件 RTDXTM； 实时分析和数据可视化软件； CCS 不仅具有一系列的调试、分析能力，还提供了实时分析和数据可视化功能，大大降低了 DSP 系统的开发难度，使开发者可以将精力集中在应用开发上。 四、参数计算 利用 MATLAB 计算滤波系数 利用MATLAB 中的FDATOOL 设计滤波器参数设置如图 4-1，频域波形如图 4-2。其中 FDATool（Filter Design an