函数程序代码： 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); 运行结果截图：
求原函数的程序代码： syms t s
Fs =sym('[2*(s-3)*(s+3)]/[(s-5)*(s^2+16)]') ft=ilaplace(Fs) 原函数：
ft = (50*cos(4*t))/41 + (32*exp(5*t))/41 + (125*sin(4*t))/82 ②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*(s+2)*(s+5)]; FFss=abs(FFs); ezmesh(FFss); ezsurf(FFss); colormap(hsv); 运行结果截图：
求原函数的程序代码： syms t s
Fs =sym('[(s+1)*(s+3)]/[s*(s+2)*(s+5)]') ft=ilaplace(Fs) 原函数：
ft = exp(-2*t)/6 + (8*exp(-5*t))/15 + 3/10
3. 已知连续时间信号f(t)?cos(2?t)??(t)??(t?4)?，请分别求出该信号的拉氏变换F(s)及其傅里叶变换F(j?)，并用MATLAB绘出F(s)的曲面图及振幅频谱F(j?)的波形，观察F(s)的曲面图在虚轴上的剖面图，并将它与信号的振幅频谱曲线进行比较，分析两者的对应关系。
(1)拉氏变换： 程序代码： syms t s
ft=sym('cos(2*pi*t)*(Heaviside(t)-Heaviside(t-4))');
Fs=laplace(ft) 运算结果：
F(s)=laplace(Heaviside(t), t, s-pi*2i)/2 + laplace(Heaviside(t), t, s+pi*2i)/2 - laplace(Heaviside(t - 4), t, s - pi*2i)/2 - laplace(Heaviside(t - 4), t, s + pi*2i)/2
(2)傅里叶变换： 程序代码： syms t w
Gt=sym('cos(2*pi*t)*(Heaviside(t)-Heaviside(t-4))'); Fw=fourier(Gt) 运算结果： F(w) = fourier(cos(2*pi*t)*Heaviside(t), t, w) - fourier(Heaviside(t - 4)*cos(2*pi*t), t, w)
四、总结报告
由于平时都是在上课过程中学习理论知识，而这次实验是在理论知识的基础上来进行实验操作，但并是不全是上课时学习的理论知识，也存在许多的新知识。所以对于这次把上课时的理论知识并结合新知识一起应用于实践操作来说是有点困难的。
信号与系统的实验不同于大物实验，一开始说可以多人合作完成的实验，到最后是一个人单独完成。在为数不多的四次实验中，我深深感受到了团队合作在实验中的重要性。在自己对自己写出的代码，运行出现错误的时候，两个人或者多个人对实验的共同理解是实验高效、误差小完成的基础。
参考文献：
《MATLAB应用大全》 《MATLAB无师自通》