信号与线性系统试验报告5

本文格式为Word版，下载可任意编辑——信号与线性系统试验报告5连续信号与系统的复频域分析及MATLAB实现

一、

试验目的

1、把握MATLAB实现连续时间信号拉普拉斯变换及逆变换的方法。2、把握MATLAB绘制拉普拉斯变换的三维曲面图，并分析复频域特性和时

移特性。

二、试验原理及知识要点

1、连续时间非周期信号的拉普拉斯变换及逆变换（laplace()及ilaplace()函数）；2、拉普拉斯变换的数值算法；

3、绘制拉普拉斯变换的三维曲面图（meshgrid()及mesh()函数）三、试验软件：MATLAB软件四、试验内容及试验记录

12.1利用MATLAB的laplace函数，求以下信号的拉普拉斯变换。（1）symst;

F=(1-exp(-0.5*t))*Heaviside(t);L=laplace(F)

运行的结果为：L=

1/s-1/(s+1/2)

12.2利用MATLAB的ilaplace函数，求以下像函数F(s)的拉普拉斯逆变换。（1）symss;

L=(s+1)/(s*(s+2)*(s+3));F=ilaplace(L)运行的结果：F=

1/6+1/2*exp(-2*t)-2/3*exp(-3*t)

12.3利用MATLAB的residue函数求12.2题中（1）小题的拉普拉斯逆变换，并与ilaplace函数的计算结果进行比较。（1）a=[11];b=[1560];

[k,p,c]=residue(a,b)运行的结果为：k=

-0.66670.50000.1667p=

-3.0000-2.00000c=[]

由上述程序的运行结果知，F(s)有三个单实极点，部分分式展开结果：

F(s)=(-2/3)/(s+3)+0.5/(s+2)+(1/6)/s则拉普拉斯逆变换：

f(t)=(-2/3e^(-3t)+0.5e^(-2t)+1/6)u(t)

用residue函数求出的结果与用ilaplace函数求出的结果是一样的。只是后者简单点。

12.4试用MATLAB绘出以下信号拉普拉斯变换的三维曲面图，并通过三维曲面图观测分析信号的幅频特性。

（4）f(t)=exp(-t)*cos(pi*t/2)*u(t)其对应的拉氏变换为：（1）symst;

F=exp(-t)*cos(pi/2*t);L=laplace(F)

L=

(s+1)/((s+1)^2+pi^2/4)曲面图及代码为：

x=-1:0.08:0.2;y=-2:0.08:2;

[x,y]=meshgrid(x,y);s=x+i*y;

F=abs(4./pi.^2.*(s+1)./(4.*(s+1).^2./pi.^2+1));mesh(x,y,F);surf(x,y,F)

colormap(hsv);

title('单边指数信号拉普拉斯变换幅度曲面图');xlabel('实轴')ylabel('虚轴')

单边指数信号拉普拉斯变换幅度曲面图

20

15

10

5

0210-1虚轴-2-1-0.5实轴00.5

由该曲面可直观看出，曲面图中有两个峰点位置P1,2=-1±pi/2,这正是单边余弦信号cos(t)*u(t)的拉普拉斯变换曲面图沿S平面实轴方向平移-1的结果。

12.6试用MATLAB分别绘出单位阶跃信号u（t）及其时间平移信号u(t-t0)的时域波形和拉普拉斯变换三维曲面图，交互式地改变t0的大小，观测分析拉普拉斯变换的时移特性。

t0=0;t=-2:0.001:60;

y=Heaviside(t);subplot(2,3,1);

plot(t,y,'b');set(gca,'color',[111]);set(gca,'XColor',[000]);

set(gca,'YColor',[000]);set(gca,'ZColor',[000]);

title('u(t)的时域波形','Color',[000]);axis([-2,60,-0.2,2]);

b=[01];a=[10];subplot(2,3,2);

lapulas1(b,a,t0);

title('拉氏变换幅度曲面图','Color',[000]);subplot(2,3,3);

lapulas2(b,a,t0);

title('u(t)拉氏变换相位曲面图','Color',[000]);holdon

t0=input('请输入信号尺度变换因子t0:')subplot(2,3,4);t=-8:0.001:60;y=Heaviside(t-t0);plot(t,y,'b');

axis([-8,60,-0.2,2]);

title('u(t-t0)的时域波形','color',[000]);holdon

subplot(2,3,5);

lapulas1(b,a,t0);title('u(t-t0)的拉氏变换曲面图','Color',[000]);holdon

subplot(2,3,6);

lapulas2(b,a,t0);title('u(t-t0)拉氏变换相位曲面图','Color',[000]);

u(t)的时域波形21.510.500204060拉氏变换幅度曲面图u(t)拉氏变换相位曲面图105050-500.5420-2-450-500.511u(t-t0)的时域波形21.510.500204060u(t-t0)的拉氏变换曲面图u(t-t0)拉氏变换相位曲面图105050-500.5420-2-450-500.511

当输入t0=6时，由上述的波形图知信号在时间轴的平移后，其幅度曲面图按指数规律改变，，而相位曲面图则发生线性改变。12.7已知连续时间信号f(t)=exp(-2*t)*u(t).(1)求f（t）的拉普拉斯变换F（s）及傅里叶变换F(jw):

symstw;

f=exp(-2*t)*sym('heaviside(t)');L=laplace(f)F=fourier(f)L=1/(s+2)F=

1/(w*i+2)

（1）用MATLAB分别绘出上述信号的拉

傅里叶变换（幅度频谱曲线）0.5普拉斯变换幅度曲面图|F(s)|和振幅频谱曲线|F(jw)|；其波形和代码如下：

0.450.40.350.30.250.20.150.10.050-20-15-10-50频率?5101520w=-20:0.1:20;Fw=1./(w*i+2);

plot(w,abs(Fw))title('傅里叶变换（幅度频谱曲线）')xlabel('频率')pause

0.5clf;

0.4a=-0:0.1:5;

b=-20:0.1:20;

0.3[a,b]=meshgrid(a,b);

c=a+i*b;0.2c=1./(c+2);

0.1c=abs(c);

mesh(a,b,c);020surf(a,b,c);

view(-60,20)axis([-0,5,-20,20,0,0.5]);

title('拉普拉斯变换幅度曲面图');colormap(hsv);

拉普拉斯变换幅度曲面图41002-10-200

（3）观测比较信号振幅频谱曲线与拉普拉斯变换幅度曲面图在虚轴上的剖面曲线的关系，

分析频域与复频域的对应关系与规律。

通过上述两图的拉氏变换的虚轴剖面与幅度频谱曲线的对比分析知，幅度频谱曲线是在拉氏变换F(s)的曲面图中当令s=jw所得，即F(jw)=F(s)|s=jw

Pole-ZeroMap12.8已知连续时间系统的系统函数H2.5（s）分别如下：2（4）H(s)=(s^2-4)/(s^2+4)1.5(1)利用MATLAB绘出系统的零、极点分1布图，判断系统的稳定性；0.5num=[10-4];0den=[104];-0.5H=tf(num,den);-1pzmap(H)-1.5[p,z]=pzmap(H)-2

运行结果为：

ImaginaryAxis(seconds-1)-2.5-2-1.5-1-0.500.5-111.52RealAxis(seconds)p=

0+2.0000i0-2.0000iz=

2.0000-2.0000

由分布图观测知，系统有位于虚轴上的一对共轭虚极点，故该系统不满足系统稳定的复频域条件，是不稳定系统。

(1)求出系统的冲击响应h(t)并绘出其时域波形。

系统的冲击响应为：symss;

H(s)=(s^2-4)/(s^2+4);f=ifourier(H(s)

运行结果为：

ImpulseResponse54321Amplitude0-1-2-3f=-4-(4*pi*heaviside(-x)*exp(2*x)-

-502*pi*exp(-2*x)*dirac(x)-2*pi*exp(2*x)*dirac(x)+(pi*exp(-2*x)*dirac(x,1))/2-

(pi*exp(2*x)*dirac(x,1))/2+4*pi*exp(-2*x)*heaviside(x))/(2*pi)其对应的时域波形及代码如下：num=[10-4];den=[104];

impulse(num,den)