实验二、信号与系统时域分析的MATLAB实现.docxVIP

实验二、信号与系统时域分析的 MATLAB实现 一、 实验目的 掌握利用Matlab求解LTI系统的冲激响应、阶跃响应和零状态响应，理解 卷积概念。 二、 实验内容 1、卷积运算的MATLAB实现： 计算连续信号卷积用 MATLAB中的函数conv,可编写连续时间信号卷积通用 函数sconv, function [f,n]=sconv(f1,f2,n1,n2,p) f=conv(f1,f2);f=f*p; n3=n1(1)+n2(1); L=length(n1)+length(n2)-1; n=n3:p:(n3+(L-1)*p); 例 2.1 f (t) f1 (t) f2(t) 在MATLAB编写： p=0.01; n1=-1:p:1;f1=ones(1,length(n1)); n2=0:p:1;f2=2*n2; [f,n]=sconv(f1,f2,n1,n2,p); subplot(3,1,1),plot(n1,f1), axis([-1.5,1.5,0,2]),grid on subplot(3,1,2),plot(n2,f2), axis([-0.1,1.2,0,3]),grid on subplot(3,1,3),plot(n,f),axis([-1.5,5,0,2]),grid on 运算结果如图例2.1所示。 计算离散信号卷积用 MATLAB中的函数conv,其调用格式为： c conv(a,b) 式中，a、b为待卷积两序列的向量表示。向量 c的长度为a、b长度之和减1。 例2.2 已知序列x[k]1,2,3,4;k 0,1,2,3,y[k] 1,1,1,1;k 0,1,2,3 ,计 例2.2 已知序列x[k] x[k] y[ k]。 解：MATLAB的程序为： x=[1,2,3,4];y=[1,1,1,1]; z=conv(x,y) N=length(z); stem(0:N-1,z) 运行结果：如图例 2.2 z = 1 3 6 10 9 7 2、连续时间系统零状态响应的求解: LTI连续时间系统的零状态响应可通过求解初始状态为零的常系数微分方程得 到。在MATLAB中控制系统工具箱中提供了一个用于求解零初始条件微分方程数值解 的函数lsim。其调用方式为： y lsim (sys, f ,t) 式中，sys表示LTI系统模型，用来表示微分方程、差分方程、状态方程， t表示计算系 统响应的抽样点向量，f是系统输入信号向量。在求解微分方程时，微分方程的 LTI系 统模型sys借助于tf函数获得，其调用方式为： sys tf (b, a) 式中，b和a分别为微分方程右端和左端各项的系数向量。例如：对 2阶微分方程 3y(t) 2y(t) 5y(t) f(t) 2f(t) 可用 a=[3 2-5], b=[-1 0 2] , sys=tf(b,a)获得 LTI 模型。 例2.3 系统的微分方程为y(t) 2y(t) 100y(t) f(t),输入信号为 f (t) 10sin2 t ,求系统的零状态响应 y(t)。 解：MATLAB的程序为： sys=tf([0 0 1],[1 2 100]); t=0:0.01:5; f=10*sin(2*pi*t); y=lsim(sys,f,t); plot(t,y) 运行结果如图例2.3所示。 3、连续时间系统冲激响应和阶跃响应的 求解 MATLABimpulse MATLAB impulse和step分别表示。其调用方式为: 中控制系统工具箱中提供函数 y impulse(sy§ t) y step(sy§t) 例2.4 求例2.3所示系统的冲激响应。 解：MATLAB的程序为： sys=tf([0 0 1],[1 2 100]); t=0:0.01:5; y=impulse(sys,t); plot(t,y) title(连续系统冲激响应) xlabel(time(sec)) ylabel(y(t)) 运行结果如图例2.4所示。 4、 离散时间系统零状态响应的求解 LTI离散系统一般用线性常系数差分方程描述： aiy[k i] bj f[k j] i 0 j 0 f[k]、y[k]分别表示系统的输入和输出， n是差分方程的阶数。 已知差分方程的n个初始 状态和输入f[k],就可以通过编程由下式迭代计算出系统的输出。 在零初始状态下，MATLAB信号处理工具中提供了一个 filter函数，计算差分方程 描述的系统的响应。其调用方式为： y filter (b, a, f) 式中，b=[b0,b1,---,bn],a=[a0,a1,---an]分别是差分方程左、右端的系数向量， f表示输入 序列，y表示输出序列，它们长度要相同。 5、 离散时间系统单位脉冲响应的求解