实验三基于SIMULINK的控制系统时域分析

实验三基于SIMULINK的控制系统时域分析
实验三基于SIMULINK的控制系统时域分析

实验三 基于SIMULINK 的控制系统时域分析(2学时)

一.实验目的:

掌握使用SIMULINK 、控制工具箱求解系统的输入和输出响应的仿真方法。

二.实验方法及预习内容:

利用SIMULINK 工具进行控制系统模型分析、系统设计与仿真的相关原理。

三.实验内容:

1.分别使用解微分方程方法、控制工具箱、Simulink 求解具有如下闭环传递函数的系统的阶跃响应。

432

10

()8364010

s s s s s φ=

++++

解微分方程方法求解: Wffc.m 文件:

function dx=wffc(t,x) u=1;

dx=[-8*x(1)-36*x(2)-40*x(3)-10*x(4)+u;x(1);x(2);x(3)];

主文件:

%-------------------------------------实验1.1 [t,x]=ode45('wffc',[0,18],[0;0;0;0]) y=10*x(:,4); figure(2) plot(t,y) grid on

title('解微分方程方法求解系统阶跃响应曲线') xlabel('时间') ylabel('输出')

结果:

控制工具箱求解:

程序:

num=[10];

den=[1 8 36 40 10];

sys=tf(num,den);

step(sys)

title('控制工具箱求解系统阶跃响应曲线') xlabel('时间')

ylabel('输出')

grid

结果:

Simulink求解设置:

单位阶跃

仿真设置

SCOPE设置

程序

figure(3)

plot(ScopeData(:,1),ScopeData(:,2)) grid on

title('Simulink求解系统阶跃响应曲线') xlabel('时间')

ylabel('输出')

结果

2. 某小功率随动系统动态结构如图所示,已知:10.01T =, 20.05T =,

01K =,1300K =,21K =,0.08c K =。若系统输入分别为()1()sr t t θ=,sr t θ=,

[1()1(1.5)]sr t θ=-,试用Simulink 分析系统的输出()sc t θ分别如何?

设置:

单位阶跃:

单位斜坡

两阶跃叠加

仿真设置

Scope设置:

程序:

figure(4)

plot(gg(:,1),gg(:,2),'*') grid on

title('不同信号下的控制系统输出响应曲线')

xlabel('时间')

ylabel('输出')

hold on

plot(gg(:,1),gg(:,3),'LineWidth',1)

hold on

plot(gg(:,1),gg(:,4),'g','LineWidth',2)

legend('单位阶跃下系统输出响应曲线','单位斜坡下系统输出响应曲线','两阶跃叠加下系统输出响应曲线')

结果:

四.实验要求:

1.熟悉Simulink法进行控制系统时域分析的基本步骤,并与微分方程法、控制工具箱法进行比较;

2.熟悉Simulink法在输入不同信号(单位阶跃、单位斜坡、两阶跃叠加)下的控制系统输出响应实验方法、图形显示。

二阶系统时域分析

1.有一位置随动系统,其结构图如下图所示,其中K = 4。求该系统的:1)自然 k 振荡角频率;2)系统的阻尼比;3)超调量和调节时间;4)如果要求 <0.707 , 值。 应怎样改变系统参数 K k 2.已知受控对象的开环传递函数为

(1)单位反馈时,计算单位脉冲响应的输出。 (2)试采用速度反馈方法,使得系统的阻尼比ζ=05.,确定速度反馈系数τ的值,并计算性能改善后的动态性能。 解 (1)单位反馈时,闭环传递函数为 其单位脉冲响应为 响应曲线为等幅振荡的,所以该系统仅作单位反馈,不能实现调节作用。 (2)增加速度反馈如图所示。 闭环传递函数为 ζωτ=,所以 阻尼比ζ=05.,则有2 n τ=?= 20.50.95 此时,系统阶跃响应的超调量为 调节时间为 3.已知速度反馈控制系统如图所示,要求系统的超调量为20%,峰值时间为1秒,试计算相应的前向增益K与速度反馈系数K 的值。如果保持K值不变,Kf为零时,计算超调量增大值。

解上述系统的闭环传递函数为 比较二阶系统的标准式有 给定的性能指标为 上述指标与系统特征参数ζ和ωn的关系为: 解得 所以: 当K=125.,Kf=0时,也就是没有速度反馈时,闭环传递函数成为: 阻尼比:

超调量增大为: 4.对下图所示系统,试求K为何值时,阻尼比ζ=0.5。并求此时系统单位阶跃响应的最大超调量和调整时间。 解:系统开环传函为: 系统闭环传函为: 最大超调量: 调整时间

5. 系统结构如图,欲使超调量бp =0. 2, 过渡过程时间t s =1秒(Δ=0.02), 试确定K 和τ的值。 答案: ()2222(2)2n n n K s s K s K s ωτζωωΦ==+++++ 0.456ζ= 8.77 n ω= 277n K ω== 0.078τ= 6. 题图所示机械系统,当受到 F =40N 力的作用时,位移量xt ()的阶跃响应如图所示,试确定机械系统的参数m ,k, f 的值。 解: 图示机械系统的传递函数为 由图所示稳态值()1c ∞=,由终值定理 得到 K=40N/m 由超调量: 峰值时间:

大作业1(机电控制系统时域频域分析)

《机电系统控制基础》大作业一 基于MATLAB的机电控制系统响应分析 哈尔滨工业大学 2013年11月4日

1 作业题目 1. 用MATLAB 绘制系统2 ()25()() 425 C s s R s s s Φ== ++的单位阶跃响应曲线、单位斜坡响应曲线。 2. 用MATLAB 求系统2 ()25 ()()425 C s s R s s s Φ==++的单位阶跃响应性能指标:上升时间、峰值时间、调节时间和超调量。 3. 数控直线运动工作平台位置控制示意图如下: X i 伺服电机原理图如下: L R (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 的取值对于系统的性能的影响。

2 题目1 单位脉冲响应曲线 单位阶跃响应曲线

源代码 t=[0:0.01:1.6]; %仿真时间区段和输入 nC=[25]; dR=[1,4,25]; fi=tf(nC,dR); %求系统模型 [y1,T]=impulse(fi,t); [y2,T]=step(fi,t); %系统响应 plot(T,y1); xlabel('t(sec)'),ylabel('x(t)'); grid on; plot(T,y2); xlabel('t(sec)'),ylabel('x(t)'); grid on; %生成图形 3 题目2 借助Matlab,可得: ans = 0.4330 0.6860 25.3826 1.0000 即

第三章控制系统的时域分析法知识点

第三章 控制系统的时域分析法 一、知识点总结 1.掌握典型输入信号(单位脉冲、单位阶跃、单位速度、单位加速度、正弦信号)的拉氏变换表达式。 2.掌握系统动态响应的概念,能够从系统的响应中分离出稳态响应分量和瞬态响应分量;掌握系统动态响应的性能评价指标的概念及计算方法(对于典型二阶系统可以直接应用公式求解,非典型二阶系统则应按定义求解)。 解释:若将系统的响应表达成拉普拉氏变换结果(即S 域表达式),将响应表达式进行部分分式展开,与系统输入信号极点相同的分式对应稳态响应;与传递函数极点相同的分式对应系统的瞬态响应。将稳态响应和瞬态响应分式分别进行拉氏逆变换即获得各自的时域表达式。 性能指标:延迟时间、上升时间、峰值时间、调节时间、超调量 3.掌握一阶系统的传递函数形式,在典型输入信号下的时域响应及其响应特征;掌握典型二阶系统的传递函数形式,掌握欠阻尼系统的阶跃响应时域表达及其性能指标的计算公式和计算方法;了解高阶系统的性能分析方法,熟悉主导极点的概念,定性了解高阶系统非主导极点和零点对系统性能的影响。 tr tp ts td

4.熟悉两种改善二阶系统性能的方法和结构形式(比例微分和测速反馈),了解两种方法改善系统性能的特点。 5.掌握系统稳定性分析方法:劳斯判据的判断系统稳定性的判据及劳斯判据表特殊情况的构建方法(首列元素出现0,首列出现无穷大,某一行全为0);掌握应用劳斯判据解决系统稳定裕度问题的方法。了解赫尔维茨稳定性判据。 6.掌握稳态误差的概念和计算方法;掌握根据系统型别和静态误差系数计算典型输入下的稳态误差的方法(可直接应用公式);了解消除稳态误差和干扰误差的方法;了解动态误差系数法。 二、相关知识点例题 例1. 已知某系统的方块图如下图1所示,若要求系统的性能指标为: δδ%=2222%,tt pp=1111,试确定K和τ的值,并计算系统单位阶跃输入下的特征响应量:tt,tt。 图1 解:系统闭环传递函数为:Φ(s)=CC(ss)RR(ss)=KK ss2+(1+KKKK)ss+KK 因此,ωnn=√KK,ζζ=1+KKKK2√KK, δ%=e?ππππ?1?ππ2?ζζ=0.46, t pp=ππωωdd=1ss?ωdd=ωnn?1?ζζ2=3.14 ?ωnn=3.54 K=ωnn2=12.53,τ=2ζζωnn?1KK=0.18 t ss=3ζζωωnn=1.84ss

控制系统时域与频域性能指标的联系

控制系统时域与频域性能指标的联系 经典控制理论中,系统分析与校正方法一般有时域法、复域法、频域法。时域响应法是一种直接法,它以传递函数为系统的数学模型,以拉氏变换为数学工具,直接可以求出变量的解析解。这种方法虽然直观,分析时域性能十分有用,但是方法的应用需要两个前提,一是必须已知控制系统的闭环传递函数,另外系统的阶次不能很高。 如果系统的开环传递函数未知,或者系统的阶次较高,就需采用频域分析法。频域分析法不仅是一种通过开环传递函数研究系统闭环传递函数性能的分析方法,而且当系统的数学模型未知时,还可以通过实验的方法建立。此外,大量丰富的图形方法使得频域分析法分析高阶系统时,分析的复杂性并不随阶次的增加而显著增加。 在进行控制系统分析时,可以根据实际情况,针对不同数学模型选用最简洁、最合适的方法,从而使用相应的分析方法,达到预期的实验目的。 系统的时域性能指标与频域性能指标有着很大的关系,研究其内在联系在工程中有着很大的意义。 一、系统的时域性能指标 延迟时间t d 阶跃响应第一次达到终值h (∞)的50%所需的时间 上升时间 t r 阶跃响应从终值的10%上升到终值的90%所需的时间;对有振荡的系 统,也可定义为从0到第一次达到终值所需的时间 峰值时间t p 阶跃响应越过终值h (∞)达到第一个峰值所需的时间 调节时间 t s 阶跃响应到达并保持在终值h (∞)的±5%误差带内所需的最短时间 超调量%σ 峰值h( t p )超出终值h (∞)的百分比,即 %σ= () ()() ∞∞-h h h t p ?100% 二、系统频率特性的性能指标 采用频域方法进行线性控制系统设计时,时域内采用的诸如超调量,调整时间等描述系统性能的指标不能直接使用,需要在频域内定义频域性能指标。

一阶系统时域分析

1.已知一单位负反馈系统的单位阶跃响应曲线如下图所示,求系统的闭环传递函数。 解答: ①max ()100100()X X %%e %X δ-∞=?=?∞ 由 2.1820.090.6082e ξ-==?= ②0.8 4.946m n t ω==?= ③2222224.4648.9222 6.01424.46 6.01424.46 n B n n W K s s s s s s ωωω=?=?=++++++ 2.已知系统如下图所示,求系统的单位阶跃响应,并判断系统的稳定性。 解答: ()() ()210 1101061010.511B s s W s s s s s +==+++++ 3.16n ω==, 260.95n ξωξ=?

( )()1sin n t c X t ξωωθ-= ,arctg θ= ()31 3.2sin 0.98718.19t e t -=-+? (5分) 系统根为 1,2632P j -= =-±,在左半平面,所以系统稳定。 3.一阶系统的结构如下图所示。试求该系统单位阶跃响应的调节时间t s ;如果要求t s (5%)≤ 0.1(秒),试问系统的反馈系数应取何值? (1)首先由系统结构图写出闭环传递函数 得 T =0.1(s ) 因此得调节时间 t s =3T =0.3(s),(取5%误差带) (2)求满足t s (5%) ≤0.1(s )的反馈系数值。 假设反馈系数K t (K t >0) ,那么同样可由结构图写出闭环传递函数 由闭环传递函数可得 T = 0.01/K t 100()10()100()0.1110.1c B r X s s W s X s s s ===++?1001/()1000.0111t B t t K s W s K s s K ==+?+

控制系统的时域分析

实验报告 实验名称:实验1:控制系统的时域分析 课程名称:自控控制原理 专业:电气工程及其自动化 班级:130037 学生姓名:施苏伟 班级学号:13003723 指导教师:杨杨 实验日期:2015 年10 月16日

一、实验目的 1.观察控制系统的时域响应; 2.记录单位阶跃响应曲线; 3.掌握时间响应分析的一般方法; 4.初步了解控制系统的调节过程。 二.实验步骤: 1.将‘实验一代码’这个文件夹拷贝到桌面上; 2.开机进入Matlab6.1 运行界面(其他版本亦可); 3.通过下面方法将当前路径设置为‘实验一代码’这个文件夹所在的路径 4.Matlab 指令窗>>后面输入指令:con_sys; 进入本次实验主界面。 5.分别双击上图中的三个按键,依次完成实验内容。

6.本次实验的相关Matlab 函数: 传递函数G=tf([num],[den])可输入一传递函数,其中num、den 分别表示分子、分母按降幂排列的系数。 三、仿真结果: (一)观察一阶系统G=1/(T+s)的时域响应: T=5s T=8s

T=13s 结果分析:一阶系统 G=1/(T+s)的,通过观察曲线发现,随着时间常数T的增大,同种响应要达到相同响应的时间增大,说明T越大,响应越慢。 (二)二阶系统的时域性能分析 (1)

结果分析:自然频率和阻尼比的适当时,通过调节相应的时间,阶跃响应可以得到稳定值。 (2)数据一:自然频率=5.96rad/sec 阻尼比=0.701

数据二:自然频率=8.2964rad/sec 阻尼比=0.701 结果分析:要达到既定范围,自然频率增大阻尼比要随之增大 (3)

第3章线性系统的时域分析习题答案

第3章 线性系统的时域分析 学习要点 1控制系统时域响应的基本概念,典型输入信号及意义; 2控制系统稳定性的概念、代数稳定判据及应用; 3控制系统的时域指标,一阶二阶系统的阶跃响应特性与时域指标计算; 4高阶系统时域分析中主导极点和主导极点法; 5 控制系统稳态误差概念、计算方法与误差系数,减小稳态误差的方法。 思考与习题祥解 题 思考与总结下述问题。 (1)画出二阶系统特征根在复平面上分布的几种情况,归纳ξ值对二阶系统特征根的影响规律。 【 (2)总结ξ和n ω对二阶系统阶跃响应特性的影响规律。 (3)总结增加一个零点对二阶系统阶跃响应特性的影响规律。 (4)分析增加一个极点可能对二阶系统阶跃响应特性有何影响 (5)系统误差与哪些因素有关试归纳减小或消除系统稳态误差的措施与方法。 (6)为减小或消除系统扰动误差,可采取在系统开环传递函数中增加积分环节的措施。请问,该积分环节应在系统结构图中如何配置,抗扰效果是否与扰动点相关 答:(1)二阶系统特征根在复平面上分布情况如图所示。 图 二阶系统特征根在复平面上的分布 当0ξ=,二阶系统特征根是一对共轭纯虚根,如图中情况①。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,变化轨迹是 以n ω为半径的圆弧,如图中情况②。 @ 当1ξ=,二阶系统特征根是一对相同的负实根,如图中情况③。 当1ξ>,二阶系统特征根是一对不等的负实根,如图中情况④。

(2)ξ和n ω是二阶系统的两个特征参量。 ξ是系统阻尼比,描述了系统的平稳性。 当0ξ=,二阶系统特征根是一对共轭纯虚根,二阶系统阶跃响应为等幅振荡特性,系统临界稳定。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,二阶系统阶跃响应为衰减振荡特性,系统稳定。ξ越小,二阶系统振荡性越强,平稳性越差; ξ越大,二阶系统振荡性越弱,平稳性越好。因此,二阶系统的时域性能指标超 调量由ξ值唯一确定,即001_ 100%2 ?=-π ξξ σe 。在工程设计中,对于恒值控制系 统,一般取 ξ=~;对于随动控制系统ξ=~。 n ω是系统无阻尼自然振荡频率,反映系统的快速性。当ξ一定,二阶系统的 时域性能指标调节时间与n ω值成反比,即34 s n t ξω≈。 (3)二阶系统增加一个零点后,增加了系统的振荡性,将使系统阶跃响应的超调量增大,上升时间和峰值时间减小。 所增加的零点越靠近虚轴,则上述影响就越大;反之,若零点距离虚轴越远,则其影响越小。 (4)二阶系统增加一个极点后,减弱了系统的振荡性,将使系统阶跃响应的超调量减小,上升时间和峰值时间减小; 所增加的极点越靠近虚轴,则上述影响就越大;反之,若极点距离虚轴越远,则其影响越小。 & (5)系统误差与系统的误差度(开环传递函数所含纯积分环节的个数或系统型别)、开环放大系数,以及作用于系统的外部输入信号有关。如果是扰动误差还与扰动作用点有关。 因此,减小或消除系统稳态误差的措施与方法有:增大开环放大系数,增加系统开环传递函数中的积分环节,引入按给定或按扰动补偿的复合控制结构。 无论采用何种措施与方法减小或消除系统稳态误差,都要注意系统须满足稳定的条件。 (6)采取在系统开环传递函数中增加积分环节的措施来减小或消除系统扰动误差时,所增加的积分环节须加在扰动作用点之前。若所增加的积分环节加在扰动作用点之后,则该积分环节无改善抗扰效果作用。这一点可以通过误差表达式分析得到。 题系统特征方程如下,试判断其稳定性。 (a )0203.002.023=+++s s s ; (b )014844122345=+++++s s s s s ; (c )025266.225.11.0234=++++s s s s ! 解:(a )稳定; (b )稳定; (c )不稳定。

控制系统的时域分析实验报告

课程名称:控制理论指导老师:成绩: 实验名称:控制系统的时域分析实验类型:冋组学生姓名: 、实验目的和要求 1用计算机辅助分析的办法,掌握系统的时域分析方法。 2. 熟悉SimUlink仿真环境。 二、实验内容和原理 (一)实验原理 系统仿真实质上就是对系统模型的求解,对控制系统来说,一般模型可转化成某个微分方程或差分方程表示,因此在仿真过程中,一般以某种数值算法从初态出发,逐步计算系统的响应,最后绘制出系统的响应曲线,进而可分析系统的性能。控制系统最常用的时域分析方法是,当输入信号为单位阶跃和单位冲激函数时,求出系统的输出响应,分别称为单位阶跃响应和单位冲激响应。在MATLAB中,提供了求取连 续系统的单位阶跃响应函数step,单位冲激响应函数impulse,零输入响应函数initial等等。 (二)实验内容 二阶系统,其状态方程模型为 U X I y = [1.9691 6.4493] +[0] U X2 1?画出系统的单位阶跃响应曲线; 2. 画出系统的冲激响应曲线; 3. 当系统的初始状态为x0=[1,0]时,画出系统的零输入响应; 4. 当系统的初始状态为零时,画出系统斜坡输入响应; (三)实验要求 1. 编制MATLAB程序,画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应; 2. 在SimUIink仿真环境中,组成系统的仿真框图,观察单位阶跃响应曲线并记录之。 三、主要仪器设备 计算机一台以及matlab软件,SimUIink仿真环境 四、操作方法与实验步骤 1、程序解决方案: 在MATLAB 中建立文件shiyu.m ,其程序如下: %时域响应函数 fun ction G1 = shiyu( A,B,C,D)

高阶系统的时域分析

题 目: 高阶系统的时域分析 初始条件:设单位系统的开环传递函数为 ) )(105() ()(2 a s s s s b s K s G ++++= 要求完成的主要任务: (包括课程设计工作量及其技术要求,以及说明书撰写等具体要 求) (1) 当K=10,a=1,b=5时用劳斯判据判断系统的稳定性。 (2) 如稳定,则求取系统的单位阶跃响应、单位斜坡响应和单位加速度响应,用 Matlab 绘制相应的曲线,并计算单位阶跃响应的动态性能指标和稳态性能指标,计算单位斜坡响应和单位加速度响应的稳态性能指标。 (3) 如不稳定,则计算系统稳定时K 、a 和b 的取值范围,在稳定范围内任取一值 重复第2个要求。 (4) 绘制稳定时系统的根轨迹(在稳定范围内任取a 、b 值)。分析K 变化对系统 性能的影响。 时间安排:

指导教师签名:年月日系主任(或责任教师)签名:年月日

目录 摘要........................................................... I 1系统稳定性分析.. (1) 2不同输入信号的时域响应曲线 (2) 2.1系统单位阶跃响应曲线 (2) 2.2系统单位斜坡函数响应曲线 (3) 2.3系统单位加速度响应曲线 (4) 3动态性能指标与稳态性能指标 (6) 3.1动态性能指标计算 (6) 3.1.1采用主导极点分析 (6) 3.1.2应用MATLAB软件进行分析 (6) 3.2稳态性能指标 (8) 4根轨迹图绘制 (9) 4.1根轨迹数据计算 (9) 4.2用MATLAB软件绘制根轨迹 (10) 5体会与总结.................................. 错误!未定义书签。 5.1总结 ........................................... 错误!未定义书签。 5.2体会 ........................................... 错误!未定义书签。本科生课程设计成绩评定表.. (13)

二阶系统时域分析

专业:电气工程及其自动化 学号:07050443 05 姓名: 实验一 二阶系统时域分析 一、 实验目的 1. 研究二阶系统的两个重要参数ξ、n ω与系统结构之间的关系。 2. 观察系统在阶跃输入作用下的响应,运用基本理论,分析系统过度过程特点及各种参数对其学习过程的影响,从而找出改善系统动态性能的方法,并在实验中加以验证。 3. 学习二阶系统阶跃响应的测试方法。 4. 掌握开环传递函数与闭环传递函数之间的对应关系,以及ξ、n ω与传递函数系数之间的关系。 二、 实验内容 选择适当的元器件建立单位负反馈二阶系统。 开环传递函数由积分环节和惯性环节构成:()() 1S T S T K S G 21+= 令T T T 21==。 1. 设1T = 改变K 值,使阻尼比ξ,分别为0、0.5、0.7、1、1.5;观察并记录在单位阶跃信号作用下,不同阻尼比时,系统输出响应曲线,并测量系统的超调量σ%、上升时间r t 、峰值时间p t 、调节时间s t 。 (1)当阻尼比ξ无限大时: (2)当阻尼比ξ=0.5时:

(3)当阻尼比ξ=0.7时: (4)当阻尼比ξ=1时: (5)当阻尼比ξ=1.5时:

2. 设定K 值 使ξ=0.707,改变时间常数T ,观察并记录在单位阶跃信号作用下,系统输出曲线,并测量系统的超调量σ%、上升时间r t 、峰值时间p t 、调节时间s t 。并与(1)的结果加以比较。 (1) 当T=0.1时: (2) 当T=1时:

(3) 当T=1.5时: 3. 改变时间常数 使1T 不等于2T ,观察并记录输出波形的变化情况。 (1) 当1T 1=,2T 2=时: (2) 当2T 1=,1T 2=时:

高阶系统的时域分析

课程设计任务书 学生姓名: 专业班级: 自动化1002班 指导教师: 肖纯 工作单位: 自动化学院 题 目: 高阶系统的时域分析 初始条件:设单位系统的开环传递函数为 ) )(105() ()(2a s s s s b s K s G ++++= 要求完成的主要任务: (包括课程设计工作量及其技术要求,以及说明书撰写等 具体要求) (1) 当K=10,a=1,b=5时用劳斯判据判断系统的稳定性。 (2) 如稳定,则求取系统的单位阶跃响应、单位斜坡响应和单位加速度响应,用 Matlab 绘制相应的曲线,并计算单位阶跃响应的动态性能指标和稳态性能指标,计算单位斜坡响应和单位加速度响应的稳态性能指标。 (3) 如不稳定,则计算系统稳定时K 、a 和b 的取值范围,在稳定范围内任取一值 重复第2个要求。 (4) 绘制稳定时系统的根轨迹(在稳定范围内任取a 、b 值)。分析K 变化对系统 性能的影响。 时间安排: 指导教师签名: 年 月 日 系主任(或责任教师)签名: 年 月 日

目录 1 高阶系统的数学模型 (1) 2 系统稳定性分析 (2) 3 高阶系统的时域分析 (5) 3.1 单位阶跃响应 (5) 3.1.1 单位阶跃响应 (5) 3.1.2 单位阶跃响应动态性能 (7) 3.1.3 单位阶跃响应稳态性能 (8) 3.2 单位斜坡响应 (9) 3.2.1 单位斜坡响应 (9) 3.2.2 单位斜坡响应稳态性能 (10) 3.3 单位加速度响应 (11) 3.3.1 单位加速度响应 (11) 3.3.2 单位加速度响应稳态性能 (12) 4 系统根轨迹 (13) 5 设计心得体会 (14) 参考文献 (14)

二阶系统的时域分析

实验三 二阶系统的时域分析 一、实验目的 1、通过考察系统的过渡过程指标,研究二阶系统的特征参数—阻尼比和自然频率对系统特性的影响,以及系统特征根的位置与过渡过程的关系。 2、学习自己设计实验,安排适当的实验参数,达到以上实验目标。 二、实验内容 根据传递函数2 22 2)(n n n s s s G ωζωω++=的单位阶跃响应,求取过渡过程的质量指标。按表1的形式整理实验数据,分析实验结果,完成实验报告。 此时,系统的特征根为j j s n n βαζωζω±=-±-=2 2,11。 1、令ζ=0.5,取三种不同的n ω,观察根在根平面上的位置,求其过渡过程和它的质量指

标,进行比较。说明当ζ相同时,过渡过程的哪些指标是相同的? 00.2 0.4 0.6 0.8 1 1.2 1.4 ωn 改变,ζ=0.5不变 Tim e (sec) A m p l i t u d e

2、固定n ω,取ζ=0、0. 3、 0.5、0.7、1,观察根在根平面上的位置,求其过渡过程和它的质量指标。总结当ζ不同时,质量指标有哪些变化? 24681012141618 00.20.40.60.811.2 1.41.61.82 Time (sec) A m p l i t u d e

通过上面两图形与表格总结可以得出: n ω影响二阶系统过渡过程中的峰值时间,过渡时间(在ζ不变的情况下,峰值时间随n ω增 大而减小,过渡时间随n ω的增大而减小) ζ影响几乎全部过渡过程指标,其中超调量,衰减比仅与ζ有关(超调量随着ζ的增大而 减小,衰减比随着ζ的增大而增大;在n ω不变的情况下,峰值时间随ζ增大而增大,过渡时间随ζ的增大而减小。) n ω,ζ对系统的稳态误差均没有影响,且均为0.

实验三基于SIMULINK控制系统时域分析

实验三基于SIMULINK的控制系统时域分析(2学时)一.实验目的: 掌握使用SIMULINK、控制工具箱求解系统的输入和输出响应的仿真方法。二.实验方法及预习内容: 利用SIMULINK工具进行控制系统模型分析、系统设计与仿真的相关原理。三.实验内容: 1.分别使用解微分方程方法、控制工具箱、Simulink求解具有如下闭环传递函数的系统的阶跃响应。 解微分方程方法求解: Wffc.m文件: function dx=wffc(t,x) u=1; dx=[-8*x(1)-36*x(2)-40*x(3)-10*x(4)+u;x(1);x(2);x(3)]; 主文件: %-------------------------------------实验1.1 [t,x]=ode45('wffc',[0,18],[0;0;0;0]) y=10*x(:,4); figure(2) plot(t,y) grid on title('解微分方程方法求解系统阶跃响应曲线') xlabel('时间') ylabel('输出') 结果: 控制工具箱求解:

程序: num=[10]; den=[1 8 36 40 10]; sys=tf(num,den); step(sys) title('控制工具箱求解系统阶跃响应曲线') xlabel('时间') ylabel('输出') grid 结果: Simulink 求解 设置: 单位阶跃 仿真设置 SCOPE 设置 程序 figure(3) plot(ScopeData(:,1),ScopeData(:,2)) grid on title('Simulink 求解系统阶跃响应曲线') xlabel('时间') ylabel('输出') 结果 2. 某小功率随动系统动态结构如图所示,已知:10.01T =, 20.05T =,01K =,1300K =,21K =,0.08c K =。若系统输入分别为()1()sr t t θ=,sr t θ=,[1()1(1.5)]sr t θ=-,试用Simulink 分析系统的输出()sc t θ分别如何? 设置: 单位阶跃:

实验七 控制系统的时域分析方法

实验七 控制系统频域分析方法 1.实验目的 (1)熟练掌握Nyquist 图和Bode 图的绘制。 (2)熟练掌握利用Nyquist 图和Bode 图分析系统的性能。 2.实验仪器 (1)Matlab6.5应用软件安装版 一套 (3)PC 机 一台 3. 实验原理 依据MA TLAB 的建模指令,利用MATLAB 对系统仿真,分析系统的频率特性。 4. 实验步骤 (1)建立系统的MATLAB 模型,绘制系统Nyquist 图和Bode 图,分析系统稳定性 (2)求系统的幅值穿越频率和相位穿越频率,分析系统的稳定性。 (3)依据系统框图建立系统模型,利用LTI Viewer 分析系统的稳定性。 (4)绘制离散系统开环传递函数的Nyquist 图和Bode 图,绘制系统单位阶跃响应图。 5. 实验报告内容(选做其中三题) 1、绘制下列各单位反馈系统开环传递函数的Bode 图和Nyquist 图,并根据其稳定裕度判断系统的稳定性。(使用subplot 指令) ) 31)(2s 1)(s 1(10)s (G 1k s +++=)( )101)(s 1(s 10)s (G 2k s ++= )( ) 2.01)(s 1.01(s 10)s (G 32k s ++=)( )101)(s 1.01(s 10)s (G 42k s ++= )( 2、设单位反馈系统的开环传递函数为)12s (s K )s (G 2k ++=n n w s w ξ,其中无阻尼固有频率 Wn=90rad/s ,阻尼比ξ=0.2,试确定是系统稳定的K 的范围。 3、设系统如图7-22所示,试用LTI Viewer 分析系统的稳定性,并求出系统的稳定裕度及单位阶跃响应峰值. 4、设闭环离散系统结构如图7-23所示,其中) 1(10s +=s s G )(,1s =)(H ,绘制T=0.01s,1s 时离散系统开环传递函数的Bode 图和Nyquist 图,以及系统的单位阶跃响应曲线..

控制系统的时域分析实验报告

一、实验目的和要求 1.用计算机辅助分析的办法,掌握系统的时域分析方法。 2.熟悉Simulink 仿真环境。 二、实验内容和原理 (一)实验原理 系统仿真实质上就是对系统模型的求解,对控制系统来说,一般模型可转化成某个微分方程或差分方程表示,因此在仿真过程中,一般以某种数值算法从初态出发,逐步计算系统的响应,最后绘制出系统的响应曲线,进而可分析系统的性能。控制系统最常用的时域分析方法是,当输入信号为单位阶跃和单位冲激函数时,求出系统的输出响应,分别称为单位阶跃响应和单位冲激响应。在MA TLAB 中,提供了求取连续系统的单位阶跃响应函数step ,单位冲激响应函数impulse ,零输入响应函数initial 等等。 (二)实验内容 二阶系统,其状态方程模型为 ?1x -0.5572 -0.7814 1x 1 = + u ?2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 1.画出系统的单位阶跃响应曲线; 2.画出系统的冲激响应曲线; 3.当系统的初始状态为x0=[1,0]时,画出系统的零输入响应; 4.当系统的初始状态为零时,画出系统斜坡输入响应; (三)实验要求 1.编制MA TLAB 程序,画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应; 2.在Simulink 仿真环境中,组成系统的仿真框图,观察单位阶跃响应曲线并记录之。 三、主要仪器设备 计算机一台以及matlab 软件,simulink 仿真环境 四、操作方法与实验步骤 1、程序解决方案:

在MATLAB命令窗口中输入下列命令:并返回系统的传递函数 其输出的曲线如下

关于典型二阶系统的时域分析10页

林美花(1班)学号:200900192029 二、1:. 在过阻尼情况下,典型二阶系统有两个相异的实数极点,其阶跃响应实际上是两个一阶系统响应的叠加。请以例【3-1】中的系统为例(ωn=5),不断增大ζ值,观察每个ζ值下两个实数极点间的距离;同时绘出两个实数极点分别对应的一阶系统响应和二阶系统的响应,观察它们间的关系。你能得出什么结论?为什么? 解:(1)根据理论推算两实数极点之间的距离为2*ωn*(ζ2-1)0.5 ,所以增大ζ值,两个实数极点间的距离随之增大。 (2)源程序如下: clc; clear; wn=5; num=wn^2; zeta=[1.1:0.1:2.0]; for i=1:10 figure(i) hold on s1=-zeta(i)*wn+wn*(zeta(i)^2-1)^0.5; s2=-zeta(i)*wn-wn*(zeta(i)^2-1)^0.5; num1=wn^2/(s1-s2); num2=-wn^2/(s1-s2); den=[1,2*zeta(i)*wn,wn^2];

step(num,den) den=[1,-s1]; step(num1,den) den=[1,-s2]; step(num2,den) hold off end title('stepresponse')

结论:在过阻尼的状态下,由图像可知其阶跃响应实际上是两个一阶系统响应的叠加。随着ζ的不断增加,一个极点不断靠近原点,另一个不断远

离。当两个极点相距较近时,对阶跃响应产生的影响都不能忽略。ζ的增大使不断远离原点的极点所产生的影响越来越小,最后趋近于零。当两个极点的绝对值之比达到某一倍数(五倍)以上时,则可以忽略离虚轴较远的极点的影响,将二阶系统近似为一阶系统来考虑。同理,在考虑高阶问题时可以找到主导极点,可以降阶处理,化简运算。 二、2:请绘制出图3-21。根据典型二阶系统的脉冲响应,可以分析出系统的哪些暂态性能指标,为什么? 解: clc; clear; wn=5; num=wn^2; zeta=[0.1:0.2:0.7,1.0]; figure(1) hold on for i=1:5 den=[1,2*zeta(i)*wn,wn^2]; impulse(num,den) end hold off title('stepresponse')

自动控制原理实验报告《线性控制系统时域分析》讲述

实验一线性控制系统时域分析 1、设控制系统如图1 所示,已知K=100,试绘制当H分别取H=0.1 ,0.2 0.5,1, 2,5,10 时,系统的阶跃响应曲线。讨论反馈强度对一阶 系统性能有何影响? 图1 答: A、绘制系统曲线程序如下: s=tf('s'); p1=(1/(0.1*s+1)); p2=(1/(0.05*s+1)); p3=(1/(0.02*s+1)); p4=(1/(0.01*s+1)); p5=(1/(0.005*s+1)); p6=(1/(0.002*s+1)); p7=(1/(0.001*s+1)); step(p1);hold on; step(p2);hold on; step(p3);hold on; step(p5);hold on; step(p6);hold on; step(p7);hold on;

B 、绘制改变H 系统阶跃响应图如下: 00.050.10.150.20.250.30.350.40.450.5 0.2 0.4 0.6 0.8 1 1.2 1.4 Step Response Time (seconds) A m p l i t u d e 结论: H 的值依次为0.1、0.2、0.5、1、2、5、10做响应曲线。matlab 曲线默认从第一条到第七条颜色依次为蓝、黄、紫、绿、红、青、黑,图中可知随着H 值得增大系统上升时间减小,调整时间减小,有更高的快速性。 2、 二阶系统闭环传函的标准形式为 22 2 ()2n n n s s s ωψξωω=++,设已知 n ω=4,试绘制当阻尼比ξ分别取0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 5 等值时,系统的单位阶跃响应曲线。求出ξ取值 0.2 ,0.5 ,0.8时的超调量,并求出ξ取值 0.2 ,0.5 ,0.8,1.5,5时的调节时间。讨论阻尼比变化对系统性能的影响。

控制系统设计

附录1:设计报告封面 辽宁工程技术大学 二○一五~二○一六学年第一学期《自动控制原理》综合训练项目二 报告书 题目:控制系统设计 班级:电网13-1 学号:1305080116 姓名:苏小平 指导教师:邱彬 二○一五年十二月

附录Ⅱ:设计报告内封——设计任务书

一题目、任务及要求 已知:某单位负反馈系统的开环传递函数 ()()()1110K G s H s s s s = ?? ++ ? ?? ,利用时 域分析法,根轨迹法和频率分析法分析系统,并用频域法设计控制器。 1、 第一阶段设计任务:时域分析法 1) 利用劳斯判据,求K 的稳定域;若要保持系统特征根都在1s =-的左侧,求K 的取值 范围。 2)计算当5K =时,将此三阶系统近似为二阶系统,并计算系统的暂态特性指标,超调量,调节时间。 3) 当输入为单位斜坡信号时,计算系统的稳态误差。 2、第二阶段设计任务:根轨迹法 1)绘制系统根轨迹图。 2)利用根轨迹图,求系统临界稳定时的K 值。 3 )应用主导极点的概念,计算ξ=时,系统的超调量和调节时间。 3、第三阶段设计任务:频域分析法 1) 当10K =时,绘制系统开环幅相特性,应用奈奎斯特稳定判据判定系统的稳定性。 2) 当10K =时,绘制系统对数频率特性曲线,并计算相角稳定裕量和幅值稳定裕量,并根据相角裕量和幅值裕量判定系统稳定性。 3) 4、第四阶段设计任务:频域法校正设计 要求系统在单位斜坡输入下稳态误差小于0.02,且相位裕量()45c γω≥ ,请利用串联 校正方法,设计控制器()c G s ,并写出控制器的实现方式。

二报告及书写内容要求 课程设计任务完成后,每位同学必须独立书写一份课程设计报告,注意:不得抄袭他人的报告(或给他人抄袭),一旦发现,成绩为零分。 (一)内容要求 1)分析过程完整,计算正确,利用MATLAB进行结果验证; 2)设计要求,说明串联校正的设计思路(滞后校正,超前校正); 3)详细设计(包括的图形有:串联校正结构图,校正前系统的Bode图,校正装置的Bode 图,校正后系统的Bode图); 4)MATLAB编程代码及运行结果(包括图形、运算结果); 5)校正实现的电路图及仿真实验结果(校正前后系统的阶跃响应图-MATLAB或SMULINK 辅助设计); 6)校正前后的系统性能指标的计算; 7)总结(包括课程中的学习体会与收获、对本次综合训练项目设计的认识等内容)。(二)格式字体要求 1) A4文档输出(保留电子文档); 2)每部分另起一行,1.5倍行距; 3)公式需用公式编辑器编辑

高阶系统的时域分析(课程设计)

课程设计任务书 学生姓名: 专业班级: 指导教师: 肖 纯 工作单位: 自动化学院 题 目: 高阶系统的时域分析 初始条件:设单位系统的开环传递函数为2 () ()(48)() p K s b G s s s s s a +=+++ 要求完成的主要任务: (包括课程设计工作量及其技术要求,以及说明书撰写等具体要求) (1) 当K=10,a=1,b=4时用劳斯判据判断系统的稳定性。 (2) 如稳定,则求取系统的单位阶跃响应、单位斜坡响应和单位加速度响应,用 Matlab 绘制相应的曲线,并计算单位阶跃响应的动态性能指标和稳态性能指标,计算单位斜坡响应和单位加速度响应的稳态性能指标。 (3) 如不稳定,则计算系统稳定时K 、a 和b 的取值范围,在稳定范围内任取一值 重复第2个要求。 (4) 绘制a=1,b=4时系统的根轨迹。 时间安排: 指导教师签名: 年 月 日 系主任(或责任教师)签名: 年 月 日

目录 1 高阶系统的数学模型 (1) 2 系统稳定性分析 (1) 3 高阶系统的时域分析 (3) 3.1 单位阶跃响应 (4) 3.1.1 求单位阶跃响应 (4) 3.1.2 单位阶跃响应动态性能 (7) 3.1.3 单位阶跃响应稳态性能 (9) 3.2 单位斜坡响应 (10) 3.2.1 求单位斜坡响应 (10) 3.2.2 单位斜坡响应稳态性能 (11) 3.3 单位加速度响应 (11) 3.3.1 求单位加速度响应 (11) 3.3.2 单位加速度响应稳态性能 (13) 4 系统根轨迹 (13) 5 设计心得体会 (15) 参考文献 (15)

高阶系统的时域分析 1 高阶系统的数学模型 一个高阶系统的闭环传递函数的一般形式为: 10111011()(),()m m m m n n n n b s b s b s b C s s m n R s a s a s a a ----++++Φ==≤++++ 对分子、分母进行因式分解,得到零极点形式: 11 () () ()() () m i i n j j K s z C s s R s s p ==-Φ= =-∏∏ (1) 式(1)中,K=b 0/a 0;z i ,p j 分别为系统闭环零、极点。 本设计给定的单位反馈系统的开环传递函数为 2 () ()(48)()p K s b G s s s s s a +=+++ (2) 则其闭环传递函数为(假设为负反馈): 2432()() ()(48)()()(4)(84)(8)K s b K s b s s s s s a K s b s a s a s a K s Kb ++Φ==++++++++++++ (3) 2 系统稳定性分析 线性系统稳定的充分必要条件为:闭环系统特征方程的所有根均具有负实部;或者说,闭环传递函数的极点均位于s 左半平面。 若求出闭环系统特征方程的所有根,就可判定系统的稳定性。但对于高阶系统来说,求特征方程根很困难,并且不易对参数进行分析。现使用一种不用求解特征根来判别系统稳定性的方法—劳斯稳定判据。 设系统的特征方程为10110()0,0n n n n D s a s a s a s a a --=++++=>,则可列出劳斯表如 表1所示。

控制系统的时域分析

实验二控制系统的时域分析 一、实验目的 学习利用MATLAB 进行控制系统时域分析,包括典型响应、判断系统稳定性和分析系统的动态特性; 二、预习要点 1、系统的典型响应有哪些? 2、如何判断系统稳定性? 3、系统的动态性能指标有哪些? 三、实验方法 (一)典型响应 1、阶跃响应: 阶跃响应常用格式: 1、)(sys step 或[y,t]=step(sys);其中sys 可以为连续系统,也可为离散系统。 2、),(Tn sys step ;表示时间范围0---Tn。 3、),(T sys step ;表示时间范围向量T 指定。 4、),(T sys step Y =;可详细了解某段时间的输入、输出情况。例:假设一连续模型为:s e s s s s s s G -+++++=10 232623102010)(234,则可以通过下面的命令直接输入系统模型,并绘制出阶跃响应曲线。 解:>>num=[0,0,0,10,20]; >>den=[10,23,26,23,10]; >>G=tf(num,den); >>G.iodelay=1; >>step(G,30)%终止时间为30。 2、脉冲响应: 脉冲响应函数常用格式:①)(sys impulse ; ②); ,(); ,(T sys impulse Tn sys impulse ③) ,(T sys impulse Y =3、任意输入响应: 任意输入响应的几种常用格式: ),,(T U sys lsim ;其中sys 可为任意模型;T 为时间向量;U 为响应时间对应的系统输入,例 如:)sin(T U =; (二)分析系统稳定性 有以下三种方法: 1、利用pzmap 绘制连续系统的零极点图;pzmap(G); 2、Pole(G)和zero(G)可以分别求出系统的极点和零点。 3、利用roots 求分母多项式的根来确定系统的极点。roots(den). (三)系统的动态特性分析 方法一:图解法 在控制理论中,介绍典型线性系统的阶跃响应分析时,常用一些指标来定量描述系统的超调

相关文档
最新文档