信号与系统-实验报告-实验五

实验五连续信号与系统的S域分析学院班级姓名学号一、实验目的1.熟悉拉普拉斯变换的原理及性质2.熟悉常见信号的拉氏变换3.了解正/反拉氏变换的MATLAB实现方法和利用MATLAB绘制三维曲面图的方法4.了解信号的零极点分布对信号拉氏变换曲面图的影响及续信号的拉氏变换与傅氏变换的关系二、实验原理拉普拉斯变换是分析连续时间信号的重要手段。对于当t∞时信号的幅值不衰减的时间信号，即在f(t)不满足绝对可积的条件时，其傅里叶变换可能不存在，但此时可以用拉氏变换法来分析它们。连续时间信号f(t)的单边拉普拉斯变换F(s)的定义为：0()()stFsftedt拉氏反变换的定义为：1()()2jstjftFsedsj显然，上式中F(s)是复变量s的复变函数，为了便于理解和分析F(s)随s的变化规律，我们将F(s)写成模及相位的形式：()()()jsFsFse。其中，|F(s)|为复信号F(s)的模，而()s为F(s)的相位。由于复变量s=σ+jω，如果以σ为横坐标（实轴），jω为纵坐标（虚轴），这样，复变量s就成为一个复平面，我们称之为s平面。从三维几何空间的角度来看，|()|Fs和()s分别对应着复平面上的两个曲面，如果绘出它们的三维曲面图，就可以直观地分析连续信号的拉氏变换F(s)随复变量s的变化情况，在MATLAB语言中有专门对信号进行正反拉氏变换的函数，并且利用MATLAB的三维绘图功能很容易画出漂亮的三维曲面图。①在MATLAB中实现拉氏变换的函数为：F=laplace(f)对f(t)进行拉氏变换，其结果为F(s)F=laplace(f,v)对f(t)进行拉氏变换，其结果为F(v)F=laplace(f,u,v)对f(u)进行拉氏变换，其结果为F(v)②拉氏反变换f=ilaplace(F)对F(s)进行拉氏反变换，其结果为f(t)f=ilaplace(F,u)对F(w)进行拉氏反变换，其结果为f(u)f=ilaplace(F,v,u)对F(v)进行拉氏反变换，其结果为f(u)注意：在调用函数laplace()及ilaplace()之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量。对laplace()中的f及ilaplace()中的F也要用符号定义符sym将其说明为符号表达式。具体方法参见第一部分第四章第三节。例①：求出连续时间信号()sin()()fttt的拉氏变换式，并画出图形求函数拉氏变换程序如下：symsts%定义符号变量ft=sym('sin(t)*Heaviside(t)');%定义时间函数f(t)的表达式Fs=laplace(ft)%求f(t)的拉氏变换式F(s)运行结果：Fs=1/(s^2+1)绘制拉氏变换三维曲面图的方法有2种：方法一：symsxyss=x+i*y;%产生复变量sFFs=1/(s^2+1);%将F(s)表示成复变函数形式FFss=abs(FFs);%求出F(s)的模ezmesh(FFss);%画出拉氏变换的网格曲面图ezsurf(FFss);%画出带阴影效果的三维曲面图colormap(hsv);%设置图形中多条曲线的颜色顺序方法二：figure(2)%打开另一个图形窗口x1=-5:0.1:5;%设置s平面的横坐标范围y1=-5:0.1:5;%设置s平面的纵坐标范围[x,y]=meshgrid(x1,y1);%产生矩阵s=x+i*y;%产生矩阵s来表示所绘制曲面图的复平面区域，%其中矩阵s包含了复平面-6σ6，-6jω6范围内%以间隔0.01取样的所有样点fs=1./(s.*s+1);%计算拉氏变换在复平面上的样点值ffs=abs(fs);%求幅值mesh(x,y,ffs);%绘制拉氏变换的三维网格曲面图surf(x,y,ffs);%绘制带阴影效果的三维曲面图axis([-5,5,-5,5,0,8]);%设置坐标显示范围colormap(hsv);%设置图形中多条曲线的颜色顺序说明：从拉普拉斯变换的三维曲面图中可以看出，曲面图上有象山峰一样突出的尖峰，这些峰值点在s平面的对应点就是信号拉氏变换的极点位置。而曲面图上的谷点则对应着拉氏变换的零点位置。因此，信号拉氏变换的零极点位置决定了其曲面图上峰点和谷点位置。例②：求出函数21()1Fss的拉氏反变换式MATLAB程序如下：symsts%定义符号变量Fs=sym('1/(1+s^2)');%定义F(s)的表达式ft=ilaplace(Fs)%求F(s)的拉氏反变换式f(t)运行结果：ft=sin(t)注意:在MATLAB中，求拉氏反变换的函数ilaplace()，在默认情况下是指拉氏右变换，其运行结果是单边函数。如例②中的运行结果为ft=sin(t)，实际上是指ft=sin(t)。三、实验内容1.求出下列函数的拉氏变换式，并用MATLAB绘制拉氏变换在s平面的三维曲面图①3()2()5()ttftetet解：symstsft=sym('(-2*exp(-t)+5*exp(-3*t))*Heaviside(t)');Fs=laplace(ft)symsxyss=x+i*y;FFs=-2/(s+1)+5/(s+3);FFss=abs(FFs);ezmesh(FFss);ezsurf(FFss);colormap(hsv②()()(2)fttt解：symstsft=sym('Heaviside(t)-Heaviside(t-2)');Fs=laplace(ft)symsxyss=x+i*y;FFs=1/s-exp(-2*s)/s;FFss=abs(FFs);ezmesh(FFss);ezsurf(FFss);colormap(hsv);③3()sin()()tftett解：symstsft=sym('exp(-3*t)*sin(t)*Heaviside(t)');Fs=laplace(ft)symsxyss=x+i*y;FFs=1/((s+3)^2+1);FFss=abs(FFs);ezmesh(FFss);ezsurf(FFss);colormap(hsv);④()sin()()(2)ftttt解：symstsft=sym('sin(pi*t)*(Heaviside(t)-Heaviside(t-2))');Fs=laplace(ft)symsxyss=x+i*y;FFs=pi/(s^2+pi^2)-exp(-2*s)*pi/(s^2+pi^2);FFss=abs(FFs);ezmesh(FFss);ezsurf(FFss);colormap(hsv);2.已知信号的拉氏变换如下，请用MATLAB画出其三维曲面图，观察其图形特点，说出函数零极点位置与其对应曲面图的关系，并且求出它们所对应的原时间函数f(t)，①22(3)(3)()(5)(16)ssFsss解：symsxyss=x+i*y;FFs=(2*(s-3)*(s+3))/((s-5)*(s^2+16));FFss=abs(FFs);ezmesh(FFss);ezsurf(FFss);colormap(hsv);②(1)(3)()(2)(5)ssFssss解：symsxyss=x+i*y;FFs=((s+1)*(s+3))/((s+2)*(s+5)*s);FFss=abs(FFs);ezmesh(FFss);ezsurf(FFss);colormap(hsv);3.已知连续时间信号()s(2)()(4)ftcottt，请分别求出该信号的拉氏变换()Fs及其傅里叶变换()Fj，并用MATLAB绘出()Fs的曲面图及振幅频谱()Fj的波形，观察()Fs的曲面图在虚轴上的剖面图，并将它与信号的振幅频谱曲线进行比较，分析两者的对应关系。解：symstsft=sym('cos(2*pi*t)*(Heaviside(t)-Heaviside(t-4))');Fs=laplace(ft)symsxyss=x+i*y;FFs=s/(s^2+4*pi^2)-exp(-4*s)*s/(s^2+4*pi^2);FFss=abs(FFs);ezmesh(FFss);ezsurf(FFss);colormap(hsv);symstwGt=sym('cos(2*pi*t)*(Heaviside(t)-Heaviside(t-4))');Fw=fourier(Gt,t,w);FFw=maple('convert',Fw,'piecewise');FFP=abs(FFw);ezplot(FFP,[-10*pi10*pi]);grid;axis([-10*pi10*pi02.2])