一、实验目的
1、学会运用MATLAB求拉普拉斯变换(LT);
2、学会运用MATLAB求拉普拉斯反变换(ILT)。
二、实验原理及实例分析
LT是分析连续信号与系统的重要方法。运用LT可以将连续LTI系统的时域模型
简便地进行变换,经求解在还原为时域解。从数学角度看,LT是求解常系数线性
微分方程的工具。由LT导出的系统函数对系统特性分析也具有重要意义。
(一)拉普拉斯变换(LT)
对于一些不满足绝对可积条件的时域信号,是不存在傅里叶变换的。为了使更多的函数存在变换,并简化某些变换形式或运算过程,引入衰减因子 =,其中,为任意实数,使得满足绝对可积条件,从而求的傅里叶变换,即把频域扩展为复频域。
连续时间信号的LT定义为:
ILT定义为:
上式构成了拉普拉斯变换对,称为的像函数,而称为的原函数。可以将拉普拉斯变换理解为广义的傅里叶变换。
考虑到实际问题,人们用物理手段和实验方法所能记录和处理的一切信号都是有起始时刻的,对于这类单边信号或因果信号,我们引入单边LT,定义为:
如果连续信号可用符号表达式表示,则可用MATLAB的符号数学工具箱中的laplace函数来实现其单边LT,其语句格式为:。式中L返回的是默认符号为自变量s的符号表达式,f则为时域符号表达式,可通过sym函数来定义。
2、基于MATLAB符号数学工具箱实现ILT
如果连续信号可用符号表达式表示,则可用MATLAB的符号数学工具箱中的 ilaplace函数来实现其ILT,其语句格式为:。
式中f返回的是默认符号为自变量t的符号表达式,L则为s域符号表达式,也可通过sym函数来定义。
三、实验内容
1、完成以下例题:
例1:用MATLAB的laplace函数求的FT。
实验代码:
syms a t s;
ft = exp(-1*t)*sin(a*t)*heaviside(t);
Fs = laplace(ft, s);
实验结果:
>> Fs
Fs = a/((s + 1)^2 + a^2)
例2:用MATLAB的ilaplace函数求的ILT。
实验代码:
syms t s;
Fs = s^2/(s^2 + 1);
ft = ilaplace(Fs, t);
实验结果:
>> ft
ft = dirac(t) - sin(t)
2、完成以下练习:
(1)利用MATLAB求信号的LT。
实验代码:
syms t s;
ft = (exp(-3*t) - exp(-5*t))/t;
Fs = laplace(ft, s);
实验结果:
>> Fs
Fs = log(s + 5) - log(s + 3)
(2)利用MATLAB求的ILT,并画出其原函数波形图
实验代码:
syms t s;
Fs = exp(-s)/(4*s*(s^2 + 1));
ft = ilaplace(Fs, t);
fplot(ft, [-2*pi, 5*pi+1]); grid on;
title('f(t) Graph');
xlabel('t');
ylabel('(u(t)*(1 - cos(t - 1)))/4');
实验结果:
>> ft
ft = -( heaviside(t - 1)*(cos(t - 1) - 1))/4
实验图像:
四、实验总结
总的说来,实验是一个发现问题并解决问题的过程,只有经历过问题的困扰,才能收获解决时的喜悦,拉普拉斯变化是解决零输入和零状态响应的又一种巧妙方法,因此在此实验基础上熟练应用拉普拉斯变换才是最大的目的。