2_用MATLAB进行控制系统的时域分析
实验5-控制系统时域分析
实验5-控制系统时域分析实验目的：1. 掌握控制系统的时域分析方法；2. 熟悉控制系统的基本概念；3. 比较不同控制系统的性能指标，并对其优化。
实验原理：控制系统是由控制器、被控对象和传感器等组成的系统。
它的主要功能是将被控对象的输出值与预期输出值（设定值）进行比较，并根据比较结果对控制器的输出信号进行调整，以实现预期的控制系统动态响应。
系统的状态可以用输入输出关系来表示，通常用系统函数表示，它是输入信号与输出信号的转换函数。
根据系统函数的性质，系统的特性可以分析出来，比如稳态误差、响应时间和阻尼等。
控制系统的时域分析方法主要包括以下内容：1. 稳态误差分析稳态误差是指当控制系统到达稳定状态时，被控对象的输出值与设定值之间的差值。
它是一个反映控制系统偏离设定状态能力的指标。
稳态误差对于不同类型的系统有不同的计算方法，常见的系统类型包括比例控制系统、积分控制系统和派生控制系统。
比例控制系统的稳态误差是：$e_{ss}= \frac {k_p}{1+k_p}, (k_p \neq 0)$派生控制系统的稳态误差是0。
2. 基本响应特性分析一个控制系统的基本响应特性主要包括死区、超调量和稳定时间等。
死区是指当控制器输出的信号在一定范围内时，被控对象的输过不会发生变化。
死区对控制系统的响应时间和稳态误差有很大影响，通常需要根据系统的特点对死区进行调整。
超调量是指被控对象的输出值在达到设定值后，超出设定值的程度。
常见的超调量有百分比超调量和绝对超调量。
3. 阻尼及其影响阻尼是指系统的阻尼比，它是表征系统阻尼程度的一个参数。
阻尼对控制系统的稳定性和性能有很大影响。
当阻尼比为1时，系统的响应最快，但容易出现震荡现象。
阻尼比小于1时，系统的响应相对较慢，但是不会出现震荡现象。
当阻尼比大于1时，系统的响应速度较慢，但相对稳定。
实验步骤：本实验采用MATLAB软件对几种常见的控制系统进行时域分析，具体步骤如下：1. 打开MATLAB软件，新建文件进行编程。
利用MATLAB进行时域分析
实验一 利用MATLAB 进行时域分析一 实验目的1 掌握利用MA TLAB 绘制控制系统时域响应方法。
2 研究二阶系统的两个重要参数阻尼比ζ和无阻尼自然频率ωn 对系统动态性能的影响。
3 掌握利用绘制的图形进行线性系统时域分析的方法二 实验内容1 初步掌握MA TLAB 时域分析中的基本命令；2 绘制系统响应曲线图并根据图形进行性能分析；三 实验步骤1 初步掌握MA TLAB 时域分析中的基本命令；A 熟悉MA TLAB 运行环境；B 控制系统模型的描述如考虑函数432543232546()34276s s s s G s s s s s s ++++=+++++ 传递函数模型描述sys=tf(num,den)执行如下程序：num=[3 2 5 4 6]; den=[1 3 4 2 7 2];sys=tf(num,den)执行结果如下：Transfer function:3 s^4 + 2 s^3 +5 s^2 + 4 s + 6------------------------------------- s^5 + 3 s^4 + 4 s^3 + 2 s^2 + 7 s + 2C 绘制系统的响应曲线考虑如下传递函数225()425G s s s =++ 绘制系统的单位阶跃响应曲线。
执行如下程序：num=[0 0 25];den=[1 4 25];impulse （num, den ）;step(num,den);grid或者num=[0 0 25];den=[1 4 25];sys2=tf(num,den);impulse(sys2, t); (或者可缺省t，即impulse(sys2) ;)step(sys2,t);grid执行结果如下D 绘制多条响应曲线求典型二阶系统222()2nn nss sωζωωΦ=++，当ωn=10时，ζ分别为0,0.25,0.5，0.7,1,2时的单位阶跃响应。
基于MATLAB的时域分析法实现
性 定 常 系 统 中 ，我 们 根 据 系 统 的 闭环 传 递 函数 分 母 的最 高 次 幂 把 系 统 分 为 一 阶 、二 阶 、高 阶
系统 ，一 阶 、二 阶 系 统 可 以求 得 时 问响 应 的 解 析 解 ，亦 能 求 得 性 能 指 标 超 调 量 ％ 、 调 节 时
中 图 分 类 号 ：Ｔ １ ＭＡ Ｐ３２ 文 献标识 码 ：Ａ
ＭＡ Ｌ Ｂ是 由 Ｃ Ｅ Ｅ博 士 和 Ｊ Ｈ Ｉ Ｉ Ｅ工 程 师 开 发 的 一 套 高 性 能 的 数 值 计 算 和 数 据 ＴＡ ＬＶ Ｏ Ｎ ＬＴ Ｌ 可 视化 软 件 ． 具 有 语 言 简单 、紧 凑 ，语 法 限制 不 严 格 ，库 函数 、运 算 符 丰 富 ，程 序 设 计 自 它 由度 大 ，绘 图功 能 强 大 的 特点 ． 另外 ，它 针 对 不 同 学科 研 制 了 功 能 强 大 的 工 具 箱 ，这 些 工 具
具 有 极 高 的 编 程 效 率 ． 用 Ｍ Ｔ Ａ 编 程 实 现 线 性 定 常 系 统 分 析 ，使 整 个 分 析 过 程 变 得 简 单 、 易 操 运 Ａ ＬＢ 作 ． 绍 了用 ＭＡ Ｌ Ｂ实 现 时 域 分 析 的 编 程 方 法 ，并 给 出相 应 的 仿 真 波 形 ， 介 ＴＡ 关 键 词 ：ＭＡ Ｌ Ｂ；根 轨 迹 法 ； 时 域 分 析 ； 线 性 定 常 系 统 ＴＡ
２ 基于 Ｍ Ｔ Ａ Ａ Ｌ Ｂ的 时域 分 析 法 实现
收 稿 日期 ：２ ｏ ～ Ｏ ～ ｌ ｏ２ ｌ Ｏ 作 者 简 介 ：宋 维 君 （９７ ） 女 ，辽 宁 抚 顺 人 ，讲 师 ， 主 要 从 事 计 算 机 语 言 研 究 １６ 一 ，
自动控制原理MATLAB仿真实验报告
实验一 MATLAB 及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点1、 系统的典型响应有哪些2、 如何判断系统稳定性3、 系统的动态性能指标有哪些 三、实验方法（一） 四种典型响应1、 阶跃响应：阶跃响应常用格式：1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。
2、),(Tn sys step ；表示时间范围0---Tn 。
3、),(T sys step ；表示时间范围向量T 指定。
4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。
2、 脉冲响应：脉冲函数在数学上的精确定义：0,0)(1)(0〉==⎰∞t x f dx x f其拉氏变换为：)()()()(1)(s G s f s G s Y s f ===所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式： ① )(sys impulse ； ②);,();,(T sys impulse Tn sys impulse③ ),(T sys impulse Y =（二） 分析系统稳定性 有以下三种方法：1、 利用pzmap 绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容 (一) 稳定性1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出253722)(2342++++++=s s s s s s s G 的极点。
%Matlab 计算程序num=[3 2 5 4 6];den=[1 3 4 2 7 2];G=tf(num,den);pzmap(G);p=roots(den)运行结果： p =+ - + -P ole-Zero MapReal AxisI m a g i n a r y A x i s-2-1.5-1-0.500.5-1.5-1-0.50.511.5图1-1 零极点分布图由计算结果可知，该系统的2个极点具有正实部，故系统不稳定。
控制系统的时域分析
2、时域分析例一、对于典型的二阶系统的传递函数为：2n22n nG(s)2s+sωξωω=+，其中固有频率ωn=3，绘制当阻尼比ξ=0.1,0.2,0.4，0.707，0.9,1.0,2.0时系统的单位阶跃响应。
解：在MATLAB中建立M文件，程序如下：wn=3;sigma=[0.1 0.2 0.4 0.707 0.9 1.0 2.0];hold on;t=linspace(0,12,120)';for k=sigmanum=wn.^2;den=[1,2*wn*k,wn.^2];step(num,den,t)endtitle('典型二阶系统不同阻尼比时的单位阶跃响应曲线 ')gtext('sigma=0.1');gtext('sigma=0.2');gtext('sigma=0.4');gtext('sigma=0.707');gtext('sigma=0.9');gtext('sigma=1.0');gtext('sigma=2.0');运行结果为：分析：由图可见，在ξ=0.4-0.9的范围内，系统上升比较快，且超调量不是很大，故一般工程系统中，ξ就选在这个范围中，其中尤其当ξ=0.707时，响应较快i，此时的超调量为4.31%，通常称其为最佳阻尼，具有最佳阻尼的系统成为二阶最佳系统。
例二、已知连续系统的开环传递函数为：(s5)G(s)50(s3)(s4)s+=++试判断系统闭环的稳定性以及绘制闭环系统的零极点图并作出该系统的单位阶跃响应曲线求出单位阶跃响应稳态误差。
解：（1）判断稳定性方法一：利用roots()函数编制M文件，程序如下：k=50;z=[-5];p=[0 -3 -4];[n1,d1]=zp2tf(z,p,k);P=n1+d1;roots(P)运行结果为：ans =-1.0760 + 7.1000i-1.0760 - 7.1000i-4.8480 + 0.0000i计算结果表明，所有特征根的实部均为负，故闭环系统是稳定的。
第3章用MATLAB进行控制系统时域分析
第3章用MATLAB进行控制系统时域分析MATLAB是一种功能强大的计算机软件，被广泛用于进行控制系统分析和设计。
在控制系统领域，时域分析是一种常用的方法，它可以用来评估和改进系统的性能及稳定性。
在MATLAB中，有几个重要的函数可以用于进行控制系统时域分析，包括step、impulse、lsim和initial等。
这些函数可以帮助我们了解系统的响应和行为。
首先，我们可以使用step函数来绘制系统的阶跃响应。
阶跃响应表示系统对输入信号的反应。
通过分析阶跃响应曲线的特性，我们可以评估系统的稳定性和性能。
例如，我们可以从阶跃响应曲线中获取系统的超调量、上升时间和峰值时间等信息。
step函数的使用方法如下：```step(sys)```其中，sys是一个控制系统的传递函数或状态空间模型。
调用step 函数后，MATLAB会自动绘制系统的阶跃响应曲线，并且返回一个包含系统响应数据的结构体。
另一个常用的时域分析函数是impulse，它可以用来绘制系统对冲击输入的响应。
冲击响应是系统对于单位冲击信号的输出。
通过分析冲击响应曲线，我们可以了解系统的固有特性，例如共振频率和阻尼比。
impulse函数的使用方法与step函数类似：```impulse(sys)```同样，调用impulse函数后，MATLAB会绘制系统的冲击响应曲线，并且返回一个包含系统响应数据的结构体。
除了阶跃响应和冲击响应外，我们还可以使用lsim函数来绘制系统对任意输入信号的响应。
lsim函数可以接受用户自定义的输入信号，并给出系统的输出响应。
通过分析系统的输出信号，我们可以更加全面地了解系统的性能和行为。
lsim函数的使用方法如下：```lsim(sys,u,t)```其中，sys是一个控制系统的传递函数或状态空间模型，u是输入信号，t是时间向量。
调用lsim函数后，MATLAB会根据输入信号和时间向量绘制系统的输出响应曲线，并返回一个包含系统响应数据的结构体。
基于MATLAB自动控制系统时域频域分析与仿真
基于MATLAB自动控制系统时域频域分析与仿真MATLAB是一款强大的数学软件，也是自动控制系统设计的常用工具。
它不仅可以进行时域分析和频域分析，还可以进行相关仿真实验。
本文将详细介绍MATLAB如何进行自动控制系统的时域和频域分析，以及如何进行仿真实验。
一、时域分析时域分析是指对系统的输入信号和输出信号进行时域上的观察和分析，以了解系统的动态特性和稳定性。
MATLAB提供了一系列的时域分析工具，如时域响应分析、稳态分析和步骤响应分析等。
1.时域响应分析通过时域响应分析，可以观察系统对于不同的输入信号的响应情况。
在MATLAB中，可以使用`lsim`函数进行系统的时域仿真。
具体步骤如下：- 利用`tf`函数或`ss`函数创建系统模型。
-定义输入信号。
- 使用`lsim`函数进行时域仿真，并绘制系统输出信号。
例如，假设我们有一个二阶传递函数模型，并且输入信号为一个单位阶跃函数，可以通过以下代码进行时域仿真：```num = [1];den = [1, 1, 1];sys = tf(num, den);t=0:0.1:10;u = ones(size(t));[y, t, x] = lsim(sys, u, t);plot(t, y)```上述代码中，`num`和`den`分别表示系统的分子和分母多项式系数，`sys`表示系统模型，`t`表示时间序列，`u`表示输入信号，`y`表示输出信号。
通过绘制输出信号与时间的关系，可以观察到系统的响应情况。
2.稳态分析稳态分析用于研究系统在稳态下的性能指标，如稳态误差和稳态标准差。
在MATLAB中，可以使用`step`函数进行稳态分析。
具体步骤如下：- 利用`tf`函数或`ss`函数创建系统模型。
- 使用`step`函数进行稳态分析，并绘制系统的阶跃响应曲线。
例如，假设我们有一个一阶传递函数模型，可以通过以下代码进行稳态分析：```num = [1];den = [1, 1];sys = tf(num, den);step(sys)```通过绘制系统的阶跃响应曲线，我们可以观察到系统的稳态特性。
哈工大机电系统控制基础大作业Matlab时域分析
《机电系统控制基础》大作业一基于MATLAB的机电控制系统响应分析哈尔滨工业大学2013年12月12日1作业题目1. 用MATLAB 绘制系统2()25()()425C s s R s s s Φ==++的单位阶跃响应曲线、单位斜坡响应曲线。
2. 用MATLAB 求系统2()25()()425C s s R s s s Φ==++的单位阶跃响应性能指标：上升时间、峰值时间、调节时间和超调量。
3. 数控直线运动工作平台位置控制示意图如下：X i伺服电机原理图如下：LR(1)假定电动机转子轴上的转动惯量为J 1，减速器输出轴上的转动惯量为J 2，减速器减速比为i ，滚珠丝杠的螺距为P ，试计算折算到电机主轴上的总的转动惯量J ；(2)假定工作台质量m ，给定环节的传递函数为K a ，放大环节的传递函数为K b ，包括检测装置在内的反馈环节传递函数为K c ，电动机的反电势常数为K d ，电动机的电磁力矩常数为K m ，试建立该数控直线工作平台的数学模型，画出其控制系统框图；(3)忽略电感L 时，令参数K a =K c =K d =R=J=1，K m =10，P/i =4π，利用MATLAB 分析kb 的取值对于系统的性能的影响。
源代码：t=[0:0.01:5];u=t;C=[25],R=[1,4,25];G=tf(C,R);[y1,T]=step(G,t);y2=lsim(G,u,t);subplot(121),plot(T,y1);xlabel('t(sec)'),ylabel('x(t)'); grid on;subplot(122),plot(t,y2);grid on;xlabel('t(sec)'),ylabel('x(t)');仿真结果及分析：源代码：t=[0:0.001:1];yss=1;dta=0.02;C=[25],R=[1,4,25];G=tf(C,R);y=step(G,t);r=1;while y(r)<yss;r=r+1;endtr=(r-1)*0.001;[ymax,tp]=max(y);tp1=(tp-1)*0.001;mp=(ymax-yss)/yss;s=1001;while y(s)>1-dta && y(s)<1+dta;s=s-1;endts=(s-1)*0.001;[tr tp1 mp ts]仿真结果及分析：C = 25ans = 0.4330 0.6860 0.2538 1.0000由输出结果知：上升时间为0.4330秒，峰值时间为0.6860秒，最大超调量为0.2538，调整时间1.0000秒。
实验三 利用MATLAB进行时域分析
实验三
一、实验目的
(1)学会使用MATLAB编程绘制控制系统的单位阶跃响应曲线；
(2)研究二阶控制系统中，、n对系统动态特性和时域指标的影响；
(3)掌握准确读取动态特性指标的方法；
lab2='zunibi = 0.25';text(0.3,1.5,lab2),
>> lab3='zunibi = 0.5';text(0.3,1.2,lab3),
lab4='zunibi = 0.75';text(0,3,1.05,lab4),
lab5='zunibi = 1';text(0,35,0.9,lab5),
den=[1,4,25];
step(num,den)
grid%绘制网格线。
title(¹Unit-Step Response of G(s)=25/(s^2+4s+25) ¹) %图像标题
图3-2MATLAB绘制的响应曲线
还可以用下面的语句来得出阶跃响应曲线
>> G=tf([0,0,25],[1,4,25]);
den = [1 2*sigma*10 100];
damp(den)
sys = tf(num,den);
i=i+1;
step(sys,2)
hold on
end
Eigenvalue Damping Freq. (rad/s)
0.00e+000 + 1.00e+001i 0.00e+000 1.00e+001
基于MATLAB的控制系统时域分析
基于MATLAB的控制系统时域分析控制系统的时域分析是通过对控制系统的输入与输出进行观察、计算和分析，研究系统的动态特性和稳态特性，以便设计出性能良好的控制系统。
MATLAB是一种功能强大的数学软件，也是控制系统分析和设计常用的工具。
基于MATLAB的控制系统时域分析包括系统的零极点分析、步跃响应分析、频域响应分析以及稳态误差分析等方面。
首先，零极点分析是控制系统时域分析的基础。
零极点是系统传递函数的根，可以通过MATLAB的roots函数来求解。
根据零极点的位置，可以判断系统的稳定性、阻尼比和共振频率等信息。
例如，根在左半平面的系统为稳定系统，零根的个数与极根的个数之差即为系统的阶数；根在右半平面的系统为不稳定系统；根在虚轴上的系统为临界稳定系统等。
其次，步跃响应分析是研究系统对输入步跃信号的响应情况。
通过MATLAB中的step函数，可以绘制系统的单位阶跃响应曲线，并可以获得系统的一些重要参数，如上升时间、峰值时间、峰值大小和稳态误差等。
其中，上升时间是指系统从初始值到达其稳定值的时间；峰值时间是指系统响应曲线达到其峰值的时间；峰值大小是指系统响应曲线在峰值时的大小；稳态误差是指系统响应与输入信号的差异。
接下来，频域响应分析是研究系统对输入频率信号的响应情况。
通过MATLAB中的freqresp函数，可以绘制系统的频率响应曲线。
频域响应曲线可以展示系统在不同频率下的增益和相位特性，用于研究系统的稳定性和频率补偿。
频域响应分析常用的工具包括Bode图和Nyquist图等。
最后，稳态误差分析是分析系统输出与输入的差异，用于评估系统的性能。
通过MATLAB中的step函数和sys2tf函数，可以计算系统的稳态误差，并绘制稳态误差曲线。
稳态误差可以分为位置误差和速度误差，用于评估系统的静态精确性。
综上所述，基于MATLAB的控制系统时域分析涵盖了零极点分析、步跃响应分析、频域响应分析和稳态误差分析等方面。
实验二典型系统的时域响应分析实验仿真报告答案分析解析
实验二典型系统的时域响应分析实验仿真报告答案分析解析Company number：【WTUT-WT88Y-W8BBGB-BWYTT-19998】实验二典型系统的时域响应分析1. 实验目的1) 通过用MATLAB 及SIMULINK 对控制系统的时域分析有感性认识。
2) 明确对于一阶系统，单位阶跃信号、单位斜坡信号以及单位脉冲信号的响应曲线图。
3) 对于二阶系统阶跃信号的响应曲线图以及不同阻尼比、不同自然角频率取值范围的二阶系统曲线比较图。
4) 利用MATLAB 软件来绘制高阶控制系统的零极点分布图，判断此系统是否有主导极点，能否用低阶系统来近似，并将高阶系统与低阶系统的阶跃响应特性进行比较5）编制简单的M文件程序。
2. 实验仪器PC计算机一台，MATLAB软件1套3. 实验内容1）一阶系统的响应（1) 一阶系统的单位阶跃响应在SIMULINK 环境下搭建图1的模型，进行仿真，得出仿真曲线图。
理论分析：C（s）=1/[s+1)]由拉氏反变换得h(t)=1-e^(-t/ (t>=0)由此得知，图形是一条单调上升的指数曲线，与理论分析相符。
（2) 一阶系统的单位斜坡响应在SIMULINK 环境下搭建图2的模型，将示波器横轴终值修改为12进行仿真，得出仿真曲线图。
理论分析：C （s ）=1/[s^2(4s+1)]可求的一阶系统的单位斜坡响应为c(t)=(t-4)+4e^(-t/4)e(t)=r(t)-c(t)=4-4e^(-t/4) 当t=0时，e(t)=0,当趋于无穷时，误差趋于常值4.3） 一阶系统的单位脉冲响应在medit 环境下，编译一个.m 文件，利用impulse （）函数可以得出仿真曲线图。
此处注意分析在SIMULINK 环境中可否得到该曲线图。
理论分析：C （s ）=5/+2)=(5/2)/+1)可求的g(t)=^(-t/,是一个单调递减的函数。
两种环境下得到的曲线图不一致。
2）二阶系统的单位阶跃响应二阶系统的闭环传递函数标准形式为 其阶跃响应可以分以下情况解出①当0=ζ时，系统阶跃响应为 )cos(1)(t t c n ω-=②当10<<ζ时，系统阶跃响应为 )sin(111)(2θωζζω+--=-t e t c d tn其中ζζθ/121-=-tg ，21ζωω-=n d③当1=ζ时，系统阶跃响应为 t n n e t t c ωω-+-=)1(1)(④当1>ζ时，系统阶跃响应为 ⎪⎪⎭⎫ ⎝⎛---=21221121)(λλζωλλt t ne e t c 其中121---=ζζλ，122-+-=ζζλ （1）自然角频率1=n ω选取不同阻尼比=ζ0,,,,,,，用MATLAB 得到二阶系统阶跃响应曲线。
实验二 离散控制系统的性能分析1
实验二离散控制系统的性能分析（时域/频域）一、实验目的1．掌握离散闭环系统的动态性能时域参数的分析与计算方法；2．掌握离散系统稳定性的频域典型参数分析与计算方法。
二、实验工具1．MATLAB 软件（6.5 以上版本）；2．每人计算机一台。
三、实验内容1．在 Matlab 语言平台上，通过给定的闭环离散系统，深刻理解时域参数的物理意义与计算方法，内容包括如下：●阻尼比参数分析：Z 平面与 S 平面的极点相互转换编程实现；分析 S/Z 两个平面域特殊特性（水平线、垂直线、斜线、圆周等）的极点轨迹相互映射方法；系统阶跃响应参数：上升时间和超调量等。
2．采用频域分析方法，通过编程计算，进一步理解离散系统的稳定性参数，包括如下：●通过幅频图，进行增益裕度分析；●通过相频图，进行相位裕度分析。
四、实验步骤1.阻尼比计算注释：Example 1 Damping ratio computationts=0.1;gp=tf(1,[1 1 0])gz=c2d(gp,ts,'zoh')kz=tf(5*[1,-0.9],[1 -0.7],ts);sys_ta=feedback(gz*kz,1,-1)p=pole(sys_ta)- 2 -radii=abs(p);angl=angle(p)damp(sys_ta)real_s=log(radii)/tsimg_s=angl/tszeta=cos(atan(-img_s./real_s))wn=sqrt(real_s.^2+img_s.^2)运行结果：2.水平 S 平面线到 z 平面的映射注释：Example 2 Mapping of horizontal s-plane line to z-planexx=[0:0.05:1]'N=length(xx)s0=-xx*35;s=s0*[1 1 1 1 1]+j*ones(N,1)*[0,0.25,0.5,0.75,1]*pi/tsplot(real(s(:,1)),imag(s(:,1)),'-o',real(s(:,2)),imag(s(:,2)),'-s',... real(s(:,3)),imag(s(:,3)),'-^',real(s(:,4)),imag(s(:,4)),'-*',...real(s(:,5)),imag(s(:,5)),'-v'),sgridz=exp(s*ts)plot(real(z(:,1)),imag(z(:,1)),'-o',real(z(:,2)),imag(z(:,2)),'-s',... real(z(:,3)),imag(z(:,3)),'-^',real(z(:,4)),imag(z(:,4)),'-*',...real(z(:,5)),imag(z(:,5)),'-v'),zgrid3.垂直 S 平面线到 z 平面的映射注释：Example 3 Mapping of vertical s-plane line to z-planes0=j*xx*pi/ts;s=ones(N,1)*[0,-5,-10,-20,-30]+s0*[1 1 1 1 1]plot(real(s(:,1)),imag(s(:,1)),'-o',real(s(:,2)),imag(s(:,2)),'-s',...real(s(:,3)),imag(s(:,3)),'-^',real(s(:,4)),imag(s(:,4)),'-*',...real(s(:,5)),imag(s(:,5)),'-v'),sgridz=exp(s*ts)plot(real(z(:,1)),imag(z(:,1)),'-o',real(z(:,2)),imag(z(:,2)),'-s',...real(z(:,3)),imag(z(:,3)),'-^',real(z(:,4)),imag(z(:,4)),'-*',...real(z(:,5)),imag(z(:,5)),'-v'),zgrid4.恒定阻尼比 S 平面线映射到 z 平面注释：Example 4 Mapping of constant damping ratio s-plane lines into z-plane s=s0*[1 1 1 1]-imag(s0)*[0,1/tan(67.5*pi/180),...1/tan(45*pi/180),1/tan(22.5*pi/180)]s=[s,real(s(:,4))];plot(real(s(:,1)),imag(s(:,1)),'-o',real(s(:,2)),imag(s(:,2)),'-s',...real(s(:,3)),imag(s(:,3)),'-^',real(s(:,4)),imag(s(:,4)),'-*',...real(s(:,5)),imag(s(:,5)),'-v'),sgridz=exp(s*ts)plot(real(z(:,1)),imag(z(:,1)),'-o',real(z(:,2)),imag(z(:,2)),'-s',...real(z(:,3)),imag(z(:,3)),'-^',real(z(:,4)),imag(z(:,4)),'-*',...real(z(:,5)),imag(z(:,5)),'-v'),zgrid5.将圆 s 平面线映射到 z 平面注释：Example 5 Mapping of circle s-plane line to z-planephi=xx*pi/2s0=(pi/ts)*(-cos(phi)+j*sin(phi))s=s0*[1,0.75,0.5,0.25,0]plot(real(s(:,1)),imag(s(:,1)),'-o',real(s(:,2)),imag(s(:,2)),'-s',... real(s(:,3)),imag(s(:,3)),'-^',real(s(:,4)),imag(s(:,4)),'-*',...real(s(:,5)),imag(s(:,5)),'-v'),sgridz=exp(s*ts)plot(real(z(:,1)),imag(z(:,1)),'-o',real(z(:,2)),imag(z(:,2)),'-s',... real(z(:,3)),imag(z(:,3)),'-^',real(z(:,4)),imag(z(:,4)),'-*',...real(z(:,5)),imag(z(:,5)),'-v'),zgrid6.阶跃响应注释：Example 6 Step response measurek=[0:1:60];step(sys_ta,k*ts);7.根轨迹注释：Example 7 Root-locus analysisrlocus(gz*kz)Amplitude；注释：Example 8 Root-locus analysis in page 56 numg=[1 0.5];deng=conv([1 -0.5 0],[1 -1 0.5]);sys_z=tf(numg,deng,-1)rlocus(sys_z)注释：Example 9 Root-locus analysis in page 57numg=[1];deng=[1 4 0];ts=0.25sys_s2=tf(numg,deng)sys_z2=c2d(sys_s2,ts,'imp')rlocus(sys_z2)8.频率响应注释：Example 10 Analysis of frequency response and roots locus in page 59 a=1.583e-7;k=[1e7,6.32e6,1.65e6];w1=-1;w2=1;ts=0.1;v=logspace(w1,w2,100);deng=[1.638 1 0];numg1=k(1,1)*a*[-1 1]numg2=k(1,2)*a*[-1 1]numg3=k(1,3)*a*[-1 1]sys_s1=tf(numg1,deng)sys_s2=tf(numg2,deng)sys_s3=tf(numg3,deng)bode(sys_s1,sys_s2,sys_s3,v),grid onnumg=1.2e-7*[1 1]deng=conv([1 -1],[1 -0.242]);sys_z2=tf(numg,deng,ts)rlocus(sys_z2),grid on五、实验思考1. S 平面与 Z 平面不同位置的映射关系分析s平面虚轴的映射s平面整个虚轴映射为z平面单位圆，左半平面任一点映射在z平面单位圆内，右半平面任一点映射在单位圆外。
基于matlab的一阶系统及二阶系统时域特性仿真实验(1)
实验一 一阶系统及二阶系统时域特性MatLab 仿真实验（2学时）一、实验目的1、使学生通过实验中的系统设计及理论分析方法，帮助学生进一步理解自动控制系统的设计与分析方法。
2、熟悉仿真分析软件。
3、利用Matlab 对一、二阶系统进行时域分析。
4、掌握一阶系统的时域特性，理解常数T 对系统性能的影响。
5、掌握二阶系统的时域特性，理解二阶系统重要参数对系统性能的影响。
二、实验设备计算机和Matlab 仿真软件。
三、实验内容1、一阶系统时域特性 一阶系统11)(+=Ts s G ，影响系统特性的参数是其时间常数T ，T 越大，系统的惯性越大，系统响应越慢。
Matlab 编程仿真T=0.4，1.2，2.0，2.8，3.6，4.4系统单位阶跃响应。
2、二阶系统时域特性a 、二阶线性系统 16416)(2++=s s s G 单位脉冲响应、单位阶跃响应、单位正弦输入响应的 Matlab 仿真。
b 、下图为具有一微分负反馈的位置随动系统框图，求出系统的闭环传递函数，根据系统瞬态性能指标的定义利用Matlab 分别计算微分反馈时间常数τ为0，0.0125，0.025时系统的上升时间、峰值时间、最大超调量和调整时间。
C 、二阶线性系统3612362++s s ξ，当ξ为0.1，0.2，0.5，0.7，1.0，2.0时，完成单位阶跃响应的Matlab 仿真，分析ξ值对系统响应性能指标的影响。
四、实验要求1、进入机房，学生要严格遵守实验室规定。
2、学生独立完成上述实验，出现问题，教师引导学生独立分析和解决问题。
3、完成相关实验内容，记录程序，观察记录响应曲线，响应曲线及性能指标进行比较，进行实验分析4、分析系统的动态特性。
5、并撰写实验报告，按时提交实验报告。
五、Matlab 编程仿真并进行实验分析一、一阶系统时域特性:clearclcnum=1for del=0.4:0.8:4.4den=[del 1];step(tf(num,den))hold onendlegend('T=0.4','T=1.2','T=2.0','T=2.8','T=3.6','T=4.4')如图为T=0.4，1.2，2.0，2.8，3.6，4.4系统单位阶跃响应。
MATLAB论文-基于MATLAB的自动控制系统及案例分析
自动控制系统是在无人直接参与的情况下，可使生产过程按照期望的运行规律去运行的控制系统。系统是完成要求和任务的部分的组合[2]。
2.2
如果控制系统的输出量对系统运行过程或结果没有控制作用，这种系统称为开环控制系统。图2-1表示了开环控制系统输入量与输出量之间的关系。
1、MATLAB图形处理系统
这是MATLAB图形系统的基础，它包括生成二维数据和三维数据可视化、图像处理、动画及演示图形和创建完整的图形用户接口的命令。
2、MATLAB数学函数库
MATLAB提供了非常丰富的数学计算函数，可以进行简单和复杂的数学计算，比如矩阵求逆、级数求和、贝塞尔函数和快速傅里叶变换等。
MATLAB has powerful functions of numerical calculation and symbolic computation, as well as powerful data visualization and human-computer intelligent interactive ability, and its development is constantly bringing forth new ideas. Today's MATLAB has a richer data type and structure, a more friendly object-oriented, faster, more beautiful graphics visualization, a wider range of mathematical and data analysis resources, and more application development tools.
仿真实验一：控制系统的时域分析
仿真实验一：控制系统的时域分析一、实验目的：1．观察控制系统的时域响应；2．记录单位阶跃响应曲线；3．掌握时间响应分析的一般方法；4．初步了解控制系统的调节过程。
二、实验步骤：1．开机进入Matlab7.1运行界面。
2．Matlab指令窗："Command Window". 运行指令：con_sys; 进入本次实验主界面。
3．分别双击上图中的三个按键，依次完成实验内容。
4．本次实验的相关Matlab函数：tf([num],[den])可输入一传递函数。
step(G,t)在时间范围t秒内，画出阶跃响应图。
三、实验内容：1.观察一阶系统G=1/(T+s) 的时域响应：取不同的时间常数T，分别观察该系统的脉冲响应、阶跃响应、斜坡响应以及单位加速度响应。
实验操作：打开一阶系统的时域性能分析软件，调节时间常数大小，分别选择脉冲响应、阶跃响应、斜坡响应、加速度响应，则会得到相应响应的曲线。
T=2.5899sT=7.6259sT=12.0863s1、二阶系统的时域性能分析：（1）调节时间滑块，使阶跃响应最终出现稳定值。
操作过程：打开二阶系统性能分析软件，调节时间滑块，使阶跃响应最终出现稳定值，得到响应响应曲线。
固定阻尼比和响应时间，调节自然频率，观察零极点分布图和阶跃响应的变化。
固定自然频率和响应时间，调节阻尼比，观察零极点分布图和阶跃响应的变化。
根据要求，调节出符合上升时间、峰值时间、调解时间、超调量的阶跃响应，记录此时的自然频率、阻尼比、响应时间。
在一定的自然频率和阻尼比的情况下，通过调节相应时间，阶跃响应最终出现稳定值。
增大响应时间，可以使调节时间、峰值时间缩短。
1）阻尼比和响应时间一定的条件下，调节自然频率，对零极点分布的阶跃响应的影响：应的影响：（2）结合系统的零极点图，观察自然频率与阻尼比对极点位置的影响。
当阻尼比和响应时间一定的条件下，改变自然频率，对极点位置有一定的影响。
由图可以观察出，增大自然频率，角度不变，但极点的实部和虚部数值都同时增大。
(完整版)利用MATLAB进行时域分析
自动控制原理与系统课程实验报告实验题目：利用MATLAB进行时域分析班级：机电1131班姓名：刘润学号：38号一、实验目的及内容时域分析法是一种直接在时间域中对系统进行分析的方法，具有直观、准确的优点，并且可以提供系统时间响应的全部信息。
在此实验中，主要介绍时域法进行系统分析，包括一阶系统、二阶系统以及高阶系统,以及系统的性能指标。
通过实验，能够快速掌握、并利用MATLAB及控制系统箱对各种复杂控制系统进行时域分析。
二、实验设备三、实验原理典型的二阶系统在不同的阻尼比的情况下，它们的阶跃响应输出特性的差异是很大的。
若阻尼比过小，则系统的振荡加剧，超调量大幅度增加；若阻尼比过大，则系统的响应过慢，又大大增加了调整时间，下面通过此实验课题分析输出响应变化规律：已知二阶振荡环节的传递函数为：G（s）=ωn*ωn/（s*s+2*ζ*ωn*s+ωn*ωn），其中ωn=0.4，ζ从0变化到2，求此系统的单位阶跃响应曲线，并分析当ζ发生变化时，二阶系统的响应有什么样的变化规律。
四、实验步骤编出程序如下图：五、实验结果画出图表如下图：六、结果分析（1）当ξ=0（无阻尼）（零阻尼）时：无阻尼时的阶跃响应为等幅振荡曲线。
如图ξ=0曲线。
（2）当0<ξ<1（欠阻尼）时：对应不同的ξ，可画出一系列阻尼振荡曲线，且ξ越小，振荡的最大振幅愈大。
如图ξ=0.4曲线。
(3)当ξ=1（临界阻尼）时：临界阻尼时的阶跃响应为单调上升曲线。
如图ξ=1曲线。
（4）当ξ>1（过阻尼）时：过阻尼时的阶跃响应也为单调上升曲线。
不过其上升的斜率较临界阻尼更慢。
如图ξ=1.6曲线七、教师评语。
MATLABsimulink稳定性分析时域分析
Re( pi ) 0,i 1,2, , n
对于离散时间系统，如果系统全部极点都位于Z平 面的单位圆内，则系统是稳定的；否则系统是不稳定 的。
pi 1,i 1,2, , n 2
直接判定方法
对于传递函数模型tf(num,den)，利用求根 函数roots(den)来求极点。对于状态空间模型 SS(A,B,C,D)利用求特征值函数eig(A)来求特征 值。这样根据极点或特征值即可直接判定系统 的稳定性。
[u,t]=gensig(type,Ta,Tf,T)
其中：产生一个类型为type的信号序列u(t)， type为 以下标识字符串之一：’sin’—正弦波；’square’— 方波；’pulse’—脉冲序列，Ta为周期，Tf为持续时 间，T为采样时间。
运行结果显示： System is Stable
8
3. 利用李雅普诺夫第二法来判断系统的稳定性
线性定常连续系统
x Ax
在平衡状态xe=0处渐近稳定的充要条件是：对任给的 一个正定对称矩阵Q，存在一个正定的对称矩阵P，
且满足李雅普诺夫方程
ATP+PA＝-Q
而标量函数V(x)=xTPx是这个系统的一个二次型李雅
roots(d) 运行结果显示：
ans =
-12.8990
-5.0000
-3.1010
4
例7-1 已知闭环系统的传递函数为
G(s) 3s4 2s3 s2 4s 2 3s5 5s4 s3 2s2 2s 1
试判断系统的稳定性，并给出不稳定极点。
解：MATLAB程序如下
10
控制系统的时域分析实验报告
控制系统的时域分析实验报告实验目的：1.了解控制系统的时域分析方法；2.学习使用MATLAB进行时域分析；3.通过实验验证时域分析的准确性。
实验原理：时域分析是控制系统研究中的一种方法，通过研究系统在时间上的响应来研究系统的动态特性和稳定性。
在时域分析中，常用的方法包括脉冲响应、阶跃响应和正弦响应等。
通过对这些响应进行观察和分析，可以得到系统的各种性能指标，如超调量、响应时间、稳态误差等。
实验步骤：1.使用MATLAB编写程序，生成一个二阶控制系统的传递函数。
2.通过给控制系统输入一定的信号，观察系统的脉冲响应，并记录脉冲响应图像。
3.给控制系统输入一个阶跃信号，观察系统的阶跃响应，并记录阶跃响应图像。
4.给控制系统输入一个正弦信号，观察系统的正弦响应，并记录正弦响应图像。
5.根据实验数据，使用MATLAB分析系统的性能指标，如超调量、响应时间和稳态误差等。
实验结果：通过实验测得的数据和MATLAB分析，得到了控制系统的各种性能指标。
例如，测得的脉冲响应图像显示了系统的初值响应特性；阶跃响应图像显示了系统的过渡过程；正弦响应图像显示了系统的频率响应特性。
通过分析这些响应图像，可以得到系统的超调量、响应时间和稳态误差等指标。
实验结论：1.通过实验和分析，了解了控制系统的时域分析方法；2.掌握了使用MATLAB进行时域分析的技巧；3.实验证明了时域分析在控制系统研究中的重要性和准确性。
实验心得：通过进行控制系统的时域分析实验，我深刻认识到了时域分析在控制系统研究中的重要性。
通过观察和分析系统的脉冲响应、阶跃响应和正弦响应，可以全面了解系统的动态特性和稳定性。
同时，学会了使用MATLAB进行控制系统的时域分析，这将在我未来的研究工作中发挥重要作用。
实验结果验证了时域分析的准确性，这对我提高对控制系统的理解和研究能力有着积极影响。
实验二 控制系统的时域分析
自动控制原理MATLAB仿真实验化工机械系12自动化应凯业学号1220301015实验二 控制系统的时域分析一、实验目的学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性。
二、预习要点1、 系统的典型响应有哪些？2、 如何判断系统稳定性？3、 系统的动态性能指标有哪些？ 三、实验方法（一） 四种典型响应1、 阶跃响应：阶跃响应常用格式：1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。
2、),(Tn sys step ；表示时间范围0---Tn 。
3、),(T sys step ；表示时间范围向量T 指定。
4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。
2、 脉冲响应：脉冲函数在数学上的精确定义：0,0)(1)(0〉==⎰∞t x f dx x f其拉氏变换为：)()()()(1)(s G s f s G s Y s f ===所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式： ① )(sys impulse ； ②);,();,(T sys impulse Tn sys impulse③ ),(T sys impulse Y =（二） 分析系统稳定性 有以下三种方法：1、 利用pzmap 绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出253722)(2342++++++=s s s s s s s G 的极点。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验三用MATLAB进行控制系统的时域分析[实验目的]1．研究线性系统在典型输入信号作用下的暂态响应；2．熟悉线性系统的暂态性能指标；3．研究二阶系统重要参数阻尼比ξ对系统动态性能的影响；4．熟悉在MATLAB下判断系统稳定性的方法；5．熟悉在MATLAB下求取稳态误差的方法。
[实验指导]MATLAB中有两类用于求解系统时域响应的方法。
其一是利用MATLAB 中的控制系统工具箱(Control System Toolbox)提供的函数（命令）；其二是Simulink仿真，它主要用于对复杂系统进行建模和仿真。
一、用MATLAB函数（命令）进行暂态响应分析1 求取线性连续系统的单位阶跃响应的函数——step基本格式为：step(sys)step(num,den)step(A,B,C,D)step(sys,t)step(sys1,sys2,…,t)y=step(sys,t)[y,t]=step(sys)[y,t,x]=step(sys)其中模型对象的类型如下：sys = tf(num,den) 多项式模型sys = zpk(z,p,k) 零点极点模型sys = ss(a,b,c,d) 状态空间模型参数无t，表示时间向量t的范围自动设定。
参数有t，表示给定时间向量t，应该有初值，时间增量，末值，如t=0:0.01:2。
前5种函数可以绘出阶跃响应曲线；后3种函数不绘阶跃响应曲线，而是返回响应变量y，时间向量t，以及状态变量x。
2 求取线性连续系统的单位脉冲响应的函数——impulse基本格式为：impulse(sys)impulse(num,den)impulse (sys,tf)impulse (sys,t)impulse (sys1,sys2,…,t)y=impulse(sys,t)[y,t]=impulse(sys)[y,t,x]=impulse (sys)3 求取线性连续系统的单位斜坡响应MATLAB没有直接求系统斜坡响应的功能函数。
在求取控制系统的斜坡响应时，通常用阶跃响应函数step()求取传递函数为G (s)/s的系统的阶跃响应，则其结果就是原系统G (s)的斜坡响应。
原因是，单位阶跃信号的拉氏变换为1/s ，而单位斜坡信号的拉氏变换为1/s2。
4．求取线性连续系统对任意输入的响应的函数——lsim其格式为y=lsim(sys,u,t)其中，t为仿真时间，u为控制系统的任意输入信号。
5．暂态响应性能指标在阶跃响应曲线窗口，使用右键弹出浮动菜单，选择其中的Characteristics子菜单，有4个子项：①Peak Response 峰值响应，点击将出现标峰值记点，单击此标记点可获得峰值幅值，超调量和峰值时间。
②Settling Time 调节时间，点击将出现调节时间标记点，单击此标记点即可获得调节时间。
③Rise Time 上升时间，点击将出现上升时间标记点，单击此标记点即可获得上升时间。
④Steady State 稳定状态，若系统稳定，点击将在稳态值处出现标记点，单击此标记点即可获得稳态值；若系统不稳定，标记点不会出现。
对于不同的系统响应类型，Characteristics菜单的内容并不相同。
虽然不同响应曲线的特性参数不相同，但是均可以使用类似的方法从系统响应曲线中获得相应的信息。
Step ResponseTime (sec)A m p l i t u d e0123450.511.5Step ResponseTime (sec)A m p l i t u d e024680.511.56、其它①hold on 命令：可以允许在已经画曲线的图形窗口上再画新曲线；hold off 命令取消该功能。
②figure （i ）命令：打开第i 个图形窗口，把曲线绘在该图形窗口。
③grid on 命令：使图上出现网格。
④subplot(m,n,p)命令；把一个画面分成m×n 个图形区域, p 代表当前的区域号，可在每个区域中分别画一个图。
⑤也可以通过主界面菜单file/new/figure 打开1个新图形窗口，系统自动为其编号。
7．举例例1：系统传递函数为22033()210s G s s s +=++，求其阶跃响应、脉冲响应、斜坡响应、 ①分3个图形窗口分别显示的程序：>> n=[20,33];d=[1,2,10];sys=tf(n,d);step(sys);figure(2);impulse(sys); n1=[20,33];d1=[1,2,10,0];sys1=tf(n1,d1);figure(3);step(sys1) 运行结果：Step ResponseTim e (sec)A m p l i t u d e1234567Time (sec)A m p l i t u d e0123456Step ResponseTime (sec)A m p l i t u d e0123456789105101520253035③在1个窗口中的3个子图形窗口显示的程序：>> figure(1);n=[20,33];d=[1,2,10];sys=tf(n,d);subplot(2,2,1);impulse(sys); subplot(2,2,2);n1=[20,33];d1=[1,2,10,0];sys1=tf(n1,d1);step(sys1); subplot(2,2,3);step(sys) 运行结果：Time (sec)A m p l i t u d eStep ResponseTime (sec)A m p l i t u d e10203040Step ResponseTime (sec)A m p l i t u de2468例2：系统传递函数为210()210G s s s =++，求输入分别是自定义的1(t)和4*1(t)时的响应。
①2个输入和2个输出波形分别在4个子图形窗口显示的程序：>> figure(7);n=[10];d=[1,2,10];sys=tf(n,d);subplot(2,2,1);u=1+0*t;t=0:0.01:10;plot(t,u); subplot(2,2,3);y=lsim(sys,u,t);plot(t,y) ;>> u1=4+0*t;subplot(2,2,2);plot(t,u1);subplot(2,2,4);y1=lsim(sys,u1,t);plot(t,y1) 运行结果：05100.511.5205100.511.5051033.544.550510246②2个输入和2个输出波形在1个图形窗口同时显示的程序：>> n=[10];d=[1,2,10];sys=tf(n,d); t=0:0.01:10;u=1+0*t;sys2=tf([1],[1]) ; lsim(sys2,u,t);hold on;lsim(sys,u,t);u1=4+0*t;plot(t,u1);lsim(sys,u1,t)Linear Simulation ResultsTime (sec)A m p l i t u d e012345678910123456例3：系统传递函数为220()210G s s s =++，求输入分别是自定义的1(t)+3*sin(t)时的响应。
程序；>> n=[20];d=[1,2,10];sys=tf(n,d); t=0:0.001:15;u=1+3*sin(t);sys2=tf([1],[1]) ;figure(16); lsim(sys2,u,t);hold on;lsim(sys,u,t)Linear Simulation ResultsTime (sec)A m p l i t u d e51015-50510例4：系统传递函数为3220()2210G s s s s =+++，求系统的阶跃响应。
程序；>> n=[20];d=[1,2,6,10];sys=tf(n,d); figure(24);step(sys)Step ResponseTime (sec)A m p l i t u d e05101520253035400.511.522.53二、用Simulink 进行暂态响应分析1．系统仿真方框图的建立方框图的建立与实验一中所述相同，不同点是不用输入点与输出点标记，输入点安置信号发生器，比如阶跃输入信号；输出点安置示波器。
需要如下操作：打开Simulink →Sources 子库，将step 模块（阶跃输入信号）复制到（拽到）模型文件窗口，放到相应位置。
（或其他输入信号模块）打开Simulink →Sinks 子库，将scope 模块（示波器）复制到（拽到）模型文件窗口，放到相应位置。
输入信号模块和示波器模块都可以进行参数设置。
2．设置仿真控制参数打开Simulation 菜单，找到Parameters 选项，可打开参数设置对话框。
它包括仿真时间范围的选择、仿真算法的选择、仿真步长的指定及仿真精度（误差）的定义等。
3．运行可选择 Simulatio n → Start 。
点击示波器，在示波器窗口中可以看到响应仿真曲线。
举例；系统方框图和阶跃输入下示波器显示的响应曲线如下：三、在MATLAB 下判断系统稳定性首先求得闭环传递函数，再使用MATLAB 函数（命令）roots(den)解出特征方程的根，即闭环极点，再根据极点位置，判断系统是否稳定。
举例：判断系统43250()28526G s s s s s =++++ 是否稳定。
程序和结果：>> d=[1 2 8 5 26]; roots(d) ans =-1.4509 + 2.1633i -1.4509 - 2.1633i 0.4509 + 1.9049i 0.4509 - 1.9049i分析判断：有2个根在s 平面右半部分，系统不稳定。
或者：>> n=[50]; d=[1 2 8 5 26];sys=tf(n,d)Transfer function: 50---------------------------------------- s^4 + 2 s^3 + 8 s^2 + 5 s + 26>> roots(sys.den{1})ans =-1.4509 + 2.1633i -1.4509 - 2.1633i 0.4509 + 1.9049i 0.4509 - 1.9049i四、在MATLAB 下求取稳态误差求取稳态误差终值的函数（命令）为dcgain( ) 调用格式为dcg=dcgain(G ) 其中 G=s ·R(S)·φe(S)R(S): 输入信号的拉氏变换； φe(S)：误差传递函数； 举例：系统前向通道传递函数为3250()2858G s s s s =+++，反馈通道传递函数为()4/25H s =，求输入为r(t)=1(t)时的稳态误差。