《北京理工大学信号与系统实验讲解电子版》由会员分享，可在线阅读，更多相关《北京理工大学信号与系统实验讲解电子版（14页珍藏版）》请在金锄头文库上搜索。
1、实验1 信号的时域描述与运算一、实验目的 掌握信号的MATLAB表示及其可视化方法。掌握信号基本时域运算的MATLAB实现方法。利用MATLAB分析常用信号，加深对信号时域特性的理解。二、实验原理与方法1. 连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信
2、号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： t=0:0.01:10; x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： x=sin(t); ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形常用的信号产生函数 函数名 功能 函数名 功能heaviside单位阶跃函数rectpuls
3、门函数sin正弦函数tripuls三角脉冲函数cos余弦函数square周期方波sincsinc函数sawtooth周期锯齿波或三角波exp指数函数2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号
4、的微分和积分。这里微分使用差分来近似求取的，由时间向量和采样值向量表示的连续时间信号，其微分可以通过下式求得 其中表示采样间隔。MATLAB中用diff函数来计算差分。 连续时间信号的定积分可以由MATLAB的qud函数实现，调用格式为 quad (function_name,a,b)其中，function_name为被积函数名，a、b为积分区间。 对于符号对象表示的连续时间信号，MATLAB提供了diff函数和quad函数分别用于求微分和积分。3.离散时间信号的MATLAB表示离散时间信号仅在一些离散时刻有定义。在MATLAB中离散时间信号需要使用两个向量来表示，其中一个向量用于表示离散的时
5、间点，另一个向量表示在这些时间点上的值。例如对于如下时间信号 采用MATLAB可以表示如下： n=-3:4; x=-3 2 -1 2 1 -1 2 3; stem(n,x,filled); xlabel(n); title(x(n);Stem函数用于绘制离散时间信号波形，为了与我们表示离散时间信号的习惯相同，在绘图时一般需要添加filled选项，以绘制实心的杆状图形。上述命令绘制的信号时域波形如图3所示。4.离散时间信号的时域运算 离散时间信号的相加相乘是将两个信号对应的时间点上的值相加或相乘，可以直接使用算术运算的运算符“+”和“*”来计算。离散时间信号的位移，则可看作是将表示时间的向量平移
6、，而表示对应时间点上的值的向量不变。离散时间信号的反转，则可以看作是将表示时间的向量和表示对应时间点上的值的向量以零点为基准点，一纵轴为对称轴反折，向量的反折可以利用MATLAB的fliplr函数实现。实验2 LTI系统的时域分析一、实验目的掌握利用MATLAB对系统进行时域分析的方法。掌握连续时间系统零状态响应、冲激响应和阶跃响应的求解方法。掌握求解离散时间系统响应、单位抽样响应的方法。加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。二、实验原理1、连续时间系统时域分析的MATLAB实现1)连续时间系统的MATLAB表示设LTI因果系统的微分方程一般式为：则在MA
7、TLAB里，可以建立系统模型如下： sys=tf(b,a);2）连续时间系统的零状态响应用lsim(sys,x,t)表示求解零状态响应。3）连续时间系统的冲激响应与阶跃响应。用impulse函数来调用。2、离散时间系统时域分析的MATLAB实现1）离散时间系统的MATLAB表示。LTI离散系统通常可以由系统差分方程描述;则在MATLAB里，可以建立系统模型如下： 2）离散时间系统对任意输入的响应。用filter（b,a,x）函数调用。3）离散时间系统的单位抽样响应。用impz函数来调用。3、卷积和与卷积积分1）离散时间序列的卷积和：调用格式为x=conv(x1,x2);2）连续时间信号的卷积积
8、分连续时间信号x1(t)和x2(t)的卷积积分x(t)定义如下实验3 信号频域分析实验目的深入理解信号频谱的概念，掌握信号的频域分析方法。观察典型周期信号和非周期信号的频谱，掌握其频谱特性。实验原理与方法1.连续周期信号的频谱分析如果周期信号满足狄里赫利条件，就可以展开为傅里叶级数形式，即 （1） （2）式中，表示基波周期，为基波频率，表示任一个基波周期内的积分。式（1）和式（2）定义为周期信号复指数形式的傅里叶级数，系数称为的傅里叶系数。周期信号的傅里叶级数还可以由三角函数的线性组合来表示，即 （3）其中： （4）式（3）中同频率的正弦项和余弦项可以合并，从而得到三角函数形式的傅里叶级数，即
9、 （5）其中： （6）可见，任何满足狄里赫利条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。一般来说周期信号表示为傅里叶级数时需要无限多项才能完全逼近原信号，但在实际应用中经常采用有限项级数来替代，所选项数越多就越逼近原信号。2.连续非周期信号的频谱分析对于非周期连续时间信号，吸纳后的傅里叶变换和傅里叶逆变换定义为 （7） （8）式（7）和式（8）把信号的时域特性和频域特性联系起来，确立了非周期信号和频谱之间的关系。采用MATLAB可以方便地求取非周期连续时间信号的傅里叶变换，这里我们介绍常用的集中方法。符号运算法MATLAB的符号数学工具箱提供了直接求解傅里叶变换和反变换
10、的函数，fourier函数和ifourier函数，基本调用格式为 X=fourier(x) X=ifourier(X)默认的时域变量为t，频域变量为w。数值积分法除了采用符号运算的方法外，我们还可以利用MATLAB的quad函数，采用数值积分的方法来进行连续信号的频谱分析，quad函数是一个用来计算数值积分的函数。利用quad函数可以计算非周期连续时间信号的频谱。Quad函数的一般调用格式为：y=quad(fun,a,b)y=quad(fun,a,b,TOL,TRACE,p1,p2,)其中fun指定被积函数，可以采用inline命令来创建，也可以通过传递函数句柄的形式来指定，a、b表示定积分的
11、下限和上限，TOL表示允许的相对或绝对积分误差，TRACE表示以被积函数的点绘图形式来跟踪该函数的返回值，如果TOL和TRACE为空矩阵，则使用缺省值，“p1,p2，”表示被积函数出时间t之外所需的其他额外输入参数。3）数值近似法我们还可以利用MATLAB的数值计算的方法近似计算连续时间傅里叶变换。傅里叶变换可以由式（9）近似计算 （9）当为时限信号，且足够小，式（9）可以演变为 （10）而式（10）中求和部分又可以表示成一个行向量和一个列向量的乘积 （11）式（11）可以很方便地利用MATLAB实现。3.离散周期时间信号的频域分析基波周期为N的周期序列可以用N个成谐波关系的复指数序列的加权和
12、表示，即 （12） 这里k=表示求和仅需包括一个周期内的N项，周期序列在一个周期内的求和与起点无关。将周期序列表示成式（12）的形式，成为离散傅里叶级数，而系数则称为离散傅里叶系数。离散傅里叶系数可以由式（13）确定。 （13）傅里叶系数也称为的频谱系数，而且可以证明是以N为周期的离散频率序列。这说明了周期的离散时间函数对应于频域为周期的离散频率。这里，我们用周期N与傅里叶系数的乘积来表示周期离散时间信号的频谱，即 （14）可以利用MATLAB提供的函数fft用来计算，调用格式为 该函数返回一个周期内的值，其中x表示一个周期内的样本值。4.离散非周期时间信号的频域分析非周期序列可以表示成一组复指数序列的连续和 （15）其中 （16）式（16）称为的离散时间傅里叶变换，式（15）和式（16）确立了非周期离散时间信号及其离散时间傅里叶变换之间的关系。是连续频率的函数，称为频谱函数，且是周期的连续频率函数，其周期为。可见，非周期离散时间函数对应于频域中是一个连续的周期的频率函数。对于有限长的时间序列，式（16）可以表示为 (17)式（17）可以方便地利用MATLAB实现。实验4 LTI系统的频域分析一、实验目的加深对LTI系统频域响应基本概念的掌握和理解
