1 请解释下列名字术语:自动控制系统、受控对象、扰动、给定值、参考输入、反馈。
解:自动控制系统:能够实现自动控制任务的系统,由控制装置与被控对象组成; 受控对象:要求实现自动控制的机器、设备或生产过程
扰动:扰动是一种对系统的输出产生不利影响的信号。如果扰动产生在系统内部称为内扰;扰动产生在系统外部,则称为外扰。外扰是系统的输入量。 给定值:受控对象的物理量在控制系统中应保持的期望值 参考输入即为给定值。
反馈:将系统的输出量馈送到参考输入端,并与参考输入进行比较的过程。
2 请说明自动控制系统的基本组成部分。
解: 作为一个完整的控制系统,应该由如下几个部分组成: ① 被控对象: 所谓被控对象就是整个控制系统的控制对象;
② 执行部件: 根据所接收到的相关信号,使得被控对象产生相应的动作;常
用的执行元件有阀、电动机、液压马达等。
③ 给定元件: 给定元件的职能就是给出与期望的被控量相对应的系统输入量(即参考量);
④ 比较元件: 把测量元件检测到的被控量的实际值与给定元件给出的参考值
进行比较,求出它们之间的偏差。常用的比较元件有差动放大器、机械差动装置和电桥等。
⑤ 测量反馈元件:该元部件的职能就是测量被控制的物理量,如果这个物理量
是非电量,一般需要将其转换成为电量。常用的测量元部件有测速发电机、热电偶、各种传感器等;
⑥ 放大元件: 将比较元件给出的偏差进行放大,用来推动执行元件去控制被
控对象。如电压偏差信号,可用电子管、晶体管、集成电路、晶闸管等组成的电压放大器和功率放大级加以放大。
⑦ 校正元件: 亦称补偿元件,它是结构或参数便于调整的元件,用串联或反
馈的方式连接在系统中,用以改善系统的性能。常用的校正元件有电阻、电容组成的无源或有源网络,它们与原系统串联或与原系统构成一个内反馈系统。
3 请说出什么是反馈控制系统,开环控制系统和闭环控制系统各有什么优缺点?
解:反馈控制系统即闭环控制系统,在一个控制系统,将系统的输出量通过某测量机构对其进行实时测量,并将该测量值与输入量进行比较,形成一个反馈通道,从而形成一个封闭的控制系统;
开环系统优点:结构简单,缺点:控制的精度较差;
闭环控制系统优点:控制精度高,缺点:结构复杂、设计分析麻烦,制造成本高。 4 请说明自动控制系统的基本性能要求。
解:(1)稳定性:对恒值系统而言,要求当系统受到扰动后,经过一定时间的调整能够回到原来的期望值。而对随动系统而言,被控制量始终跟踪参考量的变化。稳定性通常由系统的结构决定的,与外界因素无关,系统的稳定性是对系统的基本要求,不稳定的系统不能实现预定任务。
(2)准确性:控制系统的准确性一般用稳态误差来表示。即系统在参考输入信号作用下,系统的输出达到稳态后的输出与参考输入所要求的期望输出之差叫做给定稳态误差。显然,这种误差越小,表示系统的输出跟随参考输入的精度越高。 (3)快速性:对过渡过程的形式和快慢的要求,一般称为控制系统的动态性能。系统的快速性主要反映系统对输入信号的变化而作出相应的快慢程度,如稳定高射炮射角随动系统,虽然炮身最终能跟踪目标,但如果目标变动迅速,而炮身行动迟缓,仍然抓不住目标。
2-1 设质量-弹簧-摩擦系统如图2-1所示,途中f为黏性摩擦系数,k为弹簧系数,系统的输入量为力p(t),系统的输出量为质量m的位移x(t)。试列出系统的
图2-1 习题2-1 质量-弹簧-摩擦系统示意图
输入输出微分方程。
解:显然,系统的摩擦力为fdx(t),弹簧力为kx(t),根据牛顿第二运动定律有 dtdx(t)d2x(t)p(t)fkx(t)m
dtdt2移项整理,得系统的微分方程为
d2x(t)dx(t)mfkx(t)p(t) 2dtdt
2-2 试列写图2-2所示机械系统的运动微分
方程。
解:由牛顿第二运动定律,不计重力时,得
k[yd2y1dy22(t)y1(t)]M1dt2k1y1f1dtF
整理得
Md2y1dy11dt2fdt(k1k2)y1(t)Fk2y2(t)
2-3 求下列函数的拉氏变换。 (1)f(t)3(1sint) (2)f(t)teat (3)f(t)cos(3t4)
解:(1)L[f(t)]L[3(1sint)]
3(L[1]L[sint])3(11ss21) 3(s2s1)s(s21)(2)f(t)teat
L[t]1s2 L[f(t)]L[teat]1(sa)2
(3)f(t)cos(3t24)2[sin(3t)cos(3t)] L[f(t)]22[sin(3t)cos(3t)]
图2-2 习题2-2 机械系统示意图
2(L[sin(3t)]L[cos(3t)])223s(22) 2s9s92s32s29
2-4 求下列函数的拉氏反变换 (1)F(s)s1
(s2)(s5)s6 s2(s3)(2)F(s)2s25s1(3)F(s)
s(s21)解:(1)F(s)s112
(s2)(s5)s2s5L1[F(s)]L1[12] s2s512]2L1[]s2s5 2e5tL1[e2t(2)F(s)s6211
s2(s3)s2ss3L1[F(s)]L1[211] s2ss311111]L[]L[]2sss3 2t1e3t2L1[2s25s11s52(3)F(s)
s(s21)ss11s5L1[F(s)]L1[2]
ss11s5L1[]L1[2]ss1 1cost5sint
2-5 试分别列写图2-3中各无源网络的微分方程(设电容C上的电压为uc(t),电容C1上的电压为uc1(t),以此类推)。
+uc(t)C-+uc1(t)C1-+uR1(t)-R1R1uiR2uouiR+uc2(t)-C2RuouiC+uc1(t)-C-uc2(t)+R2uo(a)(b)
(c)
图2-3 习题2-5 无源网络示意图
解:(a)设电容C上电压为uc(t),由基尔霍夫定律可写出回路方程为
uc(t)ui(t)uo(t)Cduc(t)uc(t)uo(t)
dtR1R2整理得输入输出关系的微分方程为
Cduo(t)du(t)u(t)11()uo(t)Cii dtR1R2dtR1(b)设电容C1、C2上电压为uc1(t),uc2(t),由基尔霍夫定律可写出回路方程为
uc1(t)ui(t)uo(t)ui(t)uc2(t)uo(t)uc2(t)du(t)C2c2
RRdtdu(t)uo(t)uc2(t)RC1c1dt整理得输入输出关系的微分方程为
d2uo(t)duo(t)uo(t)d2ui(t)dui(t)ui(t)RC1C2(2CC)RCC2C 12121dt2dtRdt2dtR(c)设电阻R2上电压为uR2(t),两电容上电压为uc1(t),uc2(t),由基尔霍夫定律可写出回路方程为
uc1(t)ui(t)uR2(t) (1) uc2(t)uo(t)uR2(t) (2)
Cduc1(t)du(t)u(t)Cc2R2 (3) dtdtR2ui(t)uo(t)du(t)Cc2 (4)
R1dt(2)代入(4)并整理得
duR2(t)duo(t)ui(t)uo(t) (5) dtdtR1C(1)、(2)代入(3)并整理得
Cdui(t)du(t)du(t)u(t)Co2CR2R2 dtdtdtR2两端取微分,并将(5)代入,整理得输入输出关系的微分方程为
d2uo(t)duo(t)uo(t)d2ui(t)11dui(t)ui(t) R2C(1)RC2dt2R1CdtR1Cdt2R1CdtR1C
2-6 求图2-4中各无源网络的传递函数。
+Uc(s)-C+Uc1(s)-C1R1+Uc1(s)--Uc2(s)+R1Ui(s)R2Uo(s)RUi(s)+Uc2(s)-C2RUo(s)Ui(s)C+UR2(s)-CR2Uo(s)(a)(b)(c) 图2-4 习题2-6示意图 解:(a)由图得
CsUC(s)UC(s)Uo(s) (1) R1R2UC(s)Ui(s)Uo(s) (2)
(2)代入(1),整理得传递函数为
U(s)Cs1oUR1R1R2CsR2i(s)Cs1
1R1R2CsR1R2R1R2(b)由图得
UC1(s)Ui(s)Uo(s) Ui(s)UC2(s)RUo(s)UC2(s)RC2sUC2(s) RC1sUC1(s)Uo(s)UC2(s)
整理得传递函数为
RC1U1so(s)RC2s2R2C1C2s22RC1s1U)1R2C i(sRC1sRC2s1C2s2Rs(2C1C2)1RC2s2(c)由图得
UC1(s)Ui(s)UR2(s) UC2(s)Uo(s)UR2(s) CsUC1(s)CsUC2(s)UR2(s)R 2Ui(s)Uo(s)RCsUC2(s) 1整理得传递函数为
1)
2)
1) 2)
3) 4)
((((((Uo(s)R1R2C2s2R1Cs 22211Ui(s)CsR1R2Cs(R12R2)Cs1R1R2R1R2CsCs1R2
图2-5 习题2-7 无源网络示意图
2-7 求图2-5中无源网络的传递函数。 解:由图得
U1(s)U2(s)1(Cs)U2(s)
R1R2Ls整理得
1R1U2(s)R2Ls U1(s)Cs11R1CLs2(R1R2CL)sR1R2R1R2Ls
2-8 试简化图2-6中所示系统结构图,并求传递函数C(s)/R(s)和C(s)/N(s)。 解:(a)
⑴求传递函数C(s)/R(s),按下列步骤简化结构图:
图2-6 习题2-8 系统结构图示意图
① 令N(s)0,利用反馈运算简化如图2-8a所示
R(S)-G11G1H1G21G2H2H3图2-8a C(S)
②串联等效如图2-8b所示
R(S)-G1G21G1H11G2H2C(S)H3 图2-8b
③根据反馈运算可得传递函数
G1G2C(s)G1G21G1H11G2H2
G2R(s)1G1(1GH)(1GH)GGH1122123H31G1H11G2H2G1G2
1G1H1G2H2G1H1G2H2G1G2H3⑵求传递函数C(s)/N(s),按下列步骤简化结构图:
①令R(s)0,重画系统结构图如图2-8c所示
H2N(S)H1--H3G1+++C(S)G2 图2-8c
② 将H3输出端的端子前移,并将反馈运算合并如图2-8d所示
N(S)-H1+G1+G21G2H2C(S)H3/H1 图2-9d
③G1和H1串联合并,并将单位比较点前移如图2-8e所示
-1/G1H1N(S)+H3/H1-G1H1G21G2H2C(S) 图2-8e
④串并联合并如图2-8f所示
N(S)11G1H1+G1G2H11G2H2H3/H1图2-8f C(S) ⑤根据反馈和串联运算,得传递函数
G1G2H1C(s)11G2H2(1)
GGHHN(s)G1H1112131G2H2H1G1H11G1G2H1 G1H11G2H2G1G2H3G2G1G2H11G2H2G1G2H3
(b)求传递函数C(s)/R(s),按下列步骤简化结构图: ①将H2的引出端前移如图2-8g所示
H2R(S)-H1G1-G21/G3C(S)-G3H3 图2-8g ②合并反馈、串联如图2-8h所示
H2/G3R(S)-H1 图2-8h G1-G2G31G3H3C(S)
③ 将H1的引出端前移如图2-8i所示
H2/G3R(S)-H1G1-G2G31G3H31G3H3G2G3C(S) 图2-8i ④ 合并反馈及串联如图2-8j所示
R(S)-G1G2G31G2H2G3H3C(S)1G3H3H1G2G3 图2-8j ⑤根据反馈运算得传递函数
G1G2G3C(s)1G2H2G3H3
G1G2G31G3H3R(s)1H11G2H2G3H3G2G3G1G2G3
1G1HG2H2G3H3G1H1G3H32-9 试简化图2-7中所示系统结构图,并求传递函数C(s)/R(s)。
习题2-4 无源网络示意图
图2-7 习题2-9 系统结构图示意图
解:求传递函数C(s)/R(s),按下列步骤简化结构图: ① 将H1的引出端前移如图2-9a所示
H1-G1-G2-1/G4R(S)G3G4C(S)H2H3 图2-9a
② 合并反馈及串联如图2-9b所示
H1/G4-G1-R(S)G2G3G41G3G4H2C(S)H3 图2-9b ③ 合并反馈、串联如图2-9c所示
R(S)-G1G2G3G41G2G3H1G3G4H2C(S)H3 图2-9c
④根据反馈运算,得传递函数
G1G2G3G4C(s)G1G2G3G41G2G3H1G3G4H2
GGGGR(s)11234H31G2G3H1G3G4H2G1G2G3G4H31G2G3H1G3G4H22-10 根据图2-6给出的系统结构图,画出该系统的信号流图,并用梅森公式求
系统传递函数C(s)/R(s)和C(s)/N(s)。
解:(a)根据结构图与信号流图的对应关系,用节点代替结构图中信号线上传递的信号,用标有传递函数的之路代替结构图中的方框,可以绘出系统对应的信号流图。如图2-10a所示。
-H1R(S)1G1N(S)1H2G21C(S)-H3 (1)令N(s)0,求系统传递函数C(s)/R(s)
由信号流图2-10a可见,从源节点R(s)到阱节点C(s)之间,有一条前向通路,其增益为
p1G1G2
有三个相互接触的单独回路,其回路增益分别为
L1G1H1,L2G2H2,L3G1G2H3 L1与L2互不接触
L12G1H1G2H2
流图特征式
1(L1L2L3)L121G1H1G2H2G1G2H3G1G2H1H2
由于前向通路与所有单独回路都接触,所以余因子式
11
根据梅森增益公式,得系统闭环传递函数为
C(s)p11G1G2R(s)1G1H1G2H2G1G2H3G1G2H1H2
(2)令R(s)0,求系统传递函数C(s)/N(s)
由信号流图2-10a可见,从源节点N(s)到阱节点C(s)之间,有两条前向通路,其增益为
p1G2,p2G1G2H1
有两个相互接触的单独回路,其回路增益分别为
L1G2H2,L2G1G2H3
没有互不接触的回路,所以流图特征式为
1(L1L2)1G2H2G1G2H3
由于前向通路与所有单独回路都接触,所以余因子式
11,21
根据梅森增益公式,得系统闭环传递函数为
C(s)12G2G1G2H1pii R(s)i11G2H2G1G2H3(b)根据结构图与信号流图的对应关系,用节点代替结构图中信号线上传递的信号,用标有传递函数的之路代替结构图中的方框,可以绘出系统对应的信号流图。如图2-10b所示。
-H2R(S)1G11G21G31C(S)-H1-H3 图2-10b
求系统传递函数C(s)/R(s)
由信号流图2-10b可见,从源节点R(s)到阱节点C(s)之间,有一条前向通路,其增益为
p1G1G2G3
有三个相互接触的单独回路,其回路增益分别为
L1G1H1,L2G2H2,L3G3H3 L1与L3互不接触
L13G1G3H1H3
流图特征式为
1(L1L2L3)L131G1H1G2H2G3H3G1G3H1H3
由于前向通路与所有单独回路都接触,所以余因子式
11
根据梅森增益公式,得系统闭环传递函数为
C(s)p11G1G2G3 R(s)1G1H1G2H2G3H3G1G3H1H32-11 根据图2-7给出的系统结构图,画出该系统的信号流图,并用梅森公式求
系统传递函数C(s)/R(s)。
解:根据结构图与信号流图的对应关系,用节点代替结构图中信号线上传递的信号,用标有传递函数的之路代替结构图中的方框,可以绘出系统对应的信号流图。如图2-11a所示
-H1R(S)1G1G2G3G41C(S)-H2-H3 图2-11a 由信号流图2-11a可见,从源节点R(s)到阱节点C(s)之间,有一条前向通路,其增益为
p1G1G2G3G4
有三个相互接触的单独回路,其回路增益分别为
L1G2G3H1,L2G3G4H2,L3G1G2G3G4H3
没有互不接触回路。因此,流图特征式
1(L1L2L3)1G2G3H1G3G4H2G1G2G3G4H3
由于前向通路与所有单独回路都接触,所以余因子式
11
根据梅森增益公式,得系统闭环传递函数为
C(s)p11G1G2G3G4 R(s)1G2G3H1G3G4H2G1G2G3G4H3
3-2 已知各系统得脉冲响应,试求系统的闭环传递函数:
(1)k(t)0.0125e1.25t; (2)k(t)5t10sin(4t45); (3)k(t)0.1(1et3)。 解:
(1)(s)L[k(t)]0.0125
s1.2510(sin4tcos4t)] 2(2)(s)L[k(t)]L[5t s52() s2s242s242 5(231422ss1)1616 2ss2(1)16111(3)(s)L[k(t)]0.1[ ]1ss10s(3s1)3
3-3 已知二阶系统的单位阶跃响应为h(t)1012.5e1.2tsin(1.6t53.1),试求系
统的超调量%,峰值时间tp和调节时间ts。 解:
h(t)1012.5e1.2tsin(1.6t53.1) =10[11.25e1.2tsin(1.6t53.1)]
由上式可知,此二阶系统的放大系数是10,但放大系数并不影响系统的动态性能指标。
由于标准的二阶系统单位阶跃响应表达式为
h(t)1112entsin(n12t)
n1.2所以有 1121.25
211.6n0.6解上述方程组,得
n2所以,此系统为欠阻尼二阶系统,其动态性能指标如下 超调量 %e峰值时间 tp12100%e0.61.25100%9.5%
n1220.81.96s
调节时间 ts
3.5n3.52.92 20.63-4 设单位负反馈系统的开环传递函数为G(s)跃输入下的动态性能。 解题过程:
由题意可得系统得闭环传递函数为
0.4s1,试求系统在单位阶
s(s0.6)2G(s)0.4s1nsa (s)2221G(s)ss1as2dnsn其中a2,n1,dn2z0.5,z2.5。这是一个比例-微分控制二阶系统。
比例-微分控制二阶系统的单位阶跃响应为 h(t)1redntsin(n1d2t) 故显然有 r2z22dnnz1d22 3n1d21d2)arctan1.686 arctan(zdnd darctan1d2d31.047
此系统得动态性能指标为 峰值时间 tpdn12d2d3.155
dtp21n超调量 %r1e16.2%
1123ln(z22nnn)lnzln(1d2)225.134 调节时间 tsdn
3-5 已知控制系统的单位阶跃响应为h(t)10.2e60t1.2e10t,试确定系统的阻
尼比和自然频率n。 解:
系统的单位脉冲响应为k(t)h(t)12e60t12e10t12(e10te60t) 系统的闭环传递函数为(s)L[k(t)]12(自然频率 n60024.5 阻尼比
3-6 已知系统特征方程为3s410s35s2s20,试用劳斯稳定判据和赫尔维
茨稳定判据确定系统的稳定性。 解:
先用劳斯稳定判据来判定系统的稳定性,列出劳斯表如下
s4 3 5 2s3 10 147s2 2
10153s1 47s0 2701.429
2600•11600 )2s10s60s10s600显然,由于表中第一列元素得符号有两次改变,所以该系统在s右半平面有两个闭环极点。因此,该系统不稳定。
再用赫尔维茨稳定判据来判定系统的稳定性。显然,特征方程的各项系数均为正,则
2a1a2a0a310531470
a12a41022 2002
a31显然,此系统不稳定。
3-7 设单位负反馈系统的开环传递函数为G(s)K,试应2(s2)(s4)(s6s25)用劳斯稳定判据确定义为多大值时,特使系统振荡,并求出振荡频率。 解:
由题得,特征方程是s412s369s2198s200K0 列劳斯表
s4 1 69 200+K s3 12 198 s2 52.5 200+K s1 7995-12K s0 200+K
由题意,令s1所在行为零得K666.25 由s2行得 52.5s2200666.250
解之得 s4.062i,所以振荡角频率为 4.062rad/s
3-8 已知单位负反馈系统的开环传递函数为G(s)系统稳定时的K值范围。 解:
由题可知系统的特征方程为
D(s)s43s34s2(2K)s2K0
K(0.5s1),试确定2s(s1)(0.5ss1)列劳斯表如下
s4 1 4 s3 3 2+K10-K 2K 3(10-K)(2+K)6K13s 10-K3s0 2K s2 由劳斯稳定判据可得
10K30[(10K)(2K)/3]6K0 (10K)/32K0
解上述方程组可得 0K1.705
3-9系统结构如图3-1所示,G(s)K,定义误差e(t)r(t)c(t),
s(Ts1)(1) 若希望图a中,系统所有的特征根位于s平面上s2的左侧,且阻尼
比为0.5,求满足条件的K,T的取值范围。 (2) 求图a系统的单位斜坡输入下的稳态误差。
(3) 为了使稳态误差为零,让斜坡输入先通过一个比例微分环节,如图b所
示,试求出合适的K0值。
(a) (b)
解:(1)闭环传递函数为(s)KTs2sKK/T
1Ks2sTT即nK111,2n,0.5n,K TTTTD(s)Ts2sK,令s's2,代入上式得,
D'(s)T(s'2)2s'2KTs'2(4T1)s'4T1/T20
列出劳斯表,
s2 T 4T+1T2 s1 1-4T s0 4T+1T2
T0,14T0,4T1/T200T1/4 或T0,14T0,4T1/T20无解
0T1/4,4K
(2) R(t)t,系统为I型系统 ∴ess1/K (3) G'(s)(K0s1)KKsKK20
s(Ts1)KTssKTs1KK01Ts2(1KK0)s E(s)R(s)C(s)R(s)[1G'(s)]222sTssKs(TssK)令esslimsE(s)lims0Ts1KK01KK00K01/K
s0Ts2sKKK0并没有改变系统的稳定性。
3-10 已知单位反馈系统的开环传递函数: (1)G(s)100;
(0.1s1)(s5)50
s(0.1s1)(s5)(2)G(s)试求输入分别为r(t)2t和r(t)22tt2时,系统的稳态误差。 解: (1)G(s)10020
(0.1s1)(s5)(0.1s1)(0.2s1)由上式可知,该系统是0型系统,且K20。
11,,。根据线性0型系统在1(t),t,t2信号作用下的稳态误差分别为:
21K叠加原理有该系统在输入为r(t)2t时的稳态误差为ess22•,该系统在输入为r(t)22tt2时的稳态误差为ess22• (2) G(s)12• 1K5010
s(0.1s1)(s5)s(0.1s1)(0.2s1)由上式可知,该系统是型系统,且K10。
11根据线性叠加型系统在1(t),t,t2信号作用下的稳态误差分别为:0,,。K21原理有该系统在输入为r(t)2t时的稳态误差为ess22•0.2,该系统在输入
K1为r(t)22tt2时的稳态误差为ess2202
K
3-11已知闭环传递函数的一般形式为
bmsmbm1sm1b1sb0G(s) (s)nn11G(s)H(s)san1sa1sa0误差定义为e(t)r(t)c(t)。试证,
(1) 系统在阶跃信号输入下,稳态误差为零的充分条件为
(s)a0
snan1sn1a1sa0(2)系统在斜坡信号输入下,稳态误差为零的充分条件为
(s)a1sa0 nn1san1sa1sa0(3)推导系统在斜坡信号输入下稳态误差为零的充分条件 (4)求出系统闭环传递函数与系统型别之间的关系 解:(1)(s)a0 nn1san1sa1sa0 E(s)R(s)C(s)R(s)[1(s)]
snan1sn1a1s1 •n
ssan1sn1a1sa0sn1an1sn2a1 n n1san1sa1sa0 满足终值定理的条件,
snan1sn1a1s e()limsE(s)limn0 n1s0s0saa1sa0n1s 即证 (2)(s)a1sa0 nn1san1sa1sa0 E(s)R(s)C(s)R(s)[1(s)]
snan1sn1a2s21 2•n n1ssan1sa1sa0sn1an1sn2a2 n n1san1sa1sa0 满足终值定理的条件,
snasn1as e()lims0sE(s)limn12s0snan1n1sa0 1sa0 即证
(3) 对于加速度输入,稳态误差为零的必要条件为
(s)a2s2a1sa0snan1 n1sa1sa0同理可证
(4)系统型别比闭环函数分子最高次幂大1次。 3-12 已知单位反馈系统的开环传递函数为:
(1)G(s)50(0.1s1)(2s1);
(2)G(s)Ks(s24s200); (3)G(s)10(2s1)(4s1)s2(s22s10)
试求位置误差系数Kp,速度误差系数Kv,加速度误差系数Ka。解:
(1) 此系统是一个0型系统,且K20。故查表可得KpK10,Ka0
(2) 根据误差系数的定义式可得
KKplims0G(s)H(s)lims0s(s24s200) KKvlims0s•G(s)H(s)lims0s•s(s24s200)K200
Klims0s2G(s)H(s)limKas0s2s(s24s200)0(3) 根据误差系数的定义式可得
Kv0,
10(2s1)(4s1)s0s0s2(s22s10)10(2s1)(4s1)Kvlims•G(s)H(s)lims•22
s0s0s(s2s10)10(2s1)(4s1)Kalims2G(s)H(s)lims2221s0s0s(s2s10)KplimG(s)H(s)lim
3-13设单位反馈系统的开环传递函数 G(s)K0•Km1•
Tfs1s(Tms1)i•Kf 输入信号为 r(t)(abt)•1(t)
其中K0, Km, Kf, i, Tf, Tm均为正数,a和b为已知正常数。如果要求闭环系统的稳态误差ess<0, 其中00, 试求系统各参数满足的条件。 解:首先系统必须是稳定的,系统的闭环特征方程为
TfTms3(TfTm)s2sK0
式中,KK0KfKm/i,为系统的开环增益,各参数满足: K0, (TfTm)KTmTf0 即稳定条件为 0KTfTm TfTm由于本例是I型系统,其Kp, KvK,故在r(t)(abt)•1(t)作用下,其稳态误差
essbb 0 必有 K0K于是,即能保证系统稳定,又满足对系统稳态误差要求的各参数之
间的条件为
b0K0KfKm/iTfTmTfTm
3-14 设单位反馈系统的开环传递函数为G(s)1Ts。试用动态误差系数法求出
当输入信号分别为r(t)t22时,系统的稳态误差。 解:
系统的误差传递函数为
e(s)E(s)1TsTs(Ts)2(Ts)3(Ts)4R(s)1G(s)1Ts
所以有E(s)e(s)•R(s)Ts•R(s)(Ts)2•R(s)(Ts)3•R(s)(Ts)4•R(s)对上式进行拉氏反变换可得
e(t)Tr(t)Tr(t)Tr(t)T4r(4)(t)•2••3••• (1)
当r(t)t22时,显然有
r(t)t r(t)1•••
r(t)r(4)(t)•••0T(tT)
将上述三式代入(1)式,可得e(t)T•tT2•1T3•0T4•0系统的稳态误差为esslime(t)limT(tT)
tt
3-15 假设可用传送函数
C(s)1描述温度计的特性,现在用温度计测量盛R(s)Ts1在容器内的水温,需要一分钟时间才能指出实际水温的98%的数值。如果给容器加热,使水温依10/min的速度线性变化,问温度计的稳态误差有多大? 解:
由题意,该一阶系统得调整时间ts1min,但ts4T,所以T0.25min。
系统输入为r(t)10t,可推得R(s)因此可得 C(s)10 s2110R(s)2 Ts1s(Ts1) c(t)10t10T10TetT
c(t)的稳态分量为css(t)10t10T
稳态误差为 ess(t)r(t)css(t)10T100.252.5 所以,稳态误差为2.5C
3-16如图3-2所示的控制系统结构图,误差E(s)在输入端定义,扰动输入
n(t)21(t).
(1) 试求K40时,系统在扰动输入下的稳态输出和稳态误差。 (2) 若K20, 其结果又如何?
1(3) 在扰动作用点之前的前向通道中引入积分环节,对其结果有何影响?
s1在扰动作用点之后的前向通道中引入积分环节,对其结果又有何影响?
s
图3-2 习题 3-16 示意图
解:令G1K1 ,G2,H2.5
0.05s1s5则C(s)G2N(s)G1G2E(s) 代入E(s)R(s)HC(s)
得 C(s)G2G1G2N(s)R(s)
1G1G2H1G1G2H令R(s)0,得扰动作用下的输出表达式: Cn(s)G2N(s)
1G1G2HG2HN(s)
1G1G2H此时的误差表达式为:En(s)R(s)HCn(s)若在s 右半平面上解析,则有 essnlimsEn(s)lims0G2HsN(s)
s01GGH12在扰动输入下的稳态输出为 Cn()limsCn(s)lims0G2sN(s)
s01GGH12代入N(s),G1,G2,H的表达式,可得
15 ,essn12.5K12.5K25(1) 当K40时,cn() ,essn10110125(2) 当K20时,cn(),essn
5151 cn() 可见,开环增益的减小将导致扰动作用下系统稳态输出的增大,且稳态误差的绝对值也增大。
(3) 若1s加在扰动之前,则 G1K1 G2 H2.5
s(0.05s1)s5 得 cn()0,essn0 若1s加在扰动之后,则 G11K G2 H2.5
(s5)0.05s1 cn() essn20.02(K40),0.04(K20) 2.5K50.05(K40),0.1(K20)
2.5K 可见在扰动作用点之前的前向通路中加入积分环节,可以消除阶跃输入引起的稳态误差。
3-17 设随动系统的微分方程为:
d3c(t)d2c(t)dc(t)TmTaTmKc(t)Kr(t) 32dtdtdt 其中,c(t)为系统输出量,r(t)为系统输入量,Tm为电动机机电时间常数,
Ta为电动机电磁时间常数,K为系统开环增益。初始条件全部为零,试讨论:
(1) Ta、Tm与K之间关系对系统稳定性的影响
(2) 当Ta0.01, Tm0.1, K500时,可否忽略Ta的影响?在什么影响下Ta的影响可以忽略?
解:(1)对系统微分方程在零初始条件下进行拉氏变换,得闭环系统特征方程
TmTas3Tms2sK0
当 TmTa K均为正值时,且有 D2Tm(1TaK)0
即 0K1Ta 时 闭环系统稳定。
(2)由于Ta0.01,因此只有当0K100
闭环系统才稳定,显然,对于K500, 闭环不稳定。此时若略去Ta, 闭环特征方程为
Tms2sK0.1s2s5000
上式中各项系数为正,从而得到得出闭环系统稳定的错误结论。如果 K100。如果K100,则略去Ta不会影响闭环稳定性。
对于本例,当K1Ta时,不能忽略Ta对稳定性的影响,否则可以忽略。 3-18 设计题
飞机的自动控制,是一个需要多变量反馈方式的例子。在该系统中,飞机的飞行姿态由三组翼面决定,分别是:升降舵,方向舵和副翼,如附图3-3(a)所示。飞行员通过操纵这三组翼面,可以使飞机按照既定的路线飞行。
这里所要讨论的自动驾驶仪是一个自动控制系统,它通过调节副翼表面来控
图3-3(a) 飞机副翼模型图
制倾角,只要使副翼表面产生一个的变形,气压在这些表面上会产生一个扭
1矩,使飞机产生侧滚。飞机副翼是由液压操纵杆来控制的,后者的传递函数为。
s测量实际的倾角,并与输入设定值进行比较,其差值被用来驱动液压操纵杆,而液压操纵杆则反过来又会引起副翼表面产生变形。
图3-3(b) 飞机控制倾角结构图
为简单化起见,这里假定飞机的侧滚运动与其他运动无关,其结构图如图3-3(b)所示,又假定K11,且角速率由速率陀螺将其值进行反馈,期望的阶跃响应的超调量%10%,调节时间(以2%的标准)ts9s,试选择合适的Ka和K2值。 解:
由于过阻尼响应缓慢,故通常不希望采用过阻尼系统,在本题中0,1欠阻尼
G(s)1KaKa2
ss1K2Kas(K2Ka1)sKa 2s(K2Ka1)sKa(s)因此,
2nKa 2nK2Ka1计算可得
nKaK2Ka1 2Ka又因,%e12100%,ts3.5n9
由题计算可得0.59,n0.659
故Ka0.4,K20.634
图4-1 习题4-1系统零极点分布图
4-1 已知系统开环零极点分布如图4-1所示,试绘制相应的根轨迹图。 解:
j0j0(b)j0(d)(e)图4-1a根轨迹图 j0(c)(a)j0(f)j0 (a)根轨迹的渐近线条数为nm0 (b)根轨迹的渐近线条数为nm0
(c)根轨迹的渐近线条数为nm3,渐近线的倾斜角为160,2180,
3240
(d)根轨迹的渐近线条数为nm0 (e)根轨迹的渐近线条数为nm0
(f)根轨迹的渐近线条数为nm1,渐近线的倾斜角为180
4-2 已知单位反馈控制系统的前向通道传递函数为:
(1) G(s)K(s1)KG(s) (2)
s(s1)(s2)(s5)s2(s2)(s4)KK(s1)G(s) (4) 22s(s4)(s4s20)s(s1)(s4s16)(3) G(s)K0,画出各系统的根轨迹图。 解:(1)按下列步骤绘制根轨迹:
① 系统开环有限零点为z11;开环有限极点为p1,20,p32,p44 ②实轴上的根轨迹区间为,4,2,1
③根轨迹的渐近线条数为nm3,渐近线的倾角为
160,2180,360
渐近线与实轴的交点为
apzii1i1nminm5
3闭环系统根轨迹如下图4-2a所示
j-4-2-10 图4-2a 闭环系统根轨迹图 (2)按下列步骤绘制根轨迹:
①系统没有开环有限零点;开环有限极点为p10,p21,p32,p45 ②实轴上的根轨迹区间为5,2,1,0
③根轨迹的渐近线条数为nm4,渐近线的倾角为
145,2135,3135,445
渐近线与实轴的交点为
a④分离点方程为
pzii1i1nminm2
11110 dd1d2d5解得分离点d14.06,d20.40 闭环系统根轨迹如下图4-2b所示
j-5-2-10 图4-2b
(3)按下列步骤绘制根轨迹:
①系统没有开环有限零点;开环有限极点为p10,p24,p3,42j4 ② 实轴上根轨迹区间为4,0
③ 根轨迹的渐近线条数为nm4,a2,a45,135,225,315 ④根轨迹的起始角:复数开环有限极点p3,42j4处,p390,p490 ⑤分离点方程为
11110 dd4d2j4d2j4解得分离点d12,d2,32j6
检查
d12时,K*
d2,32j6时,K*100
d1,d2,d3皆为闭环系统根轨迹的分离点。
⑥确定根轨迹与虚轴的交点:系统闭环特征方程为
D(s)s48s336s280sK*0
列写劳斯表
s4 1 36 K*s3 8 80 s2 26 K* 80268K*s 26s0 0 1当K*260时,劳斯表出现全零行,辅助方程为
A(s)26s22600
解得根轨迹与虚轴交点为10。 根轨迹如下图4-2c所示:
j42-6-4-20-2-4 图4-2c
(4)按下列步骤绘制根轨迹:
①系统开环有限零点为z11;开环有限极点为p10,p21,
p3,4-2j23
②实轴上根轨迹区间为(,1],0,1
2③根轨迹的渐近线条数为nm3,a,a60,180,60
3④分离点方程为
11111 dd1d2j23d2j23d1解得分离点d12.26,d20.45
根轨迹如下图4-2d所示:
j-101 图4-2d
4-3 给定系统如图4-2所示,K0,试画出系统的根轨迹,并分析增益对系统
图4-2 习题4-3系统零极点分布图
阻尼特性的影响。
解:(1)作系统的根轨迹。开环传递函数为
G(s)F(s)K(s2)(s3)
s(s1)①开环极点为0和1,开环零点为2和3。 ②所以实轴上的根轨迹区间为[3,2]和[1,0]。 ③分离点方程
1111 dd1d2d3得分离点d12.366,d20.634 检查
d12.366时,K*s(s1)(s2)(s3)s(s1)(s2)(s3)s2.3660.0718
d20.634时,K*s0.63413.93
可得到根轨迹如下图4-3a所示
j-3-2-10 图4-3a
(2)分析增益对阻尼特性的影响。
从根轨迹图可以看出,对于任意K0,闭环系统都是稳定的,但阻尼状况不同。
增益较小时(0K0.0718)系统过阻尼; 增益很大时(K13.93),系统过阻尼; 增益中等时(0.0718K13.93),系统欠阻尼。
4-4 给定控制系统如图4-3所示, K0,试用系统的根轨迹图确定,速度反
馈增益K为何值时能使闭环系统极点阻尼比等于0.7。
解:(1)求系统的闭环特征方程并划成标准形式。通过方块图变换或代数运算可以求得单位反馈系统的开环传递函数
G(s)10/(s1)110
110k/(s1)ss(s110k)因为可变参数K不是分子多项式的相乘因子,所以先求系统的闭环特征方程
s(s110k)10s2s10ks100
改写为
110ks0 2ss10即,上述闭环特征方程也相当于开环传递函数为
G'(s)的系统的闭环特征方程。
Ks0,K10k
s2s10(2)根据G'(s)作出根轨迹图。
G'(s)有两个极点0.5j3.1225,一个零点0,所以负实轴是根轨迹,而且
其上有分离点。将闭环特征方程改写为
s2s10K
s由dK/ds0可以求得s10,其中s10在根轨迹上,对应增益为
K5.32460,故s10是实轴上的分离点。根轨迹如图4-4a所示。
jζ=0.70 图4-4a (3)求反馈增益k。首先要确定闭环极点。设途中虚线代表0.7,则闭
环极点为根轨迹和该虚线的交点,由0.7可得arccos45.57。设
s1njn120.7nj0.51n
列出该点对应的辐角条件
argG'(s)args
(s10.5j3.1225)(s10.5j3.1225)arctan0.51n3.12250.51n3.12250.51arctanarctan0.70.7n0.50.7n0.5
180(2k1)经整理得
arctan0.51n3.12250.51n3.12250.51 arctan180(2k1)arctan0.7n0.50.7n0.50.7两边同取正切,整理得
21.0202n10.20200
解得,n3.1623。所以该闭环极点为s12.2136j2.2583。再由
s2s10Ks3.4272
s-2.2136j2.2583得速度反馈增益为kK/100.3427。
4-5 已知单位反馈系统的开环传递函数为:G(s)K。要求系统的
s(s1)(0.5s1)闭环极点有一对共轭复数极点,其阻尼比为0.5。试确定开环增益K,并近似分析系统的时域性能。
解:根据绘制常规根轨迹的基本法则,作系统的概略根轨迹如图4-5a所示。
0.5 j21-2-10-1-2 图4-5a
欲确定K,需先确定共轭复极点。设复极点为s1,2xjy 根据阻尼比的要求,应保证
yxtg(180arccos)1.732x
在图上作0.5的阻尼线,并得到初始试探点的横坐标x0.3,由此求得纵坐标y0.52。在s0.3j0.52处检查相角条件
G(s)-173.6
不满足相角条件;修正x0.32,则y0.5,点s0.32j0.5处的相角为
177.4;再取x0.33,则y0.572,点s0.33j0.572处的相角为180。
因此共轭复极点s1,20.33j0.572。由模值条件求得
K1*K2s0.33j0.5720.513
运用综合除法求得另一闭环极点为s32.34。共轭复极点的实部与实极点的实部之比为0.14,因此可视共轭复极点为系统的主导极点,系统的闭环传递函数可近似表示为
(s)0.436
s20.665s0.436并可近似地用典型二阶系统估算系统的时域性能
ts3.5n10.6
%e/12100%16.3%
4-6 已知单位反馈系统的开环传递函数为:
K(s22s4)G(s),K0 2s(s4)(s6)(s1.4s1)试画出系统的根轨迹图,并分析系统的稳定时K的取值范围。 解:由题得
开环极点:0,4,6和0.7j0.714 开环零点:1j1.7321
分离、会合点:从s平面的零点、极点分布可知在区间内(4,0)可能有分离、会合点。 记
A(s)s(s4)(s6)(s21.4s1)s511.4s439s343.6s224sB(s)s2s4由A'(s)B(s)A(s)B'(s),可得
2
(5s445.6s3117s287.2s24)(s22s4)(s11.4s39s43.6s24s)(2s2)经整理后得到
32
3s630.8s5127.4s4338.4s3531.2s2348.8s960
用试探法或程序算得区间(4,0)内的一个根为2.3557,它就是实轴上的分离点。 根轨迹自复数极点的出射角:.88 根轨迹趋向复数零点的入射角:102.52 根轨迹与虚轴的交点:闭环特征方程为
f(s)s511.4s439s3(43.6k)s2(242K)s4K0
4211.4(43.6K)4K0令sj,可得 43[39(242K)]0由第二式得K0.5419.5212,代入第一式,得 620.2492.92960 解得11.2115,22.15,33.7537 根据以上数据画根轨迹图,如图4-6a所示。
图4-6a 根轨迹图
再分析系统得稳定情况:根轨迹与虚轴第一个交点的频率为 11.2115, 利用幅值条件可以计算出对应的增益
s(s4)(s6)(s21.4s1)K|sj1.211515.56
s22s4同样可以算得与22.15和33.7537对应的增益K2.74,K3163.43 参看根轨迹图可知:系统稳定时K的取值范围为:K15.或.74K163.43 4-7 已知单位反馈系统的开环传递函数为:
G(s)K
s(s1)(s2)K的变化范围是0,试画出系统的根轨迹图。 解:按下列步骤绘制根轨迹:
①系统没有开环有限零点;开环有限极点为p10,p21,p32 ②实轴上的根轨迹区间为,2,1,0
③根轨迹的渐近线条数为nm3,渐近线的倾角为
160,2180,360
渐近线与实轴的交点为 a④分离点方程为 解得分离点d0.42
pzii1i1nminm1
1110 dd1d2闭环系统根轨迹如下图4-7a所示
j21-2-10-1-2 图4-7a
4-8 已知反馈控制系统的开环传递函数为:
G(s)K,H(s)1,K0,a0
s(sa)试画出K和a同时变化的根轨迹簇。
解:(1)列写闭环特征方程。闭环特征方程为s2asK0
(2)画a0,K从0到的根轨迹。a0时闭环特征方程为s2K0。
这相当于一个开环传递函数为
G1(s)H1(s)K 2s的系统。它的根轨迹是与虚轴重合的直线。见图4-8a中由圆圈构成的根轨迹。 (3)画K为常数,a从0到的根轨迹。给定K,则闭环特征方程为
as0 2sKas它相当于一个开环传递函数为G2(s)H2(s)2的系统,该系统的开环极点为
sK1jK,开环零点为0。图4-8a中不带圆圈的根轨迹是K1,4,9,16时的根轨迹。
j43210-1-2-3-4 1(sa)44-9 已知单位反馈系统的开环传递函数为:G(s)2 s(s1)a的变化范围是[0, ],试画出系统的闭环根轨迹。
解:系统闭环特征方程为
11 D(s)s3s2sa0
44即有11s3s2s41a40
K*等效开环传递函数为G1(s)
12s(s)2K*1a,变化范围为0, 4按照绘制常规根轨迹的基本法则确定根轨迹的各项参数:
1(1)等效系统无开环有限零点;开环有限极点为:p10,p2p3
2(2)实轴上的根轨迹区间为,0
1(3)根轨迹有3条渐近线,且a,a60,180,300
31K*(3s22s)d40 (4)根轨迹的分离点:由分离点方程G1(s)1dss2(s)4211 解得d1,d2
26(5)根轨迹与虚轴的交点:根据闭环特征方程列写劳斯表如下:
1s3 1 4a s2 1
41as1 - 44当a1时,劳斯表的s1行元素全为零,辅助方程为A(s)s2解得s1,2j10 41 2绘制系统参数根轨迹如图4-9a所示
j121216012 图4-9a 4-10 已知反馈控制系统中,其开环传递函数为:
K(s22s4)G(s) 2s(s4)(s6)(s1.4s1)s4时的闭环根轨迹概略图; ss1.05(2) 绘制H(s)时的闭环根轨迹概略图;
s(1) 绘制H(s)(3) 比较开环零点变化对根轨迹形状的影响。 解:(1)开环传递函数
K(s22s4)G(s)H(s)2
s(s6)(s21.4s1)按下列步骤绘制根轨迹:
①系统开环有限零点为z1,21j1.732;开环有限极点为p1,20,p36,
p4,50.7j0.714
②实轴上的根轨迹区间为,6
③根轨迹的渐近线条数为nm3,渐近线的倾角为
160,2180,360
渐近线与实轴的交点为
apzii1i1nminm1.8
闭环系统根轨迹如下图4-10a所示
j-60 图4-10a根轨迹图
(2)开环传递函数
K(s1.05)(s22s4)G(s)H(s)2
s(s4)(s6)(s21.4s1)按下列步骤绘制根轨迹:
①系统开环有限零点为z11.05,z2,31j1.732;开环有限极点为
p1,20,p34,p46,p5,60.7j0.714
②实轴上的根轨迹区间为,6 和4,1.05 ③根轨迹的渐近线条数为nm3,渐近线的倾角为
160,2180,360
渐近线与实轴的交点为
apzii1i1nminm2.78
闭环系统根轨迹如下图4-10b所示
j-6-40 图4-10b根轨迹图 4-11 给定控制系统的开环传递函数为:
G(s)sa,a0
s(2sa)试作出以a为参变量的根轨迹,并利用根轨迹分a取何值时闭环系统稳定。 解:(1)求系统的闭环特征方程并化成标准的形式。因为可变参数a不是分子多项式的相乘因子,所以先求系统的闭环特征方程
2s2assa0
可改写为
1a(s1)0
s(2s1)则开环传递函数为
G'(s)a(s1)K(s1),Ka0
s(2s1)s(2s1)(2)根据G'(s)作系统的根轨迹。G'(s)中的增益为负值,所以要作系统的补根轨迹。开环极点为0.5和0,开环零点为1。按照补根轨迹的作图规则,实轴上的根轨迹区间为0.5,0和1,。在0.5,0 区间有会合点,在1,有分离点。为求分离、会合点,将闭环特征方程改写为
Ks(2s1)
(s1)由dK/ds0,得s24s10,解得s12.2247,s20.2247,分别对应的增益为K9.90和K0.1010,所以是分离、会合点。可以证明,不在
ja1(K1)a0a01a(K)12-10-1a(K) 图4-11a 实轴上的根轨迹是一个圆,圆心在1,0,半径为1.2227。以Ka为参变量的根轨迹如图4-11a所示,图中箭头表示a从0到的方向,也即K从0到
的方向。
(3)求a使闭环系统稳定的取值范围。首先求根轨迹与虚轴的交点。由闭环特征方程
2s2(1K)sK0
可知,K1时系统处于临界稳定状态,这相当于a1,所以使闭环系统稳定的范围为0a1。
4-12 实系参数多项式函数为:
A(s)s35s2(6a)sa
欲使A(s)0的根均为实数,试确定参数a的范围。 解:对A(s)0作等效变换得
1等效开环函数为
a(s1)0 32s5s6sG1(s)H1(s)a(s1)
s(s2)(s3)当a0时,需绘制常规根轨迹:
系统开环有限零点为z11;开环有限极点为p10,p22,p33 实轴上的根轨迹区间为3,2和1,0 根轨迹有2条渐近线,且
a2;a90,90
由分离点方程
1111 d1dd2d3在实轴区间3,2内用试探法求得d2.47。绘制根轨迹图,如图4-12a所示。
当a0时,需绘制零度根轨迹。实轴上,零度根轨迹区间为(-∞,-3],[-2,-1]和[0,+∞]。作零度根轨迹图,如图4-12b所示。
当多项式有根2.47时,根据模值条件得
add2d30.419
d1根据常规根轨迹图,知当0a0.419时,多项式的根皆为实数;根据零度
根轨迹图,知当a0时,多项式的根亦全为实数。因此所求参数a的范围为
a0.419。
21-3-20-1-221-3-20-1-2-1-1 图4-12a 常规根轨迹 图4-12b 零度根轨迹
4-13 设系统开环传递函数为:
G(s)K
s(s1)(s10)(1) 大致画出系统的根轨迹图;
(2) 用文字说明当K0时,如何求系统单位阶跃响应的超调量%,峰值时间tp及调节时间ts。 解:(1)绘根轨迹图
1103.67;a60,180 3111分离点:由0,得d0.487
dd1d10渐近线:a相应的根轨迹增益Kd2.377
根轨迹与虚轴交点:闭环特征方程s311s210sK0 列劳斯表
s3 1 10 s2 11 K 110-Ks 0 10s0 K 1
当K110时,劳斯表出现全零行,由辅助方程11s21100 得根轨迹与虚轴交点处为K110,3.16 根轨迹图如下图4-13a所示:
j321-100-1-2-3-1 图4-13a
(2)求动态性能指标
当0K2.377时,系统%0,tp0,闭环有两个实主导极点1和2,且
12,因此求得调节时间如下:
ts4.751, 12 13ln(1)2ts, 121当2.377K100时,闭环系统有一对共轭复极点,则
A3lnt%D tp,%100e1p,tsD1A1由于Ddn12,An,1n,ln()ln()
2D1因此
tp1A, %100e12%n12, ts=3.5
, 0.8 n
4-14 设单位负反馈系统的开环传递函数为:
G(s)K(s4)
s(s2)试画出系统根轨迹图,并求出系统具有最小阻尼比时的闭环极点和对应的增益K。
解:系统在实轴上的根轨迹区域为0,2和4, 在这两段区域内,均存在分离点。为了求出分离点,令
111 dd2d4求出
d14221.172d24226.828
因而复数根轨迹是以(4,j0)为圆心,22为半径的一个圆,如图4-14a所示
ζ=0.707jβ-7-6-5-4-3-2-10 图4-14a
在图上,过原点作圆得切线,得最小阻尼比线。由根轨迹图知,对于等腰直角三角形,必有
45,故最小阻尼比cos0.707 响应的闭环极点s1,22j2
由根轨迹模值条件,可求出相应的增益为
K
2j2•2j222
2j244-15 已知单位负反馈系统的开环传递函数为:
K* G(s)2s(s4)(s4s5)试按照步骤作出K*0时的根轨迹图。 解:开环极点:p10,p24,p3,42j 根轨迹在实轴上的区间0,4
根轨迹的渐近线nm4,a2,a45,135 分离点:
11110 dd4d2jd2j 即
112(d2)20 dd4d4d5整理得 d36d210.5d50
为了求取分离点方程的根,将上式表示为110.5(d0.476)0
d2(d6)K1*(d0.476)令等效开环传递函数为G1(d)
d2(d6)其中K1*10.5。若令K1*从0变到,其根轨迹如图4-15a所示。图中,渐近线
j2d1-2-10-1-21d2-6-5-4-3 图4-15a a2.78,a90;分离点d1'1.088,d2'2.627。 在图上,试探d2,检验模值条件 K1*22410.5
1.524故符合要求,故d2为分离点方程的一个根。利用综合除法,有
d36d210.5d5(d2)(d24d2.5)(d2)(d0.775)(d3.225)0
求得分离点d10.775,d22,d33.225 分离角为90 根轨迹的起始角
p180(311(pp))180arctanarctan1809090 3j22j1(j3)4p90
4根轨迹与虚轴的交点:闭环特征方程为s48s321s220sK*0 列劳斯表
s4 1 21 K* s3 8 20 s2 18.5 K* 370-8K*s 0 18.5s0 K* 1
显然,当K*46.25时,根轨迹和虚轴相交,由辅助方程
18.5s246.250
求得交点处K*46.25,1.58
根据以上步骤,绘制系统根轨迹图4-15b
图4-15b 根轨迹图
4-16 设某单位负反馈系统的开环传递函数为:
G(s)4K(1s)
s[(K1)s4](1) 绘制K从0时系统的根轨迹图;
(2) 求系统阶跃响应中含有eatcos(t)时的K值范围,其中a0, 0; (3) 求系统有一个闭环极点为2时的闭环传递函数。 解:绘制根轨迹图 闭环特征方程为
(s24s)K(s24s4)0 写成根轨迹方程形式为:
K(s2)20 1s(s4)令等效开环传递函数为
K(s2)2G1(s)
s(s4)实轴上根轨迹:4,0 分离点:由
112求得d1 dd4d2与虚轴交点:列劳斯表
s2 K+1 4K s1 4-4K 0 s0 4K 显然,当时系统处于临界稳定,由辅助方程并代入K1,解出交点处
K1,2 分离点处根轨迹增益:由模值条件得:Kd绘出系统根轨迹如图4-16a所示
j1131 333K1,2-4-3-2-1-1012 图4-16a (2)求K值范围
当系统阶跃响应含有etcos(t)分量时,系统处于欠阻尼状态,系统有一对
1具有负实部的共轭极点,K值范围为K1
3(3)求闭环传递函数
当系统具有s12闭环极点时,由模值条件,其对应的K值为K于是 G(s)221 444(1s)
5s(s4)4G(s)0.8(1s)
1G(s)(s0.4)(s2)闭环传递函数为(s)
5-1 设系统闭环稳定,闭环传递函数为(s),试根据频率特性的定义证明,输
入为余弦函数r(t)Acos(t)时,系统的稳态输出为
css(t)A|(j)|cos[t(j)]
解:
由题目可得r(t)Acos(t)AcostcosAsintsin 对等式两边同时进行拉氏变换可得
R(s)Ascossinscossin AAs22s22s22由于系统闭环稳定,所以(s)不存在正实部的极点。假设(s)可表示为如下表达式
(s)M(s)
(ss1)(ss2)(ssn)由以上分析可得,系统的闭环传递函数为
C(s)(s)R(s)M(s)scossinA
(ss1)(ss2)(ssn)s22对上述闭环传递函数作如下分解
C(s)i1nDiB1B2 ssisjsj对上式等式两边进行拉氏反变换可得
c(t)DiesitB1ejtB2ejt
i1n由系统稳态输出的定义可得
css(t)limc(t)B1ejtB2ejt
t利用留数法确定待定的系数B1,B2
B1limA(s)sjscossincossinA(j)ej(j)()
sj22jB2limA(s)sjscossincossinA(j)ej(j)()
sj22j所以可得
css(t)A(j)(cosjt(j)sinjt(j)ee 22j cosjt(j)sinjt(j)ee) 22j A(j)coscost(j)sinsint(j) A(j)cost(j)
5-2 若系统阶跃响应为:
h(t)11.8e4t0.8e9t
试确定系统频率特性 解:
单位阶跃输入信号的拉氏变换为R(s)系统单位阶跃响应的拉氏变换为
11.80.836 H(s)
ss4s9s(s4)(s9)1 s系统的闭环传递函数为G(s)H(s)36 R(s)(s4)(s9)将sj代入传递函数G(s)可得 G(j)
5-3 设系统结构图如图5-1所示,试确定输入信号
36
(j4)(j9)r(t)sin(t300)cos(2t450)
作用下,系统的稳态误差ess(t)。
图5-1 习题5-3控制系统结构图
解:
如图5-1所示,系统的误差传递函数为
11(s)1ss1e(s)111s2 s11: 备注:为什么稳态误差? 其幅频特性和相频特性分别为
(j)2124,e(j)arctanarctane2 当r(t)sin(t300)cos(2t450)sin(t300)sin(2t450)时
e1ss(t)114sin(t30arctan1arctan12)4144sin(2t45arctan2arctan1) 105sin(t304526.57)104sin(2t4563.4345) 0.63sin(t48.43)0.79sin(2t63.43)
5-4已知系统开环传递函数
G(s)H(s)K(s1)s2(Ts1); K,,T0
试分析并绘制T和T情况下的概略幅相曲线。 解:
由题可知,系统的频率特性如下 G(j)H(j)K(j1)K(j1) 22(j)(Tj1)(Tj1)由于系统2,所以开环幅相曲线要用虚线补画180的半径为无穷大的圆弧 当0时,G(j0)H(j0) 当时,G(j)H(j)K(j01),(0)180 20(Tj01)K(j1)0,()180 2(Tj1)K(j1)K(T21)jK(T)又由于G(j)H(j)2,所以有
(Tj1)2(T221)当T时,开环幅相曲线始终处于第三象限,如图5-4a所示; 当T时,开环幅相曲线始终处于第二象限,如图5-4b所示。
图5-4a开环幅相曲线 图5-4b开环幅相曲线
5-5 已知系统开环传递函数
G(s)H(s)1
s(s1)(s2)试分别绘制1,2,3,4时系统的概略开环幅相曲线。 解:
由题目可知,系统的频率特性如下 G()H()1
(j)(j1)(j2)当1时,开环幅相曲线要用虚线补画90的半径为无穷大的圆弧。 若0,则G(j0)H(j0)1,(0)90
(0)1(01)(02)若,则G(j)H(j)10,()270
()1(1)(2)由以上分析可知,系统概略开环幅相曲线如图5-5a所示。
当2时,开环幅相曲线要用虚线补画180的半径为无穷大的圆弧。 若0,则G(j0)H(j0)1,(0)180 2(0)(01)(02)若,则G(j)H(j)10,()360 2()(1)(2)由以上分析可知,系统概略开环幅相曲线如图5-5a所示。
当3时,开环幅相曲线要用虚线补画270的半径为无穷大的圆弧。 若0,则G(j0)H(j0)1,(0)270
(0)3(01)(02)若,则G(j)H(j)10,()420
()3(1)(2)由以上分析可知,系统概略开环幅相曲线如图5-5a所示。
当4时,开环幅相曲线要用虚线补画360的半径为无穷大的圆弧。 若0,则G(j0)H(j0)1,(0)360 4(0)(01)(02)若,则G(j)H(j)10,()0 4()(1)(2)由以上分析可知,系统概略开环幅相曲线如图5-5a所示。
图5-5a系统开环幅相曲线
5-6已知系统开环传递函数
G(s)H(s)10
s(2s1)(s20.5s1)试分别计算0.5和2时,开环频率特性的幅值A()和相位()。 解:
系统的开环频率特性表达式如下
1010(242.52)10j(23)G(j)H(j)
j(2j1)(20.5j1)(242.52)2(23)2当0.5时
10(242.52)10j(23)G(j)H(j)(242.52)2(23)2168j
0.5此时A()1628217.,()arctan当2时
826.57 1610(242.52)10j(23)G(j)H(j)(242.52)2(23)220.32j0.21
此时A()0.3220.2120.38,()arctan
0.2133.27 0.325-7 绘制下列传递函数的对数幅频渐进特性曲线
a.G(s)1
(10.5s)(12s)
图2-7a对数幅频渐进特性曲线
b.G(s)(10.5s) 2s
图2-7b对数幅频渐进特性曲线
c.G(s)s10
s26s10
图2-7c对数幅频渐进特性曲线
d.G(s)30(s8)
s(s2)(s4)
图2-7d对数幅频渐进特性曲线
5-8 已知系统开环传递函数
G(s)K
s(s51)(s2001)试绘制K10的对数频率特性曲线,并算出截止频率c。 解:由题可得G(j)则
G(j)10(11K
j(j51)(j2001)225)(11222002)12
()2arctan5arctan200
因此 20lnG(j)2020ln10ln(1对数频率特性曲线如图5-8a所示
225)10ln(122002)
图5-8a 对数频率特性曲线
又20lnG(j)0,可得G(j)1,即10(11225)(11222002)121
计算可得c50rad/s
5-9 已知系统开环传递函数为:
G(s)H(s)50(s2) 2s11s10a.计算截止频率c。
b.确定对数幅频渐进特性曲线的低频渐进线的斜率。 c.绘制对数幅频特性曲线。
解:G(s)H(s)50(s2)50(s2)
s211s10(s1)(s10)50(j2)
(j1)(j10)122212212G(j)H(j) G(j)H(j)50(4)(1)(100)
20lnG(j)H(j)20ln(50(4)(1)(100))
221221212 20ln5010ln(42)10ln(12)10ln(1002)0 计算可得c49rad/s 当1时,斜率为0;
当12时,斜率为-20dB/d; 当210时,斜率为0; 当10时,斜率为-20dB/d;
绘制对数幅频特性曲线,如图5-9a所示。
图5-9a 对数幅频特性曲线
5-10 利用奈氏判据分别判断题5-4,5-5系统的闭环稳定性。 解: (1) 性。
当T时,系统的开环幅相曲线如图5-4a所示,由图可知,系统的开环幅相曲线不包围(1,j0),根据奈奎斯特判据可得N0 又由系统得开环传递函数可知P0
即ZP2N0,闭环系统在s右半平面无极点,T时闭环系统稳定。 当T时,系统的开环幅相曲线如图5-4b所示,由图可知,N1 又由系统得开环传递函数可知P0
即ZP2N2,闭环系统在s右半平面有2个极点, T时闭环系统不稳定。(2) 对于题5-5的系统,其开环幅相曲线如图所示,由图5-5a可知 当1时,N0,又由系统得开环传递函数可知P0
即ZP2N0,闭环系统在s右半平面无极点,1时闭环系统稳定。 当2,3,4时,N1,又由系统得开环传递函数可知P0
即ZP2N2,闭环系统在s右半平面有2个极点,2,3,4时闭环系统不稳定。
5-11 用劳斯判断据验证题5-10的结果。
解:
(1)对于题5-4的系统,由题得闭环系统特征方程为
对于题5-4的系统,分T和T的两种情况来讨论系统的闭环稳定
Ts3s2KsK0
列劳斯表
s3 T K s2 1 K s T-KTs0 K则当T时,KKT0,即第一列各值为正,即闭环系统稳定; 当T时,KKT0,即第一列各值不全为正,即闭环系统不稳定。 (2)对于题5-5的系统,由题得闭环系统特征方程为
sr(s1)(s2)10,即sr23sr12sr10
1
当r1时,列劳斯表
s3 1 2 s2 3 1 s 4 0 s0 1第一列各值为正,即闭环系统稳定; 当r2时,列劳斯表
1
s4 1 2 1s3 3 0 s2 2 32s0 0 s1 -第一列各值不全为正,即闭环系统不稳定;
当r3,4时,情况与r2相同,即闭环系统不稳定。
5-12 已知三个系统的开环传递函数为
G1(s)K(T2s1),
s2(T1s1)K(T2s1),
s2(T1s1)(T10,T20,T30,T40)
G2(s)G(s)K(T2s1)(T4s1), 3s(T1s1)(T3s1)又知它们的奈奎斯特曲线如图5-2(a)(b)(c)所示。找出各个传递函数分别对应的奈奎斯特曲线,并判断单位反馈下闭环系统的稳定性
图5-2 习题5-12控制系统乃奎斯特曲线图
解:三个传递函数对应的奈奎斯特曲线分别为b,c,a 对G1(s)K(T2s1)式,P0,N0 2s(T1s1)则ZP2N0,故系统稳定; 对G2(s)K(T2s1)式,P0,N0
s2(T1s1)则ZP2N0,故系统稳定; 对G(s)K(T2s1)(T4s1)式,P0,N0
s3(T1s1)(T3s1)则ZP2N0,故系统稳定;
5-13 已知系统开环传递函数
G(s)K; K,T0
s(Ts1)(s1)试根据奈氏判据,确定其闭环稳定条件: a.T2时,K值的范围; b.K10时,T值的范围; c.K,T值的范围。 解:
由系统的开环传递函数可知,系统的开环曲线图如图5-13a所示
图5-13a 系统开环曲线
由于P0,故想要闭环系统稳定,必有N0,即幅相曲线不包围点(1,j0)。 系统的频率特性表达式如下
KK2(T1)jK(T21)G(j) 24222j(Tj1)(j1)(T1)(T1)a、T2时,对于开环幅相曲线与实轴的交点有
K(T21)K(221)0 242224222(T1)(T1)9(21)由上式可得2,则交点的实轴坐标为 2K2(T1)3K21
(T1)242(T21)2942(221)2由上式可得 0K3 2b、K10时,对于开环幅相曲线与实轴的交点有
K(T21)10(T21)0
(T1)242(T21)2(T1)242(T21)2由上式可得1,则交点的实轴坐标为 T1(T1)K(T1)T1 (T1)242(T21)2(T1)211(T11)2T2TT1由上式可得 0T
9210c、对于开环幅相曲线与实轴的交点有
K(T21)K(T21)0
(T1)242(T21)2(T1)242(T21)2由上式可得1,则交点的实轴坐标为 T1(T1)K(T1)T1 (T1)242(T21)2(T1)211(T11)2T2TTT11由上式可得 0K ,0TTK12K
5-14 某系统的开环传递函数为
Q(s)K(T2s1)
s2(T1s1)要求画出以下4种情况下的奈奎斯特曲线,并判断闭环系统的稳定性: a.T20; b.0T2T1; c.0T2T1;
d.0T1T2。 解:
a. 当T20时,Q(s)K, 2s(T1s1)其开环幅相曲线如图5-14a所示,P0,N1
则ZP2N2,故在s平面右半平面有2个闭环极点,闭环系统不稳定;
K(jT21)K(1T1T22)K(T2T1)b.当0T2T1时,Q(j) 2(1jT1)2(1T122)若0,则|Q(j0)|,(0)180 若,则|Q(j0)|0,(0)180
其开环幅相曲线如图5-14b所示,P0,N1 则ZP2N2,故系统不稳定; c. 当0T2T1时,Q(s)K 2s若0,则|Q(j0)|,(0)180 若,则|Q(j0)|0,(0)180
1其开环幅相曲线如图5-14c所示,P0,N
2则ZP2N1,故系统不稳定;
K(jT21)K(1T1T22)K(T2T1)d.当0T1T2时,Q(j) 2(1jT1)2(1T122)由0T1T2可得Re[Q(j)]0,Im[Q(j)]0
故可得其开环幅相曲线如图5-14d所示,P0,N0
图5-14a 开环幅相曲线 图5-14b 开环幅相曲线
则ZP2N0,故系统稳定。
图5-14c 开环幅相曲线 图5-14d开环幅相曲线
5-15 已知反馈控制系统的开环传递函数为
Q(s)K,(T1,T1,T1,T10)
s(T1s1)(T2s1)(T3s1)(T4s1)如果闭环系统不稳定,闭环传递函数会有几个极点在复数平面的右半平面? 解:
|Q(j)|K(T1)(T1)(T1)(T1)
121212222122321224212(j)2arctanT1arctanT2arctanT3arctanT4
当0时,|Q(j)|,()90 当时,|Q(j)|0,()450
由于系统不稳定,故可得其开环幅相曲线如图5-15a所示
由图可得P0,N1
则ZP2N2,故闭环传递函数有2个极点在复数平面的右半平面。
图5-15a 开环幅相曲线
5-16 设控制系统的结构图如图5-3所示。
a.求出开环传递函数; b.画出对数相频特性曲线;
c.求出临界开环比例Kc和截止频率c;
d.用奈氏判据判断该系统是否稳定,如果稳定再分别求出当输入信号
u(t)1(t)和u(t)t的情况下系统的静态误差。
图5-3 习题5-16控制系统结构图
解:
(a)系统开环传递函数为G(s)0.1K1010••
0.1s13s1ss(13s)(10.1s)(b)()090arctan3arctan0.1 0,()0900090 ,()0909090270
图5-16a
(c)A()K19•10.0122
0,A(0) ,A()0 系统开环频率特性为
K[(T1T2)j(1T1T22)] G(j)H(j)
(1T122)(1T222)与实轴的交点 x10.131.826
G(jx)H(jx)Re[G(jx)H(jx)]故幅相曲线为
K30.10.3K
30.13.1
图5-16b
当G(jx)H(jx)1时,系统临界稳定,得Kc10,cx1.826 当0.3K1时,ZP2N000,系统稳定 3.1 图5-16c
当0.3K1时,ZP2N0(01)1,系统不稳定 3.1当u(t)1(t)时,ess0, 当u(t)t时,essRK0.1
5-17 已知某最小相位系统的开环对数幅频特性如图5-4所示。
a.写出其开环传递函数;
b.画出其相频特性草图,并从图上求出和标明相角裕度和幅值裕度; c.求出该系统达到临界稳定时的开环比例系数值K;
d.在复数平面上画出其奈奎斯特曲线,并标明点1j0的位置。
图5-4 习题5-17控制系统结构图
解:
(1)确定系统积分或微分环节的个数。因对数幅频渐近特性曲线的低频渐近线的斜率为20dB/dec,由图,低频渐近斜率为40dB/dec,故2,系统含有2个积分环节。
(2)确定系统传递函数结构形式。由于对数幅频渐近特性曲线为分段折线,其各转折点对应的频率为所含一阶或二阶环节的交接频率,每个交接频率处斜率的变化取决于环节的种类。
10.5处,斜率变化20dB/dec,对应微分环节;
22处,斜率变化20dB/dec,对应惯性环节; 25处,斜率变化20dB/dec,对应惯性环节。 因此,所测系统具有下述传递函数
K(1T1s) G(s)2s(1T2s)(1T3s)其中K待定。
(3)低频渐近线方程为
K(1s2(1111s)1
2s)(13s)La()20lnK20lnK20ln20lnK40ln
由给定点(,La())(0.5,8),得K10 故所测系统传递函数为G(s)
5-18设单位反馈控制系统的开环传递函数
10(12s) 2s(10.5s)(10.2s)G(s)as1 s2试确定相角裕度为450时的参数值。 解:
系统的频率特性表达式为 G(j)ja1 2设系统的截止频率为c,则由相角裕度的定义可得
180(c)180arctanac180arctanac45 即ac1 又由于G(j(c))(ac)2121
c2c2由上式得 c1.19 所以 a
10.84
c5-19 若高阶系统的时域指标为18%%25%,0.1ts0.2,试根据经验公
式确定系统的截止频率和相角裕度的范围。 解:根据经验公式,
0.160.4(tsK011)sin
c11)25% sin 根据题意有, 18%0.160.4(K0 0.1c0.2
可求得.772.2, 39.02c65.35
5-20 典型二阶系统的开环传递函数
2n G(s)s(s2n)若已知10%%30%,试确定相角裕度的范围;若给定n10,试确定系统带宽b的范围。 解:由于0.160.4(11)且10%%30%, sin可解得39.158.6
而根据题意 180G(jc)arctan(2(412))
2212又有bn((12)(12)1),且n10
2212故计算可得:11.61b14.11
5-21 设二阶系统如图5-5(a)所示。若分别加入测速反馈校正,0.1Kt1.5(图
5-5(b))和比例-微分校正,0.1Kd1.5(图5-5(c)),并设n1,0.2,试确定各种情况下相角裕度的范围,并加以比较。
图5-5 习题5-21控制系统结构图
解:(a)由题意可知系统开环频率特性
2n G(j)H(j)j(j2n)A()121242n.1.121n, ,()0arctan2n22设c为截止频率,当c时,则有
A(c)和
12c1c242n.1.12n1 2(c)0arctan12nc2
把n1,0.2代入上式,得:
c0.9608,(c)157.40 1800(c)1800157.4022.60 (b)由题意可知系统开环传递函数为
2n G(s)2s(s2nKtn)其开环频率特性为
2n G(j)H(j)2j(j2nKtn)A()1.122112K2.ntn.12n2Ktn()0arctan1 222nKtn设c为截止频率,当c时,则有
A(c)1c.12112K2.cntn2.12n2Ktn1
和
(c)0arctan1. c222nKtn把n1,0.2,设Kt0.6,代入上式,得:
c0.7862,(Kt0.6),(c)128.170 1800(c)51.830
(c)由题意可知系统开环传递函数为
G(s)K(Kds1),其中Kn2
s(s2n1)其开环频率特性为
G(j)H(j)K(jKd1) 2j(j2n1)A()22K1Kd1224n212,()0arctan12n.arctanKd2
设c为截止频率,当c时,则有
A(c)22K1Kdcc142c22n121
和
(c)0arctan12n.carctanKdc2
把n1,0.2,设Kd1,代入上式,得:
c1.226,(Kd1),(c)111.10 1800(c)68.90
5-22 已知单位反馈系统的开环幅相特性曲线如图5-6所示。当K50时,系统
幅值裕度h1,穿越频率c1,求输入为r(t)t25sinct,幅裕度为下述值时,系统的稳态误a.h0.5 b.h3
图5-6 习题5-22控制系统结构图
试值差。
解:设系统开环传递函数为:
G(s)K(1T2s)
s(1T1s)开环系统幅频特性为:
A()K1T2221T212
系统的开环频率特性为:
K[(T1T2)j(1T1T22)]G(j)
(1T22)解得x1 T1T2501T2211 当g1有x1,
2T1T21T1得T20.02,T150
则系统开环传递函数可写成
G(s)K(10.02s)
s(150s)系统与实轴的交点为G(jx)KT20.02K
当h0.5时, K100,ess11000.01 当h3时,K16.7,ess116.70.05
5-23 设单位反馈系统如图5-7所示。其中,截止频率c5,K10;T0.1时,
若要求c不变,问K与T如何变化才能使系统相角裕度提高至450? 解:开环系统幅频特性为:
A()K(1T22).A1()
1212相频特性为:
()0arctanTarctan1() 当c时,
(c)0arctanTcarctanc1(c),把K10;T0.1;c5代入得: (c)0026.57078.6901(c)
1800(c)
若要求相角提高450,即要求(c)提高450,设调整后的系统相频特性为:
(c)00(26.570450)78.6901(c)
调整后的T值为:[tan(26.570450)]c0.6,K值不做调整。
图5-7 习题5-23单位反馈系统结构图
5-24 已知单位反馈控制系统的开环传递函数为
Ke0.1sG0(s)
s(s1)试绘制系统的对数频率特性曲线,并据此确定: a.求K1时的相角裕度; b.求K20时的幅值裕度; (1)解: 开环系统幅频特性为:
A()K12
令A()1,当K1时,得c0.75 开环系统相频特性为:
()02arctan,当c时,有 arctan0.75126.90
(c)021800(c)53.10
(2) 解:开环系统的频率特性为:
G(j)20[(sin0.1cos0.1)j(cos0.1sin0.1)]令其虚部为零,即
(12)cos0.1sin0.10
tan0.11得x25
G(jx)20[(sin0.1xxcos0.1x)]0.03
x(1x2)h133.3
G(jx)5-25 若单位反馈系统的开环传递函数
Ke0.8sG(s)
s1试确定使系统稳定的K值。 解:
Ke0.8j系统的频率特性表达式为 G(j)
j1由上式可得,系统的幅频特性和相频特性分别为 G(j)K,()arctan0.8
12系统临界稳定时开环幅相曲线穿过点(1,j0),此时
G(j)K1,()arctan0.8180
12由上式可得, 2.45,K2.45212.65
显然,当0K2.65时,由奈奎斯特稳定判据可得系统闭环稳定。 故K的取值范围为 0K2.65
5-26 设单位反馈系统的开环传递函数
5s2esG(s)
(s1)4试确定闭环系统稳定时,延迟时间的范围。 解:
52ej系统的频率特性表达式为 G(j)
(j1)4由上式可得,系统的幅频特性和相频特性分别为
G(j)5221,()1804arctan2180
系统临界稳定时开环幅相曲线穿过点(1,j0),此时
G(j)521,()1804arctan180180
212由幅频特性可得 215
解之可得 11.618,20.618(舍去) 又(1)1804arctan1.6181.618180180 即1.37
显然,当01.37时,由奈奎斯特稳定判据可得系统闭环稳定。故的取值范围为
01.37
6-1 设单位反馈系统开环传递函数为:
G(s)200
s(0.1s1)试设计一无源校正网络,使已校正系统的相角裕度不小于450,截止频率不低于50。 解:
作待校正系统对数幅频特性L'(),如图6-1a所示,得c'40,'14,故应选择超前网络。
图6-1a系统特性
取c''m55,量得L'()6dB,由10lgaL'(m),求得
a4,T10.009
ma取无源超前网络aGc(s)1aTs10.036s 1Ts10.009s将放大增益提高4倍,作校正后系统L''(),见图6-1a,得满足设计要求得如下指标:
c55, (c'')47.2 6-2 设单位反馈系统的开环传递函数:
G(s)K s(s1)试设计一串联超前校正装置,使系统满足如下指标: (1) 相角裕度450;
(2) 在单位斜坡输入作用下的稳态误差ess(3) 截止频率c7.5rads。 解:首先确定开环增益ess则开环传递函数为G(s)1rad; 1511得K15,取K15 k1515 由G(jc)1,得未校正前的截止频率
s(s1)c3.81
对应得相角裕度'900arctanc'14.710不符合要求,进行串级超前校正。 取mc''7.5 计算L'(7.5)20lgG(j7.5)11.56 由10lgaL'(7.5)11.56,得a14,所以T所以设计得超前网络传递函数为G(s)10.036
ma0.504s1
0.036s1150.504s1
s(s1)0.036s1最终校正系统的开环传递函数为Gc(s)G(s)验算''900arctanc''arctan0.504c''arctan0.036c''67.745 满足性能指标要求,设计合理。
6-3 已知单位反馈系统的开环传递函数为:
G(s)4K
s(s2)试设计串联校正装置,使校正后系统的相位裕度500,幅值裕度h10dB,静态速度误差系数K20s1。
解:给定系统的稳定裕量时宜采用频率响应校正设计方法。
确定期望的开环增益K。因为KvlimGp(s)2K,所以取K10。
s0分析增益校正后的系统。图6-2中的虚线为Gp(j)40的对数幅
j(j2)频特性和相频特性。图6-3a中的对数幅频特性采用的是渐近线,渐近线的拐点处的分贝数用数字表示,相频特性为示意图。
从图6-3a虚线所示的对数幅频特性可以测算出增益穿越频率
c6.325rad/s
相位裕量17.550。校正的任务是增加相位裕量。由图可以看出,采用超前校正,可以提高相位裕量。因为增益已经确定,所以超前校正装置采用
Gc(s)1Ts,a1的形式。在0时,Gc(j)1,因此校正装置不会影响低1aTs频增益,故而不会改变已获得的静态误差系数。
图6-3a系统校正前后的伯德图
由500可得5517.5533,并进而取m538。
(1sinm)/(1sinm)0.24
超前校正装置的最大相角频率为m1/(aT),而且在该频率的增益为
Gc(j)m1/a。要使增益穿越频率等于m,Gp(j)Gc(j)曲线必须在m处穿过轴,即
20lg|Gp(jc)Gc(jc)|20lg|Gp(jc)|20lg|Gc(jc)| 20lg|Gp(jc)|10lg0所以 20lg|Gp(jc)|10lg6.20dB
由图6-2可以算出c9.04rad/s。进而取m9.04可得 T10.2258,114.4286,T0.019,18.4528 TTm故校正装置的传递函数为 Gc(s)
6-4 设系统开环传递函数为:
G(s)K
s(10.1s)10.0226ss4.429 4.16710.02ss18.453试用比例—微分装置进行校正,使系统K200,(c'')500,并确定校正参数。
解:首先确定开环增益,取KK200 所以未校正开环传递函数为G0(s)计算校正前截止频率为c44.7
计算相角裕度为90arctan4.4712.6
相角裕度低于性能指标,可用比例微分装置进行校正。设比例微分校正装置传递函数为
Gc(s)Ts1
200
s(0.1s1)需要补偿德超前角为05012.637.4 取38, 又因为arctancTarctan44.7T38 可得比例微分装置的时间常数T0.0175
所以比例微分校正装置的传递函数为Gc(s)0.0175s1 校正后系统的开环传递函数为G(s)200200(0.0175s1)•(0.0175s1)
s(0.1s1)s(0.1s1)验算90arctan0.1carctan0.0175c50.650 符合要求。
6-5 设单位反馈系统的开环传递函数为:
G(s)K
s2(10.2s)试设计串联校正装置Gc(s),使系统的Ka10,(c'')350。
解:取KKa10,绘待校正系统L'(),'(),如6-5a图,由图6-5a查得
c'3.16,(c')32
采用超前网络,其最大超前相角应为max3235976 由于max较大,应采用两级超前校正,每级'max38
1s/1超前网络传递函数为Gc(s),2a1
1s/22依据sin'maxa1,算得a4.204,取a5,故251 a110(1s/1)2校正后系统开环传递函数为Gc(s)G(s)2 2s(10.2s)(1s/2)10(c''/1)2当c''时,因为c''是1和2得几何中点,因此可得A()''21
c0.2c''''c从而求得c''5012,代入c''a151, 有12.8,214,于是c''6.26
图6-5a系统的对数频率特性
经放大补偿后,
(1s/2.8)2Gc(s) 2(1s/14)验算:(c'')180''(c'')180147.832.235
故将a增大,取a6,c''a16.86,2a116.8,算得(c'')37.3
(10.357s)2满足设计要求。因此 Gc(s),放大器增益需提高6倍。 2(10.06s)
6-6 设单位反馈系统开环传递函数为:
G(s)7
11ss1s126试设计一串联滞后校正网络,使已校正系统的相角裕度为40020,幅值裕度不低于10dB,开环增益保持不变,截止频率不低于“1”。 解:画校正前系统开环对数频率特性曲线,由图6-6a得
c'3.8,'(c')4.4,20lgh1dB
表明待校正系统不稳定,由于c'大于要求得c'',故可采用串联滞后校正。
图6-6a系统的对数频率特性曲线
11由90arctanc''arctanc''画(c'')曲线
26根据题目意思,估计c(c'')6,而''40,因此(c'')''c(c'')46 由(c'')曲线查得c''1.25 ,满足c''1的要求。 当c''1.25时,b(dB)L'(c'')15dB,故b0.18
令1bT0.1c'',求得T44.4。于是串联滞后网络传递函数为
Gc(s)1bTs18s 1Ts144.4s校正后系统开环传递函数为
Gc(s)G(s)1(18s)
11s(1s)(1s)(144.4s)26验算:c''1.25, 41.2, 20lgh10dB,满足要求。
6-7 对于题6-4试用比例—积分装置进行串联校正。
解:加PI控制器后,系统成为Ⅱ型,有Kv,必满足稳态性要求。因此K可取任何满足要求得任意值。因待校正系统(c'')90arctan0.1c''。今要求
(c'')50,再考虑校正元件产生的滞后相位,可以选择c'',使arctan0.1c''30。 设取c''6,由题可得K6,所以未校正开环传递函数为G(s)6
s(0.1s1)选PI校正装置传递函数Gc(s)1Ts Ts根据1T0.1c'',取1T0.6,则校正后系统的传递函数为
G(s)Gc(s)3.6(11.67s)
s2(10.1s)校正后的相角裕度(c'')arctanc''0.6arctanc''1053.3满足设计要求。
6-8 已知单位反馈系统的开环传递函数为:
G(s)K
s(1s)(10.5s)试设计串联校正装置使系统具有相位裕度400,幅值裕度h10dB,静态速度误差系数K5s1。
解:由题意 KvlimGp(s)K可得K5
s0 画Gp(j)的伯德图,从图6-8a中可以看出,只要将对数幅频渐近线德中段下降,即可满足相位裕量要求。所以,采用滞后校正,传递函数取为
Gc(s)1Ts,1
1Ts令401252,它对应的相角为18052128 对应于该期望相位裕量的频率为
argGp(j)arg190arctanjarctanj0.5128j(1j)(10.5j)解得0.5rad/s,即为新的增益穿越频率c
图6-8a系统校正前后的伯德图
令1Tc/50.1,得T10
因此,20lgGp(jc)Gc(jc)0dB,即20lg20,解得10。 滞后校正装置的传递函数为:Gc(s)1Ts1s0.1•
1Ts10s0.015(10s1)
s(1s)(10.5s)(1100s)校正后系统的开环传递函数为Gp(s)Gc(s)验算c0.5rad/s,39,h11dB,K5s1,符合要求。
6-9 设单位反馈系统开环传递函数为:
G(s)126
11ss1s11060要求设计—串联校正装置,使系统满足: (1) 输入速度为1时,稳态速度误差不大于1126; (2) 许可的放大器增益不变;
(3) 相角裕度不小于300,截止频率为20。
解:绘待校正系统L'(),由图6-9a得c'36,算出'15.5。表明待校正系统不稳定,且c'c''要求,宜采用串联滞后-超前校正。
图6-9a系统对数幅频特性
由图知,b10,则Tb1/b0.1,于是在60时,L''()的斜率均为-20。
由于要求c''20,故可得16dB 因此,已校正系统开环频率特性为
Gc(j)G(j)126(1ja)
j(1j/60)(1j/a)(1j/63)令''20,30,由
''90arctan20/aarctan20/60arctan126/aarctan20/63 可以求出a8.9,于是校正网络为
Gc(s)(10.112s)(10.1s)
(10.71s)(10.016s)验算:''33.730,满足指标要求。
6-10 已知单位反馈系统的开环传递函数为:
G(s)K
s(10.1s)(10.01s)试设计串联滞后超前校正装置使校正后系统具有相角裕度400,增益穿越频率c20rads,静态速度误差系数K100s1。 解:由题意 KvlimsGp(s)K可得K100
s0 所以增益校正后的开环传递函数为Gp(s)100
s(10.1s)(10.01s)由题意,采用超前和滞后分别设计的滞后超前装置,即
G(s)Gc1(s)Gc2(s), Gc1(s)1T1s1T2s, Gc2(s)
1T1s1T2s其中1,1,且不一定等于1。
设计超前部分:根据题目要求,超前部分至少应再提供超前角25,故取超前装置的最大超前角为m40,由此可算的,由1T1c0.2174209.33,故T10.107 因此,超前部分的传递函数为Gc1(s)1T1s10.107s
1T1s10.023s1sin400.2174
1sin40由于它的零点和对象的一个极点十分接近,故该取T10.1,0.23 所以Gc1(s)10.1ss10 4.34810.023ss43.48设计滞后部分:要使20rads成为增益穿越频率,必须满足
可解得20lgGc2(j20)20lg14dB,即5 Gp(j20)Gc1(j20)Gc2(j20)1,令1T2c/54,得T20.25,所以 滞后部分的传递函数为Gc2(s)1T2s10.25s1s4
1T2s11.25s5s0.8从而可得,超前滞后装置的传递函数为
Gc(s)(10.1s)(10.25s)(s10)(s4)0.87
(10.023s)(11.25s)(s43.48)(s0.8)校正后的开环传递函数为
Gp(s)Gc(s)100(10.25s)
s(10.01s)(10.023s)(11.25s)验算,c20rad/s,44.97,K100s1,符合要求。
6-11 已知系统开环传递函数为:
G(s)K
s(10.5s)(10.1s)试设计PID校正装置,使系统K10,(c'')500且c''4。 解:令KKv10,作待校正系统L'(),如图6-11a。 由图6-7知,c'4.47,(c')0
图6-11a 系统特性曲线
设PID校正装置传递函数为Gc(s)则校正后系统频率特性为
(11s)(12s)(1s1)(1s2)
1ss1G(j)Gc(j)K1(1j1)(1j2), K1K1
(j)2(1j2)(1j10)由于校正后为Ⅱ型系统,故Kv的要求肯定满足,系统开环增益可任选,由其他条件而定。
初选c''4。为降低系统阶次,选22,并选10.1,c''0.4,此时
G(j)Gc(j)''cK1(1j0.4)
(j)2(1j10)其对数幅频特性应通过截止频率4,故由近似式得K11.6,从而K4。
验算:(c'')arctanc''0.4arctanc''1062.5 满足设计要求。
K1c''/0.4''2c1
6-12 设系统结构图如图6-1所示,图中,G1(s)K1,
T0s1G2(s)K2K,G3(s)3,其中,K10~6000可调,K212,
(T1s1)(T2s1)s
图6-1 习题6-12示意图
K31400,T00.014,T10.1,T20.02。试设计反馈校正装置H(s),
使系统满足K150,%40%,ts1。
解:令K15000,画待校正系统G0(s)150的对数幅
s(10.014s)(10.02s)(10.1s)频特性曲线,如图6-12a。由图6-9得c'38.7 由对数幅频特性可得G(s)150(10.25s),因此
s(10.014s)(12.86s)(10.0133s)G2(s)H(s)2.86s
(10.25s)(10.1s)(10.02s)图6-12a 系统对数幅频特性曲线
当475时
(4)18090arctan0.254arctan0.14arctan0.02444.3,因而小闭环稳定。
由于G2(s)已知,故有H(s)0.238s0.25s 0.9510.25s10.25s验算,K150,,Mr1.27,%27%,ts0.63,符合设计要求。
6-13 设系统结构图如图6-2所示,待校正系统的开环传递函数为:
40
s(10.003s)G0(s)试用三阶最佳工程设计法设计校正装置Gc(s)。
图6-2 习题6-13示意图
解:选Gc(s)作为PI调节器,即Gc(s)s1Ts
40(1s)/T
s2(10.003s)校正后系统开环传递函数为G(s)Gc(s)G0(s)确定校正装置参数和T,40.0030.012,40/T1/(80.0032),即
T0.0029
于是PI调节器为Gc(s)10.012s14.14(1)
0.0029s0.012s本题如采用最小Mr设计法,取H5,则可得PI调节器为
Gc(s)10.015s15(1)
0.003s0.015s6-14 设复合控制系统如图6-3所示,图中,Gn(s)为顺馈装置传递函数;
G2(s)1s2。Gc(s)Kts,G1(s)K1;为测速发电机及分压器的传递函数;
试确定K1,Gn(s)及Gc(s),使系统输出量完全不受扰动n(t)的影响,且单位阶跃响应超调量%25%,峰值时间tp2。
(a)
(b)
解:
(1G1G2GC)GnG2s2K1GC(s)Gn(s) N(s)21G1G2GCG1G2sK1GC(s)K1要使系统输出完全不受扰动影响,应使N(s)0,于是s2K1GC(s)Gn(s)0 即 Gn(s)s2K1GC(s)s2K1Kts
系统对输入的开环传递函数为
2nG1G2K1 G(s)221G1G2GCsK1Ktss2ns按题意要求 %e120.25,tpn122
从而解得 ln4(ln4)220.404,n211.7172
22.984,Kt因此 K1n2n0.471 K1于是 Gc(s)Kts0.471s,Gn(s)s2K1Ktss21.387s
1.计算机控制系统的组成部分有哪些?它的特点是什么?
答:计算机控制系统由硬件部分和软件部分组成,其中硬件部分主要由主机、外部设备、过程输入输出设备和广义被控对象组成,软件部分包括系统软件和应用软件。计算机控制系统的功能强大而且安全可靠。归纳起来,主要有如下几个特点:
(1)可以同时实现模拟变送器、控制器、指示器、手操器以及记录仪等多种模拟仪表的功能,并且便于集中监视和操作;――利用了计算机的存储、数字运算和显示功能;
(2)一台计算机可以同时控制多个回路,并且还可以同时实现DDC、顺序控制、监督控制等多种控制功能――利用了计算机的快速运算功能;
(3)可以实现模拟控制难以实现的复杂控制规律,如最优控制、自适应控制、多变量控制等等――利用了计算机强大的信息处理能力;
(4)计算机控制系统的调试、整定灵活方便,系统控制方案、控制策略以及控制算法只需要修改软件即可实现;
(5)利用网络分布结构可实现计算机控制管理集成系统; (6)计算机控制系统中同时存在连续性和离散性两类信号。 2.计算机控制系统与经典自动控制系统在信号上有什么不同?
答:经典自动控制系统的信号是模拟信号,而计算机控制系统的信号是数字信号。 3.什么是采样定理?采样周期的一般选择原则是什么?
答:一个连续时间信号f(t),设其频率带宽是有限的,其最高频率为m,如果在等间隔点上对该信号f(t)进行连续采样,为了使采样后的离散信号f*(t)能包含原信号f(t)的全部信息量。则采样角频率(s)只有满足表达式s2m,采
样后的信号f*(t)才能够无失真地复现f(t),否则不能从f*(t)中恢复f(t),这个定理就称为采样定理。
采样周期的一般选择原则是:TS4.试推导零阶保持器的传递函数? 解:零阶保持器的输入输出关系式为
fh(nTt)fh(nT) 0≤Δt gh(t)1(t)1(tT) 对脉冲过渡函数gh(t)取拉氏变换,可得零阶保持器的传递函数: 1eTs1eTsGh(s) sss5.求下列函数的Z变换 (1) f(t)1eat 6 s(s2) 1(2)f(t) 4(4)F(s)s2 (s1)(s3)t(3)F(s) 解:(1)F(z)Z[1eat] Z[1(t)]Z[eat] zz aTz1zez(1eaT)2 z(1eaT)zeaT11(2)F(z)Z[]zn1(4Tz)1(4Tz)2(4Tz)n 4n04tnT在上式中,若(4Tz)11,则无穷级数是收敛的,利用等比级数求和公式,可得 14TzF(z)T T11(4z)4z1(3)将F(s)展开成部分分式: F(s)633 s(s2)ss2查Z变换表,可得: F(z)3zz 3z1ze2T3z(1e2T)2 2T2Tz(1e)ze(4)将F(s)展开成部分分式: F(s)s21111 (s1)(s3)2s12s3查Z变换表,可得: F(z)1z1z T3T2ze2ze1z(eTe3T)22T z(ee3T)ze4T6.求下列函数的初值和终值 10z1 (1) F(z) (1z1)2 14z13z2(2)F(z) 12z16z22.5z3z2(z2z1)(4)F(z)2 2(z0.8z1)(zz0.8)z5(3)F(z)2 z4z3解:(1)由初值定理得 f(0)limf*(t)limF(z) t0z10z1lim10 z(1z1)2f()limf*(t)lim(z1)F(z) tz110z1 lim(z1) 12z1(1z)(2)由初值定理得 f(0)limf*(t)limF(z) t0z14z13z2lim1 z12z16z22.5z3f()limf*(t)lim(z1)F(z) tz114z13z2lim(z1)0 z112z16z22.5z3(3)由初值定理得 f(0)limf*(t)limF(z) t0zlimtz50 zz24z3z1f()limf*(t)lim(z1)F(z) lim(z1)z1z50 2z4z3(4)由初值定理得 f(0)limf*(t)limF(z) t0zz2(z2z1)lim21 z(z0.8z1)(z2z0.8)由于F(z)的四个极点都在单位圆内,故可以用终值定理求解。 f()limf*(t)lim(z1)F(z) tz1z2(z2z1)lim(z1)2 z1(z0.8z1)(z2z0.8)0 7. 求下列函数的Z反变换 z(1) F(z) z0.5 z2 (2)F(z) (z0.8)(z0.1)0.5z2(3)F(z) (z1)(z0.5)(4)F(z)tTz (z1)(z2)解:(1)查Z变换表,可得f(t)2(2)将F(z)展开,得 8z1z F(z)7z0.87z0.1查Z变换表,可得 tt81TTf(t)0.80.1 77(3)将F(z)展开,得 F(z)z0.5z z1z0.5查Z变换表,可得 f(t)10.50.510.5tTt1T (4)将F(z)展开,得 F(z)12 z1z2查Z变换表,可得 f(t)122t1T21 tT8. 求下列系统的Z传递函数 (1) y(k)2y(k2)3y(k4)u(k)u(k1) (2) G(s)2 s(0.1s1)解:(1)对方程两边取Z变换得: (12z23z4)Y(z)(1z1)U(z) Y(z)1z1 G(z) U(z)12z23z4(2)将G(s)展开成部分分式 G(s)222 s(0.1s1)ss10查Z变换表,可得: G(z)2z2z 10Tz1ze2z(1e10T)2 z(1e10T)ze10T 图7-10 习题9系统示意图 9.已知一离散系统如下图所示,试分析要使系统稳定,K的取值范围。 解:从上图可以得出系统的闭环传递函数为: W(z)KT z1KT根据系统稳定性的充要条件可得: 当1KT1 0K2时,系统稳定,否则系统不稳定 T10.已知一离散系统的闭环特征方程为: D(z)z31.03z20.43z0.000 试用Routh判据判断系统的稳定性。 解:针对上式做W变换,即将zw1代入得: 1wD(w)2.45w33.62w21.52w0.40 作Routh表如下: w3 w2 w1 w0 2.45 3.62 1.25 0.4 0 1.52 0.4 0 由于Routh表第一列均大于零,故系统稳定。 因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- oldu.cn 版权所有 浙ICP备2024123271号-1
违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务