《信号与线性系统实验报告.docx》由会员分享，可在线阅读，更多相关《信号与线性系统实验报告.docx（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 信号与线性系统实验报告5 连续信号与系统的复频域分析及MATLAB实现 一、 试验目的 1、 把握MATLAB 实现连续时间信号拉普拉斯变换及逆变换的方法。 2、 把握MATLAB 绘制拉普拉斯变换的三维曲面图，并分析复频域特性和时 移特性。 二、 试验原理及学问要点 1、 连续时间非周期信号的拉普拉斯变换及逆变换（laplace( )及ilaplace( )函数）； 2、 拉普拉斯变换的数值算法； 3、 绘制拉普拉斯变换的三维曲面图（meshgrid()及mesh()函数） 三、 试验软件： MATLAB软件 四、 试验内容及试验记录 12.1 利用MATLAB的laplace函数，求以下
2、信号的拉普拉斯变换。 （1） syms t; 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）syms s; 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）
3、a=1 1; b=1 5 6 0; k,p,c=residue(a,b) 运行的结果为： k = -0.6667 0.5000 0.1667 p = -3.0000 -2.0000 0 c = 由上述程序的运行结果知，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、观看分析信号的幅频特性。 （4） f(t)=exp(-t)*cos(pi*t/2)*u(t) 其对应的拉氏变换为： （1） syms t; F=exp(-t)*cos(pi/2*t); L=laplace(F) L = (s + 1)/(s + 1)2 + pi2/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(”单边指数信号拉普
5、拉斯变换幅度曲面图”); xlabel(”实轴”) ylabel(”虚轴”) 单边指数信号拉普拉斯变换幅度曲面图 20 15 10 5 0 210-1虚轴-2-1-0.5实轴00.5 由该曲面可直观看出，曲面图中有两个峰点位置P1,2 =-1pi/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); s
6、ubplot(2,3,1); plot(t,y,”b”); set(gca,”color”,1 1 1); set(gca,”XColor”,0 0 0); set(gca,”YColor”,0 0 0); set(gca,”ZColor”,0 0 0); title(”u(t)的时域波形 ”,”Color”,0 0 0); axis(-2,60,-0.2,2); b=0 1; a=1 0; subplot(2,3,2); lapulas1(b,a,t0); title(”拉氏变换幅度曲面图”,”Color”,0 0 0); subplot(2,3,3); lapulas2(b,a,t0);
7、title(”u(t)拉氏变换相位曲面图”,”Color”,0 0 0); hold on 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”,0 0 0); hold on subplot(2,3,5); lapulas1(b,a,t0); title(”u(t-t0)的拉氏变换曲面图”,”Color”,0 0 0); hold on subplot(2,3,6)
8、; lapulas2(b,a,t0); title(”u(t-t0)拉氏变换相位曲面图”,”Color”,0 0 0); u(t)的时域波形 21.510.500204060拉氏变换幅度曲面图u(t)拉氏变换相位曲面图105050-500.5420-2-450-500.511 u(t-t0)的时域波形21.510.500204060u(t-t0)的拉氏变换曲面图u(t-t0)拉氏变换相位曲面图105050-500.5420-2-450-500.511 当输入t0=6时，由上述的波形图知信号在时间轴的平移后，其幅度曲面图按指数规律转变，而相位曲面图则发生线性转变。 12.7 已知连续时间信号f(
9、t)=exp(-2*t)*u(t). (1)求f（t）的拉普拉斯变换F（s）及傅里叶变换F(jw): syms t w; 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; F
10、w=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.3a,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 （
11、3）观看比拟信号振幅频谱曲线与拉普拉斯变换幅度曲面图在虚轴上的剖面曲线的关系， 分析频域与复频域的对应关系与规律。 通过上述两图的拉氏变换的虚轴剖面与幅度频谱曲线的比照分析知，幅度频谱曲线是在拉氏变换F(s)的曲面图中当令s=jw所得，即F(jw)=F(s)|s=jw Pole-Zero Map12.8 已知连续时间系统的系统函数H2.5（s）分别如下： 2（4）H(s)=(s2-4)/(s2+4) 1.5(1)利用MATLAB绘出系统的零、极点分1布图，推断系统的稳定性； 0.5num=1 0 -4; 0den=1 0 4; -0.5H=tf(num,den); -1pzmap(H) -1.
12、5p,z=pzmap(H) -2 运行结果为： Imaginary Axis (seconds-1)-2.5-2-1.5-1-0.500.5-111.52Real Axis (seconds)p = 0 + 2.0000i 0 - 2.0000i z = 2.0000 -2.0000 由分布图观看知，系统有位于虚轴上的一对共轭虚极点，故该系统不满意系统稳定的复频域条件，是不稳定系统。 (1) 求出系统的冲击响应h(t)并绘出其时域波形。 系统的冲击响应为： syms s; H(s)=(s2-4)/(s2+4); f=ifourier(H(s) 运行结果为： Impulse Response54
13、321Amplitude0-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=1 0 -4; den=1 0 4; impulse(num,den); 50100150202350Time (seconds) 五、试验分析及试验小结 1.编程时遇到函数源代码不存在导致无法显示结果。 解决方法：从书上找到源代码手动输入保存为脚本文件。 2.编程时还有图形参数设置不当，导致图形显示有问题。 解决方法：不断调试范围，直到显示正常。2、试验参数对显示图形结果有重要影响，若参数不适宜，图形结果不能表达出函特点，甚至不正确。 3.试验结果与理论计算值相吻合。