实验六 连续LTI系统时域分析.pptVIP

实验六 连续LTI系统的时域分析 实验目的 了解连续时间信号的卷积积分、连续系统的阶跃响应和冲激响应以及零输入响应的Matlab实现方法。 了解函数conv()、impulse()、 step() 的调用格式及作用。 通过该实验，掌握连续LTI系统的时域基本分析方法及编程思想。 实验内容 1、验证性实验 采用Matlab语言编程，实现连续时间信号的卷积积分运算（或系统的零状态响应）； 采用 step()、impulse()两函数编程，实现连续LTI系统的单位阶跃响应和单位冲激响应； 采用Matlab语言编程，求连续LTI系统的零输入响应。 2、程序设计实验。 函数impulse()调用格式 impulse(b,a,t) 此调用格式将绘出由向量a和b定义的连续系统在0~t时间范围内的冲击响应的时域波形。 impulse(b,a,t1:p:t2) 此调用格式将绘出由向量a和b定义的连续系统在t1~t2时间范围内，以时间间隔p均匀取样的冲击响应的时域波形。 y= impulse(b,a,t1:p:t2) 此调用格式将绘出由向量a和b定义的连续系统在t1~t2时间范围内，以时间间隔p均匀取样的冲击响应的数值解。 step() 函数有相同的调用格式 step(b,a) step(b,a,t) step(b,a,t1:p:t2) y=step(b,a,t1:p:t2) 实验要求 在计算机中输入程序，验证实验结果，并将实验结果存入指定存储区域。 对于程序设计实验，要求通过对验证性实验的练习，自行编制完整实验程序，并得出实验结果。 在实验报告中写出完整的自编程序，并给出实验结果。 实验步骤 一）验证性实验 1．连续时间信号的卷积积分（零状态响应） 实现连续信号卷积的通用函数M文件如下 function [f,k]=lxjuanji(f1,f2,k1,k2,p) % f: 卷积积分f(t) % k:卷积积分f(t)对应的时间向量 % f1:信号f1(t)的向量值 % f2:信号f2(t)的向量值 % k1:f1(t)对应的时间向量 % k2:f2(t)对应的时间向量 % p:取样间隔 f=conv(f1,f2) f=f*p; k0=k1(1)+k2(1); % 计算f(t)的起点位置 k3=length(f1)+length(f2)-2 % 计算f(t)的宽度 k=k0:p:k3*p % 计算f(t)的对应的时间向量 subplot(2,2,1) plot(k1,f1) title(f1(t)) xlabel(t) ylabel(f1(t)) subplot(2,2,2) plot(k2,f2) title(f2(t)) xlabel(t) ylabel(f2(t)) subplot(2,2,3) plot(k,f) title(f(t)) xlabel(t) ylabel(f(t)) 1．1已知两连续信号如下所示，试采用以上通用函数M文件，画出卷积结果的波形。 1．2 已知连续LTI系统，激励信号和系统的单位冲激响应如下所示，试直接采用Matlab语言编程，画出该系统零状态响应的波形。 程序如下 t1=1 t2=2 t3=0 t4=1 T=0.01 t=0:T:t2+t4 f1=ones(size(t)).*((tt1)-(tt2)) f2=ones(size(t)).*((tt3)-(tt4)) yf=conv(f1,f2)*T subplot(3,1,1) plot(t,f1) title(f1(t)) subplot(3,1,2) plot(t,f2) title(f2(t)) subplot(3,1,3) plot(t,yf(1:(t2+t4)/T+1)) title(yf(t)) 1．3已知连续LTI系统，激励信号和系统的单位冲激响应分别如下所示，试采用符号函数编程的方法，确定该系统零状态响应的数学表达式。 程序如下 syms x t=sym(t,positive) % 设定限制性符号变量 f=exp(-t) h=t*exp(-t/3) fh_x=subs(f,t,x)*subs(h,t,t-x) % 形成被积函数 yf=int(fh_x,x,0,t) % 积分，结果输出 输出结果表达式为 yf = 9/4*exp(-t)+3/2*exp(-1/3*t)*t-9/4*exp(-1/3*t) 此结果仅当时成立。 1．4已知连续LTI系统，激励信号及系统的单位冲激响应分别如下所示，试绘出该系统零状态响应的波形。 2．连续系统的单位阶跃响应和单位冲激响应