信号与系统实验三【DOC精选】.doc

信号与系统实验报告 姓名：穆佳美 学号：096040156 实验七 拉普拉斯变换 1、利用mesh函数画出信号f(t)=sin(t)u(t)的拉普拉斯变换的曲面图。 实验代码： 实验图像： function shiyan7_1() %y=sin(t)u(t)laplace变化的曲面图 clc;clear all; x=-3:0.1:3; y=-3:0.1:3; [x y]=meshgrid(x,y); z=x+i*y; z=1./(1+z.^2); mesh(x,y,z); surf(x,y,z);axis([-3 3 -3 3 -4 4]) 2、利用meshgrid、mesh、surf函数画出信号f(t)= u(t)-u(t-2)的拉普拉斯变换的曲面图，观察曲面图在虚轴剖面上的曲线，并将其与信号傅里叶变换绘制的振幅频谱进行比较。 实验代码： function shiyan7_2() %y=u(t)-u(t-2)laplace变化的曲面图,与其fourier变换的比较 clc;clear all; x=0:0.1:5; y=-20:0.1:20; [x y]=meshgrid(x,y); z=x+i*y; z=(1-exp(-2*z))./z; subplot(211);mesh(x,y,z);surf(x,y,z);view(-70,10);title(laplace变换); %fourier变换 w=-20:0.1:20; Fw=(2*sin(w).*exp(-j*w))./w; subplot(212);plot(w,(Fw));title(fourier变换); 实验图像： 3、画出的曲面图，观察拉普拉斯变换的零极点。 实验代码： 实验图像： function shiyan7_3() %H(s)=2(s-3)(s+3)/(s-5)/(s^2+10)的曲面图 clc;clear all; x=-5:0.1:8; y=-5:0.1:5; [x y]=meshgrid(x,y); z=x+i*y; z=2.*(z-3).*(z+3)./(z-5)./(z.^2+10); mesh(x,y,z);surf(x,y,z);view(-60,10); 利用roots函数画出a.和的零极点图。 实验代码： 实验图像： function shiyan7_4() %画系统零极点图 A1=[1 2 -3 2 1];B1=[1 0 -4]; ljdt(A1,B1);title(连续系统零极点图1); A2=[1 5 16 30];B2=[5 20 25 0]; figure;ljdt(A2,B2);title(连续系统零极点图1); function [p q]=ljdt(A,B) %用于画系统零极点图 %p 系统极点位置行向量 %q 系统零点位置行向量 %A极点B零点 p=roots(A); q=roots(B); x=max(abs([p q])); x=x+0.1;y=x; hold on; axis([-x x -y y]);%确定坐标轴位置 axis(square); plot(real(p),imag(p),x);%画极点 plot(real(q),imag(q),o);%画零点 alpha=0:pi/30:2*pi;%角度[0,2*pi] x=cos(alpha); y=sin(alpha); plot(x,y,-.); 5、已知拉普拉斯变换，利用residue函数求其拉普拉斯逆变换。 6、已知系统函数为，利用residue函数求该系统的冲击响应h(t)，并利用impulse函数画出其时域波形，判断系统的稳定性。 实验代码： function shiyan7_5() clc; a1=[2 4];b1=[1 0 4 0]; bffs(a1,b1) a2=[1 4];b2=[1 3 2 0]; bffs(a2,b2) impulse(a2,b2); function y=bffs(a,b) %部分分式 [r,p,k]=residue(a,b); n=length(r); syms t ; for i=1:n y(i)=r(i)*exp(p(i)*t); end if(length(k)~=0) y(n+1)=k(1)*dirac(t); end 实验结果： ans = [ (-1/2-1/2*i)*exp(2*i*t), (-1/2+1/2*i)*exp(-2*i*t), 1] ans = [ exp(-2*t), -3*