实验四 LTI系统的时域分析.pptVIP

实验五 离散LTI系统时域分析 离散时间序列卷积和 连续系统的阶跃响应和冲击响应 离散系统的单位响应 实验目的 了解离散时间序列卷积和、离散系统单位序列响应的Matlab实现方法； 了解函数conv()、impz( )的调用格式及作用； 了解差分方程迭代解法的Matlab实现方法； 通过该实验，掌握离散LTI系统的时域基本分析方法及编程思想。 实验内容 1、验证性实验 采用函数conv()编程，实现离散时间序列的卷积和运算（或系统的零状态响应）； 采用差分方程的迭代解法，求离散时间系统的全响应； 采用函数impz( )编程，求离散时间系统的单位序列响应； 2、程序设计实验 预备知识 离散序列卷积和序列号的确定 采用函数conv()，可以快速求出两个离散时间序列的卷积和，但是此函数不需要给出两序列的对应的时间序列号，也不返回卷积和序列对应的序列号。因此需要讨论卷积和序列对应的序列号的确定问题。 函数impz()调用格式 impz(b,a) 该调用格式以默认方式绘出向量a，b定义的离散系统的单位序列响应的离散时间波形。 impz(b,a,n) 该调用格式将绘出由向量a,b定义的离散系统在0~n（n必须为整数）离散时间范围内的单位序列响应的时域波形。 impz(b,a,n1:n2) 该调用格式将绘出由向量a,b定义的离散系统在n1~n2（n1,n2 必须为整数,且n1n2）离散时间范围内的单位序列响应的时域波形。 y=impz(b,a,n1:n2) 该调用格式并不绘出系统单位序列响应的时域波形，而是求出向量a，b定义的离散系统在n1~n2（n1,n2必须为整数,且n1n2）离散时间范围内的单位序列响应的数值解。 实验要求 在计算机中输入程序，验证实验结果，并将实验结果存入指定存储区域； 对于程序设计实验，要求通过对验证性实验的练习，自行编制完整实验程序，并得出实验结果。 在实验报告中写出完整的自编程序，并给出实验结果。 实验步骤 一）验证性实验1．离散时间序列卷积和（零状态响应） 计算两离散序列卷积和的实用程序 function[f,k]=lsjuanji (f1,f2,k1,k2) %The function of compute f=f1*f2 % f:卷积和序列f(k)对应的值 % k:序列f(k)对应的序列号 % f1:序列f1(k) % f2:序列f2(k) % k1:序列f1(k)对应的序列号 % k2:序列f2(k)对应的序列号 f=conv(f1,f2) k0=k1(1)+k2(1); %序列f(k)对应的起始序列号 k3=length(f1)+length(f2)-2; %序列f(k)的宽度 k=k0:k0+k3; %序列f(k)对应的序列号向量 subplot(2,2,1) stem(k1,f1) title(f1(k)) xlabel(k) ylabel(f1(k)) subplot(2,2,2) stem(k2,f2) title(f2(k)) xlabel(k) ylabel(f2(k)) subplot(2,2,3) stem(k,f) title(f(k)) xlabel(k) ylabel(f(k)) y0=0 % 初值y(0)=0 y(1)=1 y(2)=5/6*y(1)-1/6*y0+cos(2*pi/4)/6 for k=3:20 y(k)=5/6*y(k-1)-1/6*y(k-2)+cos(k*pi/4)/6 y(k) end yy=[y0 y(1:20)] % 取y(k)从y(0)到y(20) k=1:21 stem(k-1,yy) grid on xlabel(k) ylabel(y(k)) 3、离散系统的单位序列响应 某离散LTI系统的差分方程为， 则对应的向量为： 试绘出此系统的单位序列响应的波形。 程序如下a=[1,-1,0.9]b=[1]impz(b,a)结果如图 二）程序设计实验 1．已知离散LTI系统，激励 单位序列响应 画出系统的零状态响应在有限区间的波形。（有限区间自行设定） 2．已知两离散序列，试画出两序列的卷积和波形。 3．描述LTI离散系统的差分方程如下，请绘出该系统在0~50时间范围内单位序列响应的波形，并求出数值解。 * 下一页 上一页 * 计算两序列的卷积和， 对应的 对应的 程序如下： f1=[1 2 1] k1=[-1,0,1] f2=ones(1,5) k2=-2:2 [f,k]=lsjuanji (f1,f2,k1,k2) 试采用上述函数M文件编程，完成两序列的卷积和