信号与系统实验 第八章.docVIP

一、实验目的 1.用MATLAB观察LTI连续时间系统的零状态响应 2.用MATLAB观察系统冲激响应，阶跃响应的数值解，波形 3.用MATLAB观察差分方程描述的离散系统 4用MATLAB观察卷积积分 二、实验内容 1.求解一定激励微分方程的零状态响应，并与用经典求法得到的结果比较 2求解微分方程在一定时间范围内的冲击响应和阶跃响应 3求解离散系统差分方程在一定激励下的零状态响应，并与用经典求法得到的结果比较 4.求解差分方程描述的了离散系统在一定时间内的单位系列响应和阶跃响应 5求解已知图形的卷积积分 三、主要算法与程序8.1题(1) a=[1 3 2]; b=[2 6]; sys=tf(b,a); p=0.1; t=0:p:10; f=1.*u(t); lsim(sys,f,t); 8.2题 8.1第（1）题方程解得如下 t=[0:0.01:10]; plot (-4*exp(-t)+exp(-2*t)+3) 8.3题（3） a=[1 4 5] b=[1 0] subplot(2,1,1) impulse(b,a) subplot(2,1,2) step(b,a) 8.4题 a=[1 2 1]; b=[ 1]; n=0:20; x=0.25.^n; y=filter(b,a,x) subplot(2,1,1) stem(n,x,*) title(输入) subplot(2,1,2) stem(n,y,*) title(输出) 8.5解得 n=-1:21; y=(1/6).*(-1).^n+(2/3).*(1/4).^n.*jyxl(n); stem(n,y,*) 8.6(3题) a=[1 1 0.25] b=[1] x=jyxl(n); subplot(2,1,2) y=filter(b,a,x) stem(n,y,*) title(单位阶跃响应) subplot(2,1,1) impz(b,a,21) title(单位系列响应) 8.7题 x1=[0,1,2,1,0,0]; n1=-2:3; x2=[0,1,1,1,1,0,0]; n2=[-1:5]; [x,n]=gghconv(x1,x2,n1,n2) 8.9题(5) t4=-3:0.01:3; f4=0.5*(t4-2).*H(t4-2)-t4.*H(t4)+0.5*(t4+2).*H(t4+2); t3=t4; f3=0.5*t3.*(H(t3)-H(t3-2)) [t,f]=gggfconv(f4,f3,t4,t3); 四、实验结果与分析 8.1 8.2 8.3 8.4 8.5 8.6 8.7 8.9 请输入取样时间间隔d: 0.01 五、实验小结 本次试验对连续时间系统和离散时间系统分别进行了观察 用经典方法求的的结果与用MATLAB做出的图形相符，进一步通过图形直观的观察了结果 对函数Lsim（系统输出信号仿真）impulse（连续系统冲击响应的时域波形），impz(离散系统单位系列响应)以及filter(对指定时间范围的输入序列产生的响应)进行了学习 利用自己编写的函数求已知图形的卷积积分