信号与系统实验四 LTI系统的时域分析
目录
【实验目的】
【实验原理】
【实验设备】
【实验内容】
1. 至少用两种方法求以下系统的单位阶跃响应并绘图表示出来。
2. 编程实现教材p106, 2-15(3),并绘图表示出来。
【实验感悟】
【实验目的】
(1)学会运用MATLAB求解连续与离散系统的冲激响应和阶跃响应。
(2)学会运用MATLAB求解LTI系统的零输人响应和零状态响应。
【实验原理】
1.连续时间系统的冲激响应和阶跃响应求解
在连续时间LTI系统中,冲激响应和阶跃响应是系统特性的描述﹐对它们的分析是线性系统中极为重要的问题。输入为单位冲激函数àt)所引起的零状态响应称为单位冲激响应,简称冲激响应﹐用h(t)表示;输人为单位阶跃函数u(t)所引起的零状态响应称为单位阶跃响应,简称为阶跃响应﹐用g(t)表示。
在MATLAB中,对于连续LTI系统的冲激响应和阶跃响应的数值解,可分别用控制系统工具箱提供的函数impulse和 step来求解。其语句格式分别为:
y =impulse(syt ,t)
y =step(sys,t)
其中,t表示计算系统响应的时间抽样点向量,sys表示LTI系统模型。
2.连续时间系统的响应求解
LTI连续系统可用线性常系数微分方程来描述,即:
该微分方程的全解(又名全响应)分为两个组成部分:一是与该方程相应的齐次方程的齐次解,记作 yh(t),另一个是满足非齐次方程的特解,记作yp(t)。系统的全响应也可以分解为零输人响应和零状态响应两部分。零输人响应是指外加激励为零时,只由初始状态作用于系统所产生的响应﹐通常用yzi(t)表示;零状态响应是指初始状态为零时,由外加激励作用于系统所产生的响应,通常用yzs(t)表示。
MATLAB符号工具箱提供了dsolve函数,可实现常系数微分方程的符号求解,其调用格式为
其中,参数eql,eq2,…表示各微分方程,它与MATIAB符号表达式的输入基本相同,微分或导数的输人是用Dy,D2y,D3y,…来表示y的一阶导数、二阶导数、三阶导数等;参数condl,cond2,…表示各初始条件或起始条件;参数v表示自变量,默认为是变量t。
【实验设备】
(1)计算机。
(2)MATLAB软件。
【实验内容】
1. 至少用两种方法求以下系统的单位阶跃响应并绘图表示出来。
方法一:
t=0:0.001:6;%设置起点、步长和终点
sys=tf([1,3,7],[1,4,6,4,1]);%应用传递函数
g=step(sys,t);%用函数命令step求解阶跃响应的数值解
plot(t,g),grid on;%绘图命令
xlabel('Time(sec)'),ylabel('g(t)')%设置横纵轴标签
title('阶跃响应')%设置图像名称
方法二:
a=[1 4 6 4 1];
b=[1,3,7];
t=0:0.01:6;%设置起点、步长和终点
f=heaviside(t);%f是以t为变量的阶跃函数
lsim(b,a,f,t);%用lsim命令画出LTI系统对f和t所描述的信号对应出的时间响应
xlabel('Time'),ylabel('g(t)')%设置横纵轴标签
title('阶跃响应')%设置图像名称
2. 编程实现教材p106, 2-15(3),并绘图表示出来。
a=[1 2 1];
b=[1 1 1];
t=0:0.001:6;% 设置起点、步长和终点
f=cos(t).*heaviside(t);%设置表达式
lsim(b,a,f,t);%lsim命令画出LTI系统对f和t所描述的信号对应出的时间响应
3. 求其零状态解,并绘图表示出来
n=(-10:10);
a=[1,-0.5];
b=1;
x=sin(n*pi/4).*heaviside(n);%设置激励的表达式
y=filter(b,a,x);%用filter函数求解在x激励时所产生的响应序列的数值解
stem(n,y);%绘图命令
title('输出序列')%设置图像名称
【实验感悟】
通过本次实验,我学会了MATLAB中对于连续LTI系统的冲激响应和阶跃响应的数值解的求解方法,掌握了控制系统工具箱提供的函数lsim和step命令,尤其是对于一个系统的响应的不同求解方法,step(sys,t)命令是使用时间矢量t来模拟。在系统时间单位中,表达式t在sys的时间单位属性中是指定的。而lsim函数是针对线性时不变模型,给定任意输入,得到任意输出。lsim函数表示任意输入函数的响应,连续系统对任意输入函数的响应可以利用lsim函数求取。
通过此次实验,我也掌握了filter函数的使用,同时我也通过官方文献了解到表达式Y = filter(b,a,X) 滤除向量X中的数据,其中b是分子系数向量,a是分母系数向量。输入X为滤波前序列,Y为滤波结果序列,整个滤波过程是通过差分方程实现的。此次实验我也复习了前面已经学习过的命令,例如heaviside命令等等。
信号与系统实验四 LTI系统的时域分析相关推荐
- LTI的频域分析matlab,matlab与信号实验连续lti系统的频域分析
matlab与信号实验连续lti系统的频域分析 上机实验 3 连续 LTI 系统的频域分析一.实验目的 (1).掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法,以及傅立叶变换的时移特性,傅立叶变换 ...
- LTI的频域分析matlab,matlab与信号实验——连续lti系统的频域分析x26lc.doc
matlab与信号实验--连续lti系统的频域分析x26lc 上机实验3 连续LTI系统的频域分析一.实验目的 (1).掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法,以及傅立叶变换的时移特性,傅 ...
- php实现汉明码,实验四 汉明码系统
<实验四 汉明码系统>由会员分享,可在线阅读,更多相关<实验四 汉明码系统(11页珍藏版)>请在人人文库网上搜索. 1.实验四 汉明码系统 一. 实验原理和电路说明这些多差错控 ...
- 实验4 LTI系统的时域分析
实验4 LTI系统的时域分析 1.冲激响应:y=impulse(syt,t) 2.阶跃响应:y=step(sys,t) t表示系统响应的时间抽样点向量,sys表示LTI系统模型 表示系统的特性参数:s ...
- LTI的频域分析matlab,MATLAB与信号实验 —— 连续LTI系统的频域分析
上机实验3 连续LTI 系统的频域分析 一.实验目的 (1).掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法,以及傅立叶 变换的时移特性,傅立叶变换的频移特性的实现方法; (2).了解傅立叶变换的 ...
- 领悟《信号与系统》之 LTI 系统的卷积积分及性质
LTI 系统的卷积积分及性质 一.卷积积分及其性质 1. 任意信号的分解 2. 任意信号作用下的零状态响应 3. 卷积的定义 4. 卷积性质 1. 代数性质 2. f(t)与奇异信号的卷积 (非常重要 ...
- 领悟《信号与系统》之 LTI 系统的时域响应
连续时间信号的时域分析法 一.LTI 连续系统的响应 1. LTI 微分方程的建立与求解 2. LTI 连续系统的零输入响应与零状态响应 1. 零输入响应 2. 零状态响应 3. 全响应 二.冲激响应 ...
- 练习系统 实验四 red and black
当前编程题:实验四 基础算法问题(16级)---Red and Black 1. 问题描述 There is a rectangular room,covered with square tiles. ...
- 信号与系统实验 01 - | 连续系统频域分析的MATLAB实现
一.实验目的 1. 掌握连续时间信号频谱特性的MATLAB分析方法: 2. 掌握连续系统的频率响应MATLAB分析方法方法. 二.实验原理 1. 连续时间信号的频谱---傅里叶变换 syms t w ...
- C# 实验四 获取系统时间、点击加一秒功能
题目要求 运行效果 有60秒进位,60分钟进位,24小时进位(清零)功能 代码 Form1.cs using System; using System.Collections.Generic; usi ...
最新文章
- Python实现一元及多元线性回归
- 代码编译突然变缓慢问题解决办法(codeblock)
- libevent源码学习-----Reactor模型
- 文件隐藏服务器版本信息,如何隐藏Apache版本号和其他敏感信息
- 字体感觉小了 引入的vant_动态字体海报如何设计与制作?
- php 返回设置时间戳,PHP-返回int时间戳而不是datetime
- jquery城市选择案例
- 陀螺仪数据转换成角度_请教怎么把用陀螺仪积分得到的角度转换到大地坐标系下?...
- vi+ctags+cscope联合使用
- 一起学ORBSLAM2(9)ORBSLAM的PNP解决方案
- 博科BROCADE交换机开启snmp服务
- wps怎么做文档分享
- 2014年8月22日 百度发布“冰桶算法”
- 诛仙服务器技能修改,诛仙玩家浅析诛仙3技能修改所带来的轰动
- 2020IT从业者如何找到高薪工作
- 玩转群晖NAS——实现IPV6访问并部署动态域名解析程序(Golang)
- mysql数据库的介绍与初始配置
- 可以分屏的软件_mac必备软件
- Django框架之美多商城项目
- C++文件操作——fopen/fseek/ftell/fread/fwrite/fclose等函数用法