信号与系统实验报告实验一-信号与系统的时域分析

实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MATLAB函数；2、掌握连续时间和离散时间信号的MATLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程；3、牢固掌握系统的单位冲激响应的概念，掌握LTI系统的卷积表达式及其物理意义，掌握卷积的计算方法、卷积的基本性质；4、掌握利用MATLAB计算卷积的编程方法，并利用所编写的MATLAB程序验证卷积的常用基本性质；掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MATLAB求解LTI系统响应，绘制相应曲线。基本要求：掌握用MATLAB描述连续时间信号和离散时间信号的方法，能够编写MATLAB程序，实现各种信号的时域变换和运算，并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MATLAB描述的方法，掌握卷积运算、线性常系数微分方程的求解编程。二、实验原理信号（Signal）一般都是随某一个或某几个独立变量的变化而变化的，例如，温度、压力、声音，还有股票市场的日收盘指数等，这些信号都是随时间的变化而变化的，还有一些信号，例如在研究地球结构时，地下某处的密度就是随着海拔高度的变化而变化的。一幅图片中的每一个象素点的位置取决于两个坐标轴，即横轴和纵轴，因此，图像信号具有两个或两个以上的独立变量。在《信号与系统》课程中，我们只关注这种只有一个独立变量（Independentvariable）的信号，并且把这个独立变量统称为时间变量（Timevariable），不管这个独立变量是否是时间变量。在自然界中，大多数信号的时间变量都是连续变化的，因此这种信号被称为连续时间信号（Continuous-TimeSignals）或模拟信号（AnalogSignals），例如前面提到的温度、压力和声音信号就是连续时间信号的例子。但是，还有一些信号的独立时间变量是离散变化的，这种信号称为离散时间信号。前面提到的股票市场的日收盘指数，由于相邻两个交易日的日收盘指数相隔24小时，这意味着日收盘指数的时间变量是不连续的，因此日收盘指数是离散时间信号。而系统则用于对信号进行运算或处理，或者从信号中提取有用的信息，或者滤出信号中某些无用的成分，如滤波，从而产生人们所希望的新的信号。系统通常是由若干部件或单元组成的一个整体（Entity）。系统可分为很多不同的类型，例如，根据系统所处理的信号的不同，系统可分为连续时间系统（Continuous-timesystem）和离散时间系统（Discrete-timesystem），根据系统所具有的不同性质，系统又可分为因果系统（Causalsystem）和非因果系统（Noncausalsystem）、稳定系统（Stablesystem）和不稳定系统（Unstablesystem）、线性系统（Linearsystem）和非线性系统（Nonlinearsystem）、时变系统（Time-variantsystem）和时不变系统（Time-invariantsystem）等等。然而，在信号与系统和数字信号处理中，我们所分析的系统只是所谓的线性时不变系统，这种系统同时满足两个重要的基本性质，那就是线性性和时不变性，通常称为线性时不变（LTI）系统。1.信号的时域表示方法1.1将信号表示成独立时间变量的函数例如x(t)=sin(ωt)和x[n]=n(0.5)nu[n]分别表示一个连续时间信号和一个离散时间信号。在MATLAB中有许多内部函数，可以直接完成信号的这种表达，例如：sin()：正弦信号cos()：余弦信号exp()：指数信号1.2用信号的波形图来描述信号用函数曲线表示一个信号，图1.1就是一个连续时间信号和一个离散时间信号的波形图。图1.1连续时间信号与离散时间信号的波形图1.3将信号用一个数据序列来表示对于离散时间信号，还可以表示成一个数的序列，例如：x[n]={....,0.1,1.1,-1.2,0,1.3,….}↑n=0在《信号与系统》和《数字信号处理》课程中，上述三种信号的描述方法是经常要使用的。2用MATLAB仿真连续时间信号和离散时间信号在MATLAB中，无论是连续时间信号还是离散时间信号，MATLAB都是用一个数字序列来表示信号，这个数字序列在MATLAB中叫做向量(vector)。通常的情况下，需要与时间变量相对应。如前所述，MATLAB有很多内部数学函数可以用来产生这样的数字序列，例如sin()、cos()、exp()等函数可以直接产生一个按照正弦、余弦或指数规律变化的数字序列。2.1连续时间信号的仿真程序Program1_1是用MATLAB对一个正弦信号进行仿真的程序，请仔细阅读该程序，并在计算机上运行，观察所得图形。%Program1_1%Thisprogramisusedtogenerateasinusoidalsignalanddrawitsplotclear,%Clearallvariablescloseall,%Closeallfigurewindowsdt=0.01;%Specifythestepoftimevariablet=-2:dt:2;%Specifytheintervaloftimex=sin(2*pi*t);%Generatethesignalplot(t,x)%Openafigurewindowanddrawtheplotofx(t)title('Sinusoidalsignalx(t)')xlabel('Timet(sec)')常用的图形控制函数axis([xmin,xmax,ymin,ymax])：图型显示区域控制函数，其中xmin为横轴的显示起点，xmax为横轴的显示终点，ymin为纵轴的显示起点，ymax为纵轴的显示终点。有时，为了使图形具有可读性，需要在所绘制的图形中，加上一些网格线来反映信号的幅度大小。MATLAB中的gridon/gridoff可以实现在你的图形中加网格线。gridon：在图形中加网格线。gridoff：取消图形中的网格线。x=input(‘Typeinsignalx(t)inclosedform:’)在《信号与系统》课程中，单位阶跃信号u(t)和单位冲激信号δ(t)是二个非常有用的信号。它们的定义如下0,0)(1)(ttdttt1.1(a)0,00,1)(tttu1.1(b)这里分别给出相应的简单的产生单位冲激信号和单位阶跃信号的扩展函数。产生单位冲激信号的扩展函数为：functiony=delta(t)dt=0.01;y=(u(t)-u(t-dt))/dt;产生单位阶跃信号的扩展函数为：%Unitstepfunctionfunctiony=u(t)y=(t=0);%y=1fort0,elsey=0请将这二个MATLAB函数分别以delta和u为文件名保存在work文件夹中，以后，就可以像教材中的方法使用单位冲激信号δ(t)和单位阶跃信号u(t)。2.2离散时间信号的仿真程序Program1_2用来产生离散时间信号x[n]=sin(0.2πn)。%Program1_2%Thisprogramisusedtogenerateadiscrete-timesinusoidalsignalanddrawitsplotclear,%Clearallvariablescloseall,%Closeallfigurewindowsn=-10:10;%Specifytheintervaloftimex=sin(0.2*pi*n);%Generatethesignalstem(n,x)%Openafigurewindowanddrawtheplotofx[n]title('Sinusoidalsignalx[n]')xlabel('Timeindexn')请仔细阅读该程序，比较程序Program1_1和Program1_2中的不同之处，以便自己编程时能够正确使用这种方法方针连续时间信号和离散时间信号。程序Program1_3用来仿真下面形式的离散时间信号：x[n]={....,0.1,1.1,-1.2,0,1.3,….}↑n=0%Program1_3%Thisprogramisusedtogenerateadiscrete-timesequence%anddrawitsplotclear,%Clearallvariablescloseall,%Closeallfigurewindowsn=-5:5;%Specifytheintervaloftime,thenumberofpointsofnis11.x=[0,0,0,0,0.1,1.1,-1.2,0,1.3,0,0];%Generatethesignalstem(n,x,'.')%Openafigurewindowanddrawtheplotofx[n]gridon,title('Adiscrete-timesequencex[n]')xlabel('Timeindexn')由于在程序的stem(n,x,'.')语句中加有'.'选项，因此绘制的图形中每根棒条线的顶端是一个实心点。如果需要在序列的前后补较多的零的话，可以利用函数zeros()，其语法为：zeros(1,N)：圆括号中的1和N表示该函数将产生一个一行N列的矩阵，矩阵中的所有元素均为零。利用这个矩阵与序列x[n]进行组合，从而得到一个长度与n相等的向量。例如，当x[n]={0.1,1.1,-1.2,0,1.3}时，为了得到程序Program1_3中的序列，↑n=0可以用这个MATLAB语句x=[zeros(1,4)xzeros(1,2)]来实现。用这种方法编写的程序如下：%Program1_4%Thisprogramisusedtogenerateadiscrete-timesinusoidalsignal%anddrawitsplotclear,%Clearallvariablescloseall,%Closeallfigurewindowsn=-5:5;%Specifytheintervaloftimex=[zeros(1,4),0.1,1.1,-1.2,0,1.3,zeros(1,2)];%Generatethesequencestem(n,x,'.')%Openafigurewindowanddrawtheplotofx[n]gridon,title('Adiscrete-timesequencex[n]')xlabel('Timeindexn')离散时间单位阶跃信号u[n]定义为0,00,1][nnnu1.2离散时间单位阶跃信号u[n]除了也可以直接用前面给出的扩展函数来产生，还可以利用MATLAB内部函数ones(1,N)来实现。这个函数类似于zeros(1,N)，所不同的是它产生的矩阵的所有元素都为1。值得注意的是，利用ones(1,N)来实现的单位阶跃序列并不是真正的单位阶跃序列，而是一个长度为N单位门(Gate)序列，也就是u[n]-u[n-N]。但是在一个有限的图形窗口中，我们看到的还是一个单位阶跃序列。在绘制信号的波形图时，有时我们需要将若干个图形绘制在图一个图形窗口中，这就需要使用MATLAB的图形分割函数subplot()，其用法是在绘图函数stem或plot之前，使用图形分割函数subplot(n1,n2,n3)，其中的参数n1，n2和n3的含义是，该函数将把一个图形窗口分割成n1xn2个子图，即将绘制的图形将绘制在第n3个子图中。2.3信号的时域变换2.3.1信号的时移信号的时移可用下面的数学表达式来描述：设一个连续时间信号为x(t)，它的时移y(t)表示为：y(t)=x(t-t0)1.3其中，t0为位移量。若t0为正数，则y(t)等于将x(t)右移t0秒之后的结果。反之，若t0为负数，则y(t)等于将x(t)左移t0秒之后的结果。在MATLAB中，时移运算与数学上习惯表达方法完全相同。