数字信号课程设计报告
数字信号处理课程设计报告.doc
数字信号处理课程设计报告淮阴工学院数字信号处理课程设计报告课题名称基于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）画出处理过程中所得的各种波形及频谱图。
数字信号处理课程设计
目录一、摘要.................................................................................................................................... - 2 -1、DSP简介 ................................................................................................................... - 2 -2、系统框图...................................................................................................................... - 3 -二、概述.................................................................................................................................... - 4 -1、系统框图描述...................................................................................................................... - 4 -2、设计目的.............................................................................................................................. - 4 -三、系统设计............................................................................................................................ - 4 -1、总体方案设计 ........................................................................................................... - 4 -2、设计原理.................................................................................................................... - 4 -四、硬件设计............................................................................................................................ - 6 -1、系统硬件框图.............................................................................................................. - 6 -2、 TMS320C5402简介....................................................................................................... - 6 -3、电源设计...................................................................................................................... - 7 -4、复位电路设计.............................................................................................................. - 7 -5、时钟电路设计.............................................................................................................. - 8 -6、 D/A数据转换通道....................................................................................................... - 8 -7、独立键盘设计.............................................................................................................. - 9 -五、软件设计.......................................................................................................................... - 10 -1、正弦波形产生原理.................................................................................................... - 10 -2、设计方案.................................................................................................................... - 10 -3、中断程序流程图........................................................................................................ - 11 -4、正弦信号发生器程序清单........................................................................................ - 12 -六、实验结果.......................................................................................................................... - 18 -1、CCS工程项目的调试 ............................................................................................... - 18 -2、仿真结果...................................................................................................................... - 18 -七、总结.................................................................................................................................. - 21 -八、参考文献.......................................................................................................................... - 21 -一、摘要1、DSP简介数字信号处理(Digital Signal Processing，简称DSP)是一门涉及许多学科而又广泛应用于许多领域的新兴学科。
数字信号处理教程第二版课程设计
数字信号处理教程第二版课程设计1. 项目背景数字信号处理（Digital Signal Processing，简称DSP）是近年来发展最为迅速的学科之一。
在现代通信、控制、音频、视频等领域得到了广泛的应用。
而数字信号处理教程则是DSP学习的入门教材。
本课程设计旨在加深同学们对数字信号处理理论知识的理解，提高同学们的分析及解决数字信号处理问题的能力。
2. 课程内容本课程设计基于数字信号处理教程第二版，其中包括了以下几个方面的内容：2.1 数字信号处理基础本章主要内容包括采样、量化、离散傅里叶变换、数字滤波器设计等基础概念，为后续章节的学习打下基础。
2.2 时域和频域分析时域分析包括了线性时不变系统的时域响应和卷积定理的讲解；频域分析则主要讲解了频谱、功率谱、预测、循环卷积等方面的知识。
2.3 离散傅里叶变换本章主要介绍离散傅里叶变换（Discrete Fourier Transform，简称DFT）的概念及其在数字信号处理领域中的作用，同时还包括了FFT 算法和频域滤波的内容。
2.4 数字滤波器设计及实现本章主要涵盖数字滤波器设计的整个流程，包括了IIR和FIR两类数字滤波器的设计及其在实际应用中的实现。
2.5 DSP处理器及应用本章主要介绍DSP处理器的基本原理和内部构造及其在音频处理、图像处理和通信领域中的应用。
3. 课程要求3.1 课堂讲授老师会通过PPT讲授以上5个章节的内容，讲解完后会留下问答的时间。
同学们应积极思考问题，提出询问，共同讨论解答问题。
3.2 课程设计设计一份数字滤波器，包括其IIR和FIR两个版本，需使用Matlab 或者Python实现。
同学们需要掌握数字滤波器基本概念、对数字信号和滤波器的理解，并能熟练使用Matlab或Python进行数据处理和仿真。
3.3 课程考核课程考核主要分为两部分：•课堂调查问卷满分10分，调查问卷将在整个课程结束后进行，主要考查同学们对数字信号处理知识的掌握和应用能力。
dsp综合设计课程设计报告
dsp综合设计课程设计报告一、教学目标本课程的教学目标是使学生掌握DSP（数字信号处理器）综合设计的基本理论和实践技能。
通过本课程的学习，学生应能够：1.知识目标：理解DSP的基本概念、原理和应用；熟悉DSP芯片的内部结构和编程方法；掌握DSP算法的设计和实现。
2.技能目标：能够使用DSP芯片进行数字信号处理的设计和实现；具备DSP程序的编写和调试能力；能够进行DSP系统的故障诊断和优化。
3.情感态度价值观目标：培养学生对DSP技术的兴趣和热情，提高学生的问题解决能力和创新意识，使学生认识到DSP技术在现代社会中的重要性和应用价值。
二、教学内容本课程的教学内容主要包括DSP的基本理论、DSP芯片的内部结构和工作原理、DSP程序的设计和调试方法、DSP应用系统的设计和实现等。
具体包括以下几个部分：1.DSP的基本概念和原理：数字信号处理的基本概念、算法和特点；DSP芯片的分类和特点。
2.DSP芯片的内部结构：了解DSP芯片的内部结构和工作原理，包括CPU、内存、接口、外设等部分。
3.DSP程序的设计和调试：学习DSP程序的设计方法，包括算法描述、程序编写和调试技巧。
4.DSP应用系统的设计和实现：掌握DSP应用系统的设计方法，包括系统架构、硬件选型、软件开发和系统测试等。
三、教学方法为了实现本课程的教学目标，我们将采用多种教学方法，包括讲授法、案例分析法、实验法等。
具体方法如下：1.讲授法：通过教师的讲解，使学生掌握DSP的基本理论和原理，引导学生理解DSP技术的核心概念。
2.案例分析法：通过分析具体的DSP应用案例，使学生了解DSP技术的实际应用，培养学生的实际操作能力。
3.实验法：通过实验操作，使学生熟悉DSP芯片的使用方法和编程技巧，提高学生的实践能力。
四、教学资源为了支持本课程的教学内容和教学方法的实施，我们将准备以下教学资源：1.教材：选择一本合适的教材，作为学生学习的基础资料，提供系统的DSP知识。
数字信号处理课程设计
数字信号处理课程设计
一、概述
本次信号处理课程设计主要对常见的数字信号处理算法进行实现。
主要内容包括数字信号滤波器、傅立叶变换和数字信号检测算法。
通过实验，学生将学习主要处理手段；同时了解数字信号处理的基本原理和应用。
二、主要内容
（1）数字信号滤波器：实现简单的数字滤波器，同时计算滤波器的频率响应；
（2）傅立叶变换：实现常用的傅立叶变换，并利用变换后的信号图像进行频率分析；
（3）数字信号检测算法：实现基本的一阶和二阶差分算法，并利用此算法进行实时信号检测；
三、实验步骤
（1）准备实验材料：将数字信号的原始信号数据以文件的形式存储，使用MATLAB等软件进行处理；
（2）实现数字滤波器：实现一阶以及多阶低通、高通和带通滤波器，
并计算响应的频谱；
（3）实现傅立叶变换：实现Fourier变换后的信号图像处理，如二维DFT等；
（4）实现数字信号检测算法：实现一阶和二阶差分算法，并利用此算法进行实时信号检测；
（5）数字信号处理综合应用实验：针对实际的数字信号，分析信号的特征，并基于实验结果进行信号处理算法的比较。
四、实验结果
完成本次实验后，可以实现对不同数字信号的处理，掌握其中滤波器、傅立叶变换等数字信号处理理论，并掌握常规的算法，学会运用算法实现实际信号处理工程。
数字信号处理课程设计
index=1; while n<=N IND(n)=index; index=index+k; index=mod(index, 2^Bits); n=n+1; end m=IND; IND=IND+1; y=LUT(IND); subplot(212) plot(t,y) figure(1); axis([0 endtime-1.2 1.2]); grid on 3、测试结果 ⑴输入 dds(10,100,10,0.5)，输出波形图如图 2
结论：实验测得输出频率为 10Hz,与理论值相符，实验正确。
题目二
数字锁相环的 MATLAB 实现 及仿真
一、 设计目的
1、理解数字锁相环的原理，熟悉实现流程。 2、学会用 MATlab 进行一些信号分析
二、 数字锁相环原理
1、数字锁相环简介 锁相环 (phase-locked loop)为无线电发射中使频率较为稳定的一种方法， 主要有 VCO（压控振荡器）和 PLL IC （锁相环集成电路）,压控振荡器给出一个 信号,一部分作为输出，另一部分通过分频与 PLL IC 所产生的本振信号作相位比 较，为了保持频率不变，就要求相位差不发生改变，如果有相位差的变化，则 PLL IC 的电压输出端的电压发生变化，去控制 VCO,直到相位差恢复，达到锁频 的目的， 是一种能使受控振荡器的频率和相位均与输入信号保持确定关系的闭环 电子电路。 锁相环的用途是在收、 发通信双方建立载波同步或位同步，分为模拟锁相环 和数字锁相环两种。 随着数字电路技术的发展，数字锁相环以其独有的优势在调 制解调、频率合成、FM 立体声解码、彩色副载波同步、图象处理等各个方面得 到了广泛的应用。 数字锁相环不仅吸收了数字电路可靠性高、体积小、价格低 等优点， 还解决了模拟锁相环的直流零点漂移、器件饱和及易受电源和环境温度 变化等缺点， 此外还具有对离散样值的实时处理能力，已成为锁相技术发展的方 向。 2、DPLL 基本模型和原理 全数字锁相环包括数字鉴相鉴频器（PFD） 、数字滤波器（LPF） 、数字振荡 器（NCO）三部分，如下图所示：
中南大学数字信号处理课程设计
中南大学数字信号处理课程设计报告专业班级: 电信1303指导老师：姓名:学号:目录一、课程设计要求二、设计过程（1）设计题目（2）设计源代码（3）设计结果（4）结果分析三、设计总结与心得体会四、课程设计指导书一、课程设计要求1、课程设计指导书①《数字信号处理（第二版）》，丁玉美等，西安电子科技大学出版社；②《MATLAB 及在电子信息课程中的应用》，陈怀琛等，电子工业出版社。
2、课程设计题目⑴、信号发生器用户根据测试需要，可任选以下两种方式之一生成测试信号：①、直接输入（或从文件读取）测试序列；②、输入由多个不同频率正弦信号叠加组合而成的模拟信号公式（如式1-1 所示）、采样频率（Hz）、采样点数，动态生成该信号的采样序列，作为测试信号。
⑵、频谱分析使用FFT 对产生的测试信号进行频谱分析并展示其幅频特性与相频特性，指定需要滤除的频带，通过选择滤波器类型（IIR / FIR），确定对应的滤波器（低通、高通）技术指标。
⑶、滤波器设计根据以上技术指标（通带截止频率、通带最大衰减、阻带截止频率、阻带最小衰减），设计数字滤波器，生成相应的滤波器系数，并画出对应的滤波器幅频特性与相频特性。
①IIR DF 设计：可选择滤波器基型（巴特沃斯或切比雪夫型）；②FIR DF 设计：使用窗口法（可选择窗口类型，并比较分析基于不同窗口、不同阶数所设计数字滤波器的特点）。
⑷、数字滤波根据设计的滤波器系数，对测试信号进行数字滤波，展示滤波后信号的幅频特性与相频特性，分析是否满足滤波要求（对同一滤波要求，对比分析各类滤波器的差异）。
①IIR DF：要求通过差分方程迭代实现滤波（未知初值置零处理）；②FIR DF：要求通过快速卷积实现滤波（对于长序列，可以选择使用重叠相加或重叠保留法进行卷积运算）。
⑸、选做内容将一段语音作为测试信号，通过频谱展示和语音播放，对比分析滤波前后语音信号的变化，进一步加深对数字信号处理的理解。
3、具体要求⑴、使用MATLAB（或其它开发工具）编程实现上述内容，写出课程设计报告。
fundamentalsofdigitalsignalprocessing课程设计
Fundamentals of Digital Signal Processing 课程设计一、课程设计的目的和意义数字信号处理是电子信息专业的重要课程之一，是掌握现代信号处理技术的必修课程。
本次课程设计旨在巩固和深化同学们对于数字信号处理的理解和实践能力。
通过课程设计，可以帮助学生更好地理解数字信号处理的基本概念、原理和方法，提高学生的实际操作能力，培养工程实践能力和团队协作能力，提高学生的综合素质。
二、课程设计的主要内容1. 实验器材1.电脑一台2.程序集成开发环境Keil uVision53.万用表和示波器一台2. 实验内容1.基本信号的时间和频域分析目的：了解基本信号的结构与特性，建立时间域和频域分析方法，为后续信号处理打下基础。
步骤：（1）构造三种基本信号：方波，三角波和正弦波。
（2）利用示波器测量信号的时间轴，（3）利用万用表测量信号的电压幅值，（4）在Keil uVision5中计算出信号的频谱分布。
2.数字滤波器去噪实验目的：利用数字滤波器对带有噪声的信号进行去噪。
步骤：（1）生成1000Hz正弦波作为原始信号。
（2）在信号中加入高斯白噪声进行干扰。
（3）利用FIR滤波器去除干扰后输出滤波后的信号。
（4）测量滤波前、滤波后的信号电压幅值，并对滤波前、滤波后信号的频谱分布进行比较。
3.信号变换实验目的：理解傅里叶变换和其在实际信号分析中的应用。
步骤：（1）利用Keil uVision5中的FFT工具对1000Hz正弦波进行傅里叶变换。
（2）观察输出结果并解释其意义。
（3）对另外两种基本波形进行傅里叶变换，并比较其与正弦波的差异。
（4）以实际信号为例，进行复杂信号的频域分析。
4.数字信号编码与解码实验目的：了解数字信号编码与解码的基本原理和方法。
步骤：（1）将一个信号利用模拟调制方式（例如AM、FM、PM等）进行调制，并输出模拟信号。
（2）利用AD转换器将模拟信号转换成数字信号。
（3）将数字信号解调还原成模拟信号。
数字信号处理课程设计
-40 -60 -80 -100 -120 -140
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
等波纹滤波器法设计FIR数字滤波器
Matlab应用的函数
[M,fo,ao,w]=remezord(f,a,dev) 与remez 配合使用，用于算出适合要求的滤波器阶次M ，fo和ao为有 2B个元素的向量，fo是频率轴分点，ao是在这些频率分点上理想的幅 频响应，w是有B个元素的向量，表示各频带的加权值。由于 remezord给出的阶次M有可能偏低，这时需要适当增加M。 f是具有2B个元素的向量，由通带边缘频率和阻带边缘频率构成，dev 若是两个值则为高通或低通，三个值为带通或带阻，由通带和阻带上 的偏差决定 Remez 设计出等波纹滤波器 freqz(B,A,N)：离散系统频响特性 B和A分别为离散系统的系统函数分 子、分母多项式的系数向量，N为正整数 (频率等分点的值 )
60
70
0.5
0
-0.5
0
10
20
30
40
50
60
70
小结
（1）在时域求系统响应的方法有两种,第一种是通过解差分 方程求得系统输出;第二种是已知系统的单位脉冲响应,通 过求输入信号和系统单位脉冲响应的线性卷积求得系统输 出。 （2）检验系统的稳定性，其方法是在输入端加入单位阶跃 序列， 观察输出波形,如果波形稳定在一个常数值上,系 统稳定,否则不稳定。
数字信号处理课程设计csdn
数字信号处理课程设计csdn一、教学目标本课程的教学目标旨在帮助学生掌握数字信号处理的基本概念、原理和方法，培养学生运用数字信号处理技术解决实际问题的能力。
具体目标如下：1.知识目标：–掌握数字信号处理的基本原理，包括采样、量化、数字滤波、傅里叶变换等；–了解数字信号处理的应用领域，如通信、音视频处理、图像处理等；–熟悉常用的数字信号处理算法和软件工具。
2.技能目标：–能够运用数字信号处理理论分析和解决实际问题；–具备使用Matlab等软件进行数字信号处理的能力；–学会阅读和分析数字信号处理的文献和资料。
3.情感态度价值观目标：–培养学生的团队合作精神和自主学习能力；–激发学生对数字信号处理的兴趣和热情，提高学生的专业素养。
二、教学内容本课程的教学内容主要包括以下几个部分：1.数字信号处理基本概念：采样、量化、数字信号表示和运算；2.数字滤波器：滤波器设计方法、常见滤波器类型及应用；3.离散傅里叶变换：DFT原理、快速傅里叶变换（FFT）算法及应用；4.数字信号处理算法和软件工具：Matlab及其在数字信号处理中的应用；5.数字信号处理应用案例：通信系统、音视频处理、图像处理等。
三、教学方法为了提高教学效果，本课程将采用多种教学方法相结合的方式，包括：1.讲授法：用于传授基本概念、原理和方法；2.案例分析法：通过实际案例，让学生更好地理解数字信号处理的应用；3.实验法：利用Matlab软件，让学生动手实践，培养实际操作能力；4.小组讨论法：鼓励学生分组讨论，培养团队合作精神和自主学习能力。
四、教学资源为了支持本课程的教学，我们将准备以下教学资源：1.教材：《数字信号处理》（或其他主流教材）；2.参考书：提供相关的学术文献和资料，供学生拓展阅读；3.多媒体资料：制作课件、视频等，辅助学生理解和掌握知识；4.实验设备：计算机、Matlab软件等，供学生进行实验和实践。
五、教学评估本课程的教学评估将采用多元化的评价方式，全面客观地评价学生的学习成果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、语音信号去噪处理1.设计要求：（1）在windows系统下的录音机录制一段1s左右的语音信号作为原声信号，在MATLAB软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数；（2）画出语音信号的时域波形，对采样后的语音进行fft变换，得到信号的频谱特性；对语音信号分别加入正弦噪声和白噪声，画出加噪信号的时域波形和频谱图；（3）根据对加噪语音信号谱分析结果，确定滤除噪声滤波器的技术指标，设计合适的数字滤波器，并画出滤波器的频域响应；（4）用所设计的滤波器对加噪的信号进行滤波，在同一个窗口画出滤波前后信号的时域图和频谱图，对滤波前后的信号进行对比，分析信号变化；（5）利用sound(x)回放语音信号，验证设计效果。
2.设计步骤：（1）找到7s的语音信号，利用函数wavread对语音信号进行信号读取；（2）计算样本时刻和频谱图的频率，并进行N+1点FFT变换；（3）加噪声为5000Hz的正弦信号正弦噪声，采用awgn函数加信噪比为10的高斯白噪声；（4）设计滤波器；（5）绘出相应的时域、频域图；（6）利用sound函数进行原始信号的语音播放，加噪声音播放，以及滤波之后的语言播放。
3.设计实现：（1）时域图与频谱图（加正弦）录入原始信号的时域图：加入正弦信号后的时域图：滤波后的时域图：录入原始信号的频域图：加入正弦信号后的频率图：滤波后的频域图：采用巴斯低通滤波器滤除正弦波：（2）具体代码实现：[x,fs,bits]=wavread('E:\mcpass.wav');%原信号n=size(x,1); %提取采样信号的长度t=(0:length(x)-1)/fs; %计算样本时刻f=fs*(0:(n+1)/2-1)/n+1; %计算频域图的频率X=fft(x,n+1); %进行N+1点FFT变换ts=0:1/fs:(size(x)-1)/fs; %将所加噪声信号的点数调整到与原始信号相同s=x+0.05*sin(2*pi*5000*ts)'; %加噪声为5000Hz的正弦信号正弦噪声S=fft(s,n+1); %加正弦噪声后的频域%正弦滤波wp=2000/fs*2*pi; %2000为通带截止频率ws=3000/fs*2*pi; %3000为阻带下限截止频率Rp=4; %通带波纹Rs=25; %阻带波纹T=1/fs;Fs=1/T; %定义采样间隔Wp=2/T*tan(wp/2); %计算对应的数字频率Ws=2/T*tan(ws/2);[N,wn]=buttord(Wp,Ws,Rp,Rs,'s'); %计算滤波器介数和截止频率[c,d]=butter(N,wn,'s'); %计算滤波器系统函数分子分母系数[B,A]=bilinear(c,d,Fs); %双线性变换得到数字滤波器系统函数分子分母系数[Hb,Wc]=freqz(B,A);sf=filter(B,A,s); %对加噪信号进行滤波Sf=fft(sf,n+1); %对滤波后进行N+1点FFT变换%绘图部分figure(3);plot(fs*Wc/(2*pi),20*log10(abs(Hb)));title('巴斯低通滤波器频域响应图');xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(3,1,1);plot(t,x);title('原信号时域')xlabel('时间(s)');ylabel('幅度');figure(2);subplot(3,1,1);plot(f,abs(X(1:(n+1)/2)));title('原信号频域')xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(3,1,2);plot(t,s);title('加正弦信号后的时域')xlabel('时间(s)');ylabel('幅度');figure(2);subplot(3,1,2);plot(f,abs(S(1:(n+1)/2)));title('加正弦信号后的频域图') xlabel('频率(Hz)');ylabel('幅度');figure(1)subplot(3,1,3);plot(t,sf);title('滤波后的时域图');xlabel('时间(s)');ylabel('幅度');figure(2)subplot(3,1,3);plot(f,abs(Sf(1:(n+1)/2))); title('滤波后的频域图'); xlabel('频率(Hz)');ylabel('幅度');sound(x);sound(s);sound(sf);（3）时域图与频域图（加白噪声）加白噪声后的时域图和滤除之后的时域图：加白噪声和滤除之后的频域图：采用blackman函数滤波：具体代码实现：[x,fs,bits]=wavread('E:\hbsong.wav');N=size(x,1);t=(0:length(x)-1)/fs;f=fs*(0:(N+1)/2-1)/N+1;X=fft(x,N+1);%加高斯白噪声z=awgn(x,20); %对信号加信噪比为10的高斯白噪声N1=size(z,1); %提取采样信号的长度t=(0:length(z)-1)/fs; %计算样本时刻f=fs*(0:(N1+1)/2-1)/N1+1;Z=fft(z,N1+1);Wp=2500/fs*2*pi;Ws=3000/fs*2*pi; %计算对应的数字频率B=Ws-Wp;n=ceil(1*pi/B);wc=(Wp+Ws)/2;b=fir1(n-1,wc/pi,'stop',blackman(n));%blackman窗函数滤波[H,w]=freqz(b,1);y=fftfilt(b,z);t1=(0:length(y)-1)/fs;Y=fft(y,N1+1);subplot(2,2,1);plot(t,z);title('加高斯白噪声后时域图');subplot(2,2,2);plot(f,abs(Z(1:(N1+1)/2)));title('滤波前信号频谱图')figureplot(fs*w/(2*pi),20*log10);title('blackman函数频域响应图');title('频率响应')xlabel('频率（Hz）');ylabel('幅度');subplot(2,2,3);plot(t1,y)title('滤波后信号时域图');xlabel('时间(s)');ylabel('幅度');subplot(2,2,4);plot(f,abs(Y(1:(N1+1)/2)));title('滤波后信号频谱图')xlabel('频率（Hz）');ylabel('幅度');sound(y,fs)二、语音信号的延时和混响1.设计要求：（1）利用Windows下的录音机或其他软件，录制一段自己的语音信号，时间控制在1s左右，并对录制的信号进行采样；（2）语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图；将信号加入延时和混响，再分析其频谱，并与原始信号频谱进行比较；（3）设计几种特殊类型的滤波器：单回声滤波器，多重回声滤波器，全通结构的混响器，并画出滤波器的频域响应；（4）用自己设计的滤波器对采集的语音信号进行滤波；（5）分析得到信号的频谱，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；（6）回放语音信号。
2. 设计步骤：（1）录入原始声音信号；（2）计算样本时刻和频谱图的频率，并进行N+1点FFT变换；（3）加入单回声；（4）设计单回声滤波器幅频响应函数；（5）绘出相应时域与频域图；（6）利用sound函数进行原始信号的语音播放，加单回声后语言播放，以及滤除之后的语言播放；3.设计实现：（1）时域图和频域图（加单回声）原始信号时域图：滤波后的时域图：原始信号频域图：滤波后频域图：单回声滤波器幅频响应：（2）具体代码实现（单回声）：[x,fs,bits]=wavread('E:\mcpass.wav');%原声音信号n=size(x,1);t=(0:length(x)-1)/fs;f=fs*(0:(n+1)/2-1)/n+1;X=fft(x,n+1);a=0.6;%单回声滤波R=fs*a;B=[1,zeros(1,R-2),a];A=[1,zeros(1,R-1)];y = filter(B,A,x);Y=fft(y,n+1);[H,W]=freqz(B,A); %求单回声滤波器幅频响应函数%绘图部分figure(2);plot(fs*W/(2*pi),20*log10(abs(H)));title('频率响应');xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(4,1,1);plot(t,x);title('原信号时域');xlabel('时间(s)');ylabel('幅度');subplot(4,1,2);plot(f,abs(X(1:(n+1)/2)));title('原信号频域');xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(4,1,3);plot(t,y);title('滤波后时域图');xlabel('时间(s)');ylabel('幅度');figure(1);subplot(4,1,4);plot(f,abs(Y(1:(n+1)/2)));title('滤波后频域图');xlabel('频率(Hz)');ylabel('幅度');sound(x,fs,bits);sound(y,fs,bits);（3）时域图与频域图（多重回声）原始信号时域图与多重信号时域图：原始信号时域图与多重信号频域图：多重回声滤波器幅频响应：（4）具体代码实现（多回声）%原声音信号[x,fs,bits]=wavread('E:\mcpass.wav');N=size(x,1);t=(0:length(x)-1)/fs;f=fs*(0:(N+1)/2-1)/N+1;%多重回声a=0.5;R=fs*a;B=[1,zeros(1,R-2)];A=[1,zeros(1,R-1),a];yd = filter(B,A,x); %滤波器函数%频率响应[Hb,Wc]=freqz(B,A); %求多重回声滤波器幅频响应函数%绘图部分figure(2);plot(fs*Wc/(2*pi),20*log10(abs(Hb)));title('多重滤波器频域响应图')xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(4,1,1);plot(t,x);title('原信号时域');xlabel('时间(s)');ylabel('幅度');X=fft(x,N+1);subplot(4,1,3);plot(f,abs(X(1:(N+1)/2)));title('原信号频域');xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(4,1,2);plot(t,yd);title('多重信号时域图')xlabel('时间(s)');ylabel('幅度');Yd=fft(yd,N+1); %多重回声滤波后信号FFT变换subplot(4,1,4);plot(f,abs(Yd(1:(N+1)/2)));title('多重信号频域图');xlabel('频率(Hz)');ylabel('幅度');sound(x);sound(yd);（5）时域图与频域图（全通滤波器）原始信号时域图与全通滤波器时域图：原始信号时域图与全通滤波器频域图：全通滤波器的频域响应图：（6）具体代码实现（全通滤波器）：%原声音信号[x,fs,bits]=wavread('E:\mcpass.wav');N=size(x,1);t=(0:length(x)-1)/fs;f=fs*(0:(N+1)/2-1)/N+1;%全通结构混响a=0.5;R=fs*a;B=[a,zeros(1,R-2),1];A=[1,zeros(1,R-2),a];yd = filter(B,A,x);%频率响应[Hb,Wc]=freqz(B,A); %求全通滤波器幅频响应函数figure(2);plot(fs*Wc/(2*pi),20*log10(abs(Hb)));title('全通结构频域响应图') xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(4,1,1);plot(t,x);title('原信号时域');xlabel('时间(s)');ylabel('幅度');X=fft(x,N+1);subplot(4,1,3);plot(f,abs(X(1:(N+1)/2)));title('原信号频域');xlabel('频率(Hz)');ylabel('幅度');figure(1);subplot(4,1,2);plot(t,yd);title('全通滤波器时域图');xlabel('时间(s)');ylabel('幅度');Yd=fft(yd,N+1);subplot(4,1,4);plot(f,abs(Yd(1:(N+1)/2)));title('全通滤波器频域图');xlabel('频率(Hz)');ylabel('幅度');sound(x);sound(yd);（7）时域图与频域图（延时与混响）原始信号时域图：延时时域图：混响时域图：原始信号频域：延时信号频域图：混响信号频域图：（8）具体代码实现（延时与混响）：[x,fs,bits]=wavread('E:\mcpass.wav');N=size(x,1);t=(0:length(x)-1)/fs;f=fs*(0:(N+1)/2-1)/N+1;xd=[zeros(300,1);x]; %信号延时td=(0:length(xd)-1)/fs;Xd=fft(xd);Xd1=fftshift(Xd);dFs =fs/length(xd);xd1=[x;zeros(300,1)];%信号混响x1=xd1+xd;Xhun=fft(x1);%混响信号fft变换Xh2=fftshift(x1);%平移，中心为0频率dFs = fs/length(x1);%绘图部分figure(1);subplot(4,1,1);plot(t,x);title('原信号时域');xlabel('时间(s)');ylabel('幅度');X=fft(x,N+1);subplot(4,1,2);plot(f,abs(X(1:(N+1)/2)));title('原信号频域');xlabel('频率(Hz)');ylabel('幅度');subplot(4,1,3);plot(td,xd);title('延时信号时域图')xlabel('时间(s)');ylabel('幅度');subplot(4,1,4);plot([-fs/2:dFs: fs/2-dFs],abs(Xd)); title('延时信号频域图')xlabel('频率(Hz)');ylabel('幅度');figure(2);subplot(2,1,1);plot(td,x1);title('混响信号时域图');xlabel('时间(s)');ylabel('幅度');figure(2);subplot(2,1,2);plot([-fs/2:dFs: fs/2-dFs],abs(Xh2)); title('混响信号频域图');xlabel('频率(Hz)');ylabel('幅度');sound(x);sound(xd);sound(x1);三、数字滤波器的设计及实现1.设计要求：（1）调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，观察st的时域波形和幅频特性曲线；（2）要求将st中的三路调幅信号分离，通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率，要求滤波器的通带最大衰减为0.1dB，阻带最小衰减为60dB；（3）编程调用MATLAB滤波器设计函数分别设计这三个数字滤波器，并绘图显示其幅频特性曲线；（4）调用滤波函数filter，用所设计的三个滤波器分别对复合信号st进行滤波，分离出st中的三路不同载波频率的调幅信号，并绘图显示滤波后信号的时域波形和频谱，观察分离效果。