信号与系统实验报告实验五.docx
《信号与系统实验报告实验五.docx》由会员分享,可在线阅读,更多相关《信号与系统实验报告实验五.docx(12页珍藏版)》请在冰点文库上搜索。
信号与系统实验报告实验五
信号与系统-实验报告-实验五
实验五 连续信号与系统的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)的模,而
为F(s)的相位。
由于复变量s=σ+jω,如果以σ为横坐标(实轴),jω为纵坐标(虚轴),这样,复变量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)
②拉氏反变换
%其中矩阵s包含了复平面-6<σ<6,-6fs=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平面的对应点就是信号拉氏变换的极点位置。
而曲面图上的谷点则对应着拉氏变换的零点位置。
因此,信号拉氏变换的零极点位置决定了其曲面图上峰点和谷点位置。
例②:
求出函数
的拉氏反变换式
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平面的三维曲面图
①
解:
symsts
ft=sym('(-2*exp(-t)+5*exp(-3*t))*Heaviside(t)');
Fs=laplace(ft)
symsxys
s=x+i*y;
FFs=-2/(s+1)+5/(s+3);
FFss=abs(FFs);
ezmesh(FFss);
ezsurf(FFss);
colormap(hsv
②
解:
symsts
ft=sym('Heaviside(t)-Heaviside(t-2)');
Fs=laplace(ft)
symsxys
s=x+i*y;
FFs=1/s-exp(-2*s)/s;
FFss=abs(FFs);
ezmesh(FFss);
ezsurf(FFss);
colormap(hsv);
③
解:
symsts
ft=sym('exp(-3*t)*sin(t)*Heaviside(t)');
Fs=laplace(ft)
symsxys
s=x+i*y;
FFs=1/((s+3)^2+1);
FFss=abs(FFs);
ezmesh(FFss);
ezsurf(FFss);
colormap(hsv);
④
解:
symsts
ft=sym('sin(pi*t)*(Heaviside(t)-Heaviside(t-2))');
Fs=laplace(ft)
symsxys
s=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),
①
解:
symsxys
s=x+i*y;
FFs=(2*(s-3)*(s+3))/((s-5)*(s^2+16));
FFss=abs(FFs);
ezmesh(FFss);
ezsurf(FFss);
colormap(hsv);
②
解:
symsxys
s=x+i*y;
FFs=((s+1)*(s+3))/((s+2)*(s+5)*s);
FFss=abs(FFs);
ezmesh(FFss);
ezsurf(FFss);
colormap(hsv);
3.已知连续时间信号
,请分别求出该信号的拉氏变换
及其傅里叶变换
,并用MATLAB绘出
的曲面图及振幅频谱
的波形,观察
的曲面图在虚轴上的剖面图,并将它与信号的振幅频谱曲线进行比较,分析两者的对应关系。
解:
symsts
ft=sym('cos(2*pi*t)*(Heaviside(t)-Heaviside(t-4))');
Fs=laplace(ft)
symsxys
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);
colormap(hsv);
symstw
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*pi10*pi]);grid;
axis([-10*pi10*pi02.2])