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

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

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

经典控制理论中,系统分析与校正方法一般有时域法、复域法、频域法。时域响应法是

一种直接法,它以传递函数为系统的数学模型,以拉氏变换为数学工具,直接可以求出变量

的解析解。这种方法虽然直观,分析时域性能十分有用,但是方法的应用需要两个前提,一是必须已知控制系统的闭环传递函数,另外系统的阶次不能很高。

如果系统的开环传递函数未知,或者系统的阶次较高,就需采用频域分析法。频域分析

法不仅是一种通过开环传递函数研究系统闭环传递函数性能的分析方法,而且当系统的数学

模型未知时,还可以通过实验的方法建立。此外,大量丰富的图形方法使得频域分析法分析

高阶系统时,分析的复杂性并不随阶次的增加而显著增加。

在进行控制系统分析时,可以根据实际情况,针对不同数学模型选用最简洁、最合适的方法,从而使用相应的分析方法,达到预期的实验目的。

系统的时域性能指标与频域性能指标有着很大的关系,研究其内在联系在工程中有着很大的意义。

一、系统的时域性能指标

延迟时间t d 阶跃响应第一次达到终值h()的50%所需的时间

上升时间t r阶跃响应从终值的10%上升到终值的90%所需的时间;对有振荡的系

统,也可定义为从0到第一次达到终值所需的时间

峰值时间t p 阶跃响应越过终值h()达到第一个峰值所需的时间

调节时间t s 阶跃响应到达并保持在终值h()的5%误差带内所需的最短时间

t p)超出终值h()的百分比,即

超调量% 峰值h(

h t p h

%= 100%

h

二、系统频率特性的性能指标

采用频域方法进行线性控制系统设计时,时域内采用的诸如超调量,调整时间等描述系

统性能的指标不能直接使用,需要在频域内定义频域性能指标。

1零频振幅比M(0):即 为0时闭环幅频特性值。它反映了系统 的稳态精度,M(0)越接

近于1,系统的精度越高。

M(0)工1时,表明系统有稳态误差。

2、谐振峰值 Mr :为幅频特性曲线的 A()的最大值。一般说来,

Mr 的大小表明闭环控制

系统相对稳定性的好坏。Mr 越大,表明系统对某个频率的正弦信号反映强烈,

有共振倾向,

系统的平稳性较差,相应阶跃响应的超调量越大。对应的 r 为谐振频率。

3、 谐振频率

r

:出现最大值 Mmax 时对应的频率。

4、 带宽 b

幅频特性下降至零频幅比的

70.7 %,或下降3dB 时对应的频率称为带宽(也成为闭环

截止频率)。带宽用于衡量控制系统的快速性,带宽越宽,表明系统复现快速变化信号的能 力越强,阶跃响应的上升时间和调节时间就越短。 带宽是控制系统及控制元件的重要性能指

标。

三、闭环频域性能指标与时域性能指标的关系

1、二阶系统的相互联系

对于二阶系统,其频域性能指标和时域性能指标之间有着严格的数学关系

(1)、谐振峰值 Mr 和时域超调量 之间的关系

幅频特性的谐振峰值 Mr

响应指标的要求。

在二阶系统 (s)=-

2

s

2

n S

2中,

()

=(

2)2+(2

n

)2

令 dM()

V

d

=0 , 得谐振频率

n 1-2 2

求得幅频特性峰值

M r =

二阶系统的超调量

-/

%=e

100%

由此可看出,谐振峰值 Mr 仅与阻尼比

有关,超调量 %也仅取决于阻尼比

越小,Mr 增加的越快,这时超调量

%也很大,超过40%,—般这样的系统不符和瞬态

当0.4< <0.707时,Mr 与 %的变化趋势基本一致,此时谐振峰值 Mr=1.2 ~ 1.5,超调量

% =20% ~30%,系统响应结果较满意。

当>0.707时,无谐振峰值,Mr 与 %的对应关系不再存在,通常设计时, 取在0.4至

0.7之间

由此可看出,当 为常数时,谐振频率

r

与峰值时间 t p 成反比, r 值愈大,t

p

愈小,表示系统时间响应愈快

(3) 、闭环谐振峰值Mr 和相角裕度 的关系

(j )=M( )e i ()

G(j )=A( )e j<)

_ j

< c) j(180°-)

G(i c

)=A( c

)e =A( c

)e =A( c

)(-cos -jsin )

般Mr 极大值发生在 c 附近。

在开环截止频率 c

附近,上述近似程度就越高。

(4 )、 和的关系

(2)、谐振频率

与峰值时间t

p 的关系

t p

=

与 r 之积为

t p r

\1-2一2

M( c )=

+ ( c ) (c

)=1800 -

G(j ) = A( c )

1+G(j ) |1-A( c )cos -jA( c )sin

dMLJ=0

dA()

A()

1 sin

Mr

sin

故Mr

1 sin

= 1800

(2 )、二阶系统

2 n

标准二阶系统的开环传递函数为

G(s)= —

S(S + 2 n )

G(j c )=

n

j c (j c +2

n

)=1

G(j

c)

n

— =1

—=4^+1-2

n

1/ 2 2

=1800+(-900 -arctg

)=90°-arctg

2

n

c

-——二arctg

1/2

得出

=arctg 2

1 ,4 4+1-2

对于二阶系统,一般要求:

30° < <70° 0.27< <0.8

2、带宽

b

与时域性能的关系

(1 )、一阶系统

一阶系统的闭环传递函数为

(s)=

1

1+TS

系统的闭环频率特性为

(j )=

1 1+Tj

系统的闭环幅频特性为

)=(j )

1

、、

1+T

可知,

=0时幅值为1,即零频振幅比 M(0) = 1,

则 L(0) = 20LgM(0) =_0 闭环截止频率

b

:

由b 的定义知 L( b )= L(0)-3=-3

20LgM ( b )=20 Lg [

2 =-3

V1+T b

M ( b)= ^=2

2-

=0.707

可解得:b =1/T

阶系统中调节时间、上升时间与带宽的关系

tr=2.2T ts=3T

tr=2.2/ b

, ts=3/

b

2

二阶系统的闭环传递函数为

(S )二飞

-

S +2

n

S+ n

闭环频率特性为 (j )=

2 2

n

n

2 2 2 2

(j ) +j2

n

+ n

n

+j2

n

-

系统的闭环幅频特性为

2 n

M( )= (j )=

——百

j ( n - ) +(2

n

)

可知,

=0时幅值为1,即零频振幅比 M(0) = 1,

则 L(0) = 20LgM(0) =_0 闭环截止频率 b 由b 的定义知 L( b )= L(0)-3=-3

b

= n 1-2

2

+ .(1-2 2 )2 +1

阻尼比不变,自然振荡频率越大,带宽越大;自然振荡频率不变,阻尼比越小,带宽越大; 可知带宽与系统响应速度成正比!

b

tS

之积为

b

ts=-^(1-2 2 + J - 2+

4

由此可看出,当阻尼比 给定后,闭环截止频率 b 与过渡过程时间t

s 成反比关系。换 言之, b 愈大(频带宽度0 - b 愈宽),系统的响应速度愈快。 (4) 、系统带宽的选择

带宽频率是一项重要指标。其选择要求要既能以所需精度跟踪输入信号,

又能拟制噪声

扰动信号。在控制系统实际运行中,输入信号一般是低频信号,而噪声信号是高频信号。

(5) 、带宽指标取决于下列因素: a)对输入信号的再现能力。大的带宽相应于小的上升时间,即相应于快速特性。粗略地说,

可解得:

b

)=

)2

=0.707

(3)、带宽

b

与调节时间

ts

的关系

3 5 调整时间

ts=2F

n

b

= nj 1-2 2+\(1-2 2)2

+1

2 2

、2

_ b ) +(2

带宽与响应速度成正比。 b )对高频噪声必要的滤波特性。

为了使系统能够精确地跟踪任意输入信号,

系统必须具有大的带宽。 但是,从噪声的观

点来看,带宽不应当太大。因此,对带宽的要求是矛盾的,好的设计通常需要折衷考虑。具 有大带宽的系统需要高性能的元件,因此,元件的成本通常随着带宽的增加而增大。

3、典型二阶系统频域指标与时域指标的关系

4、高阶系统频域指标与时域指标

1

谐振峰值 Mr = ----------

sin

超调量 =0.16+0.4(Mr-1) 1 Mr 1.8 ts=— K=2+1.5(Mr-1)+2.5(Mr-1)2 1 Mr 1.8

c

=J-2

b

= . (1-2 22-4 2+4

=arctg

2 ..4 2+1-2

.4 2+1-2

1/2

2

闭环阶跃响应时域指标:

%=e

100% tp= / d = /( n. 1- 2) tr =( - )/ d =(

ts=^.5

( =0.05,0< <0.9)

n

因此,若知道频域指标中的任两个,就可解算出 n

,从而求出时域指标。反之,给出

时域指标的任两个,

就可确定闭环频域指标。

%=e

100%

Mr

,ts ,相对稳定性好,超调小,振荡次数少。

不变时,n , c b

,tr ,tp

,系统灵敏度下降。

不变时,r

b

,系统灵敏,速度快。

调节时间

闭环频域指标:

c

二阶系统时域分析

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

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

实验名称: 控制系统的频域分析 实验类型:________________同组学生姓名:__________ 一、实验目的和要求 用计算机辅助分析的方法,掌握频率分析法的三种方法,即Bode 图、Nyquist 曲线、Nichols 图。 二、实验内容和原理 (一)实验原理 1.Bode(波特)图 设已知系统的传递函数模型: 1 1211121)(+-+-+???+++???++=n n n m m m a s a s a b s b s b s H 则系统的频率响应可直接求出: 1 1211121)()()()()(+-+-+???+++???++=n n n m m m a j a j a b j b j b j H ωωωωω MATLAB 中,可利用bode 和dbode 绘制连续和离散系统的Bode 图。 2.Nyquist(奈奎斯特)曲线 Nyquist 曲线是根据开环频率特性在复平面上绘制幅相轨迹,根据开环的Nyquist 线,可判断闭环系统的稳定性。 反馈控制系统稳定的充要条件是,Nyquist 曲线按逆时针包围临界点(-1,j0)p 圈,为开环传递函数位于右半s 一平面的极点数。在MATLAB 中,可利用函数nyquist 和dnyquist 绘出连续和离散系统的乃氏曲线。 3.Nicho1s(尼柯尔斯)图 根据闭环频率特性的幅值和相位可作出Nichols 图,从而可直接得到闭环系统的频率特性。在 MATLAB 中,可利用函数nichols 和dnichols 绘出连续和离散系统的Nichols 图。 (二)实验内容 1.一系统开环传递函数为 ) 2)(5)(1(50)(-++=s s s s H 绘制系统的bode 图,判断闭环系统的稳定性,并画出闭环系统的单位冲击响应。 2.一多环系统 ) 10625.0)(125.0)(185.0(7.16)(+++=s s s s s G 其结构如图所示 试绘制Nyquist 频率曲线和Nichols 图,并判断稳定性。 (三)实验要求

2018年自动控制原理期末考试题[附答案解析]

. 2017 年自动控制原理期末考试

卷与答案 一、填空题(每空1分,共20分) 1、对自动控制系统的基本要求可以概括为三个方面,即:稳定性、快速性和准确性。 2、控制系统的输出拉氏变换与输入拉氏变换在零初始条件下的比值称为传递函数。 3、在经典控制理论中,可采用劳斯判据(或:时域分析法)、根轨迹法或奈奎斯特判据( 或:频域分析 法) 等方法判断线性控制系统稳定性。 4、控制系统的数学模型,取决于系统参数, 与外作用及初始条件无关。和结构 A( )L( )lg) 或:,横坐标为( 、线性系统的对数幅频特性,纵坐标取值为5。20lg

,其中P 是指开环传函中具有正实部的极点的个数,6、奈奎斯特稳定判据中,Z = P - RZ是指闭环 传函中具有正实部的极点的个数,R 指奈氏曲线逆时针方向包围(-1, j0 )整圈数。 定义为调整时间。%是超调量。、在二阶系统的单位阶跃响应图中,7t s K)A(22 (T)1K),则其开环幅频特性为(T8、设系统的开环传递函数为1,相12 s(Ts 1)(T s 1)21110) (T ) tg (Ttg 。频特性为()9021 9、反馈控制又称偏差控制,其控制作用是通过给定值与反馈量的差值进行的。 0.5t0.2 t,则该系统的传递函数G(s) 为、若某系统的单位脉冲响应为10。5e g (t) 10e 510 0.2 sss 0.5s

11、自动控制系统有两种基本控制方式,当控制装置与受控对象之间只有顺向作用而无反向联系时,称 为开环控制系统;当控制装置与受控对象之间不但有顺向作用而且还有反向联系时,称为闭环控制系 统;含有测速发电机的电动机速度控制系统,属于闭环控制系统。 、根轨迹起始于开环极点,终止于开环零点。12 、稳定是对控制系统最基本的要求,若一个控制系统的响应曲线为衰减振荡,则该系统13稳定。判断 页脚

一阶系统时域分析

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 ==+?+

自动控制原理-线性系统的频域分析实验报告

自动调节系统频域分析 班级11081801 学号1108180135 姓名王佳炜 日期2014.1.5

线性系统的频域分析 一、实验目的 1.掌握用MATLAB 语句绘制各种频域曲线。 2.掌握控制系统的频域分析方法。 二、实验内容 1.典型二阶系统 2 2 22)(n n n s s s G ωζωω++= 绘制出6=n ω,1.0=ζ,0.3,0.5,0.8,2的bode 图,记录并分析ζ对系统bode 图的影响。 解: 程序如下: num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36]; den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36]; w=logspace(-2,3,100); bode(num,den1,w) grid hold bode(num,den2,w) bode(num,den3,w) bode(num,den4,w) bode(num,den5,w)

-100-80-60-40-200 20M a g n i t u d e (d B )10 -2 10 -1 10 10 1 10 2 10 3 P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 分析:随着.0=ζ的增大 ,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大. 2.系统的开环传递函数为 ) 5)(15(10 )(2 +-= s s s s G ) 106)(15() 1(8)(22++++= s s s s s s G ) 11.0)(105.0)(102.0() 13/(4)(++++= s s s s s s G 绘制系统的Nyquist 曲线、Bode 图,说明系统的稳定性,并通过绘制阶跃响应曲线验证。 解: 程序如下 奈氏曲线: (1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); w=logspace(-1,1,100); nyquist(num1,den1,w)

控制系统的时域分析

实验报告 实验名称:实验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 )不稳定。

自动控制原理线性系统的频域分析实验四

武汉工程大学实验报告专业电气自动化班号指导教师 姓名同组者无

M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Frequency (rad/sec) 当3.0=ζ时,程序如下: num=[0 0 36];den=[1 3.6 36];w=logspace(-2,3,100);bode(num,den,w) grid M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Bode Diagram Frequency (rad/sec) 当5.0=ζ时,程序如下: num=[0 0 36];den=[1 6 36];w=logspace(-2,3,100);bode(num,den,w) grid

M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g )Frequency (rad/sec) 当8.0=ζ时,程序如下: num=[0 0 36];den=[1 9.6 36];w=logspace(-2,3,100);bode(num,den,w) grid M a g n i t u d e (d B ) 10 1010101010P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 当2=ζ时,程序如下: num=[0 0 36];den=[1 24 36];w=logspace(-2,3,100);bode(num,den,w) grid

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

课程名称:控制理论指导老师:成绩: 实验名称:控制系统的时域分析实验类型:冋组学生姓名: 、实验目的和要求 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)

第5章_用MATLAB进行控制系统频域分析

第5章 用MATLAB 进行控制系统频域分析 一、基于MATLAB 的线性系统的频域分析基本知识 (1)频率特性函数)(ωj G 。 设线性系统传递函数为: n n n n m m m m a s a s a s a b s b s b s b s G ++???++++???++=---1101110)( 则频率特性函数为: n n n n m m m m a j a j a j a b j b j b j b jw G ++???++++???++=---)()()()()()()(1101110ωωωωωω 由下面的MATLAB 语句可直接求出G(jw)。 i=sqrt(-1) % 求取-1的平方根 GW=polyval(num ,i*w)./polyval(den ,i*w) 其中(num ,den )为系统的传递函数模型。而w 为频率点构成的向量,点右除(./)运算符表示操作元素点对点的运算。从数值运算的角度来看,上述算法在系统的极点附近精度不会很理想,甚至出现无穷大值,运算结果是一系列复数返回到变量GW 中。 (2)用MATLAB 作奈魁斯特图。 控制系统工具箱中提供了一个MATLAB 函数nyquist( ),该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时,nyquist ()函数仅在屏幕上产生奈氏图,命令调用格式为: nyquist(num,den) nyquist(num,den,w) 或者 nyquist(G) nyquist(G,w) 该命令将画出下列开环系统传递函数的奈氏曲线: ) () ()(s den s num s G = 如果用户给出频率向量w,则w 包含了要分析的以弧度/秒表示的诸频率点。在这些频率点上,将对系统的频率响应进行计算,若没有指定的w 向量,则该函数自动选择频率向量进行计算。 w 包含了用户要分析的以弧度/秒表示的诸频率点,MATLAB 会自动计算这些点的频率响应。 当命令中包含了左端的返回变量时,即: [re,im,w]=nyquist(G) 或

自动控制原理实验六线性系统的频域分析

实验六 线性系统的频域分析 一. 实验目的 (1)熟练掌握使用MA TLAB 命令绘制控制系统Nyquist 图的方法; (2)能够分析控制系统Nyquist 图的基本规律; (3)加深理解控制系统乃奎斯特稳定性判据的实际应用; (4)学会利用奈氏图设计控制系统; (5)熟练掌握运用MA TLAB 命令绘制控制系统伯德图的方法; (6)了解系统伯德图的一般规律及其频域指标的获取方法; (7)熟练掌握运用伯德图分析控制系统稳定性的方法; (8)设计超前校正环节并绘制Bode 图; (9)设计滞后校正环节并绘制Bode 图。 二. 实验原理及内容 1、频率特性函数)(ωj G 。 频率特性函数为: n n n n m m m m a j a j a j a b j b j b j b jw G ++???++++???++= ---)()()()()()()(1101110ωωωωωω 由下面的MATLAB 语句可直接求出G(jw)。 i=sqrt(-1) % 求取-1的平方根 GW=polyval(num ,i*w)./polyval(den ,i*w) 2、用MATLAB 作奈魁斯特图。 控制系统工具箱中提供了一个MA TLAB 函数nyquist( ),该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时,nyquist ()函数仅在屏幕上产生奈氏图,命令调用格式为: nyquist(num,den) ; 作Nyquist 图, nyquist(num,den,w); 作开环系统的奈氏曲线, 3、奈奎斯特稳定性判据(又称奈氏判据) 反馈控制系统稳定的充分必要条件是当ω从-∞变到∞时,开环系统的奈氏曲线不穿过点(-1,j0)且逆时针包围临界点(-1,j0)点的圈数R 等于开环传递函数的正实部极点数。 4、用MATLAB 作伯德图 控制系统工具箱里提供的bode()函数可以直接求取、绘制给定线性系统的伯德图。 命令的调用格式为: [mag,phase,w]=bode(num,den) [mag,phase,w]=bode(num,den,w) 由于伯德图是半对数坐标图且幅频图和相频图要同时在一个绘图窗口中绘制,因此,要用到半对数坐标绘图函数和子图命令。 (1) 对数坐标绘图函数 利用工作空间中的向量x ,y 绘图,要调用plot 函数,若要绘制对数或半对数坐标图,只需要用相应函数名取代plot 即可,其余参数应用与plot 完全一致。 (2) 子图命令

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

控制系统时域与频域性能指标的联系 经典控制理论中,系统分析与校正方法一般有时域法、复域法、频域法。时域响应法是一种直接法,它以传递函数为系统的数学模型,以拉氏变换为数学工具,直接可以求出变量的解析解。这种方法虽然直观,分析时域性能十分有用,但是方法的应用需要两个前提,一是必须已知控制系统的闭环传递函数,另外系统的阶次不能很高。 如果系统的开环传递函数未知,或者系统的阶次较高,就需采用频域分析法。频域分析法不仅是一种通过开环传递函数研究系统闭环传递函数性能的分析方法,而且当系统的数学模型未知时,还可以通过实验的方法建立。此外,大量丰富的图形方法使得频域分析法分析高阶系统时,分析的复杂性并不随阶次的增加而显著增加。 在进行控制系统分析时,可以根据实际情况,针对不同数学模型选用最简洁、最合适的方法,从而使用相应的分析方法,达到预期的实验目的。 系统的时域性能指标与频域性能指标有着很大的关系,研究其内在联系在工程中有着很大的意义。 一、系统的时域性能指标 延迟时间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% 二、系统频率特性的性能指标 采用频域方法进行线性控制系统设计时,时域内采用的诸如超调量,调整时间等描述系统性能的指标不能直接使用,需要在频域内定义频域性能指标。

自动控制原理线性系统的频域分析实验报告

实验四 专业 自动化 班号 03班 指导教师 陈艳飞 姓名 胡波 实验名称 线性系统的频域分析 实验日期 第 次实验 一、实验目的 1.掌握用MATLAB 语句绘制各种频域曲线。 2.掌握控制系统的频域分析方法。 二、实验内容 1.典型二阶系统 2 2 22)(n n n s s s G ωζωω++= 绘制出6=n ω,1.0=ζ,0.3,0.5,0.8,2的bode 图,记录并分析ζ对系统bode 图的影响。 解: 程序如下: num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36]; den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36]; w=logspace(-2,3,100); bode(num,den1,w) grid hold bode(num,den2,w) bode(num,den3,w) bode(num,den4,w) bode(num,den5,w)

-100-80-60-40-200 20M a g n i t u d e (d B )10 -2 10 -1 10 10 1 10 2 10 3 P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 分析:随着.0=ζ的增大 ,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大. 2.系统的开环传递函数为 ) 5)(15(10 )(2+-= s s s s G ) 106)(15() 1(8)(22++++= s s s s s s G ) 11.0)(105.0)(102.0() 13/(4)(++++= s s s s s s G 绘制系统的Nyquist 曲线、Bode 图和Nichols 图,说明系统的稳定性,并通过绘制阶跃响应曲线验证。 解: 程序如下 奈氏曲线: (1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); w=logspace(-1,1,100); nyquist(num1,den1,w)

二阶系统时域分析

专业:电气工程及其自动化 学号: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=时:

自动控制原理实验报告线性系统的频域分析讲述

武汉工程大学实验报告 专业 自动化 班号 组别 指导教师 姓名 同组者 实验名称 线性系统的频域分析 实验日期 2016/4/4 第 5 次实验 一、实验目的 1.掌握用MATLAB 语句绘制各种频域曲线。 2.掌握控制系统的频域分析方法。 二、实验内容 1.典型二阶系统 2 2 22)(n n n s s s G ωζωω++= 绘制出6=n ω,1.0=ζ,0.3,0.5,0.8,2的bode 图,记录并分析ζ对系统bode 图的影响。 解: 程序如下: num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36]; den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36]; w=logspace(-2,3,100); bode(num,den1,w) grid hold bode(num,den2,w)

bode(num,den3,w) bode(num,den4,w) bode(num,den5,w) -100-80-60-40-200 20M a g n i t u d e (d B )10 10 10 10 10 10 P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 分析:随着.0=ζ的增大 ,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大. 2.系统的开环传递函数为 ) 5)(15(10 )(2+-= s s s s G ) 106)(15() 1(8)(2 2++++= s s s s s s G ) 11.0)(105.0)(102.0() 13/(4)(++++= s s s s s s G 绘制系统的Nyquist 曲线、Bode 图,说明系统的稳定性,并通过绘制阶跃响应曲线验证。 解: 程序如下 奈氏曲线: (1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); w=logspace(-1,1,100);

二阶系统的时域分析

实验三 二阶系统的时域分析 一、实验目的 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.

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

实验七 控制系统频域分析方法 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')

相关文档
最新文档