用matlab解出冲激响应的表达式,实验三系统的时域分析

实验三系统的时域分析

重庆三峡学院 《信号与系统分析》实验 重庆三峡学院 Chongqing Three Gorges University 实 验 报 告 课程名称 信 号 与 系 统 分 析 实验名称 系 统 的 时 域 分 析 实验类型 验 证 学院名称 电子与信息工程学院 专 业 电子与信息工程学院(嵌入式) 年级班级 2014级电信1班 开出学期 2015~2016 下 学生姓名 王 英 吉 学 号 201407014119 实验教师 蒋 行 达 成 绩 2016年5月14日 实验三、系统的时域分析 1实验目的 1)学会利用MATLAB求解连续系统的零状态响应、冲激响应和阶跃响应； 2)学会利用MATLAB求解离散系统的单位取样响应； 3)学会利用MATLAB求解离散系统的卷积和。 2实验原理及实例分析 (实验原理见教材的第二章和第三章。) 2.1 连续系统零状态响应的数值求解 1：已知某LTI系统的微分方程为，其中，。试用MATLAB命令绘出范围内系统零状态响应的波形图。 解：程序如下： Clc %命令窗口清屏 close all clear all t = 0:0.01:5; sys = tf([6],[1 5 6]); %用传输函数形式表示系统 f = 10 * sin(2*pi*t) .* uCT(t); y = lsim(sys,f,t); %对输入信号模拟仿真 plot(t,y, Linewidth ,2);grid; xlabel( t(sec) );title( y(t) ); 产生的图形如图1所示。 图1 程序产生的图形 2.2连续系统的冲激响应和阶跃响应的数值求解 2：已知某LTI系统的微分方程为，试用MATLAB命令绘出范围内系统的冲激响应和阶跃响应。 解：MATLAB程序如下： clc close all clear all t = 0:0.01:5; sys = tf([1 16],[1 2 32]); h = impulse(sys,t); %计算系统的冲激响应 g = step(sys,t); %计算系统的阶跃响应 subplot(211); plot(t,h, Linewidth ,2);grid; xlabel( t(sec) );title( Impulse response --- h(t) ); subplot(212); plot(t,g, Linewidth ,2);grid; xlabel( t(sec) );title( Step response --- g(t) ); 图2 程序产生的图形 2.3离散系统的响应 3：已知系统的差分方程为，试用MATLAB命令绘出当激励信号为时，该系统的零状态响应。 解：MATLAB程序为： clc close all clear all a = [3 -4 2]; b = [1 2]; n = 0:30; x = (1/2) .^ n; y = filter(b,a,x); stem(n,y, fill );grid; xlabel( n );title( y(n) ); 程序产生的图形如图3所示。 图3程序产生的结果 2.4离散系统的单位取样响应 4：已知系统的差分方程为，利用MATLAB命令绘出该系统的单位取样响应。 解：MATLAB程序如下： clc close all clear all a = [3 -4 2]; b = [1 2]; impz(b,a,30);grid; 程序产生的图形如图4所示。 图4程序产生的图形 2.5离散卷积和的计算 5：已知系统的单位取样响应为，试用MATLAB计算当输入信号为时，系统的零状态响应。 解：MATLAB程序如下： clc close all clear all nx = -1:5; nh = -2:10; x = uCT(nx) - uCT(nx-4); h = 0.8 .^ nh .* (uCT(nh) - uCT(nh-8)); y = conv(x,h); ny = (min(nx)+min(nh)) : (max(nx) + max(nh)); subplot(311); stem(nx,x, fill );grid; xlabel( n );title( x(n) ); axis([-4 16 0 3]); subplot(312); stem(nh,h, fill );grid; xlabel( n );title( h(n) ); axis([-4 16 0 3]); subplot(313); stem(ny,y, fill );grid; xlabel( n );title( y(n)=x(n)*h(n) ); axis([-4 16 0 3]); 程序产生的图形如图5所示。 图5程序产生的图形 3实验报告与要求 在实验报告中，请简要说明连续系统的零状态响应、冲激响应和阶跃响应等概念，说明离散系统的单位取样响应、卷积和的计算方法。 答：零状态响应：系统的初始状态为零时，仅由输入信号引起的响应。 冲激响应：当初始状态为零时，输入为单位冲激函数所引起的响应。 阶跃响应：当初始状态为零时，输入为单位阶跃函数所引起的响应。 离散系统的单位取样响应定义为系统在δ(n)激励下系统的零状态响应，用h(n)表示。MATLAB求解单位取样响应可利用函数filter，并将激励设为前面所定义的impDTS函数。 离散系统的零状态响应是激励与系统的单位取样响应的卷积，因此卷积运算在离散时间信号处理领域被广泛应用。离散时间信号的卷积定义为 可见，离散时间信号的卷积运算是求和运算，因而常称为“卷积和”。 思考并完成下列各题： 1) 已知系统的微分方程为，如何利用MATLAB命令绘出该系统的冲激响应和阶跃响应的时域波形图。理论上计算出系统的单位冲激响应/单位函数响应、阶跃响应、零状态响应、全响应的表达式，并写出解题过程。 解：理论计算: 由阶跃响应定义得：g’’(t)+2g’(t)+2g(t)=δ(t)…………………① 令g’’(t)=aδ’(t)+bδ(t)+R0(t) 对上式