数字信号处理三级项目.doc

燕 山 大 学 多频窄带数字信号处理仿真系统设计 学生所在学院：信息学院 学生所在班级：电子信息工程2班 课程名称：数字信号处理 学生姓名： 姜升殿、李林、刘鹏、 左园叶、陈敏、黄品程 指导教师：陈书贞 日期：2014.12.20 2014年12月 多频窄带数字信号处理仿真系统设计 姜升殿、李林、黄品程、刘鹏、左园叶、陈敏 摘 要 本文主要是实现在TC2.0的环境下利用C语言实现矩形窗、哈明窗、汉宁窗、布莱克曼窗。最后本文对设计的系统程序进一步改进，使之可以完成信号模拟到数字的转换功能、根据数字滤波器指标自动完成滤波器设计的功能、四种滤波功能以及个步骤的时域频域图形绘制。 前言 本次数字信号三级项目是教改后用于帮助同学理解数字信号的理论和提高处理解决问题的能力而设立的，主要是利用数字信号的理论结合编程实现一个 项目组分工： 组长姜升殿，负责程序的编写和调试以及报告和汇报ppt的部分制作； 组员李林，负责报告主要编写； 组员陈敏，负责汇报ppt的主要制作； 组员黄品程，负责部分程序 组员左园叶，负责报告的理论部分整理； 组员刘鹏，负责报告的理论部分整理。 研究报告正文 1系统整体设计思想 本报告主要是用C语言在tc2.0环境下实现 由此可知，该报告设计需包括以下系统 A 模拟信号采集模块：完成对模拟信号的采样、量化，使信号变成数字信号。 B 显示模块：完成对信号的时域和频域显示，以便于对信号的频谱和时域分布进行分析。 C 滤波器模块：根据模拟信号生成的数字信号的频域特点确定要完成的滤波指标，进而按照一定的原则选择合适的滤波器进行最后要完成的滤波过程。 D 滤波模块：系统的最后一步， 具体的系统整体设计流程见图1系统设计框图 图1系统设计框图 各模块设计 2.1 模拟信号采集模块 本项目的 为例，采样后的数字信号为,其中。故编程实现即为： { int i; for ( i = 0; i n; i++) { ARRAYI[i] = 0.0; ARRAYR[i] = (float)cos(2 * PI*i*f1 / fs) + (float)cos(2 * PI*i*f2 / fs); } } 程序中 2.2 显示模块 为了实现信号的采样滤波过程的2.0编译环境实现对信号的时域和频域显示。考虑到信号包含实部和虚部，因此程序编写时为了方便，选择将一个信号用实部和虚部表示。因此显示的时候就出现了两种情况，时域显示和频域显示，信号不同域波形需要不同的求模函数实现对信号的强度的表示，具体程序函数声明为： void GetModu1(float a[],float b[],int ll) /*取模函数，主要用来对频域信号的取模*/ void GetModu2(float a[], float b[],int NNN) 在具体的DOS环境下的图形显示，考虑系统功能实现问题，设计两个绘图函数，用于专门显示窗函数 /*画图程序，用以画窗函数的时域图*/ void GRAPH(float xw[MAX],int nn,int kk)/*xw表示窗函数，nn表示描点点数，kk表示滤波器类型*/ /*画图程序，用以画函数的时域和频域的图型*/ void GRAPH1(float a[],int n,int ii,int high)/*n表示描点点数，ii表示步长，high表示所控制的画图函数的高度*/ 2.3 滤波器模块 本项目选择使用FIR滤波器设计方法。表示理想的选频滤波器，它在通带上具有单位增益和线性相位，在阻带上具有零响应。一个带宽的低通滤波器由下式给定： 为了从得到一个FIR滤波器，必须同时在两边截取。而要得到一个因果的线性相位滤波器，它的长度为，必须有： 这种操作叫做加窗，可以看作是与窗函数的乘积： 其中 根据的不同定义，可以得到不同的窗结构。 在频域中，因果FIR滤波器响应由和窗响应的周期卷积得到，即： ?? 常用的窗函数有矩形窗、巴特利特（BARTLETT）窗、汉宁（HANNING）窗、海明（HAMMING）窗、布莱克曼（BLACKMAN）窗、凯泽（KAISER）窗等。 窗函数设计的具体步骤见流程图2-1： 图2-1 窗函数设计流程图 图2-2 滤波器指标选择窗和理想响应流程 本项目实现提供四个窗函数选择，分别是矩形窗、哈明窗、汉宁窗、布莱克曼窗，通过给定的数字滤波器指标可以求出wn(n)窗函数，wn(n)窗函数得到程序实现上见附录程序中JXC(),HNC(),HMC(),BLKMC()四个窗函数子函数。具