北理工信号与系统实验报告.docxVIP

实验 2 LTI系统的时域分析 一、实验目的 ①掌握利用MATLAB对系统进行时域分析的方法。 ②掌握连续时间系统零状态响应、冲激响应和阶跃响应的求解方法。 ③掌握求解离散时间系统响应、单位抽样响应的方法。 ④加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。 二、实验原理 1、连续时间系统时域分析的MATLAB实现 1)连续时间系统的MATLAB表示 设LTI因果系统的微分方程一般式为： 则在MATLAB里，可以建立系统模型如下： b=[; a=[; sys=tf(b,a); 2）连续时间系统的零状态响应 用lsim(sys,x,t)表示求解零状态响应。 3）连续时间系统的冲激响应与阶跃响应。 用impulse函数来调用。 2、离散时间系统时域分析的MATLAB实现 1）离散时间系统的MATLAB表示。 LTI离散系统通常可以由系统差分方程描述; 则在MATLAB里，可以建立系统模型如下： b=[; a=[; 2）离散时间系统对任意输入的响应。 用filter（b,a,x）函数调用。 3）离散时间系统的单位抽样响应。 用impz函数来调用。 3、卷积和与卷积积分 1）离散时间序列的卷积和：调用格式为x=conv(x1,x2); 2）连续时间信号的卷积积分 连续时间信号x1(t)和x2(t)的卷积积分x(t)定义如下 三、实验内容 采用MATLAB绘出各系统的单位冲激响应和单位阶跃响应波形。 程序如下： b=[1]; a=[1 sqrt(2) 1]; sys=tf(b,a); subplot(121); impulse(sys); subplot(122); step(sys); 程序运行结果如下： 程序如下： b=[1 0 0]; a=[1 sqrt(2) 1]; sys=tf(b,a); subplot(121); impulse(sys); subplot(122); step(sys); 程序运行结果如下： 程序如下： b=[1 0]; a=[1 1 1]; sys=tf(b,a); subplot(121); impulse(sys); subplot(122); step(sys); 程序运行结果如下： 程序如下： b=[1 0 1]; a=[1 1 1]; sys=tf(b,a); subplot(121); impulse(sys); subplot(122); step(sys); 程序运行结果如下： 已知某系统可以由如下微分方程描述 利用MATLAB绘出该系统冲激响应和阶跃响应的时域波形。 程序如下： b=[1]; a=[1 1 6]; sys=tf(b,a); subplot(121); impulse(sys); subplot(122); step(sys); 程序运行结果如下： 根据冲激响应的时域波形分析系统的稳定性 观察上面的波形，t趋于无穷时，y(t)均趋于0，所以该系统是稳定的。 如果系统的输入为,求系统的零状态响应。 程序如下： b=[1]; a=[1 1 6]; sys=tf(b,a); t=0:0.01:10; x=exp(-1*t); lsim(sys,x,t); 程序运行结果如下： 已知描述离散系统的微分方程如下，用MATLAB绘出各系统的单位抽样响应，根据单位抽样响应的时域波形分析系统的稳定性。 程序如下： b=[1]; a=[1 3 2]; impz(b,a,0:10); 程序运行结果如下： 系统的稳定性 由图像可知，n趋于无穷时，y(n)越来越大，所以该系统是不稳定的。 程序如下： b=[1 -3]; a=[1 -0.5 0.8]; impz(b,a,0:10); 程序运行结果如下： 系统的稳定性 由图像可知，n趋于无穷时，y(n)越来越小且是有界的，所以该系统是稳定的。 已知系统可以由如下差分方程描述 程序如下： b=[1]; a=[1 1 0.25]; subplot(121); impz(b,a,0:10); subplot(122); stepz(b,a); 程序运行结果如下： 用MATLAB计算如下两个序列的卷积，并绘出图形。 程序如下： x1=[1 2 1 1]; x2=[1 1 1 1 1]; t=[-3 -2 -1 0 1 2 3 4]; x=conv(x1,x2); stem(t,x); 程序运行结果如下： 已知某LTI离散系统，其单位抽样响应，系统的输入为，计算当n=0,1,2，…，40时系统的零状态响应y(n),绘出x(n),h(n)和y(n)时域波形。 程序如下： dn=1 n1=0:dn:20 x1=sin(0.2*n1) subplot(221) st