实验2信号观察与测量

实验三信号的观察与测量

一、实验目的

1．学习使用示波器、函数信号发生器和交流电压表。

2．掌握定量测量电信号的幅值、周期和相位的方法。

二、必备知识

1.直流、正弦交流、方波信号是电路实验中常用的电源信号，可由直流稳压电源、函数信号发生器提供。这些信号的波形可用示波器进行观察。

测试信号幅度的常用仪器有万用表、交流电压表和示波器。直流电压一般用万用表测量比较方便，如需要也可用示波器测量。应注意用万用表和交流电压表测量交流电时，读数都是有效值，它们的测试对象仅限于正弦交流电；示波器可测量各种信号波形，它的读数为峰-峰值。

本实验室示波器、交流电压表和万用表测试正弦电压的性能比较如表.3.1：

2．示波器测量方法：

示波器定量测量时，垂直电压分度“VOL TS/div”旋钮和扫描时间旋钮“TIME/div”的“微调”旋钮应置于校准位置（即示波器屏幕显示区10、17和2处无“>”符号）。

a．直流电压的测量

首先使屏幕显示一水平扫描线。输入耦合方式置于“GND”，此时显示的扫描线为零电平的参考基准线，再将输入耦合方式置于“DC”位置。输入端加上被测信号，此时，“VOL TS/div”档位所指的数值与信号在垂直方向位移的格数相乘，即为测得的直流电压值。高于或低于零电平的电压分别为正值和负值。

例：被测点距基准电平为1.8格，如“VDL TS/div”档位置于5V/div，则直流电压

为：

U =1.8×5=9V 。 b ．交流电压的测量

如图2.4.1 ，如果“VDL TS/div ” 档位置于2V/div （此数值在示波器屏幕显示区位置11或18处显示），屏幕上显示被测信号峰-峰之间的高度为4格，计算方法为：

U 峰-峰=2×4=8V

电压最大值： U M =21×U 峰-峰电压有效值： U =

M U

如果用光标测量，可将两条水平光标线移到波形两测试点，直接从示波器屏幕

下方文字显示读出电压值。

图 3.1 电压和周期的测量

c ．时间、频率的测量

如图3.1，屏幕显示被测信号一个周期所占格数为6格，如果扫描时间“TIME/div ”的档位在0.5ms/div （此数值在示波器屏幕显示区位置3处显示），则周期T =6×0.5=3mS ，频率f =1/T =333H Z 。

时间、频率也可用光标测量，只需将两条竖直光标线调到信号一个周期波形位置，直接从示波器屏幕下方文字显示读出周期和频率值。 d ．同频率两信号之间相位差的测量

将两被测信号送入CH1和CH2通道，读出信号一个周期所占的格数为A ，两信

号相应点所占的格数为B ，则相位差为

??=360A

B ? 如图3.2，波形一个周期占6格，两信号相应点所占的格数为1格，则相位差

?=??=603606

1? 图3.2 相位差的测量峰

-A

e．示波器测量电流

用示波器不能直接进行电流测量，可采用间接测量法，即测量被测支路中已知电阻上的电压，然后用欧姆定律计算出电流。若被测支路中无电阻元件，需串接一个小阻值的电阻，此电阻也称取样电阻。

三、预习要求

预习电路实验光盘或本实验指导1.4~1.6中交流电压表、功率函数信号发生器和示波器的使用方法。

四、实验仪器

1．SS-7802A型双踪示波器

2．SG1645型功率函数信号发生器

3．AS2294A型交流电压表

4．十进制电容箱（RX7-O型0~1.111μF）

5．旋转式电阻箱（ZX21 0~99999.9Ω）

五、实验任务

1．观察示波器“标准信号”波形

将CH1或CH2测试线（红色夹子）接到示波器“CAL”输出端。改变触发源或调节触发电平数值，观察波形稳定情况。波形稳定后，用示波器测出该“标准信号”的峰峰值与周期，并与给定的标准值进行比较。

测试值记录：f= V峰-峰=

2．信号发生器输出电压幅值的测量

将信号发生器输出频率调为f = 1KHz，波形选择正弦波。由小到大调节输出幅值，用示波器和交流电压表分别测量，选取3个不同电压值记入表2.4.2，其中最后一次调为信号发生器最大输出电压值。由示波器测量结果计算出有效值并与交流电压表测量结果进行比较，选取一组数据画出波形图。

3．示波器测量信号的频率

将示波器接入信号发生器输出端，信号发生器输出调为U峰-峰= 4V，波形选择方波，频率分别调为200Hz、1650Hz、5KHz（由信号发生器频率计读出），用示波器测出该信号的频率（采用两种测量方法），结果记入表3.3，选取一组数据画出波形图。

测试时注意观察垂直耦合方式（DC/AC）改变对波形的影响，用文字叙述变化过程。

4．示波器测量信号的相位

按图2.4.3接线。信号发生器输出频率f=1KHz、峰-峰值U P-P = 4V的正弦波，用示波器同时观察信号源输出电压与电容电压的波形，调节R或C，观察波形的变化。记录R=2KΩ,C=0.2μF时观察到的波形，并测出它们的相位差。

图3.3 RC测试电路

注意：信号发生器和示波器的接地端需接到一起，否则可能造成电路局部短路。示波器CH1与CH2通道的接地端内部是相连的。

六、实验报告要求

1．按照实验任务的要求，用坐标纸画出所观察的波形，并标明示波器垂直电压分度值和扫描时间档位。

2．归纳示波器使用要点。

信号与系统实验2

实验报告实验二连续时间系统的时域分析一、实验目的： 1、掌握用Matlab进行卷积运算的数值方法和解析方法,加深对卷积积分的理解。 2、学习利用Matlab实现LTI系统的冲激响应、阶跃响应和零状态响应。二、实验内容及步骤实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。

1、编写程序Q2_1，完成)(1t f 与)(2t f 两函数的卷积运算。 2、编写程序Q2_2，完成)(1t f 与)(2t f 两函数的卷积运算。 3、编写程序Q2_3。利用程序Q2_1，验证卷积的相关性质。 (a) 验证性质：)()(*)(t x t t x =δ (b) 验证性质： )()(*)(00t t x t t t x -=-δ 4、编写程序Q2_4。某线性时不变系统的方程为 )(8)(2)(6)(5)(t f t f t y t y t y +'=+'+''， (a)系统的冲激响应和阶跃响应。 (b)输入()()t f t e u t -=，求系统的零状态响应)(t y zs 。三. 实验结果一： dt=0.01 t1=0:dt:2 f1=0.5*t1 t2=0:dt:2 f2=0.5*t2 f=dt*conv(f1,f2) t=0:0.01:4 plot(t,f);axis([-1 5 0 0.8])

二： dt=0.01 t=-3:dt:3 t1=-6:dt:6 ft1=2*rectpuls(t,2) ft2=rectpuls(t,4) y=dt*conv(ft1,ft2) plot(t1,y) axis([-4 4 0 5]) 以上两题出现错误点：（1）最开始模仿例1的写法用function [f,k]=sconv，总提示出现错误（2）t0+t2 ≤ t ≤ t1+t3 不大能理解的运用个特点，在编写的时候总是被忽略。导致t和t1设置的长度总出错。三：（a） dt=0.01 t=0:dt:2 t0=0 t1=0:dt:2t2=0:dt:2

信号与系统实验

《信号与系统及MATLAB实现》实验指导书

前言长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固和理解教学内容，虽然手工演算训练了计算能力和思维方法，但是由于本课程数学公式推导较多，概念抽象，常需画各种波形，作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB，借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MA TLAB的功能非常强大，我们此处仅用到它的一部分，在后续课程中我们还会用到它，在未来地科学研究和工程设计中有可能继续用它，所以有兴趣的同学，可以对MA TLAB 再多了解一些。 MA TLAB究竟有那些特点呢？ 1．高效的数值计算和符号计算功能，使我们从繁杂的数学运算分析中解脱出来； 2．完备的图形处理功能，实现计算结果和编程的可视化； 3．友好的用户界面及接近数学表达式的自然化语言，易于学习和掌握； 4．功能丰富的应用工具箱，为我们提供了大量方便实用的处理工具； MA TLAB的这些特点，深受大家欢迎，由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景，我们编写了这本《信号与系统及MA TLAB实现》指导书，内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》的同时，掌握MA TLAB的基本应用，学会应用MA TLAB的数值计算和符号计算功能，摆脱烦琐的数学运算，从而更注重于信号与系统的基本分析方法和应用的理解与思考，将课程的重点、难点及部分习题用MA TLAB进行形象、直观的可视化计算机模拟与仿真实现，加深对信号与系统的基本原理、方法及应用的理解，为学习后续课程打好基础。另外同学们在进行实验时，最好事先预习一些MA TLAB的有关知识，以便更好地完成实验，同时实验中也可利用MA TLAB的help命令了解具体语句以及指令的使用方法。

语音信号处理实验指导书

语音信号处理实验指导书实验一语音信号采集与简单处理一、实验目的、要求（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

信号与系统实验三

信号与系统实验实验三：信号的卷积小组成员：黄涛13084220 胡焰焰13084219 洪燕东13084217

一、实验目的 1. 理解卷积的物理意义； 2. 掌握运用计算机进行卷积运算的原理和方法； 3. 熟悉卷积运算函数conv的应用；二、预习内容 1. 卷积的定义及物理意义； 2. 卷积计算的图解法； 3. 卷积的应用三、实验原理说明 1.卷积的定义连续时间和离散时间卷积的定义分别如下所示： 2.卷积的计算由于计算机技术的发展，通过编程的方法来计算卷积积分和卷积和已经不再是冗繁的工作，并可以获得足够的精度，因此信号的时域卷积分析法在系统分析中得到了广泛的应用。卷积积分的数值运算可以应用信号的分段求和来实现，即：数值运算只求当时的信号值，则由上式可以得到：上式中实际上就是连续信号等间隔均匀抽样的离散序列的卷积和，当足够小的时候就是信号卷积积分的数值近似。因此，在利用计算机计算两信号卷积积分时，实质上是先将其转化为离散序列，再利用离散卷积和计算原理来计算。 3.卷积的应用 3.1 求解系统响应卷积是信号与系统时域分析的基本手段，主要应用于求解系统响应，已知一LTI系统的单位冲激响应和系统激励信号则系统响应为激励与单位冲激响应的卷积。需要注意的是利用卷积分析方法求得的系统响应为零状态响应。 3.2 相关性分析相关函数是描述两个信号相似程度的量。两信号之间的相关函数一般称之为互相关函数或者互关函数，定义如下：若是同一信号，此时相关函数称为自相关函数或者自关函数：对于相关函数与卷积运算有着密切的联系，由卷积公式与相关函数比较得：可见，由第二个信号反转再与第一个信号卷积即得到两信号的相关函数。 4.涉及的Matlab函数 4.1 conv函数格式w = conv(u,v)，可以实现两个有限长输入序列u，v的卷积运算，得到有限冲激响应系统的输出序列。输出序列长度为两个输入序列长度和减一。四、实验内容给定如下因果线性时不变系统： y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3] (1)不用impz函数，使用filter命令，求出以上系统的单位冲激响应h[n]的前20个样本; clc; N = 0:19;

信号与系统实验2

信号与系统实验信息学院通信工程 20101060163 荣华杰 8.1已知描述连续时间系统的微分方程和激烈信号f （t ）为y''(t)+4y'(t)+2y(t)=f''(t)+3f(t),f(t)=) (t u e t -试用MATLAB 的lsim 函数求出上述系统在0-10秒时间零状态响应y （t ）的值，并绘出系统零状态响应的时域仿真波形。 y''(t)+4y'(t)+2y(t)=f''(t)+3f(t),f(t)=)(t u e t - a=[1 4 2]; b=[1 0 3]; sys=tf(b,a); t=0:1:10; f=exp(-(t)); y=lsim(sys,f,t) y = 1.0000 0.1721 0.3629 0.3138 0.2162 0.1356 0.0811 0.0472 0.0271 0.0153 0.0086 波形图： a=[1 4 2]; b=[1 3]; sys=tf(b,a); t=0:0.01:10; f=exp(-2*t); lsim(sys,f,t) 8.3已知描述系统的微分方程如下，试用MATLAB 求系统在0-10秒时间范围内冲激响应和阶跃响应的数值解，并绘出系统冲击响应和阶跃响应的时域波形。（3）y''(t)+4y'(t)+5y(t)=f'(t) 冲激响应数值解： a=[1 4 5]; b=[1 0]; y=impulse(b,a,0:1:10) y =1.0000

-0.1546 -0.0409 -0.0032 0.0003 0.0001 0.0000 -0.0000 -0.0000 -0.0000 0.0000 阶跃响应数值解： a=[1 4 5]; b=[1 0]; y=step(b,a,0:1:10) y = 0 0.1139 0.0167 0.0003 -0.0003 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 冲激响应、阶跃响应的时域波形： a=[1 4 5]; b=[1 0]; subplot(1,2,1) step(b,a,10) subplot(1,2,2) impulse(b,a,10) 8.4已知描述离散系统的差分方程和输入系列x（n）如下，试用MATLAB的filter函数求出上述系统在0-20时间采样点范围内零状态响应y（n）的系列样值，并绘出系统零状态响应的时域波形。（2）y(n)+(1/2)y(n-1)=x(n)+2x(n-1),x(n)=2cos(n*pi/3)*u(n) 零状态响应y（n）的系列样值: a=[1 1/2]; b=[1 2];

信号与系统实验指导书

实验一常用信号分类与观察一、实验目的 1、了解单片机产生低频信号源； 2、观察常用信号的波形特点及产生方法； 3、学会使用示波器对常用波形参数的测量。二、实验内容 1、信号的种类相当的多，这里列出了几种典型的信号，便于观察。 2、这些信号可以应用到后面的“基本运算单元”和“无失真传输系统分析”中。三、实验原理对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、抽样信号、钟形信号、脉冲信号等。 1、正弦信号：其表达式为)sin()(θω+=t K t f ，其信号的参数：振幅K 、角频率ω、与初始相位θ。其波形如下图所示：图 1-5-1 正弦信号 2、指数信号：指数信号可表示为at Ke t f =)(。对于不同的a 取值，其波形表现为不同的形式，如下图所示：

图 1-5-2 指数信号 3、指数衰减正弦信号：其表达式为 ?? ? ??><=-)0()sin()0(0)(t t Ke t t f at ω 其波形如下图：图 1-5-3 指数衰减正弦信号 4、抽样信号：其表达式为： sin ()t Sa t t = 。)(t Sa 是一个偶函数，t = ±π,±2π,…,±n π时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

信号与系统实验二

实验二常用信号分类与观察一、实验目的 1、观察常用信号的波形特点及产生方法。 2、学会使用示波器对常用波形参数的测量。二、实验内容 1、信号的种类相当的多，这里列出了几种典型的信号，便于观察。 2、这些信号可以应用到后面的“基本运算单元”和“无失真传输系统分析”中。三、实验仪器 1、信号与系统实验箱一台（主板）。 2、20MHz 双踪示波器一台。四、实验原理对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、抽样信号、钟形信号、脉冲信号等。 1、正弦信号：其表达式为)sin()(θω+=t K t f ，其信号的参数：振幅K 、角频率ω、与初始相位θ。其波形如下图所示：图 1-5-1 正弦信号 2、指数信号：指数信号可表示为at Ke t f =)(。对于不同的a 取值，其波形表现为不同的形式，如下图所示：

图 1-5-2 指数信号 3、指数衰减正弦信号：其表达式为 ?? ???><=-)0()sin()0(0)(t t Ke t t f at ω 其波形如下图：图 1-5-3 指数衰减正弦信号 4、抽样信号：其表达式为： sin ()t Sa t t = 。)(t Sa 是一个偶函数，t = ±π,±2π,…,±n π时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

信号与系统实验四

信号与系统实验实验四：周期信号的傅里叶级数小组成员：黄涛13084220 胡焰焰13084219 洪燕东13084217

一、实验目的 1、分析典型的矩形脉冲信号，了解矩形脉冲信号谐波分量的构成。 2、观察矩形脉冲信号通过多个数字滤波器后，分解出各谐波分量的情况。 3、掌握用傅里叶级数进行谐波分析的方法。 4、观察矩形脉冲信号分解出的各谐波分量可以通过叠加合成出原矩形脉冲信号。二、预习内容 1、周期信号的傅里叶级数分解及其物理意义。 2、典型信号傅里叶级数计算方法。三、实验原理 1. 信号的时间特性与频率特性信号可以表示为随时间变化的物理量，比如电压)(t u 和电流)(t i 等，其特性主要表现为随时间的变化，波形幅值的大小、持续时间的长短、变化速率的快慢、波动的速度及重复周期的大小等变化，信号的这些特性称为时间特性。信号还可以分解为一个直流分量和许多不同频率的正弦分量之和。主要表现在各频率正弦分量所占比重的大小不同；主要频率分量所占的频率范围也不同，信号的这些特性称为信号的频率特性。无论是信号的时间特性还是频率特性都包含了信号的全部信息量。 2. 信号的频谱信号的时间特性和频率特性是对信号的两种不同的描述方式。根据傅里叶级数原理，任意一个时域的周期信号)t (f ，只要满足狄利克莱(Dirichlet)条件，就可以将其展开成三角形式或指数形式的傅里叶级数。例如，对于一个周期为T 的时域周期信号)t (f ，可以用三角形式的傅里叶级数求出它的各次分量，在区间),(11T t t +内表示为 ()∑∞ =Ω+Ω+=10sin cos )(n n n t n b t n a a t f 即将信号分解成直流分量及许多余弦分量和正弦分量，研究其频谱分布情况。 A 0t A n 0A 0t （a）（b） Ω（c）ωΩ 5Ω3Ω Ω3Ω5 3. 信号的时间特性与频率特性关系信号的时域特性与频域特性之间有着密切的内在联系，这种联系可以用图4-1来形象地表示。其中图4-1(a)是信号在幅度--时间--频率三维坐标系统中的图形；图4-1(b)是信号在幅度--时间坐标系统中的图形即波形图；把周期信号分解得到的各次谐波分量按频率的高低排列，就可以得到频谱图。反映各频率分量幅度的频谱称为振幅频谱。图4-1(c)是信号在幅度--频率坐标系统中的图形即振幅频谱图。反映各分量相位的频谱称为相位频谱。 4. 信号频谱的测量在本实验中只研究信号振幅频谱。周期信号的振幅频谱有三个性质：离散性、谐波性、收敛

语音信号处理实验报告

语音信号处理实验班级：学号：姓名：实验一基于MATLAB的语音信号时域特征分析（2学时）

1）短时能量（1）加矩形窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32; for i=2:6 h=linspace(1,1,2.^(i-2)*N);%形成一个矩形窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if(i==2) ,legend('N=32'); elseif(i==3), legend('N=64'); elseif(i==4) ,legend('N=128'); elseif(i==5) ,legend('N=256'); elseif(i==6) ,legend('N=512'); end end

00.51 1.52 2.5 3 x 10 4 -1 1 x 10 4 024 x 10 4 05 x 10 4 0510 x 10 4 01020 x 10 4 02040 （2）加汉明窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32; for i=2:6 h=hanning(2.^(i-2)*N);%形成一个汉明窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if(i==2), legend('N=32'); elseif(i==3), legend('N=64'); elseif(i==4) ,legend('N=128');

信号与系统实验报告_1(常用信号的分类与观察)

实验一：信号的时域分析一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台三、实验原理对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示：图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：

图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图：图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示：图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)

北京理工大学信号与系统实验实验报告

实验1 信号的时域描述与运算一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MA TLAB实现方法。 3. 利用MA TLAB分析常用信号，加深对信号时域特性的理解。二、实验原理与方法 1. 连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： >> t=0:0.01:10; >> x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： >> x=sin(t); >> ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 Time(seconds) 图1 利用向量表示连续时间信号

t 图 2 利用符号对象表示连续时间信号 sin(t) 2.连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量[N t t t ,,,21?]和采样值向量[N x x x ,,,21?]表示的连续时间信号，其微分可以通过下式求得 1,,2,1,|)('1-?=?-≈ +=N k t x x t x k k t t k 其中t ?表示采样间隔。MA TLAB 中用diff 函数来计算差分 k k x x -+1。连续时间信号的定积分可以由MATLAB 的qud 函数实现，调用格式为 quad ('function_name',a,b) 其中，function_name 为被积函数名，a 、b 为积分区间。

语音信号处理实验报告

通信与信息工程学院信息处理综合实验报告班级：电子信息工程1502班指导教师：设计时间：2018/10/22-2018/11/23 评语：通信与信息工程学院二〇一八年实验题目：语音信号分析与处理一、实验内容 1. 设计内容利用MATLAB对采集的原始语音信号及加入人为干扰后的信号进行频谱分析，使用窗函数法设计滤波器滤除噪声、并恢复信号。 2．设计任务与要求 1. 基本部分

（1）录制语音信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。（2）对所录制的语音信号加入干扰噪声，并对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。（3）分别利用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman 窗几种函数设计数字滤波器滤除噪声，并画出各种函数所设计的滤波器的频率响应。（4）画出使用几种滤波器滤波后信号时域波形和频谱，对滤波前后的信号、几种滤波器滤波后的信号进行对比，分析信号处理前后及使用不同滤波器的变化；回放语音信号。 2. 提高部分（5）录制一段音乐信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。（6）利用MATLAB产生一个不同于以上频段的信号；画出信号频谱图。（7）将上述两段信号叠加，并加入干扰噪声，尝试多次逐渐加大噪声功率，对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。（8）选用一种合适的窗函数设计数字滤波器，画出滤波后音乐信号时域波形和频谱，对滤波前后的信号进行对比，回放音乐信号。二、实验原理 1.设计原理分析本设计主要是对语音信号的时频进行分析，并对语音信号加噪后设计滤波器对其进行滤波处理，对语音信号加噪声前后的频谱进行比较分析，对合成语音信号滤波前后进行频谱的分析比较。首先用PC机WINDOWS下的录音机录制一段语音信号，并保存入MATLAB软件的根目录下，再运行MATLAB仿真软件把录制好的语音信号用audioread函数加载入MATLAB仿真软件的工作环境中，输入命令对语音信号进行时域，频谱变换。对该段合成的语音信号，分别用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman窗几种函数在MATLAB中设计滤波器对其进行滤波处理，滤波后用命令可以绘制出其频谱图，回放语音信号。对原始语音信号、合成的语音信号和经过滤波器处理的语音信号进行频谱的比较分析。 2.语音信号的时域频域分析在Matlab软件平台下可以利用函数audioread对语音信号进行采样,得到了声音数据变量y,同时把y的采样频率Fs=44100Hz放进了MATALB的工作空间。

信号与系统实验二的题目及答案

第二个信号实验题目 1（1）用数值法求门函数4()G t 的傅里叶变换，并给出门函数的幅频特性曲线和相频特性曲线。（2）用符号法给出函数5()2()()3 t f t e u t -=的傅里叶变换。（3）已知系统函数为34 2 1()3 s s H s s s ++=++，画出该系统的零极点图。 2 （1）用数值法给出函数5(2)2()(2)3 t f t e u t --=-幅频特性曲线和相频特性曲线。（2）对函数5(2)2()(2)3 t f t e u t --=-进行采样，采样间隔为0.01。（3）已知输入信号为()sin(100)f t t =，载波频率为1000Hz ，采样频率为5000 Hz ，试产生输入信号的调幅信号。 3（1）用符号法实现函数4()G t 的傅里叶变换，并给出门函数的幅频特性曲线和相频特性曲线。（2）已知系统函数为34 2 1()3 s s H s s s ++= ++，输入信号为()sin(100)f t t =，求该系统的稳态响应。（3）已知输入信号为()sin(100)f t t =，载波频率为100Hz ，采样频率为400 Hz ，试产生输入信号的调频信号。 4（1）已知系统函数为23 1()3 s s H s s s ++= ++，画出该系统的零极点图。（2）已知函数5()2()()3 t f t e u t -=用数值法给出函数(3)f t 的幅频特性曲线和相频特性曲线。（3）实现系统函数3421 ()3 s s H s s s ++= ++的频率响应。（4）已知输入信号为()cos(100)f t t =，载波频率为100Hz ，采样频率为400 Hz ，试产生输入信号的调相信号。 5（1）用数值法给出函数5(2)2 ()(2)3 t f t e u t -+=+幅频特性曲线和相频特性曲线。（2）用符号法实现函数 2 2i ω +的傅里叶逆变换。（3）已知输入信号为()5sin(200)f t t =，载波频率为1000Hz ，采样频率为5000 Hz ，试产生输入信号的调频信号。

信号与系统实验(新)

信号与系统实验实验1 阶跃响应与冲激响应一、实验目的 1、观察和测量RLC串联电路的阶跃响应与冲激响应的波形和有关参数，并研究其电路元件参数变化对响应状态的影响； 2、掌握有关信号时域的测量方法。二、实验原理说明实验如图1－1所示RLC串联电路的阶跃响应与冲激响应的电路连接图，图1

用周期方波通过微分电路后得到的尖顶脉冲代替冲激信号。三、实验内容 1、阶跃响应波形观察与参数测量设激励信号为方波，其幅度为1.5V 峰峰值，频率为500Hz 。实验电路连接图如图1-1（a ）所示。 ① 连接如图1-1所示 ② 调整激励源信号为方波，调节频率旋钮，使f=500Hz ，调节幅度旋钮，使信号幅度为1.5V 。（注意：实验中，在调整信号源的输出信号的参数时，需连接上负载后调节） ③ 示波器CH1接于TP909，调节滑动变阻器，使电路分别工作于欠阻尼、临界和过阻尼三种状态，并将实验数据填入表格1-1中。 ④ TP908为输入信号波形的测量点，可把示波器的CH ·接于TP908上，便于波形比较。表1-1 注:描绘波形要使三状态的X 轴坐标(扫描时间)一致。 2、冲激响应的波形观察冲激信号是由阶跃信号经过微分电路而得到。实验电路如图1—1（b ）所示。参数测量波形观察欠阻尼状态临界状态过阻尼状态状态参数测量 R< Tr= Ts= δ= R= Tr= R>

①将信号输入接于P905。（频率与幅度不变）； ②将示波器的CH1接于TP906，观察经微分后响应波形（等效为冲激激励信号）； ③连接如图1-1（b）所示 ④将示波器的CH2接于TP909，调整滑动变阻器，使电路分别工作于欠阻尼、临界和过阻尼三种状态 ④观察TP909端三种状态波形，并填于表1-2中。表1-2 表中的激励波形为在测量点TP906观察到的波形（冲激激励信号）。四、实验报告要求 1、描绘同样时间轴阶跃响应与冲激响应的输入、输出电压波形时，要标明信号幅度A、周期T、方波脉宽T1以及微分电路的τ值。 2、分析实验结果，说明电路参数变化对状态的影响。五、实验设备双踪示波器 1 台信号系统实验箱 1台上升时间t r :y(t)从0.1到第一次达到0.9所需时间。峰值时间t p :y(t)从0上升y max 所需的时间。调节时间t s :y(t)的振荡包络线进入到稳态值的% 5 误差范围所需的时间。激励波形响应波形欠阻尼状态临界状态过阻尼状态

大学本科语音信号处理实验讲义8学时

语音信号处理实验讲义时间：2011-12

目录实验一语音信号生成模型分析 (3) 实验二语音信号时域特征分析 (7) 实验三语音信号频域特征分析 (12) 实验四语音信号的同态处理和倒谱分析 (16)

实验一语音信号生成模型分析一、实验目的 1、了解语音信号的生成机理，了解由声门产生的激励函数、由声道产生的调制函数和由嘴唇产生的辐射函数。 2、编程实现声门激励波函数波形及频谱，与理论值进行比较。 3、编程实现已知语音信号的语谱图，区分浊音信号和清音信号在语谱图上的差别。二、实验原理语音生成系统包含三部分：由声门产生的激励函数()G z 、由声道产生的调制函数()V z 和由嘴唇产生的辐射函数()R z 。语音生成系统的传递函数由这三个函数级联而成，即 ()()()()H z G z V z R z = 1、激励模型发浊音时，由于声门不断开启和关闭，产生间隙的脉冲。经仪器测试它类似于斜三角波的脉冲。也就是说，这时的激励波是一个以基音周期为周期的斜三角脉冲串。单个斜三角波的频谱表现出一个低通滤波器的特性。可以把它表示成z 变换的全极点形式 12 1()(1) cT G z e z --= -? 这里c 是一个常数，T 是脉冲持续时间。周期的三角波脉冲还得跟单位脉冲串的z 变换相乘： 112 1 ()()()1(1)v cT A U z E z G z z e z ---=?= ?--? 这就是整个激励模型，v A 是一个幅值因子。 2、声道模型当声波通过声道时，受到声腔共振的影响，在某些频率附近形成谐振。反映在信号频谱图上，在谐振频率处其谱线包络产生峰值，把它称为共振峰。一个二阶谐振器的传输函数可以写成 12 ()1i i i i A V z B z C z --= -- 实践表明，用前3个共振峰代表一个元音足够了。对于较复杂的辅音或鼻音共振峰要到5个以上。多个()i V z 叠加可以得到声道的共振峰模型 12 1 11 ()()11R r r M M i r i N k i i i i k k b z A V z V z B z C z a z -=---======---∑∑∑ ∑ 3、辐射模型从声道模型输出的是速度波，而语音信号是声压波。二者倒比称为辐射阻抗，它表征了

信号与系统实验一

信号与系统实验报告实验一信号与系统的时域分析学院：信息工程学院班级：2012级电子信息工程三班姓名：学号：2012550711 指导老师：苏永新

一、实验目的 1、熟悉和掌握常用的用于信号与系统时域仿真分析的MATLAB函数； 2、掌握连续时间和离散时间信号的MATLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程； 3、牢固掌握系统的单位冲激响应的概念，掌握LTI系统的卷积表达式及其物理意义，掌握卷积的计算方法、卷积的基本性质； 4、掌握利用MATLAB计算卷积的编程方法，并利用所编写的MATLAB程序验证卷积的常用基本性质；掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MATLAB求解LTI系统响应，绘制相应曲线。基本要求：掌握用MATLAB描述连续时间信号和离散时间信号的方法，能够编写MATLAB 程序，实现各种信号的时域变换和运算，并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MATLAB描述的方法，掌握卷积运算、线性常系数微分方程的求解编程。二、实验内容及步骤实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。 Q1-1：修改程序Program1_1，将dt改为0.2，再执行该程序，保存图形，看看所得图形的效果如何？ dt = 0.01时的信号波形 dt = 0.2时的信号波形此处粘贴图形此处粘贴图形% Program1_1 % This program is used to generate a sinusoidal signal and draw its plot clear, % Clear all variables

信号与系统实验指导全部实验答案

00.51 1.52 2.53-2-1.5-1-0.50 0.5 1 1.52正弦信号 -20-15-10-505101520 -0.4 -0.200.20.40.60.81抽样信号 -0.500.51 1.52 2.53 00.5 11.5 2矩形脉冲信号 -1 012345 -0.5 0.5 1 1.5 单位跃阶信号实验一连续时间信号的MATLAB 表示实验目的 1．掌握MATLAB 语言的基本操作，学习基本的编程功能； 2．掌握MATLAB 产生常用连续时间信号的编程方法； 3．观察并熟悉常用连续时间信号的波形和特性。实验原理： 1. 连续信号MA TLAB 实现原理从严格意义上讲，MATLAB 数值计算的方法并不能处理连续时间信号。然而，可用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能够被MATLAB 处理，并且能较好地近似表示连续信号。 MATLAB 提供了大量生成基本信号的函数。比如常用的指数信号、正余弦信号等都是MATLAB 的内部函数。为了表示连续时间信号，需定义某一时间或自变量的范围和取样时间间隔，然后调用该函数计算这些点的函数值，最后画出其波形图。实验内容：

0123-2-1 1 2实部012 3 -1 12 虚部0123012取模0123 -50 5 相角00.2 0.40.6 0.81 -1-0.50 0.51 方波信号实验编程：（1）t=0:0.01:3; K=2;a=-1.5;w=10; ft=K*exp((a+i*w)*t); A=real(ft); B=imag(ft); C=abs(ft); D=angle(ft); subplot(2,2,1),plot(t,A),grid on;title('实部'); subplot(2,2,2),plot(t,B),grid on;title('虚部'); subplot(2,2,3),plot(t,C),grid on;title('取模'); subplot(2,2,4),plot(t,D),grid on;title('相角'); （2） t=0:0.001:3; y=square(2*pi*10*t,30); plot(t,y); axis([0,1,-1,1]); title('方波信号');

信号与系统实验

序列号：＿＿信号与系统实验报告课程名称信号与系统学院信息工程学院年级班别电子信息工程1班学号 3116002166 学生姓名陈俊杰指导教师黄国宏 2018年6月15日

目录实验二LTI系统的响应 (1) 一、实验目的 (1) 二、实验原理 (1) 三、实验内容 (3) 四、程序清单及实验结果 (4) 五、实验总结 (13) 实验三连续时间信号的频域分析一、实验目的 (14) 二、实验原理 (14) 三、实验内容 (17) 四、程序清单及实验结果 (17) 五、实验总结 (25) 实验五连续信号与系统的S域分析一、实验目的 (26) 二、实验原理 (26) 三、实验内容 (27) 四、程序清单及实验结果 (28) 五、实验总结 (36)

实验二 LTI 系统的响应一、实验目的 1. 熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法 2. 熟悉连续（离散）时间系统在任意信号激励下响应的求解方法 3. 熟悉应用MATLAB 实现求解系统响应的方法二、实验原理 1.连续时间系统对于连续的LTI 系统，当系统输入为f (t )，输出为y (t )，则输入与输出之间满足如下的线性常系数微分方程：() ()00()()n m i j i j i j a y t b f t ===∑∑，当系统输入为单位冲激信号δ(t )时产生的零状态响应称为系统的单位冲激响应，用h(t)表示。若输入为单位阶跃信号ε(t )时，系统产生的零状态响应则称为系统的单位阶跃响应，记为g(t)，如下图所示。系统的单位冲激响应h (t )包含了系统的固有特性，它是由系统本身的结构及参数所决定的，与系统的输入无关。我们只要知道了系统的冲激响应，即可求得系统在不同激励下产生的响应。因此，求解系统的冲激响应h(t )对我们进行连续系统的分析具有非常重要的意义。在MATLAB 中有专门用于求解连续系统冲激响应和阶跃响应，并绘制其时域波形的函数impulse( ) 和step( )。如果系统输入为f (t )，冲激响应为h(t)，系统的零状态响应为y (t )，则有：()()()y t h t f t =*。若已知系统的输入信号及初始状态，我们便可以用微分方程的经典时域求解方法，求出系统的响应。但是对于高阶系统，手工计算这一问题的过程非常困难和繁琐。在MATLAB 中，应用lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真，求出系统在任意激励信号作用下的响应。lsim( )函数不仅能

语音信号处理实验报告实验二

通信工程学院12级1班罗恒 2012101032 实验二基于MATLAB 的语音信号频域特征分析一、实验要求要求根据已有语音信号，自己设计程序，给出其倒谱、语谱图的分析结果，并根据频域分析方法检测所分析语音信号的基音周期或共振峰。二、实验目的信号的傅立叶表示在信号的分析与处理中起着重要的作用。因为对于线性系统来说，可以很方便地确定其对正弦或复指数和的响应，所以傅立叶分析方法能完善地解决许多信号分析和处理问题。另外，傅立叶表示使信号的某些特性变得更明显，因此，它能更深入地说明信号的各项红物理现象。由于语音信号是随着时间变化的，通常认为，语音是一个受准周期脉冲或随机噪声源激励的线性系统的输出。输出频谱是声道系统频率响应与激励源频谱的乘积。声道系统的频率响应及激励源都是随时间变化的，因此一般标准的傅立叶表示虽然适用于周期及平稳随机信号的表示，但不能直接用于语音信号。由于语音信号可以认为在短时间内，近似不变，因而可以采用短时分析法。三、实验设备 1.PC 机； 2.MATLAB 软件环境；四、实验内容 1.上机前用Matlab 语言完成程序编写工作。 2.程序应具有加窗（分帧）、绘制曲线等功能。 3.上机实验时先调试程序，通过后进行信号处理。 4.对录入的语音数据进行处理，并显示运行结果。 5.依次给出其倒谱、语谱图的分析结果。 6. 根据频域分析方法检测所分析语音信号的基音周期或共振峰。五、实验原理及方法 1、短时傅立叶变换由于语音信号是短时平稳的随机信号，某一语音信号帧的短时傅立叶变换的定义为：其中w(n -m)是实窗口函数序列，n 表示某一语音信号帧。令n -m=k'，则得到 ()()()jw jwm n m X e x m w n m e ∞-=-∞= -∑