数字信号处理三级项目.

《数字信号处理三级项目.》由会员分享，可在线阅读，更多相关《数字信号处理三级项目.（26页珍藏版）》请在文档大全上搜索。

1、燕山大学信息科学与工程学院燕 山 大 学多频窄带数字信号处理仿真系统设计 学生所在学院：信息学院 学生所在班级：电子信息工程2班 课程名称：数字信号处理 学生姓名： 姜升殿、李林、刘鹏、左园叶、陈敏、黄品程指导教师：陈书贞 日期：2014.12.202014年12月多频窄带数字信号处理仿真系统设计姜升殿、李林、黄品程、刘鹏、左园叶、陈敏摘 要本文主要是实现在TC2.0的环境下利用C语言实现多频窄带数字信号处理仿真系统。系统功能包括信号采集、信号时域及频域分析、FIR滤波器设计、数字信号滤波，在本文中主要设计四种滤波系统：低通滤波系统、高通滤波系统、带通滤波系统、带阻滤波系统。在FIR滤波器设计

2、时使用的是窗函数设计的方法，涉及四种窗函数，分别是矩形窗、哈明窗、汉宁窗、布莱克曼窗。最后本文对设计的系统程序进一步改进，使之可以完成信号模拟到数字的转换功能、根据数字滤波器指标自动完成滤波器设计的功能、四种滤波功能以及个步骤的时域频域图形绘制。关键词：多频窄带数字信号处理仿真系统、FIR滤波器、C语言、TC2.0前言: 本次数字信号三级项目是教改后用于帮助同学理解数字信号的理论和提高处理解决问题的能力而设立的，主要是利用数字信号的理论结合编程实现一个多频窄带数字信号处理仿真系统设计。目前的软件解决数字信号处理的有很多，比较好的是MATLAB，其强大的数字图像处理能力为我们提供了很好的处理环境

3、，但是由于其处理都是集成的，同时给我们的使用带来一定的不便性，对于我们理解信号处理过程和定制我们需要的信号处理软件没有参考价值。故我们选择利用tc2.0环境下的C语言进行系统设计与编写，预期实现信号的采集、信号时域及频域分析、FIR滤波器设计、数字信号滤波等功能的多频窄带数字信号处理软件仿真系统。项目组分工：组长姜升殿，负责程序的编写和调试以及报告和汇报ppt的部分制作；组员李林，负责报告主要编写；组员陈敏，负责汇报ppt的主要制作；组员黄品程，负责部分程序的编写；组员左园叶，负责报告的理论部分整理；组员刘鹏，负责报告的理论部分整理。研究报告正文1系统整体设计思想本报告主要是用C语言在tc2.

4、0环境下实现信号采集、信号时域及频域分析、FIR滤波器设计、数字信号滤波等功能的多频窄带数字信号处理软件仿真系统。由此可知，该报告设计需包括以下系统:A 模拟信号采集模块：完成对模拟信号的采样、量化，使信号变成数字信号。B 显示模块：完成对信号的时域和频域显示，以便于对信号的频谱和时域分布进行分析。C 滤波器模块：根据模拟信号生成的数字信号的频域特点确定要完成的滤波指标，进而按照一定的原则选择合适的滤波器进行最后要完成的滤波过程。D 滤波模块：系统的最后一步，完成对要求滤除的信号频率的滤除和有用频率分量的保留。具体的系统整体设计流程见图1-1系统设计框图:图1-1系统设计框图2 各模块设计2.

5、1 模拟信号采集模块本项目的研究的原始信号是模拟信号，而计算机只能处理二进制编码形式的数字信号，故模拟信号与数字信号之间需有一个误差比较小或者无误差的转换过程，而这个过程就是通过我们的模拟信号采集模块实现的。由理论知识可知，这个转变过程至少需要满足两个条件，一是数字信号可以无失真的还原原模拟信号，二是数字信号的频域相对于模拟信号的频域发生的改变。无失真的还原原模拟信号就要满足采样频率fs大于等于模拟信号的最高频率的二倍。以为例，采样后的数字信号为,其中。故编程实现即为：void Signal_low_high(int n)/*高通和低通信号滤波测试信号函数*/ int i;for ( i =

6、0; i <n; i+)ARRAYIi = 0.0;ARRAYRi = (float)cos(2 * PI*i*f1 / fs) + (float)cos(2 * PI*i*f2 / fs); 程序中ARRAYR,ARRAYI数组表示采集数字信号的实部和虚部。2.2 显示模块为了实现信号的采样滤波过程的图形显示，本次项目选择用TC2.0编译环境实现对信号的时域和频域显示。考虑到信号包含实部和虚部，因此程序编写时为了方便，选择将一个信号用实部和虚部表示。因此显示的时候就出现了两种情况，时域显示和频域显示，信号不同域波形需要不同的求模函数实现对信号的强度的表示，具体程序函数声明为：/*取模函

7、数，主要用来对时域信号的取模*/void GetModu1(float a,float b,int ll)/*取模函数，主要用来对频域信号的取模*/void GetModu2(float a, float b,int NNN)在具体的DOS环境下的图形显示，考虑系统功能实现问题，设计两个绘图函数，GRAPH用于专门显示窗函数的衰减曲线，GRAPH1通用绘图函数为显示步骤中函数的时域和频域，具体见程序：/*画图程序，用以画窗函数的时域图*/void GRAPH(float xwMAX,int nn,int kk)/*xw表示窗函数，nn表示描点点数，kk表示滤波器类型*/*画图程序，用以画函数的

8、时域和频域的图型*/void GRAPH1(float a,int n,int ii,int high)/*n表示描点点数，ii表示步长，high表示所控制的画图函数的高度*/2.3 滤波器模块本项目选择使用FIR滤波器设计方法。FIR滤波器通常采用窗函数方法来设计。窗设计的基本思想是，首先选择一个适当的理想选频滤波器（它总是具有一个非因果，无限持续时间脉冲响应），然后截取（加窗）它的脉冲响应得到线性相位和因果FIR滤波器。因此这种方法的重点是选择一个合适的窗函数和理想滤波器。我们用表示理想的选频滤波器，它在通带上具有单位增益和线性相位，在阻带上具有零响应。一个带宽的低通滤波器由下式给定：为了

9、从得到一个FIR滤波器，必须同时在两边截取。而要得到一个因果的线性相位滤波器，它的长度为，必须有：这种操作叫做加窗，可以看作是与窗函数的乘积：其中 根据的不同定义，可以得到不同的窗结构。在频域中，因果FIR滤波器响应由和窗响应的周期卷积得到，即：    常用的窗函数有矩形窗、巴特利特（BARTLETT）窗、汉宁（HANNING）窗、海明（HAMMING）窗、布莱克曼（BLACKMAN）窗、凯泽（KAISER）窗等。窗函数设计的具体步骤见流程图2-1：图2-1 窗函数设计流程图图2-2 滤波器指标选择窗和理想响应流程本项目实现提供四个窗函数选择，分别是矩形窗、哈明窗、汉宁窗、

10、布莱克曼窗，通过给定的数字滤波器指标可以求出wn(n)窗函数，wn(n)窗函数得到程序实现上见附录程序中JXC(),HNC(),HMC(),BLKMC()四个窗函数子函数。具体选择的窗函数类型需根据as阻带最小衰减，ws阻带截止频率，wp通带截止频率,wsl、wsh、wpl、wph带通带阻滤波器的上下边带阻带通带截止频率确定，由指标选择窗函数的子程序见choose_window()子程序，具体见流程图2-2。理想单位脉冲函数Hd(n)获取上，因为我们设计的滤波器为低通、高通、带通、带阻四种，故需要四个能产生理想逼近单位响应函数的子函数。查阅资料可知，对于高通滤波器函数来说，其理想频率响应和理想

11、单位脉冲响应函数分别为：对于带通滤波器函数来说，其理想频率响应和理想单位脉冲响应函数分别为：对于带阻滤波器函数来说，其理想频率响应和理想单位脉冲响应函数分别为：所以，由公式可以编写具体的函数产生程序，具体函数的产生见附录子函数HDN(),HDN1(),HDN2(),HDN3()程序。滤波器的单位脉冲响应h(n)得到是窗函数设计滤波器的最后一步，又流程图3-1可知，h(n)=Wn(n)*Hd(n)，通过一个简单的for循环即可实现最后的相乘工作，即程序：for(i=0;i<n;i+)/*时域加窗处理*/ Hri=Hdi*wni; Hii=0.0; 2.4 滤波模块滤波模块是设计的多频窄带数