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

实验五 连续信号与系统的S域分析

学院 班级 姓名 学号 一、实验目的

1. 熟悉拉普拉斯变换的原理及性质 2. 熟悉常见信号的拉氏变换

3. 了解正/反拉氏变换的MATLAB实现方法和利用MATLAB绘制三维曲面图的方法

4. 了解信号的零极点分布对信号拉氏变换曲面图的影响及续信号的拉氏变换与傅氏变换的关系

二、 实验原理

拉普拉斯变换是分析连续时间信号的重要手段。对于当t ∞时信号的幅值不衰减的时间信号，即在f(t)不满足绝对可积的条件时，其傅里叶变换可能不存在，但此时可以用拉氏变换法来分析它们。连续时间信号f(t)的单边拉普拉斯变换F(s)的定义为： 拉氏反变换的定义为：

显然，上式中F(s)是复变量s的复变函数，为了便于理解和分析F(s)随s的变化规律，我们将F(s)写成模及相位的形式：F(s)F(s)ej(s)。其中，|F(s)|为复信号F(s)的模，而(s)为F(s)的相位。由于复变量s=σ+jω，如果以σ为横坐标（实轴），jω为纵坐标（虚轴），这样，复变量s就成为一个复平面，我们称之为s平面。从三维几何空间的角度来看，|F(s)|和(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将其说明为符号表达式。具体方法参见第一部分第四章第三节。 例①：求出连续时间信号 f(t)sin(t)(t)的拉氏变换式，并画出图形 求函数拉氏变换程序如下：

syms t s %定义符号变量

ft=sym('sin(t)*Heaviside(t)'); %定义时间函数f(t)的表达式

1word版本可编辑.欢迎下载支持.

文档从网络中收集，已重新整理排版.word版本可编辑.欢迎下载支持.

Fs=laplace(ft) %求f(t)的拉氏变换式F(s)

运行结果：Fs = 1/(s^2+1)

绘制拉氏变换三维曲面图的方法有2种： 方法一：syms x y s

s=x+i*y; %产生复变量s

FFs=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，-6mesh(x,y,ffs); %绘制拉氏变换的三维网格曲面图 surf(x,y,ffs); %绘制带阴影效果的三维曲面图 axis([-5,5,-5,5,0,8]); %设置坐标显示范围

colormap(hsv); %设置图形中多条曲线的颜色顺序

说明：从拉普拉斯变换的三维曲面图中可以看出，曲面图上有象山峰一样突出的尖峰，这些峰值点在s平面的对应点就是信号拉氏变换的极点位置。而曲面图上的谷点则对应着拉氏变换的零点位置。因此，信号拉氏变换的零极点位置决定了其曲面图上峰点和谷点位置。 例②：求出函数F(s)1的拉氏反变换式 s21MATLAB程序如下：

syms t s %定义符号变量

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平面的三维曲面图 ① f(t)2e(t)5e(t)

解： syms t s

ft=sym('(-2*exp(-t)+5*exp(-3*t))*Heaviside(t)');

2word版本可编辑.欢迎下载支持.

t3t文档从网络中收集，已重新整理排版.word版本可编辑.欢迎下载支持.

Fs=laplace(ft) syms x y s s=x+i*y;

FFs=-2/(s+1)+5/(s+3); FFss=abs(FFs); ezmesh(FFss); ezsurf(FFss); colormap(hsv

② f(t)(t)(t2)

解： syms t s

ft=sym('Heaviside(t)-Heaviside(t-2)'); Fs=laplace(ft) syms x y s s=x+i*y;

FFs=1/s-exp(-2*s)/s; FFss=abs(FFs); ezmesh(FFss); ezsurf(FFss); colormap(hsv); ③ f(t)e3tsin(t)(t)

解： syms t s

ft=sym('exp(-3*t)*sin(t)*Heaviside(t)'); Fs=laplace(ft) syms x y s s=x+i*y;

FFs=1/((s+3)^2+1); FFss=abs(FFs); ezmesh(FFss); ezsurf(FFss); colormap(hsv);

④ f(t)sin(t)(t)(t2)

解： syms t s

ft=sym('sin(pi*t)*(Heaviside(t)-Heaviside(t-2))'); Fs=laplace(ft) syms x y s s=x+i*y;

FFs=pi/(s^2+pi^2)-exp(-2*s)*pi/(s^2+pi^2); FFss=abs(FFs);

3word版本可编辑.欢迎下载支持.

文档从网络中收集，已重新整理排版.word版本可编辑.欢迎下载支持.

ezmesh(FFss); ezsurf(FFss); colormap(hsv);

2. 已知信号的拉氏变换如下，请用MATLAB画出其三维曲面图，观察其图形特点，说出函数零极点位置与其对应曲面图的关系，并且求出它们所对应的原时间函数f (t)， ①F(s)2(s3)(s3) 2(s5)(s16)解：

syms x y s s=x+i*y;

FFs=(2*(s-3)*(s+3))/((s-5)*(s^2+16)); FFss=abs(FFs); ezmesh(FFss); ezsurf(FFss); colormap(hsv); ②F(s)(s1)(s3)

s(s2)(s5)解：

syms x y s s=x+i*y;

FFs=((s+1)*(s+3))/((s+2)*(s+5)*s); FFss=abs(FFs); ezmesh(FFss); ezsurf(FFss); colormap(hsv);

3. 已知连续时间信号f(t)cos(2t)(t)(t4)，请分别求出该信号的拉氏变换

F(s)及其傅里叶变换F(j)，并用MATLAB绘出F(s)的曲面图及振幅频谱F(j)的波

形，观察F(s)的曲面图在虚轴上的剖面图，并将它与信号的振幅频谱曲线进行比较，分析两者的对应关系。 解： syms t s

ft=sym('cos(2*pi*t)*(Heaviside(t)-Heaviside(t-4))'); Fs=laplace(ft) syms x y s s=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);

4word版本可编辑.欢迎下载支持.

文档从网络中收集，已重新整理排版.word版本可编辑.欢迎下载支持.

colormap(hsv); syms t w

Gt=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*pi 10*pi]);grid; axis([-10*pi 10*pi 0 2.2])

5word版本可编辑.欢迎下载支持.