数字信号处理课程设计报告.doc
数字信号处理课程设计报告
数字信号处理课程设计实验报告一、课程设计内容要求1、课程设计题目设计并实现一个流程如图所示的信号处理演示系统，该系统包含信号发生器、频谱分析、滤波器设计、数字滤波和输出信号分析5个主要模块，各模块的具体功能要求如下：1)信号发生器根据信号类型不同可分为两大类：（1）静态型：直接输入测试信号系列。
（2）动态型：输入如下式所示的由多个不同频率正弦信号叠加组合而成的模拟信号公式，指定采样频率和采样点数，动态生成该信号的采样序列，作为测试信号。
100sin（2pif1t）+100sin(2pif2t)+…+100sin(2pifnt)2)频谱分析是用FFT对产生的测试信号进行频域变换，展示其幅频、相频特性，指定需要滤出或保留的频带，通过选择滤波器类型（IIR或FIR），确定对应的滤波器技术指标（低通、高通、带通、带阻）。
3)滤波器设计根据IIR/FIR数字滤波器技术指标设计滤波器，生成相应的滤波器系数，并展示对应的滤波器幅频、相频特性。
(1)IIR DF设计：使用双线性变换法，可选择滤波器类型（巴特沃斯/切比雪夫型）；(2)FIR DF 设计：使用窗口法，可选择窗口类型。
4)数字滤波根据设计的滤波器系数，对测试信号进行滤波，得到滤波后信号。
（1） IIR DF：要求通过差分方程迭代实现滤波，未知初值置0处理；（2） FIR DF：要求通过快速卷积实现滤波，可以选择使用重叠相加或重叠保留法进行卷积运算，并动态展示卷积运算的详细过程。
5)输出信号分析展示滤波后信号的幅频和相频特性，分析是否满足滤波要求。
对同一滤波要求，根据输出信号频谱，对比分析各类滤波器的差异。
2、设计题目要求使用MATLAB编程实现上述信号处理演示系统，具体要求如下：（1）系统应使用图形用户界面（GUI）；（2）系统功能至少包括非语音信号的低通和高通滤波；（3）滤波器设计模块应避免使用MATLAB工具箱函数；（4）IIR DF设计必须可选基于巴特沃斯或切比雪夫1型；（5）FIR DF设计必须可选择各类窗口，且FIR滤波可选长序列卷积方法。
数字信号处理课程设计报告
目录一、课程设计要求二、设计过程（1）设计题目（2）设计源代码（3）设计结果（4）结果分析三、设计总结与心得体会四、课程设计指导书一、课程设计要求1、课程设计指导书①《数字信号处理（第二版）》，丁玉美等，西安电子科技大学出版社；②《MATLAB 及在电子信息课程中的应用》，陈怀琛等，电子工业出版社。
2、课程设计内容：⑴语音信号去噪处理主要要求：1)在windows系统下的录音机录制一段1s左右的语音信号作为原声信号，在MATLAB软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数；2)画出语音信号的时域波形，对采样后的语音进行fft变换，得到信号的频谱特性；对语音信号分别加入正弦噪声和白噪声，画出加噪信号的时域波形和频谱图；3)根据对加噪语音信号谱分析结果，确定滤除噪声滤波器的技术指标，设计合适的数字滤波器，并画出滤波器的频域响应；4)用所设计的滤波器对加噪的信号进行滤波，在同一个窗口画出滤波前后信号的时域图和频谱图，对滤波前后的信号进行对比，分析信号变化；5)利用sound(x)回放语音信号，验证设计效果。
⑵语音信号的延时和混响主要要求：1)利用Windows下的录音机或其他软件，录制一段自己的语音信号，时间控制在1s左右，并对录制的信号进行采样；2)语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图；3)将信号加入延时和混响，再分析其频谱，并与原始信号频谱进行比较；4)设计几种特殊类型的滤波器：单回声滤波器，多重回声滤波器，全通结构的混响器，并画出滤波器的频域响应；5)用自己设计的滤波器对采集的语音信号进行滤波；6)分析得到信号的频谱，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；7)回放语音信号。
⑶数字滤波器的设计及实现主要要求：1)调用信号产生函数mstg产生三路抑制载波调幅信号相加构成的复合信号st，观察st的时域波形和幅频特性曲线；2)由要求将st中的三路调幅信号分离，通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率，要求滤波器的通带最大衰减为0.1dB，阻带最小衰减为60dB；3)编程调用MATLAB滤波器设计函数分别设计这三个数字滤波器，并绘图显示其幅频特性曲线；4)调用滤波函数filter，用所设计的三个滤波器分别对复合信号st进行滤波，分离出st中的三路不同载波频率的调幅信号，并绘图显示滤波后信号的时域波形和频谱，观察分离效果。
数字信号处理DSP 课程设计报告
语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波摘要本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。
关键词课程设计；滤波去噪；巴特沃斯滤波器；脉冲响应不变法；MATLAB1.课程设计目的和要求1.1 课程设计目的《数字信号处理》课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。
本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。
这一点与验证性的基本实验有本质性的区别。
开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。
1.2课程设计的要求（1）滤波器指标必须符合工程设计。
（2）设计完后应检查其频率响应曲线是否满足指标。
（3）处理结果和分析结论应该一致，而且应符合理论。
（4）独立完成课程设计并按要求编写课程设计报告。
2 .设计原理用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR 滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。
2.1 IIR 滤波器从离散时间来看，若系统的单位抽样(冲激)响应延伸到无穷长，称之为“无限长单位冲激响应系统”，简称为IIR 系统。
无限长单位冲激响应（IIR ）滤波器有以下几个特点:(1)系统的单位冲激响应h(n)是无限长； (2)系统函数H(z)在有限z 平面（0<z <∞）； (3) 结构上存在着输出到输入的反馈，也就是结构上是递归型的。
IIR 滤波器采用递归型结构，即结构上带有反馈环路。
数字信号处理课程设计报告
课程设计报告课程名称数字信号处理课题名称数字滤波器设计及在语音信号分析中的应用专业通信工程班级通信工程1101学号************姓名皮锋指导教师胡瑛张细政2013年9月8日湖南工程学院课程设计任务书课程名称数字信号处理课题数字滤波器设计及在语音信号分析中的应用专业班级通信工程1101学生姓名皮锋学号************指导老师胡瑛张细政审批任务书下达日期2013 年9月 1 日任务完成日期2013 年9月8日《数字信号处理》课程设计任务书一、设计目的综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用MATLAB 作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。
二、设计要求1、MATLAB 的使用，掌握MATLAB 的程序设计方法。
2、Windows 环境下语音信号采集的方法。
3、数字信号处理的基本概念、基本理论和基本方法。
4、TLAB 设计FIR 和nR 数字滤波器的方法。
5、用 MATLAB 对信号进行分析和处理6、计报告4000以上，含程序设计说明，用户使用说明，源程序清单及程序框图。
7、机演示。
8、有详细的文档。
文档中包括设计思路、设计仿真程序、仿真结果及相应的分析与结论。
三、进度安排第一周星期一：课题讲解，查阅资料星期二: 总体设计，详细设计星期三：编程，上机调试、修改程序星期四：上机调试、完善程序星期五：答辩星期六-星期天：撰写课程设计报告附：课程设计报告装订顺序：封面、任务书、目录、正文、评分、附件（A4大小的图纸及程序清单）。
正文的格式:一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。
正文的内容:一、课题的主要功能；二、课题的功能模块的划分；三、主要功能的实现；四、程序调试；五、总结；六、附件（所有程序的原代码，要求对程序写出必要的注释）；七、评分表。
目录1. 课题的主要功能 (1)1.1问题描述 (1)1.2功能要求 (1)2. 课题的功能模块的划分 (2)2.1语音信号采集 (2)2.3 含噪语音信号合成 (4)2.4 数字滤波器设计及滤波 (4)2.5 回放语音信号 (4)3. 滤波器实现的原理 (5)3.1 带阻滤波器 (6)3.2 低通滤波器 (8)3.3 高通滤波 (10)4. 程序调试 (12)5. 总结 (13)6. 附件 (14)6.1 带阻滤波 (14)6.2 低通滤波 (15)6.3 高通滤波 (16)7. 评分表 (18)1.课题的主要功能1.1问题描述录制一段课程设计学生的语音信号并保存为文件，要求长度不小于10秒，并对录制的信号进行采样，使用MATLAB绘出采样后的语音信号的时域波形和频谱图。
数字信号处理课程设计报告
数字信号与处理课程设计报告课程设计题目：按频率抽选（DIF）的基-2FFT算法的matlab实现学号：********学生姓名：***班级：081132专业：信息工程指导教师：***2011年6月19 日目录一．数字信号处理的简介...........................................................................................................- 2 -二．课程设计的目的与要求.......................................................................................................- 2 -三．时间抽取法FFT的运算特点 ...............................................................................................- 2 -3.1 蝶形运算................................................ - 2 -3.2 原位计算................................................ - 3 -3.3 序数重排................................................ - 3 -3.4 蝶形类型随迭代次数成倍增加.............................. - 4 -四．基2按频率抽取（DIF）的FFT算法..................................................................................- 4 -五．程序代码及在MATLAB的实现结果 ...................................................................................- 6 -六．课程设计的心得...................................................................................................................- 9 -一．数字信号处理的简介广义来说，数字信号处理是研究用数字方法对信号进行分析、变换、滤波、检测、调制、解调以及快速算法的一门技术学科。
数字信号处理课程设计报告
数字信号处理课程设计报告《应用Matlab对信号进行频谱分析及滤波》专业：班级：姓名：指导老师：二0 0五年一月一日目录设计过程步骤（）2.1 语音信号的采集（）2.2 语音信号的频谱分析（）2.3 设计数字滤波器和画出其频谱响应（）2.4 用滤波器对信号进行滤波（）2.5滤波器分析后的语音信号的波形及频谱（）●心得和经验（）设计过程步骤2.1 语音信号的采集我们利用Windows下的录音机，录制了一段开枪发出的声音，时间在1 s内。
接着在C盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。
通过wavread函数和sound的使用，我们完成了本次课程设计的第一步。
其程序如下：[x,fs,bite]=wavread('c:\alsndmgr.wav',[1000 20000]);sound(x,fs,bite);2.2 语音信号的频谱分析首先我们画出语音信号的时域波形；然后对语音信号进行频谱分析，在Matlab中，我们利用函数fft对信号进行快速傅里叶变换，得到信号的频谱特性性。
到此，我们完成了课程实际的第二部。
其程序如下：n=1024;subplot(2,1,1);y=plot(x(50:n/4));grid on ;title('时域信号')X=fft(x,256);subplot(2,1,2);plot(abs(fft(X)));grid on ;title('频域信号');运行程序得到的图形：2.3 设计数字滤波器和画出其频谱响应紧接着着我们做了一个数字滤波器：采样频率10Hz，通带截止频率fp=3Hz，阻带截止频率fs=4Hz通带衰减小于1dB，阻带衰减大于20dB，我们主要使用双线性变换法由模拟滤波器原型设计数字滤波器。
数字信号处理课程设计
数字信号处理课程设计报告题目3运用matlab的数字音效处理系统设计系（院）：电子信息工程学院专业：通信工程班级：13通信学号：xxx姓名：xxx2015年12月19日一、摘要数字信号处理(Digital Signal Processing)技术，从20世纪60年代以来，随着计算机科学和信息科学发展，数字处理技术应运而生并得以快速发展。
语言是人们进行信息沟通的主要方式之一，它具有直接、自然、方便等优点。
语音则是语言的物理层表达方式。
语音处理主要是对语音进行机器处理，以达到传输、自动识别、机器理解等目的。
进行了语音处理过程的滤波、采样、傅立叶变换和谱包络提取的算法实现研究，讨论了在算法的DSP实现方法，Matlab 语言是一种广泛应用于工程计算及数值分析领域的新型高级语言, Matlab 功能强大、简单易学、编程效率高。
特别是Matlab 还具有信号分析工具箱, 不需具备很强的编程能力, 就可以很方便地进行信号分析、处理和设计.数字信号处理主要是研究用数字或符号序列表示和处理信号。
处理的目的可以是削弱信号中的多余内容，滤除混杂的噪声和干扰，或者是将信号变换为容易分析和识别的形式，便于估计和选择它的特征参数。
声音信号是一维连续信号,而计算机只能处理离散信号。
为了从离散信号还原连续信号,根据采样定理,可以确定采样频率的最小值。
wav文件是一种数字声音文件格式,本课程设计基于Matlab分析了wav声音文件频谱与声音的关系。
通过采集个人的一段声音进行频谱分析等处理，然后设计数字滤波器处理这个原始声音的wav文件,并比较滤波以后输出声音信号与原声音信号的异同。
基于MATLAB 的语音信号特技处理, 其实质是针对于不同的滤波器性能的分析与比较。
对于滤波器的设计可以使用MATLAB 中自带的工具, 也可以通过编程或直接利用MATLAB中的函数库来实现。
二、课程设计题目描述和要求2.1设计内容1、录制采集一段自己的语音信号；2、设计滤波器，实现单回声、多回声、混响效果和和声效果；2.2设计要求1、熟悉离散信号和系统的时域特性。
数字信号处理课设报告
数字信号处理课程设计题目无限冲激响应(IIR)数字滤波器的直接型专业电子信息工程姓名卢想学号 20130210410117时间 2016年 6 月 22日目 录一、课程设计目的 ....................................................... 1 3、掌握滤波器频响特性的绘制方法 .......................................... 2 二、课程设计原理与计算方法 ............................................... 2 三、课程设计内容 ......................................................... 3 四、程序与调试 .. (4)4.1激励信号1.15 ,2,5.0 )(cos )(001====-N a n R n en x N anπωω ............. 4 4.2.激励信号2 )()(2n n x δ= .. (5)4.3.扩展函数impseq ...................................................... 6 五、课程设计总结 (6)参考文献 (7)一、课程设计目的1、掌握IIR 滤波器的直接II 型、级联型和并联型三种结构的基本原理和特点。
2、掌握利用Matlab 实现IIR 滤波器的三种结构的程序设计方法，并能够进行三者之间的相互转换。
3、掌握滤波器频响特性的绘制方法 二、课程设计原理与计算方法按照结构划分数字滤波器，有递归式和非递归式两种。
递归式数字滤波器的差分方程为∑∑==-=-+Mi i Nk ki n x b k n y an y 11)()()((1)其中至少有一个0≠ka .非递归式数字滤波器的差分方程为∑=-=M i i i n x b n y 1)()( (2)可以看出递归式数字滤波器的响应)(n y 不仅与激励)(n x 有关，而且与以前的输出信号)(k n y -有关；而非递归式数字滤波器的响应)(n y 仅只与激励)(n x 有关。
数字信号处理课程设计实验报告
数字信号处理课程设计实验报告通信与信息工程学院数字信号处理课程设计班级：姓名：学号：指导教师：设计时间：成绩：评语：通信与信息工程学院二〇一四年题目一：采样定理的验证1.课程设计目的及要求：1). 掌握利用MATLAB分析系统频率响应的方法，增加对仿真软件MATLAB的感性认识，学会该软件的操作和使用方法。
2). 掌握利用MATLAB实现连续信号采用与重构的方法，加深理解采样与重构的概念。
3 ). 初步掌握线性系统的设计方法，培养独立工作能力。
4). 学习MATLAB中信号表示的基本方法及绘图函数的调用，实现对常用连续时间信号的可视化表示，加深对各种电信号的理解。
5). 验证信号与系统的基本概念、基本理论，掌握信号与系统的分析方法。
6). 加深对采样定理的理解和掌握，以及对信号恢复的必要性；掌握对连续信号在时域的采样与重构的方法。
2.详细设计过程及调试结果：1）．设()(0.8)()n，利用filter函数求出()*()x n u nx n x n的源程序：n=0:49xn=(0.8).^nB=1A=[1,-0.8]yn=filter(B,A,xn)stem(n,yn)xlabel('n');ylabel('yn')本题验结果及分析：2）：模拟信号()sin(20),010.01a s x t t t T s π=≤≤=，在，0.050.1s s 和间隔采样得到()x n ：a.每一个s T 画出()x n 的源程序：t=0:0.01:1;T0=0.1; T1=0.01;n1=0:100; T2=0.05;n2=0:20; T3=0.1;n3=0:10; xt=sin(20*pi*t); xn1=sin(20*pi*n1*T1); xn2=sin(20*pi*n2*T2);xn3=sin(20*pi*n3*T3);subplot(4,1,1);plot(t,xt);title('模拟信号xt图');xlabel('t');ylabel('xt');grid on; subplot(4,1,2);stem(n1,xn1,'.');title('0.01s采样图');xlabel('n');ylabel('xn1');grid on; subplot(4,1,3);stem(n2,xn2,'.');title('0.05s采样图');xlabel('n');ylabel('xn2');grid on; subplot(4,1,4);stem(n3,xn3,'.');title('0.1s采样图');xlabel('n');ylabel('xn3');grid on;调试结果分析：b.采用sin c内插从样本()y t的源程序：x n重建模拟信号()at=0:0.01:1;T0=0.1;xt=sin(20*pi*t);T1=0.01;n1=0:100;T2=0.05;n2=0:20;T3=0.1;n3=0:10;xt=sin(20*pi*t);subplot(4,1,1);plot(t,xt);title('原信号xt模拟图');xlabel('t');ylabel('xt');grid on;xn1=sin(20*pi*n1*T1);xn2=sin(20*pi*n2*T2);xn3=sin(20*pi*n3*T3);t1=0:T1:1;t2=0:T2:1;t3=0:T3:1;tn1=ones(length(n1),1)*t1-n1'*T1*ones(1,length(t1));tn2=ones(length(n2),1)*t2-n2'*T2*ones(1,length(t2));tn3=ones(length(n3),1)*t3-n3'*T3*ones(1,length(t3));yt1=xn1*sinc(tn1*pi/T1);subplot(4,1,2);plot(t1,yt1);axis([ 0 1 -1 1]); title('sinc内插0.01恢复的xt1图');xlabel('n');ylabel('xt1');grid on; yt2=xn2*sinc(tn2*pi/T2);subplot(4,1,3);plot(t2,yt2);axis([ 0 1 -1 1]); title('sinc内插0.05恢复的xt2图');xlabel('n');ylabel('xt2');grid on; yt3=xn3*sinc(tn3*pi/T3);subplot(4,1,4);plot(t3,yt3);axis([ 0 1 -1 1]); title('sinc内插0.1恢复的xt3图');xlabel('n');ylabel('xt3');grid on;调试结果分析：c.采用三次样条内插从样本()x n重建模拟信号源程序：t=0:0.01:1; xt=sin(20*pi*t);T1=0.01;n1=0:100;T2=0.05;n2=0:20;T3=0.1;n3=0:10;T1=0:T1:1;T2=0:T2:1;T3=0:T3:1;xt=sin(20*pi*t);xn1=sin(20*pi*T1);xn2=sin(20*pi*T2);xn3=sin(20*pi*T3);yt1=spline(T1,xn1,t);yt2=spline(T2,xn2,t);yt3=spline(T3,xn3,t);subplot(4,1,1);plot(t,xt);title('原信号xt模拟图');xlabel('t');ylabel('xt');grid on;subplot(4,1,2);plot(t,yt1);axis([ 0 1 -1 1]);title('三次样条0.01恢复的xt1图');xlabel('n');ylabel('xt1');grid on; subplot(4,1,3);plot(t,yt2);axis([ 0 1 -1 1]);title('三次样条0.05恢复的xt2图');xlabel('n');ylabel('xt2');grid on; subplot(4,1,4);plot(t,yt3);axis([ 0 1 -1 1]);title('三次样条0.1恢复的xt3图');xlabel('n');ylabel('xt3');grid on; 调试结果分析：3.总结体会：连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。
数字信号处理课程设计
数字信号处理 课程设计一、课程目标知识目标：1. 理解数字信号处理的基本概念、原理和方法，掌握其数学表达和物理意义；2. 掌握数字信号处理中的关键算法，如傅里叶变换、快速傅里叶变换、滤波器设计等；3. 了解数字信号处理技术在通信、语音、图像等领域的应用。
技能目标：1. 能够运用所学知识分析数字信号处理问题，提出合理的解决方案；2. 能够运用编程工具（如MATLAB）实现基本的数字信号处理算法，解决实际问题；3. 能够对数字信号处理系统的性能进行分析和优化。
情感态度价值观目标：1. 培养学生对数字信号处理学科的兴趣，激发其探索精神和创新意识；2. 培养学生严谨的科学态度和良好的团队协作精神，提高沟通与表达能力；3. 增强学生对我国在数字信号处理领域取得成就的自豪感，树立为国家和民族发展贡献力量的信心。
课程性质：本课程为专业选修课，旨在使学生掌握数字信号处理的基本理论和方法，培养其解决实际问题的能力。
学生特点：学生具备一定的数学基础和编程能力，对数字信号处理有一定了解，但缺乏系统学习和实践经验。
教学要求：结合学生特点和课程性质，注重理论与实践相结合，采用案例教学、互动讨论等教学方法，提高学生的参与度和实践能力。
通过本课程的学习，使学生能够达到上述课程目标，为后续相关课程和实际工作打下坚实基础。
二、教学内容1. 数字信号处理基础：包括数字信号、离散时间信号与系统、信号的采样与恢复等基本概念，使学生建立数字信号处理的基本理论框架。
教材章节：第一章 数字信号处理概述2. 傅里叶变换及其应用：介绍傅里叶变换的原理、性质和应用，以及快速傅里叶变换算法。
教材章节：第二章 傅里叶变换及其应用3. 数字滤波器设计：讲解数字滤波器的基本原理、设计方法和性能评价，包括IIR和FIR滤波器。
教材章节：第三章 数字滤波器设计4. 数字信号处理应用案例分析：通过通信、语音、图像等领域的实际案例，使学生了解数字信号处理技术的应用。
数字信号处理课程设计IIR.
吉林建筑工程学院电气与电子信息工程学院数字信号处理课程设计报告设计题目：IIR数字滤波器的设计专业班级：信工102学生姓名：丁航学号：10210211指导教师：杨佳吴贺君设计时间：2013.01.07－2013.01.11IIR 数字滤波器设计报告一、设计的作用、目的目的：课程设计是理论学习的延伸，是掌握所学知识的一种重要手段，对于贯彻理论联系实际、提高学习质量、塑造自身能力等于有特殊作用。
本次课程设计一方面通过MATLAB 仿真设计内容，使我们加深对理论知识的理解，同时增强其逻辑思维能力，另一方面对课堂所学理论知识作一个总结和补充。
作用：加深对脉冲响应不变法和双线性变换法设计IIR 滤波器数字滤波器基本方法的了解，熟悉这一设计的计算机编程。
观察用脉冲响应不变法和双线性变换法设计的数字滤波器和响应模拟滤波器的时域特性和频域特性，比较所涉及的数字滤波器和响应的模拟滤波器的频域特性，观察脉冲响应不变法设计中产生的频域混淆现象。
学会MATLAB 的使用，掌握运用MATLAB 设计IIR 低通滤波器。
熟悉Butterworth 滤波器、Chebyshev 滤波器和椭圆滤波器的频率特性。
二、设计任务及要求通过课程设计各环节的实践，应使学生达到如下要求：1.掌握双线性变换法及脉冲响应不变法设计IIR 数字滤波器以及窗函数法设计FIR 数字滤波器的原理、具体方法及计算机编程2.观察双线性变换法、脉冲响应不变法及窗函数法设计的滤波器的频域特性，了解各种方法的特点3.用MATLAB 画出三种方法设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查结果是否满足要求。
三、设计内容已知通带截止频率kHz f p 2.0=,通带最大衰减dB P 1=α,阻带截止频率kHz f s 3.0=，阻带最小衰减dB s 25=α，T=1ms ，按照以上技术要求，用脉冲响应不变法和双线性变换法设计巴特沃斯数字低通滤波器，并观察所设计数字滤波器的幅频特性曲线。
数字信号处理课程设计报告_4
目录一、课程设计的目的 (2)二、数字滤波器的设计步骤 (2)2.2、IIR数字滤波器与FIR数字滤波器的区别 (2)三、IIR数字滤波器 (3)3.1、IIR数字滤波器的特点 (3)3.1.2、IIR数字滤波器采用递归型结构 (3)3.1.3、借助成熟的模拟滤波器的成果 (3)3.1.4、需加相位校准网络 (3)3.2、用双线性法设计IIR数字滤波器 (3)3.3、巴特沃斯低通滤波器的设计 (4)3.4、巴特沃斯高通滤波器的设计 (5)3.4.1、巴特沃斯高通滤波器各参数图形 (5)3.4.2、巴特沃斯高通滤波器滤波效果图 (5)四、FIR数字滤波器 (5)4.1、FIR滤波器的特点 (5)4.2、窗函数法设计FIR数字滤波器 (6)五、程序实例源码 (8)六、问题分析 (12)七、心得体会 (13)八、参考文献 (13)一、课程设计的目的数字滤波是数字信号分析中最重要的组成部分之一，与模拟滤波相比，它具有精度和稳定性高、系统函数容易改变、灵活性强、便于大规模集成和可实现多维滤波等优点。
在信号的过滤、检测和参数的估计等方面，经典数字滤波器是使用最广泛的一种线性系统。
本次课程设计是通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用双线性变换法设计 IIR 数字滤波器和用窗函数法设计FIR数字滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。
二、数字滤波器的设计步骤2.1、不论是IIR滤波器还是FIR滤波器的设计都包括三个步骤：(1)按照实际任务的要求，确定滤波器的性能指标。
(2)用一个因果、稳定的离散线性时不变系统的系统函数去逼近这一性能指标。
根据不同的要求可以用IIR系统函数，也可以用FIR系统函数去逼近。
(3)利用有限精度算法实现系统函数，包括结构选择，字长选择等。
2.2、IIR数字滤波器与FIR数字滤波器的区别2.2.1、单位响应IIR数字滤波器单位响应为无限脉冲序列，而FIR数字滤波器单位响应为有限的；FIR滤波器，也就是“非递归滤波器”，没有引入反馈。
数字信号处理课程设计报告
1 绪论1.1设计目的通过本次数字信号处理课程设计，综合运用数字信号处理课程和其他有关课程的理论和生产实际知识去分析和解决具体问题，并使所学知识得到进一步巩固、深化和发展。
学习设计滤波器的一般方法，通过课程设计树立正确的设计思想，提高分析问题、解决问题的能力。
进行设计基本技能的训练，如查阅设计资料和手册、程序的设计、调试等。
1.2设计要求(1)掌握数字信号处理的基本概念，基本理论和基本方法。
(2)熟悉离散信号和系统的时域特性。
(3)掌握序列快速傅里叶变换方法。
(4)学会MATLAB的使用，掌握MATLAB的程序设计方法。
(5)掌握利用MATLAB对语音信号进行频谱分析。
(6)掌握滤波器的网络结构。
1.3 设计内容1.3.1预习题1、设计卷积运算的演示程序：（1）可输入任意两个序列x1(n)、x2(n)，指定x1(n)为自己的学号，例如x1(n)={2, 0, 1, 0, 5, 7, 0, 5, 0, 2, 0, 1}。
x2(n)的内容和长度自选。
例如x2(n)={0.5, 2.43, 6.17, 12.93, 22.17, 32.25, 40.88, 45.87, 45.87, 40.88, 32.25, 22.17, 12.93, 6.17, 2.43,1.007}。
（2）分别动态演示两个序列进行线性卷积x1(n)*x2(n)和圆周卷积x1(n)⊙x2(n)的过程；要求分别动态演示翻转、移位、乘积、求和的过程；（3）圆周卷积默认使用两个序列中的最大长度，但卷积前可以指定卷积长度N用以进行混叠分析；（4）改变圆周卷积长度N，根据实验结果分析两类卷积的关系；(5)在计算机操作系统中选一段声音文件(XP系统在“C:\WINDOWS\Media”)，读取文件取 10ms的声音数据产生时域序列x1(n)，序列内容自定义。
利用x2(n)={0.5, 2.43, 6.17, 12.93, 22.17, 32.25, 40.88, 45.87, 45.87, 40.88, 32.25, 22.17, 12.93, 6.17, 2.43, 1.007}。
数字信号处理课程设计报告
通信与信息工程学院数字信号处理课程设计班级：电子信息科学与技术姓名：学号：指导教师：设计时间：成绩：评通信与信息工程学院二〇一三年题目一 采样定理的验证1． 设()(0.8)()n x n u n =，利用filter 函数求出()*()x n x n （其中*为线性卷积）的前50个样本，并作图。
（1） n=0:49; xn=0.8.^n; B=1;A=[1,-0.8];yn=filter(B,A,xn); stem(n,yn) grid on; xlabel('n');ylabel('yn=x(n)*x(n)')ny n =x (n )*x (n )2．考虑模拟信号()sin(20),010.01a s x t t t T s π=≤≤=，在，0.050.1s s 和间隔采样得到()x n 。
（1）对每一个s T 画出()x n 。
（2）采用sin c内插从样本()y t，并从图中求出x n重建模拟信号()a在()y t中的频率（不管末端效果）。
a（3）采用三次样条内插从样本()x n重建模拟信号，并从图中求出在()y t中的频率。
a（1）clc;t=0:0.001:1;Ts1=0.01;N1=1/Ts1;n1=0:N1;xn1=sin(20*pi*n1*Ts1);subplot(3,1,1);stem(n1,xn1,'.');axis([0 100 -1 1]);xlabel('n');ylabel('x_1(n)');title('Ts=0.01');t=0:0.001:1;Ts2=0.05;N2=1/Ts2;n2=0:N2;xn2=sin(20*pi*n2*Ts2);subplot(3,1,2);stem(n2,xn2,'.');axis([0 20 -1 1]);xlabel('n');ylabel('x_2(n)');title('Ts=0.05');t=0:0.001:1;Ts3=0.1;N3=1/Ts3;n3=0:N3;xn3=sin(20*pi*n3*Ts3);subplot(3,1,3);stem(n3,xn3,'.');axis([0 10 -1 1]);xlabel('n');ylabel('x_3(n)');title('Ts=0.1');102030405060708090100n x 1(n )Ts=0.012468101214161820n x 2(n )Ts=0.0512345678910nx 3(n )Ts=0.1（2）t=0:0.001:1;T=0.1;xt=sin(20*pi*t);subplot(4,1,1);plot(t,xt); title('xt 原图');xlabel('t');ylabel('xt');grid on;Ts1=0.01;n1=0:100;t1=0:Ts1:1;xn1=sin(20*pi*n1*Ts1);tn1=ones(length(n1),1)*t1-n1'*Ts1*ones(1,length(t1));%ones （M,N)：生成M*N 大小的全1矩阵yt1=xn1*sinc(tn1*pi/Ts1);subplot(4,1,2);plot(t1,yt1);axis([0 1 -1 1]); title('sinc 内插0.01恢复的xt1图');xlabel('n');ylabel('xt1');grid on;Ts2=0.05;n2=0:20;t2=0:Ts2:1;xn2=sin(20*pi*n2*Ts2); tn2=ones(length(n2),1)*t2-n2'*Ts2*ones(1,length(t2));yt2=xn2*sinc(tn2*pi/Ts2);subplot(4,1,3);plot(t2,yt2);axis([0 1 -1 1]); title('sinc 内插0.05恢复的xt2图');xlabel('n');ylabel('xt2');grid on;Ts3=0.1;n3=0:10;t3=0:Ts3:1;xn3=sin(20*pi*n3*Ts3); tn3=ones(length(n3),1)*t3-n3'*Ts3*ones(1,length(t3));yt3=xn3*sinc(tn3*pi/Ts3);subplot(4,1,4);plot(t3,yt3);axis([0 1 -1 1]); title('sinc 内插0.1恢复的xt3图');xlabel('n');ylabel('xt3');grid on;00.10.20.30.40.50.60.70.80.91-11xt 原图tx t00.10.20.30.40.50.60.70.80.91-11sinc 内插0.01恢复的xt1图nx t 100.10.20.30.40.50.60.70.80.91-11sinc 内插0.05恢复的xt2图nx t 200.10.20.30.40.50.60.70.80.91-11sinc 内插0.1恢复的xt3图nx t 3（3）Ts1=0.01;Fs=1/Ts1;xa1=spline(Ts1*n1,xn1,t);subplot(3,1,1);plot(t,xa1);axis([0,1,-1,1]); ylabel('x_a(t)');title('Ts=0.01'); Ts2=0.05;Fs=1/Ts2;xa2=spline(Ts2*n2,xn2,t);subplot(3,1,2);plot(t,xa2);axis([0,1,-1,1]); ylabel('x_a(t)');title('Ts=0.05'); Ts3=0.1;Fs=1/Ts3;xa3=spline(Ts3*n3,xn3,t);subplot(3,1,3);plot(t,xa3);axis([0,1,-1,1]); ylabel('x_a(t)');title('Ts=0.1');00.10.20.30.40.50.60.70.80.91-101x a (t )Ts=0.0100.10.20.30.40.50.60.70.80.91-101x a (t )Ts=0.0500.10.20.30.40.50.60.70.80.91-11x a (t )Ts=0.1结果分析与体会：采用sin c 内插从样本()x n 重建模拟信号()a y t 时，0.01s 和0.05s 依据奈奎斯特采样定理应该能恢复原来的波形，但因为采样点太少，间距分散导致0.05s 时无法恢复。
数字信号处理课程设计报告
设计一信号的基本运算一、设计目的熟悉信号的基本运算，通过运用Matlab进行仿真，加深对信号基本运算的理解。
通过对数据的处理，加深对Matlab中数据存取，数值运算，矩阵运算的方式及工作原理的了解。
二、设计原理Matlab是以矩阵为基础的一种软件，其集成了数值运算、矩阵运算、信号处理和图形等众多功能。
其中，对数据的存取都是以矩阵的方式进行的。
Matlab工具箱中提供了很多已经编写好的函数，我们这用些函数的时候只需要从工具箱中调用就可以了，这些函数都十分的方便。
如其中的wavread( )函数，我们可以用来从音频文件中获取数据，然后对这些数据进行运算，然后通过sound( )函数对音频文件进行回放；还有一些特殊矩阵的生成函数，如用函数zeros生成全0矩阵：格式B=zeros(m,n)生成m×n的全0阵；用函数ones生成全1矩阵：格式B=ones(m,n)生成m×n的全1阵；用函数rand生成随机矩阵：格式B=rand(m,n)生成m×n的随机矩阵；用函数eye生成单位阵：格式B=eye(m,n)生成m×n矩阵，其中对角线元素全为1，其他元素为0。
通过类似这样的操作，我们就可以方便的对信号进行相应的处理。
本次实验中，我们对一段音频信号，进行回音的模仿，然后经过上采样和下采样，反转的处理，并演示处理后的效果。
三、设计内容本次实验，我们通过采样得到一段以采样频率为8192Hz的语音信号x（k），然后通过编写Matlab程序对这段语音信号进行回音模仿，采用函数x(k)=x(k)+a*x(k-d),期中d为时延，a为时延信号的衰减幅度。
然后对语音信号进行下采样x(k/2)、上采样x(2k)、反转x(-k)。
下采样即在得到的语音信号的基础上，隔一个k值取一个函数值；上采样，即在得到的信号的基础上进行每两个k值之间进行插值；反转即把得到的信号的k变为-k。
通过这样的处理后，回放语音信号，观察效果，再看处理后的信号的时域波形。
数字信号处理课程设计报告
数字信号处理课程设计报告姓名：班级：学号：指导老师：华南农业大学电子工程学院一、设计方案分析题目一：①录制一段自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；②给定滤波器的性能指标（技术指标：通带截止频率1khz 带内最大衰减-3db 阻带截止频率4khz 阻带内最小最小衰减-30db），采用窗函数法或双线性变换设计滤波器，并画出滤波器的频率响应；③用自己设计的滤波器对采集的语音信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化，回放语音信号。
设计方案分析：用手机录音功能录下一段语音信号，再使用迅雷转换器将MP3音频转换成.WAV格式，使用audioread('xxx.wav')函数读取音频信号，用fft函数对语音信号进行DFT变换，得到时域和频域图；对于滤波器的性能指标通带截止频率1khz 带内最大衰减-3db 阻带截止频率4khz 阻带内最小最小衰减-30db，利用双线性变换法设计该低通滤波器，根据模拟滤波器技术指标设计为响应的模拟低通滤波器；[n,Wn]=buttord(wp,ws,Ap,As,’s’)[z,p,k]=buttap(n);设计模拟滤波器[b,a]=zp2tf(z,p,k); 将系统函数的零极点转化为系统函数一般形式的系数[bn,an]=lp2lp(b,a,Wn);其中wp=2pif1/fs，ws2pif1/fs，Ap=3，As=30再利用[bz,az]= bilinear (bn,an,1/Ts)函数把模拟滤波器转换为数字滤波器；用设计的滤波器对语音信号进行滤波，使用filter(bz,az,x)函数。
题目二：老式唱机播放音乐时，常会听到“噼啪”类噪声。
“noisysound.wav”为带有这类噪声的音频文件，试设计算法滤除噪声，并从时域和频域比较处理前后的差异。
设计方案分析：通过观察所给音频文件的时域及频域图，得知噪声大多是为高频成分，所以要滤除噪声，可采用低通滤波器对语音信号进行滤波处理。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数字信号处理课程设计报告淮阴工学院数字信号处理课程设计报告课题名称基于MATLAB的语音信号分析和处理-低通部分姓名xx 学号1234567890 班级电子0000 专业电子信息工程归口系部电子与电气工程学院起迄日期2014年12月22日2014年12月26日设计地点13306 指导教师数字信号处理课程设计指导小组提交报告日期2014年12月29日一、设计目的与任务通过数字信号处理的课程设计，使学生对信号的采集，处理，传输，显示，存储和分析等有一个系统的掌握和理解。
巩固和运用数字信号处理课程中的理论知识和实验技能，掌握最基本的数字信号处理的理论和方法，培养学生发现问题，分析问题和解决问题的能力。
二、基本要求1．基本要求设计一个语音信号分析和处理系统，要求学生对所采集的语音信号在MATLAB软件平台下进行频谱分析和处理；2．提高要求对所采集的语音信号叠加干扰噪声进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。
3．基本教学要求每组一台电脑（附话筒和耳机），电脑安装MATLAB软件。
三、设计要求本次课程设计最终要求提交设计说明书，由以下各部分组成1．理论依据根据设计要求分析系统功能，掌握设计中所需的理论（采样频率，采样位数的概念，采样定理；时域信号的DFT，FFT及频谱分析；数字滤波器的设计原理和方法，各种不同性能的滤波器的性能比较），阐述设计原理。
2．信号采集采集语音信号，画出信号的时域波形图和频谱图。
3．数字滤波器设计根据语音信号的特点，设计IIR低通数字滤波器；FIR低通数字滤波器；画出各种数字滤波器的频率响应图。
4．信号处理1）利用设计的IIR数字滤波器分别对采集的信号进行滤波处理；2）利用设计的FIR数字滤波器分别对采集的信号进行滤波处理；3）在原始的语音信号3000HZ 以上频段叠加白噪声，选择所设计的一种对叠加白噪声后的语音信号处理，滤除白噪声；4）画出处理过程中所得的各种波形及频谱图。
5．信号分析对语音信号进行回放，感觉滤波前后声音的变化。
比较滤波前后语音信号的波形及频谱，对所得结果和滤波器性能进行频谱分析，阐明原因得出结论。
4、理论依据 4.1语音信号的采集和频谱分析利用WINDOWS下的录音机，用单声道录制一段自己的声音，时间在1S内。
然后在MATLAB软件平台下，利用函数WA VREAD对语音信号进行采样，记住采样频率和采样点数。
对语音信号进行快速傅立叶变换，在一个窗口同时画出信号的时域波形图和频谱图，分析语音信号的频谱特点。
4.2 设计数字滤波器低通滤波器性能指标fp2670HZ,fs3300HZ,Ap1dB,As30dB,采样频率为10000hz。
4.3 用滤波器对信号进行滤波1）利用自己设计的IIR 数字滤波器对采集的信号进行滤波处理，在一个窗口同时画出滤波前后语音信号的波形及频谱。
2）利用自己设计的FIR数字滤波器对采集的信号进行滤波处理，在一个窗口同时画出滤波前后语音信号的波形及频谱。
在MATLAB中，FIR滤波器利用fftfilt对信号进行滤波，IIR滤波器利用filter对信号进行滤波。
比较滤波前后语音信号的波形及频谱，对所得结果和滤波器性能进行频谱分析。
4.4回放语音信号对语音信号进行回放，感觉滤波前后语音信号的变化。
4.5对语音信号加入噪声后的频谱分析和滤波在MATLAB软件平台下，在3000HZ以上频段产生一个白噪声，叠加到采集的原始语音信号中，在一个窗口同时画出白噪声的时域波形及频谱，在另一个窗口同时画出叠加噪声前后语音信号的波形及频谱；设计一种滤波器叠加噪声后的语音信号进行处理，滤除叠加的白噪声，在一个窗口同时画出滤波前后语音信号的波形及频谱；比较滤波前后语音信号的波形及频谱，对所得结果进行频谱分析；对语音信号进行回放，感觉滤波前后声音的变化。
5、基于MATLAB的语音信号分析和处理5.1信号采集我们利用函数wavread对语音信号进行采集，原始语音信号及其频谱如下图，信号采集程序见附录。
图1 原始语音信号及其频谱5.2根据性能指标进行数字滤波器设计低通滤波器性能指标fp2670HZ,fs3300HZ,Ap1dB,As30dB,采样频率为10000hz。
1）FIR低通滤波器的设计FIR低通滤波器的幅频特性如下图，设计程序见附录。
图2 FIR低通滤波器的幅频特性2）IIR低通滤波器的设计IIR低通滤波器的幅频特性如下图，设计程序见附录。
图3 IIR低通滤波器的幅频特性5.3用滤波器对信号进行滤波1）FIR低通滤波图4 FIR低通滤波前后的频谱2）IIR低通滤波图5 IIR低通滤波前后的频谱5.4回放语音信号经过加噪声处理后，可在Matlab中用函数sound对声音进行回放。
其调用格式soundy,Fs，soundy和soundy,Fs,bits。
可以察觉滤波前后的声音有明显的变化。
5.5对语音信号加入噪声后的频谱分析和滤波加入噪声后的波形和频谱如下图图6 加噪声后波形和频谱6、课程设计总结数字滤波器是一种对数字信号进行滤波处理以得到期望的响应特性的离散时间系统。
它能够有效的对数字信号进行过滤，过滤数字信号中的干扰和多余成分，保证通信过程的正常实现。
在这次的数字信号处理课程设计里，比较重要的函数有wavread函数、sound函数等，通过上网和翻阅相关书籍，我了解了它们的基本用法。
关于滤波也是很重要的一部分，有FIR滤波器和IIR滤波器的设计。
通过自主学习，我很快地掌握了编程方法。
在设计过程中，我也遇到了理论知识的忘记、编程过少导致的不熟练、软件运用不自如等问题，但是经过思考、请教同学和老师、查阅相关资料，最终这些问题都得到了解决，最终设计出了符合要求的滤波器并通过相关信号进行了验证。
总的来说，通过这次的课程设计，我对MATLAB这个软件有了更深的认识，也进一步熟悉了MATLAB的使用，感觉到MATLAB软件功能的强大，大大增加了我对MATLAB的兴趣。
经过这次的课程设计，使我发现了自己的一些不足，为了在今后的学习和生活中注意改进自己的不足，我要更加努力，学习好专业知识并注重运用到实际，同时，这次课程设计还提高了自己的动手能力、思维能力和学习的积极性，为我今后的学习和工作都提供了不小的帮助，让我受益匪浅。
7、主要参考文献[1]数字信号处理季秀霞国防工业出版社2013 [2]数字信号处理与应用李正周清华大学出版社2008 [3]数字信号处理教程MATLAB释义与实现陈怀琛电子工业出版社2004. [4]数字信号处理原理及MATLAB实现丛玉良电子工业出版社2009 [5]百度文库附录MATLAB程序 1. 语音采集部分[y,fs,nbits]wavread D\UserData\lenovo\Desktop\cc\cc.wav ; soundy,fs,nbits; N length y ; Yffty,N; subplot2,1,1; ploty;title 原始信号波形; subplot2,1,2; plotabsY; title 原始信号频谱 2. FIR滤波器设计部分Ft10000; Fp2670; Fs3300; wp2*Fp/Ft; ws2*Fs/Ft; rp1; rs40; p1-10.-rp/5; s10.-rs/15; fpts[wp ws]; mag[1 0]; dev[p s]; [n,wn,beta,ftype]kaiserordfpts,mag,dev;bfir1n,wn,Kaisern1,beta; [h,w]freqzb,1; plotw/pi,absh; title FIR 低通滤波器; grid; 3. FIR滤波部分ws2*Fs/Ft; rp1; rs30; p1-10.-rp/5; s10.-rs/15; fpts[wp ws]; mag[1 0]; dev[p s]; [n,wn,beta,ftype]kaiserordfpts,mag,dev;bfir1n,wn,Kaisern1,beta; [h,w]freqzb,1; [y,fs,nbits]wavread D\UserData\lenovo\Desktop\cc\cc.wav ; n length y ; noise0.01*randnn,2; synoise; Sffts,n; zfftfiltb,s; soundz; mfftz; subplot2,2,1; plotabsS, g ; title 滤波前信号的频谱; grid; subplot2,2,2; plotabsm, r ; title 滤波后信号的频谱; grid; subplot2,2,3; plots; title 滤波前信号的波形; grid; subplot2,2,4; plotz; title 滤波后的信号波形; grid; 4．IIR滤波器设计部分Fp2670;Fs3300;As30;Ap1;Ft10000; wp2*pi*Fp/Ft; ws2*pi*Fs/Ft; fp2*Ft*tanwp/2; fs2*Fs*tanwp/2; [n11,wn11]buttordwp,ws,1,50, s ; [b11,a11]buttern11,wn11, s ; [num11,den11]bilinearb11,a11,0.5; [h,w]freqznum11,den11; plotw*8000*0.5/pi,absh; legend 用butter设计; 5.IIR滤波部分Fp2670;Fs3300;As30;Ap1;Ft10000; wp2*pi*Fp/Ft; ws2*pi*Fs/Ft; fp2*Ft*tanwp/2; fs2*Fs*tanwp/2; [n,wn]buttordwp,ws,1,40, s ; [b,a]buttern,wn, s ; [num,den]bilinearb,a,0.5; [h,w]freqznum,den; [y,fs,nbits]wavread D\UserData\lenovo\Desktop\cc\cc.wav ; n length y ; noise0.01*randnn,2; synoise; Sffts,n; zfilternum,den,s; soundz mfftz; subplot2,2,1; plotabsS, g ; title 滤波前信号的频谱; grid; subplot2,2,2; plotabsm, r ; title 滤波后信号的频谱; grid; subplot2,2,3; plots; title 滤波前信号的波形; grid; subplot2,2,4; plotz; title 滤波后的信号波形; grid; 6. 加噪声部分[y, fs, nbits] wavread D\UserData\lenovo\Desktop\cc\cc.wav ; nlengthy; noise0.01*randnn,2; zynoise; soundz; Yffty; Noisefftnoise; subplot2,1,1;ploty;hold on; plotnoise, g ;title 加噪声后波形图; hold off; subplot2,1,2;plotabsY;hold on; plotabsNoise, g ;title 加噪声后频谱;。