软件实训总结Matlab,南昌大学信号与系统matlab软件实验报告.docx

南昌大学信号与系统matlab软件实验报告.docx

1信号与系统软件实验报告班级姓名学号指导教师实验12一．实验要求绘出门函数的波形2FTT二．实验原理在MATLAB中，有一个专门用于表示单位阶跃信号的函数，即STEPFUN函数，它是用数值计算法表示的单位阶跃函数。其调用格式为TSTEPFUNT,T0,其中，T是以向量形式表示的变量，T0表示信号发生突变的时刻，在T0以前，函数值小于零，T0以后函数值大于零。有趣的是它同时还可以表示单位阶跃序列，这只要将自变量以及取样间隔设定为整数即可。有关单位阶跃序K列的表示方法，我们后面有专门论述，下面通过一个例子来说明如何调用STEPFUN函数来表示单位阶跃函数。三．MATLAB程序T40014定义时间样本向量T12指定信号发生突变的时刻U1STEPFUNT,T1产生左移位的阶跃信号T2T22指定信号发生突变的时刻U2STEPFUNT,T2产生右移位的阶跃信号T2GU1U2表示门函数PLOTT,G绘制门函数的波形AXIS4,4,05,15设定坐标轴范围4X4，05Y15四．实验结果3实验2一．实验要求若某连续系统的输入为ET，输出为RT，系统的微分方程为 5 63 YTTYTFF①求该系统的单位冲激响应HT及其单位阶跃响应GT。②若求出系统的零状态响应YT2TFE二．实验原理对于连续的LTI系统，当系统输入为FT，输出为YT，则输入与输出之间满足如下的线性常系数微分方程，当系统输入为00NMIJIJAYB单位冲激信号ΔT时产生的零状态响应称为系统的单位冲激响应，用HT表示。若输入为单位阶跃信号ΕT时，系统产生的零状态响应则称为系统的单位阶跃响应，记为GT，如下图所示。系统的单位冲激响应HT包含了系统的固有特性，它是由系统本身的结构及参数所决定的，与系统的输入无关。我们只要知道了系统的冲激响应，即可求得系统在不同激励下产生的响应。因此，求解系统的冲激响应HT对我们进行连续系统的分析具有非常重要的意义。在MATLAB中有专门用于求解连续系统冲激响应和阶跃响应，并绘制其时域波形的函数IMPULSE和STEP。如果系统输入为FT，冲激响应为HT，系统的零状态响应为YT，则有。YTHFT若已知系统的输入信号及初始状态，我们便可以用微分方程的经典时域求解方法，求出系统的响应。但是对于高阶系统，手工计算这一问题的过程非常困难和繁琐。在MATLAB中，应用LSIM函数很容易就能对上述微分方程所描述的系统的响应进行仿真，求出系统在任意激励信号作用下的响应。LSIM函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解，而且还能同时绘制出系统响应的时域波形图。以上各函数的调用格式如下⑴IMPULSE函数4函数IMPULSE将绘制出由向量A和B所表示的连续系统在指定时间范围内的单位冲激响应HT的时域波形图，并能求出指定时间范围内冲激响应的数值解。IMPULSEB,A以默认方式绘出由向量A和B所定义的连续系统的冲激响应的时域波形。IMPULSEB,A,T0绘出由向量A和B所定义的连续系统在0T0时间范围内冲激响应的时域波形。IMPULSEB,A,T1PT2绘出由向量A和B所定义的连续系统在T1T2时间范围内，并且以时间间隔P均匀取样的冲激响应的时域波形。YIMPULSEB,A,T1PT2只求出由向量A和B所定义的连续系统在T1T2时间范围内，并且以时间间隔P均匀取样的冲激响应的数值解，但不绘出其相应波形。⑵STEP函数函数STEP将绘制出由向量A和B所表示的连续系统的阶跃响应，在指定的时间范围内的波形图，并且求出数值解。和IMPULSE函数一样，STEP也有如下四种调用格式STEPB,ASTEPB,A,T0STEPB,A,T1PT2YSTEPB,A,T1PT2上述调用格式的功能和IMPULSE函数完全相同，所不同只是所绘制(求解)的是系统的阶跃响应GT，而不是冲激响应HT。⑶LSIM函数根据系统有无初始状态，LSIM函数有如下两种调用格式①系统无初态时，调用LSIM函数可求出系统的零状态响应，其格式如下LSIMB,A,X,T绘出由向量A和B所定义的连续系统在输入为X和T所定义的信号时，系统零状态响应的时域仿真波形，且时间范围与输入信号相同。其中X和T是表示输入信号的行向量，T为表示输入信号时间范围的向量，X则是输入信号对应于向量T所定义的时间点上的取样值。YLSIMB,A,X,T与前面的IMPULSE和STEP函数类似，该调用格式并不绘制出系统的零状态响应曲线，而只是求出与向量T定义的时间范围相一致的系统零状态响应的数值解。②系统有初始状态时，调用LSIM函数可求出系统的全响应，格式如下LSIMA,B,C,D,E,T,X0绘出由系数矩阵A,B,C,D所定义的连续时间系统在输入为E和T所定义的信号时，系统输出函数的全响应的时域仿真波形。T为表示输入信号时间范围的向量，E则是输入信号ET对应于向量T所定义的时间点上的取样值，X0表示系统状态变量XX1,X2,XN 在T0时刻的初值。Y,XLSIMA,B,C,D,E,T,X0不绘出全响应波形，而只是求出与向量T定义的时间范围相一致的系统输出向量Y的全响应以及状态变量X的数值解。显然，函数LSIM对系统响应进行仿真的效果取决于向量T的时间间隔的密集程度，T的取样时间间隔越小则响应曲线越光滑，仿真效果也越好。5三．MATLAB程序1A156B32SUBPLOT2,1,1,IMPULSEB,A,4SUBPLOT2,1,2,STEPB,A,42A156B32P1001定义取样时间间隔为001T10P15定义时间范围X1EXP2T1定义输入信号LSIMB,A,X1,T1,对取样间隔为001时系统响应进行仿真HOLDON保持图形窗口以便能在同一窗口中绘制多条曲线P205定义取样间隔为05T20P25定义时间范围X2EXP2T2定义输入信号LSIMB,A,X2,T2,HOLDOFF对取样间隔为05时系统响应进行仿真并解除保持四．实验结果(1)(2)6实验3一．实验要求用数值计算法实现门函数的傅里叶变换，并画出幅度频谱图1FTT二．实验原理对于连续时间信号FT，其傅里叶变换为FJ0LIMJTJNNFEDFE其中Τ为取样间隔，如果FT是时限信号，或者当|T|大于某个给定值时，FT的值已经衰减得很厉害，可以近似地看成是时限信号，则上式中的N取值就是有限的，假定为N，有FJ10NJNNFE若对频率变量Ω进行取样，得KFJ10KNJNFEM通常取，其中是要取的频率范围，或信号的频带宽度。采用2M0MATLAB实现上式时，其要点是要生成FT的N个样本值的向量，以及向量FN，两向量的内积(即两矩阵的乘积)，结果即完成上式的傅里叶变换的数值KJNE计算。三．MATLAB程序R002取样间隔Τ002T2R2T为从2到2，间隔为002的行向量,有201个样本点FTZEROS1,50,ONES1,101,ZEROS1,50产生FT的样值矩阵(即FT的样本值组成的行向量)W110PI取要计算的频率范围7M500K0MWKW1/M频域采样数为M,W为频率正半轴的采样点FWFTEXP