现代信号处理课程设计实验报告实验课题：现代信号处理专业班级：学生姓名：学生学号：指导老师：完成时间：目录一．前言-------------------------------------------------2二．课程设计内容要求及题目-------------------------3 三．设计思想和系统功能结构及功能说明-----------4四．关键部分的详细描述和介绍，流程图描述关键模块和设计思想--------------------------------------------------7五．问题分析及心得体会--------------------------20 六．参考文献------------------------------------------21 七．附录：程序源代码清单------------------------21一、前言数字滤波在通信、图像编码、语音编码、雷达等许多领域中有着十分广泛的应用。
目前，数字信号滤波器的设计在图像处理、数据压缩等方面的应用取得了令人瞩目的进展和成就。
它是数字信号处理理论的一部分。
数字信号处理主要是研究用数字或符号的序列来表示信号波形，并用数字的方式去处理这些序列，以便估计信号的特征参量，或削弱信号中的多余分量和增强信号中的有用分量。
具体来说，凡是用数字方式对信号进行滤波、变换、调制、解调、均衡、增强、压缩、固定、识别、产生等加工处理，都可纳入数字信号处理领域。
数字信号处理学科的一项重大进展是关于数字滤波器设计方法的研究。
关于数字滤波器，早在上世纪40年代末期就有人讨论设计它的可能性问题，在50年代也有人讨论过数字滤波器，但直到60年代中期，才开始形成关于数字滤波器的一整套完整的正规理论。
在这一时期，提出了各种各样的数字滤波器结构，有的以运算误差最小为特点，有的则以运算速度高见长，而有的则二者兼而有之。
出现了数字滤波器的各种实现方法，对递归和非递归两类滤波器作了全面的比较，统一了数字滤波器的基本概念和理论。
数字滤波器与模拟滤波器相比，具有精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配以及能实现模拟滤波器无法进行的特殊滤波等优点。
上学期学习了《数字信号处理》这门课，这学期的课程设计使我更加形象具体的掌握这门课程，并且可以熟练的运用MATLAB进行编程，经过此次课程设计达到了其设计内容的基本要求，实现了基本功能，并对程序进行了优化，以达到最佳效果。
在设计前制定一个设计方案，明确设计目标，方法和过程，参阅相关文献，咨询老师和同学，严格按照方案进行课程设计，最终解决问题，完成任务要求。
二、课程设计内容要求及题目1、信号发生器用户根据测试需要，可任选以下两种方式之一生成测试信号：（1）直接输入（或从文件读取）测试序列；（2）输入由多个不同频率正弦信号叠加组合而成的模拟信号公式（如式所示）、采样频率（Hz）、采样点数，动态生成该信号的采样序列，作为测试信号。
x=100*sin(2*pi*15*t)+100*sin(2*pi*20*t)+100*sin(2*pi*30*t) 2、频谱分析使用FFT对产生的测试信号进行频谱分析并展示其幅频特性与相频特性，指定需要滤除的频带，通过选择滤波器类型（IIR/FIR），确定对应的滤波器（低通、高通）技术指标。
3、滤波器设计根据以上技术指标（通带截止频率、通带最大衰减、阻带截止频率、阻带最小衰减），设计数字滤波器，生成相应的滤波器系数，并画出对应的滤波器幅频特性与相频特性。
（1）IIR DF设计：可选择滤波器基型（巴特沃斯或切比雪夫型）；（2）FIR DF设计：使用窗口法（可选择窗口类型，并比较分析基于不同窗口、不同阶数所设计数字滤波器的特点）。
4、数字滤波根据设计的滤波器系数，对测试信号进行数字滤波，展示滤波后信号的幅频特性与相频特性，分析是否满足滤波要求（对同一滤波要求，对比分析各类滤波器的差异）。
（1）IIRDF：要求通过差分方程迭代实现滤波（未知初值置零处理）；（2）FIRDF：要求通过快速卷积实现滤波（对于长序列，可以选择使用重叠相加或重叠保留法进行卷积运算）。
5、选做内容将一段语音作为测试信号，通过频谱展示和语音播放，对比分析滤波前后语音信号的变化，进一步加深对数字信号处理的理解。
三、设计思想和系统功能结构及功能说明1、设计思想首先设定由多个不同频率正弦信号叠加组合而成的模拟信号函数、决定采样频率（Hz）和采样点数，然后使用FFT对产生的测试信号进行频谱分析并展示其幅频特性与相频特性，指定需要滤除的频带，接下来使用汉宁窗低通滤波器与巴特沃斯高通滤波器进行滤波，分析滤波后的信号的幅频特性以及相频特性，最后设计语音滤波，检查是否满足滤波要求。
2、系统结构功能说明（1）信号发生器定义一个由4个不同频率叠加的正弦函数，输出该函数，观察其相频以及幅频特性。
（2）频谱分析使用FFT对产生的测试信号进行频谱分析并展示其幅频特性与相频特性，指定需要滤除的频带，通过选择滤波器类型汉宁窗和巴特沃斯滤波器，确定对应的滤波器（低通、高通）技术指标Fp、Fc、Rp、Rs。
然后对其滤波后信号进行频谱分析。
（3）数字滤波器设计1、FIR 数字滤波器设计原理FIR 数字滤波器可以用下面的差分方程来描述[∑-== -=1)(*)()()()(Nmnhnxmnxmhny其中x（n）是输入序列，y（n）是输出序列，h（n）是系统的单位脉冲响应。
由于其采用卷积的方法得到输出，故可以采用快速算法FFT进行卷积计算，提高信号处理速度。
FIR 数字滤波器满足线性相位的条件是其乘法器的系数关于中点对称，又因为其结构无反馈的回路，所以是无条件的稳定系统。
2、窗函数法设计及MATLAB 实现FIR 滤波器目前常用的设计方法有窗函数法和频率采样法，窗函数法是从时域进行设计，而频率采样法是从频域进行设计。
窗函数法由于简单、物理意义清晰，因而得到了较为广泛的应用。
窗函数法设计的基本思想是：首先根据技术指标要求，选取合适的阶数N和窗函数的类型)(nω，使其幅频特性逼近理想滤波器幅频特性。
其次，因为理想滤波器的)(nhd 是无限长的，所以需要对)(nhd进行截断，数学上称这种方法为窗函数法[1]。
)()(1jw d d e H FT n h -=)()()(n n h n h d ω=其中)(n h d 是希望逼近的滤波器的单位脉冲响应，)(ωj d e H 是希望逼近的滤波器的幅频特性，)(n h 是所设计的滤波器的单位脉冲响应，)(n ω是窗函数。
在截断的过程中要保证FIR 滤波器的稳定性与线性相位的特性。
由于截断效应使所设计滤波器的幅频特性存在误差，称之为吉布斯效应。
最后，验证所设计的滤波器是否满足要求，若不满足，重新设计。
常用的窗函数有矩形窗、三角窗、汉宁窗和凯塞窗等。
窗函数的选择一般要满足以下三个条件[：（1）具有较低的旁瓣幅度，尤其是第一旁瓣幅度。
（2）旁瓣幅度下降速度要快，以利于增加阻带衰减。
（3）主瓣宽度要窄，以获得较陡的过渡带。
FIR 滤波器的设计任务是选择有限长度的 h(n), 使函数H(ej ω)满足技术要求。
我选择的是汉宁窗低通滤波。
3、 IIR 数字滤波器的设计IIR 数字滤波器冲激响应无限宽，与模拟滤波器相匹配。
因此可采取在模拟滤波器设计的基础上进一步变换的方法设计IIR 滤波器，即通过已知的模拟滤波器系统的系统函数H(s)来设计数字滤波器的系统函数H(z )，主要是通过脉冲响应不变法或双线性变换法完成s 平面到Z 平面的变换。
通过典型的模拟滤波器（如巴特沃斯滤波器、切比雪夫滤波器等）可实现一定参数要求的数字滤波器。
IIR 数字滤波器的设计具体步骤：<1> 按照既定规则将数字滤波器技术指标转为模拟低通滤波器技术指标。
<2> 以转换后的技术指标为依据设计模拟低通滤波器H(s)（传递函数) 。
<3> 再按照一定规则将H(s)转换成H(z)（数字滤波器的传递函数）。
若设计的是高通、带通或带阻滤波器，还需进行以下步骤：将高通、带通或带阻数字滤波器的技术指标转换为低通模拟滤波器的技术指标, 然后设计出低通H(s), 再将H(s)转换为H(z)。
我选择的是巴特沃斯高通滤波器。
（4）语音滤波1、语音信号的采集MATLAB软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。
通过wavread函数的使用，我们很快理解了采样频率、采样位数等概念。
这里我直接采用了一段现成的wav格式的语音信号。
2、语音信号的频谱分析首先画出语音信号的时域波形；然后对语音号进行快速傅里叶变换，得到信号的频谱特性。
四、关键部分的详细描述和介绍，流程图描述关键模块和设计思想1、信号发生器（1）输入信号程序如下：N=512;T=0.01;Fs=1/T;t=0:pi/256:2*pi;f1=7;f2=37;f3=39;f4=43;x=100*sin(2*pi*f1*t)+100*sin(2*pi*f2*t)+100*sin(2 *pi*f3*t)+100*sin(2*pi*f4*t);plot(t,x,'g');grid;title('原始输入信号');xlabel('时间/t');ylabel('振幅');X=fft(x,N);plot(X);grid;title('输入信号频谱');设计输入信号并对其进行傅里叶变换。
输入信号如下图（2）输入信号频谱分析幅频特性Fs=100;len=512;f=Fs*(0:len/2-1)/len;plot(f,abs(X(1:len/2)));grid title('原始信号幅频特性曲线');相频特性Fs=100;len=512;f=Fs*(0:len/2-1)/len;plot(f,angle(X(1:len/2)));grid; title('相频特性曲线');幅频特性相频特性2、巴特沃斯高通滤波器高通滤波器采用Butterworth滤波器, 因为它具有较平坦的幅频特性, 而且其元件值比较合乎实际、对精度要求稍宽。
用[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s')语句实现Butterworth滤波器阶数和边缘频率的确定；[b,a]=butter(N,Wn,'s')语句来实现Butterworth滤波器的设计。
其中Wp为通带截止频率；Rp为通带截止频率处的衰减；Ws为阻带截止频率；Rs为阻带截止频率处的衰减；[ b,a]是滤波器传递函数多项式的系数；n为滤波器的阶数；Wn为截止频率。