信号分析与处理综合设计实践指导书

《信号分析与处理综合设计实践》指导

一、课程设计目的：

综合运用本课程的理论知识进行频谱分析以及滤波器设计，通过理论推导得出相应结论，并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。

二、课程设计时间安排：

第一天：布置设计任务，讲解设计要求，提示设计要点。

第二～四天：查阅资料（在图书馆或上网），弄清题目要求，提出解决方案。) 第五～十一天：根据题目要求，将理论推导与编程实现相结合，写出设计报告。

第十二天：答辩。

三、课程设计内容：

1. 设计题目一: 语音信号的处理与滤波（难度系数：0.8）

1．熟悉并掌握MATLAB中有关声音（wave）录制、播放、存储和读取的函数。

2．在MATLAB环境中，使用声音相关函数录制2秒左右自己的声音，抽样率是8000Hz/s。（考虑如何解决一个实际问题：录制刚开始时，常会出现实际发出声音落后录制动作半拍的现象，如何排除对这些无效点的采样？）3．分别取8000个和16000个数据进行频谱分析，得到幅度和相位谱，比较二者异同并分析原因。

4．针对电话信道（最高3500Hz），设计一个FIR或IIR滤波器进行滤波，把抽样率转变为7000Hz/s，并进行频谱分析，得到幅度和相位谱。

5．把处理后的所有数据储存为声音文件，与原始声音进行比较。

2. 设计题目二：编程实现任意确定信号的频谱分析算法。（难度系数：0.9）

（1）对给定的CEG和弦音音频文件取合适长度的采样记录点，然后进行频谱分析（信号的时域及幅频特性曲线要画出）。

（2）分析CEG和弦音频谱特点，对该信号频谱能量相对较为集中的频带（分低、中、高频）实现滤波（分别使用低通，带通及高通），显示滤波后信号的时域和频域曲线，并对滤波后的信号与原信号的音频进行声音回放比较。

（3）在低、中、高三个频带中，各滤出三个能量最集中的频簇，显示滤波后信号的时域和频域曲线。

（4）任意选择几个滤出的频带（或频簇）进行时域信号重建（合成），与原信号的音频进行声音回放比较。

讨论：根据上述结果，分析什么是和弦音。

3. 设计题目三：设计并实现FIR数字滤波器的窗口设计算法（难度系数：0.7）

要求：输入数字滤波器指标，包括滤波器类型（低通、带通、带阻或高通等），通带截止频率，通带最大波动，阻带开始频率，阻带衰减，设计得到FIR滤波器，并画出设计得到的滤波器的增益曲线图（要有坐标标度）。为了使编制的程序操作方便，设计处理系统的用户界面：在所设计的系统界面上可以选择滤波器的类型，输入滤波器的参数，显示滤波器的频率响应，选择输入信号并显示相应的输出信号等。

讨论：对多种窗口、不同阶数设计结果进行比较。

4. 设计题目四：语音信号的处理与滤波（难度系数：1）

内容：录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标，采用窗函数法和双线性变换法设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；回放语音信号；换一个与你性别相异的人录制同样一段语音内容，分析两段内容相同的语音信号频谱之间有什么特点；再录制一段同样长时间的背景噪声叠加到你的语音信号中，分析叠加前后信号频谱的变化，设计一个合适的滤波器，能够把该噪声滤除；最后，设计一个信号处理系统界面。

基本步骤：

1．语音信号的采集

熟悉并掌握MATLAB中有关声音（wave）录制、播放、存储和读取的函数，在MATLAB环境中，使用声音相关函数录制3秒左右自己的声音。（考虑如何解决一个实际问题：录制刚开始时，常会出现实际发出声音落后录制动作半拍的现象，如何排除对这些无效点的采样？）

2．语音信号的频谱分析

要求首先画出语音信号的时域波形；然后对语音号进行快速傅里叶变换，得到信号的频谱特性，从而加深对频谱特性的理解。

3．设计数字滤波器并画出其频率响应

给出各滤波器的性能指标：

(1)低通滤波器性能指标：f

b ＝1000 Hz，f

c

＝1200 Hz，A

s

＝100dB，A

p

＝1dB。

(2)高通滤波器性能指标：f

c ＝4800 Hz，f

b

＝5000 Hz A

s

＝100dB，A

p

＝1dB。

(3)带通滤波器性能指标：f

p1＝1200 Hz，f

p2

＝3000 Hz，f

sc1

＝1000 Hz，f

sc2

＝3200Hz，A

s ＝100dB，A

p

＝1dB。

要求用窗函数法和双线性变换法设计上面要求的3种滤波器。在Matlab中，使用函数butter，cheby1等函数设计IIR滤波器；利用Matlab中的函数freqz 画出各滤波器的频率响应。

4．用滤波器对信号进行滤波

要求用自己设计的各种滤波器分别对采集的信号进行滤波，在Matlab中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。

5．比较滤波前后语音信号的波形及频谱

要求在一个窗口同时画出滤波前后的波形及频谱。

6．回放语音信号

在Matlab中，使用相关声音函数对声音进行回放。辨别滤波前后的声音有何变化。

7．男女声语音信号频谱特点分析

换一个与你性别相异的人录制同样一段语音内容，分析两段内容相同的语音

信号频谱之间有什么特点。

8. 背景噪声滤除

再录制一段同样长时间的背景噪声（频谱能量集中在某个小范围内）叠加到你的语音信号中，分析叠加前后信号频谱的变化，设计一个合适的滤波器，能够把该噪声滤除。

9．设计系统界面（用MATLAB中GUIDE工具，要求自学）

为了使编制的程序操作方便，要求有能力的同学，设计处理系统的用户界面。在所设计的系统界面上可以选择滤波器的类型，输入滤波器的参数，显示滤波器的频率响应，选择信号等。

四、课程设计要求：

1．设计说明书（约3000～4000字），包括：

（1）封面；

（2）目录；

（3）摘要；

（4）正文。

①设计内容（简述本设计的任务和要求，可参照任务书和指导书）；

②设计原理（简述设计过程中涉及到的基本理论知识）；

③设计过程（按设计步骤详细介绍设计过程，即任务书和指导书中指定

的各项任务）

I．程序源代码：给出完整源程序清单；

II．调试分析过程描述：包括测试数据、测试输出结果，以及对程序调试过程中存在问题的思考（列出主要问题的出错现象、出错

原因、解决方法及效果等）；

III．结果分析：对程序结果进行分析，并与理论分析进行比较。

（5）结论：包括课程设计过程中的学习体会与收获、对Matlab语言和本次课程设计的认识以及自己的建议等内容；

（6）参考文献。

2．附件：（可以将设计中得出的波形图和频谱图作为附件，在说明书中涉及

相应图形时，注明相应图形在附件中位置即可；也可不要附件，所有内容全部包含在设计说明书中。所有的实验结果图形都必须有横纵坐标标注，必须有图序和图题。）

3．课程设计报告提交电子版和打印版各一份。

注意：

鼓励互相讨论交流，但不得抄袭他人的报告（或给他人抄袭），一旦发现，抄袭者和被抄袭者成绩均为零分！！！

语音信号处理实验指导书

语音信号处理实验指导书 实验一 语音信号采集与简单处理 一、 实验目的、要求 （1）掌握语音信号采集的方法 （2）掌握一种语音信号基音周期提取方法 （3）掌握短时过零率计算方法 （4）了解Matlab 的编程方法 二、 实验原理 基本概念： （a ）短时过零率： 短时内，信号跨越横轴的情况，对于连续信号，观察语音时域波形通过横轴的情况；对于离散信号，相邻的采样值具有不同的代数符号，也就是样点改变符号的次数。 对于语音信号，是宽带非平稳信号，应考察其短时平均过零率。 其中sgn[.]为符号函数 ?? ?? ?<=>=0 x(n)-1sgn(x(n))0 x(n)1sgn(x(n)) 短时平均过零的作用 1.区分清/浊音： 浊音平均过零率低，集中在低频端； 清音平均过零率高，集中在高频端。 2.从背景噪声中找出是否有语音，以及语音的起点。 （b ）基音周期 基音是发浊音时声带震动所引起的周期性，而基音周期是指声带震动频率的倒数。基音周期是语音信号的重要的参数之一，它描述语音激励源的一个重要特征，基音周期信息在多个领域有着广泛的应用，如语音识别、说话人识别、语音分析与综合以及低码率语音编码，发音系统疾病诊断、听觉残障者的语音指导等。因为汉语是一种有调语言，基音的变化模式称为声调，它携带着非常重要的具有辨意作用的信息，有区别意义的功能，所以，基音的提取和估计对汉语更是一个十分重要的问题。 ∑--= -=1 )]1(sgn[)](sgn[21N m n n n m x m x Z

由于人的声道的易变性及其声道持征的因人而异，而基音周期的范围又很宽，而同—个人在不同情态下发音的基音周期也不同，加之基音周期还受到单词发音音调的影响，因而基音周期的精确检测实际上是一件比较困难的事情。基音提取的主要困难反映在：①声门激励信号并不是一个完全周期的序列，在语音的头、尾部并不具有声带振动那样的周期性，有些清音和浊音的过渡帧是很难准确地判断是周期性还是非周期性的。②声道共振峰有时会严重影响激励信号的谐波结构，所以，从语音信号中直接取出仅和声带振动有关的激励信号的信息并不容 易。③语音信号本身是准周期性的(即音调是有变化的)，而且其波形的峰值点或过零点受共振峰的结构、噪声等的影响。④基音周期变化范围大，从老年男性的50Hz 到儿童和女性的450Hz ，接近三个倍频程，给基音检测带来了一定的困难。由于这些困难，所以迄今为止尚未找到一个完善的方法可以对于各类人群(包括男、女、儿童及不向语种)、各类应用领域和各种环境条件情况下都能获得满意的检测结果。 尽管基音检测有许多困难，但因为它的重要性，基音的检测提取一直是一个研究的课题，为此提出了各种各样的基音检测算法，如自相关函数(ACF)法、峰值提取算法(PPA)、平均幅度差函数(AMDF)法、并行处理技术、倒谱法、SIFT 、谱图法、小波法等等。 三、使用仪器、材料 微机（带声卡）、耳机，话筒。 四、 实验步骤 （1）语音信号的采集 利用Windows 语音采集工具采集语音信号，将数据保存wav 格式。 采集一组浊音信号和一组清音信号，信号的长度大于3s 。 （2）采用短时相关函数计算语音信号浊音基音周期，考虑窗长度对基音周期计算的影响。采用倒谱法求语音信号基音周期。 （3）计算短时过零率，清音和浊音的短时过零率有何区别。 五、实验过程原始记录（数据，图表，计算） 短时过零率 短时相关函数 P j j n s n s j R N j n n n n ,,1) ()()(1 =-=∑-= ∑--=-=10 )]1(sgn[)](sgn[21N m n n n m x m x Z

数字信号处理课程设计报告

抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于2倍信号最高频率。抽样频率小于2倍频谱最高频率时，信号的频谱有混叠。抽样频率大于2倍频谱最高频率时，信号的频谱无混叠。 语音信号处理是研究用数字信号处理技术和语音学知识对语音 信号进行处理的新兴学科，是目前发展最为迅速的学科之一，通过语音传递信息是人类最重要，最有效，最常用和最方便的交换信息手段，所以对其的研究更显得尤为重要。 Matlab语言是一种数据分析和处理功能十分强大的计算机应用 软件，它可以将声音文件变换成离散的数据文件，然后用起强大的矩阵运算能力处理数据。这为我们的本次设计提供了强大并良好的环境！ 本设计要求通过利用matlab对模拟信号和语音信号进行抽样，通过傅里叶变换转换到频域，观察波形并进行分析。 关键词：抽样Matlab

目录 一、设计目的： (2) 二、设计原理： (2) 1、抽样定理 (2) 2、MATLAB简介 (2) 3、语音信号 (3) 4、Stem函数绘图 (3) 三、设计内容： (4) 1、已知g1(t)=cos(6πt),g2(t)=cos(14πt),g3(t)=cos(26πt),以抽样频率 fsam=10Hz对上述三个信号进行抽样。在同一张图上画出g1(t)，g2(t)，g3(t)及其抽样点,对所得结果进行讨论。 (4) 2、选取三段不同的语音信号，并选取适合的同一抽样频率对其进 行抽样，画出抽样前后的图形，并进行比较，播放抽样前后的语音。 (6) 3、选取合适的点数，对抽样后的三段语音信号分别做DFT，画图 并比较。 (10) 四、总结 (12) 五、参考文献 (13)

五款信号完整性仿真工具介绍

现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 (一)Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，Ansoft公司的仿真工具能够从三维场求解的角度出发，对PCB设计的信号完整性问题进行动态仿真。 Ansoft的信号完整性工具采用一个仿真可解决全部设计问题： SIwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何数量的过孔和信号引线条构成。仿真结果采用先进的3D图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿真器。 (二)SPECCTRAQuest Cadence的工具采用Sun的电源层分析模块： Cadence Design Systems的SpecctraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在PCB详细设计前使用此工具，对互连线的不同情况进行仿真，把仿真结果存为拓扑结构模板，在后期详细设计中应用这些模板进行设计。 （2）DF/Signoise工具是信号仿真分析工具，可提供复杂的信号延时和信号畸变分析、IBIS 模型库的设置开发功能。SigNoise是SPECCTRAQUEST SI Expert和SQ Signal Explorer Expert进行分析仿真的仿真引擎，利用SigNoise可以进行反射、串扰、SSN、EMI、源同步及系统级的仿真。 （3）DF/EMC工具——EMC分析控制工具。 （4）DF/Thermax——热分析控制工具。 SPECCTRAQuest中的理想高速PCB设计流程： 由上所示，通过模型的验证、预布局布线的space分析、通过floorplan制定拓朴规则、由规

实验五 数字信号处理综合设计

实验五数字信号处理综合设计 一、实验目的 1．掌握在Windows环境下语音信号采集的方法； 2．掌握MATLAB设计FIR和IIR数字滤波器的方法； 3．学会用MATLAB 对信号进行分析和处理。 二、实验内容 1．语音信号的采集 要求利用windows下的录音机或其他软件，录制一段自己的话音，时间定为10秒。然后在MATLAB软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。通过wavread函数的使用，要求理解采样频率、采样位数等概念。wavread函数调用格式： y=wavread(file)，读取file所规定的wav文件，返回采样值放在向量y中。 [y,fs,nbits]=wavread(file)，采样值放在向量y 中，fs 表示采样频率（Hz），nbits 表示采样位数。 y=wavread(file,N)，读取前N 点的采样值放在向量y中。 y=wavread(file,[N1,N2])，读取从N1点到N2点的采样值放在向量y中。 例如：x1=wavread(h:\课程设计2\shuzi.wav); %读取语音信号的数据，赋给变量x1 2．语音信号的频谱分析 要求首先画出语音信号的时域波形,然后对语音信号进行频谱分析，在MATLAB中，可以利用函数fft对信号进行快速付立叶变换，得到信号的频谱特性；从而加深对频谱特性的理解。 解析： c lear;clc;clf; %语音信号的频谱分析 y=wavread('2.wav'); [y,Fs,nbits]=wavread('2. wav'); N=2048;Y=fft(y,N); Y1=fftshift(Y); plot(abs(Y));title('语音信号的幅度谱');

信号处理实验指导

目录 绪论 (1) 1离散时间信号和系统分析 1.1 离散时间信号产生与运算 (2) 1.2 离散时间系统的时域分析 (9) 1.3 离散时间系统的频域分析 (13) 1.4 离散时间系统频响的零极点确定 (14) 2快速傅立叶变换的应用 2.1 FFT的计算 (17) 2.2 利用FFT进行谱分析 (18) 2.3利用FFT实现快速卷积 (19) 3数字滤波器的设计 3.1数字滤波器的结构 (23) 3.2无限冲激响应(IIR)数字滤波器的设计 (25) 3.3有限冲激响应(FIR)数字滤波器的设计 (27) 4综合应用举例 4.1 语音信号处理 (32) 4.2 电话拨号音的合成与识别 (32)

绪论 数字信号处理主要研究如何对信号进行分析、变换、综合、估计与识别等加工处理的基本理论和方法。随着计算机技术和大规模集成电路技术的发展,数字信号处理以其方便、灵活等特点引起人们越来越多的重视。在40多年的发展过程中,这门学科基本形成了一套完整的理论体系，其中也包括各种快速、优良的算法,而且数字信号处理的理论和技术也在不断、快速地丰富和完善，新理论和新技术也层出不穷。学习这门课程的过程中，容易使人感到数字信号处理的概念抽象难懂，其中的分析方法与基本理论不容易很好地理解与掌握。因此，如何理解与掌握课程中的基本概念、基本原理、基本分析方法以及综合应用所学知识解决实际问题的能力，是本课程学习中所要解决的关键问题。 Matlab是一种面向科学和工程的高级语言，现已成为国际上公认的优秀的科技界应用软件，在世界范围内广为流行和使用。在欧美高等院校里,Matlab已成为大专院校学生、教师的必要基本技能，广泛应用于科学研究、工程计算、教学等。上世纪90年代末和本世纪初Matlab在我国也被越来越多地应用于科研和教学工作中。Matlab是一套功能强大的工程计算及数据处理软件，在工业，电子，医疗和建筑等众多领域均被广泛运用。它是一种面向对象的，交互式程序设计语言，其结构完整又具有优良的可移植性。它在矩阵运算，数字信号处理方面有强大的功能。另外，Matlab提供了方便的绘图功能，便于用户直观地输出处理结果。 本文通过Matlab系列仿真，旨在掌握基本的数字信号处理的理论和方法，提高综合运用所学知识，提高Matlab计算机编程的能力。进一步加强独立分析问题、解决问题的能力、综合设计及创新能力的培养，同时注意培养实事求是、严肃认真的科学作风和良好的实验习惯。

数字信号处理课设+语音信号的数字滤波

语音信号的数字滤波 ——利用双线性变换法实现IIR数字滤波器的设计一．课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 二．设计方案论证 1.IIR数字滤波器设计方法 IIR数字滤波器是一种离散时间系统，其系统函数为 假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和，它 是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。 2.用双线性变换法设计IIR数字滤波器 脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用z=e sT转换 平面的-π/T～π到Z平面上。也就是说，第一步先将整个S平面压缩映射到S 1 /T一条横带里；第二步再通过标准变换关系z=e s1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图1所示。 图1双线性变换的映射关系 为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现

matlab数字信号处理实验指导

电工电子实验中心实验指导书 数字信号处理 实验教程 二○○九年三月

高等学校电工电子实验系列 数字信号处理实验教程 主编石海霞周玉荣 攀枝花学院电气信息工程学院 电工电子实验中心

内容简介 数字信号处理是一门理论与实践紧密联系的课程，适当的上机实验有助于深入理解和巩固验证基本理论知识，了解并体会数字信号处理的CAD手段和方法，锻炼初学者用计算机和MATLAB语言及其工具箱函数解决数字信号处理算法的仿真和滤波器设计问题的能力。 本实验指导书结合数字信号处理的基本理论和基本内容设计了八个上机实验，每个实验对应一个主题内容，包括常见离散信号的MATLAB产生和图形显示、离散时间系统的时域分析、离散时间信号的DTFT、离散时间信号的Z变换、离散傅立叶变换DFT、快速傅立叶变换FFT及其应用、基于MATLAB的IIR和FIR数字滤波器设计等。此外，在附录中，还简单介绍了MATLAB的基本用法。每个实验中，均给出了实验方法和步骤，还有部分的MATLAB程序，通过实验可以使学生掌握数字信号处理的基本原理和方法。

目录 绪论 (1) 实验一常见离散信号的MATLAB产生和图形显示 (2) 实验二离散时间系统的时域分析 (6) 实验三离散时间信号的DTFT (9) 实验四离散时间信号的Z变换 (14) 实验五离散傅立叶变换DFT (18) 实验六快速傅立叶变换FFT及其应用 (24) 实验七基于MATLAB的IIR数字滤波器设计 (30) 实验八基于MATLAB的FIR数字滤波器设计 (33) 附录 (37) 参考文献 (40)

绪论 绪论 随着电子技术迅速地向数字化发展，《数字信号处理》越来越成为广大理工科，特别是IT领域的学生和技术人员的必修内容。 数字信号处理是把信号用数字或符号表示成序列，通过计算机或通用（专用）信号处理设备，用数值计算方法进行各种处理，达到提取有用信息便于应用的目的。数字信号处理的理论和技术一出现就受到人们的极大关注，发展非常迅速。而且随着各种电子技术及计算机技术的飞速发展，数字信号处理的理论和技术还在不断丰富和完善，新的理论和技术层出不穷。目前数字信号处理已广泛地应用在语音、雷达、声纳、地震、图象、通信、控制、生物医学、遥感遥测、地质勘探、航空航天、故障检测、自动化仪表等领域。 数字信号处理是一门理论和实践、原理和应用结合紧密的课程，由于信号处理涉及大量的运算，可以说离开了计算机及相应的软件，就不可能解决任何稍微复杂的实际应用问题。Matlab是1984年美国Math Works公司的产品，MATLAB 语言具备高效、可视化及推理能力强等特点，它的推出得到了各个领域专家学者的广泛关注，其强大的扩展功能为各个领域的应用提供了基础，是目前工程界流行最广的科学计算语言。早在20世纪90年代中期，MATLAB就己成为国际公认的信号处理的标准软件和开发平台。从1996年后，美国新出版的信号处理教材就没有一本是不用MATLAB的。 本实验指导书结合数字信号处理的基本理论和基本内容，用科学计算语言MATLAB实现数字信号处理的方法和实践，通过实验用所学理论来分析解释程序的运行结果，进一步验证、理解和巩固学到的理论知识，从而达到掌握数字信号处理的基本原理和方法的目的。

数字信号处理课程规划报告

数字信号处理课程设计报告《应用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('频域信号'); 运行程序得到的图形：

数字信号处理实验五

实验五：FIR数字滤波器设计与软件实现 信息学院 10电本2班王楚炘 2010304224 10.5.1 实验指导 1．实验目的 （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 （3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。 2．实验内容及步骤 （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示； 图10.5.1 具有加性噪声的信号x(t)及其频谱如图（3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。 （4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，

调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 （4）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示：MATLAB函数fir1和fftfilt的功能及其调用格式请查阅本书 第7章和第？章； 采样频率Fs=1000Hz，采样周期T=1/Fs； 根据图10.6.1(b)和实验要求，可选择滤波器指标参数：通带截止频率fp=120Hz，阻带截至频率fs=150Hz，换算成数字频率，通带截止频率，通带最大衰为0.1dB，阻带截至频率，阻带最小衰为60dB。]实验程序框图如图10.5.2所示，供读者参考。 Fs=1000，T=1/Fs xt=xtg 产生信号xt, 并显示xt及其频谱 用窗函数法或等波纹最佳逼近法 设计FIR滤波器hn 对信号xt滤波：yt=fftfilt(hn,xt) 1、计算并绘图显示滤波器损耗函数 2、绘图显示滤波器输出信号yt End 图10.5.2 实验程序框图 4.思考题 (1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?请写出设计步骤. 答：用窗函数法设计线性相位低通滤波器的设计步骤: a.根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口的长度N； b.构造希望逼近的频率响应函数； c.计算h d(n)； d.加窗得到设计结果h(n)=h d(n)w(n)。 (2)如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为和，阻带上、下截止频率为和，试求理想带通滤波器的截止频率。 答：希望逼近的理想带通滤波器的截止频率分别为：

数字信号处理课设共18页文档

数字信号处理课程设计 姓名：刘倩 学号：201014407 专业：信息与计算科学 实验一：常见离散信号产生和实现 一、实验目的： 1、加深对常用离散信号的理解； 2、掌握matlab 中一些基本函数的建立方法。 二、实验原理： 1.单位抽样序列 在MATLAB 中可以利用zeros()函数实现。 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： 2．单位阶越序列 在MATLAB 中可以利用ones()函数实现。 3．正弦序列 在MATLAB 中 4．复指数序列 在MATLAB 中 5．指数序列 在MATLAB 中

实验内容：由周期为10的正弦函数生成周期为20的余弦函数。 实验代码： n=0:30; y=sin(0.2*pi*n+pi/2); y1=sin(0.1*pi*n+pi/2); subplot(121) stem(n,y); xlabel ('时间序列n');ylabel('振幅');title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); subplot(122) stem(n,y1); xlabel ('时间序列n');ylabel('振幅'); title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); 实验结果： 实验二：离散系统的时域分析 实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。实验原理：离散系统 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号， 记系统单位冲激响应 则系统响应为如下的卷积计算式：

当N k d k ,...2,1,0==时，h[n]是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。 实验内容：用MATLAB 计算全解 当n>=0时，求用系数差分方程y[n]+y[n-1]-6y[n-2]=x[n]描述的一个离散时间系统对阶跃输入x[n]=8μ[n]的全解。 实验代码： n=0:7; >> [y,sf]=filter(1,[1 1 -6],8*ones(1,8),[-7 6]); >> y1(n+1)=-1.8*(-3).^n+4.8*(2).^n-2; >> subplot(121) >> stem(n,y); >> title('由fliter 函数计算结果'); >> subplot(122) >> stem(n,y1); >> title('准确结果'); 实验结果： 结果分析：有图可得由fliter 函数得出的结果与计算出的准确结果完全一致。 实验三FFT 算法的应用

数字信号处理实验指导书2016-通信

数字信号处理实验 徐俊 2016年8月

目录 实验一离散时间系统及系统响应 (2) 实验二离散傅立叶变换及其特性验证 (8) 实验三时域采样与频域采样 (17) 实验四冲激响应不变法IIR数字滤波器设计 (24)

实验一离散时间系统及系统响应 一、实验目的 1、掌握求解离散时间系统脉冲响应和阶跃响应的方法； 2、掌握用线性卷积求解离散时间系统响应的基本方法。 二、实验原理与设计方法 1、用impz和dstep函数求解离散系统的单位脉冲响应和阶跃响应 【例1-1】已知某因果系统的差分方程为 系统为零状态，求系统的脉冲响应和阶跃响应。 解：该系统是一个2阶系统，列出b m和a k系数为 a0=1，a1=0.5，a2=0， b0=1，b1=0，b2=2 MALAB程序如下（取16点作图）： a=[1,0.5,0]; b=[1,0,2]; n=16; hn=impz(b,a,n); %脉冲响应 gn=dstep(b,a,n); %阶跃响应 subplot(1,2,1),stem(hn,'k'); title('系统的单位脉冲响应'); ylabel('h(n)');xlabel('n'); axis([0,n,1.1*min(hn),1.1*max(hn)]); subplot(1,2,2),stem(gn,'k'); title('系统的单位阶跃响应'); ylabel('g(n)');xlabel('n'); axis([0,n,1.1*min(gn),1.1*max(gn)]); 结果如下图所示：

2、用conv函数进行卷积计算求系统响应 【例1-2】某离散时间系统的脉冲响应为 h b(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) 激励信号为 设A=444.128，α=50，。 试求该系统在输入信号激励下的响应。 解：MATLAB程序如下： n=1:50; %定义序列的长度是50 hb=zeros(1,50); %注意：MA TLAB中数组下标从1开始hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1; close all; subplot(3,1,1); stem(hb);title('系统h[n]'); m=1:50;T=0.001; %定义序列的长度和采样率 A=444.128;a=50*sqrt(2.0)*pi; %设置信号有关的参数 w0=50*sqrt(2.0)*pi; x=A*exp(-a*m*T).*sin(w0*m*T); subplot(3,1,2);stem(x);title('输入信号x[n]'); y=conv(x,hb); subplot(3,1,3);stem(y);title('输出信号y[n]'); 结果如下图所示：

五款信号完整性仿真分析工具

SI 五款信号完整性仿真工具介绍 （一）Ansoft公司的仿真工具 现在的高速电路设计已经达到GHz的水平，高速PCB设计要求从三维设计理论出发对过孔、封装和布线进行综合设计来解决信号完整性问题。高速PCB 设计要求中国工程师必须具备电磁场的理论基础，必须懂得利用麦克斯韦尔方程来分析PCB设计过程中遇到的电磁场问题。目前，An soft公司的仿真工具能够从三维场求解的角度出发，对PCB 设计的信号完整性问题进行动态仿真。 Ansoft 的信号完整性工具采用一个仿真可解决全部设计问题： Slwave是一种创新的工具，它尤其适于解决现在高速PCB和复杂IC封装中普遍存在的电源输送和信号完整性问题。 该工具采用基于混合、全波及有限元技术的新颖方法，它允许工程师们特性化同步开关噪声、电源散射和地散射、谐振、反射以及引线条和电源/地平面之间的耦合。该工具采用一个仿真方案解决整个设计问题，缩短了设计时间。 它可分析复杂的线路设计，该设计由多重、任意形状的电源和接地层，以及任何 数量的过孔和信号引线条构成。仿真结果采用先进的3D 图形方式显示，它还可产生等效电路模型，使商业用户能够长期采用全波技术，而不必一定使用专有仿 （二）SPECCTRAQuest Cade nee的工具采用Sun的电源层分析模块: Cade nee Design System 的SpeeetraQuest PCB信号完整性套件中的电源完整性模块据称能让工程师在高速PCB设计中更好地控制电源层分析和共模EMI 。 该产品是由一份与Sun Microsystems公司签署的开发协议而来的，Sun最初研制该项技术是为了解决母板上的电源问题。 有了这种新模块，用户就可根据系统要求来算出电源层的目标阻抗；然后基于板上的器件考虑去耦合要求，Shah表示，向导程序能帮助用户确定其设计所要求的去耦合电容的数目和类型；选择一组去耦合电容并放置在板上之后，用户就可运行一个仿真程序，通过分析结果来发现问题所在。 SPECCTRAQuest是CADENCE公司提供的高速系统板级设计工具，通过它可以控制与PCB layout相应的限制条件。在SPECCTRAQuest菜单下集成了一下工具： （1）SigXplorer 可以进行走线拓扑结构的编辑。可在工具中定义和控制延时、特性阻抗、驱动和负载的类型和数量、拓扑结构以及终端负载的类型等等。可在

数字信号处理实验指导手册

数字信号处理实验指导手册 西安文理学院 机械电子工程系

目录 实验一离散时间信号 (2) 实验二时域采样定理 (7) 实验三离散时间系统 (10) 实验四线性卷积与圆周卷积 (13) 实验五用FFT作谱分析 (16) 实验六用双线性变换法设计IIR数字滤波器 (18) 实验七 FIR滤波器设计 (20)

实验一 离散时间信号 【实验目的】 用MATLAB 实现离散时间信号的表示和运算，掌握MATLAB 的基本命令和编程方法，为后续实验打基础。 【实验原理】 在数字信号处理中，所有的信号都是离散时间信号，因此应首先解决在MATLAB 中如何表示离散信号。 设一模拟信号经A/D 变换后，得到序列信号 }),1(),0(),1(,{)}({)( x x x n x n x -== 由于MATLAB 对下标的约定为从1开始递增，因此要表示)(n x ，一般应采用两个矢量，如： ]5,4,3,2,1,0,1,2,3[---=n ]1,2,5,4,0,2,3,1,1[-=y 这表示了一个含9个采样点的矢量： )}5(,),1(),2(),3({)(x x x x n y ---= 【实验内容】 熟悉下面序列（信号）的产生方法及相关运算 1、 单位采样序列 2、 单位阶跃序列 3、 信号翻转 4、 信号相加 5、 信号折叠 6、 信号移位 【参考程序】 单位采样序列 1、impluse1.m （图1-1） n=10; x=zeros(1,n); x(1)=1;

plot(x,'*'); 2、 impluse2.m （图1-2） n=-5:5; x=[n==0]; stem(x,'*'); 3、impluse3.m （图1-3） n=1:10; n0=3; x=[(n-n0)==1]; plot(x,'*'); 单位阶跃序列 1、steps1.m （图1-4） n=10; x=ones(1,n); plot(x,'*'); 2、steps2.m （图1-5） n=10; x=ones(1,n); x(1)=0; x(2)=0; 图1-1 单位采样序列1 图1-2 单位采样序列2 图1-3 单位采样序列3

数字信号处理课程设计

数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：

摘要 滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR 滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB 窗函数法巴特沃斯

目录 摘要 (1) 1 引言 (1) 1.1课程设计目的 (1) 1.2 课程设计内容及要求 (1) 1.3课程设计设备及平台 (1) 1.3.1 数字滤波器的简介及发展 (1) 1.3.2 MATLAB软件简介 (2) 2 课程设计原理及流程 (4) 3.课程设计原理过程 (4) 3.1 语音信号的采集 (4) 3.2 语音信号的时频分析 (5) 3.3合成后语音加噪声处理 (7) 3.3.1 噪声信号的时频分析 (7) 3.3.2 混合信号的时频分析 (8) 3.4滤波器设计及消噪处理 (10) 3.4.1 设计IIR和FIR数字滤波器 (10) 3.4.2 合成后语音信号的消噪处理 (13) 3.4.3 比较滤波前后语音信号的波形及频谱 (13) 3.4.4回放语音信号 (15) 3.5结果分析 (15) 4 结束语 (15) 5 参考文献 (16)

数字信号处理实验指导书(2016)

数字信号处理实验指导书 颜华刘笑楠刘斌主编 二○○八年三

前言

目录 前言 ....................................................................................................................................................................... - 1 - 目录 ....................................................................................................................................................................... - 1 - 实验一信号、系统及系统响应 ......................................................................................................................... - 2 - 1.1实验目的及要求 ..................................................................................................................................... - 2 - 1.2实验内容 ................................................................................................................................................. - 2 - 1.3实验步骤 ................................................................................................................................................. - 2 - 1.4思考题 ..................................................................................................................................................... - 4 - 1.5选做题 ..................................................................................................................................................... - 4 - 1.6实验报告要求 ......................................................................................................................................... - 4 - 实验二信号的谱分析 ......................................................................................................................................... - 5 - 2.1实验目的及要求 ..................................................................................................................................... - 5 - 2.2实验内容 ................................................................................................................................................. - 5 - 2.3实验步骤、实验原理及思考题 ............................................................................................................. - 5 - 2.4选做题 ..................................................................................................................................................... - 7 - 2.5实验报告要求 ......................................................................................................................................... - 7 - 实验三IIR滤波器设计 ...................................................................................................................................... - 8 - 3.1实验目的及要求 ..................................................................................................................................... - 8 - 3.2实验内容 ................................................................................................................................................. - 8 - 3.3实验步骤及实验报告要求 ..................................................................................................................... - 8 - 3.4 选做题 .................................................................................................................................................... - 9 - 3.5实验报告要求 ....................................................................................................................................... - 10 - 实验四FIR滤波器设计 ................................................................................................................................... - 11 - 4.1实验目的及要求 ................................................................................................................................... - 11 - 4.2实验原理与方法 ................................................................................................................................... - 11 - 4.3实验步骤及内容 ................................................................................................................................... - 12 - 4.4思考题： ............................................................................................................................................... - 13 - 4.5选做题 ................................................................................................................................................... - 13 - 4.6 实验报告要求 ...................................................................................................................................... - 13 -