当前位置:文档之家› MATLAB的根轨迹分析法及重点习题

MATLAB的根轨迹分析法及重点习题

4.1某系统的结构如题4-1图所示,试求单位阶跃响应的调节时间t s ,若要求t s =0.1秒,系统的反馈系数应调整为多少?解:(1)由系统结构图可知系统闭环传递函数为:100()100()1001()()1001*G s s s G s H s s aa sΦ===+++ 在单位阶跃函数作用下系统输出为:12100()()()(100)100k k C s R s s s s a s s a=Φ==+++为求系统单位阶跃响应,对C(s)进行拉斯反变换:1021001001001001lim ()lim1001001lim (100)()lim 11()(100)1()(1)s s s as a at k sC s s a ak s a C s s aC s as a s a c t e a→→→-→--===+=+==-=-+=-根据定义调节时间等于响应曲线进入5%误差带,并保持在此误差带内所需要的最短时间,且根据响应系统单位阶跃响应的函数表达式可以看出系统单位阶跃响应的稳态值为1a,因此: 10010011()(1)0.950.051ln 201001=0.1ln 20=0.3s 10s s at s at s s c t e a a e t a a t --=-=⇒=⇒==因为题中,所以(2)若要求t s =0.1秒,则有:1ln 20=0.1100=0.3s t aa =⇒ 即:若要求调节时间缩小为0.1秒,则需将反馈环节的反馈系数调整为0.3。

4.2已知二阶系统的阶跃响应曲线如题4.2图所示,该系统为单位负反馈系统,试确定其开环传递函数。

解:根据系统阶跃响应曲线可以看出: 峰值时间=0.1s p t ,超调量 1.3-1%=100%30%1σ⨯=; 根据课本中对典型二阶系统222()2nn ns s s ωζωωΦ=++暂态性能指标的推导计算可知:%p t e σ-==结合本题已知阶跃响应曲线可知:0.1(1)%30%(2)p t e σ-====由式(2)可知:0.3ln 0.30.3832cot =0.3832=arccot 0.3832=69.0332=cos =0.3578eζϕζϕζϕ-=⇒-=⇒==即:将ζ带入式(1)中可得:0.1p n t ω==回顾题意对于典型二阶系统其闭环传递函数为222()2nn ns s s ωζωωΦ=++,且系统为单位负反馈系统,所以系统开环传递函数和闭环传递函数之间满足如下关系:222222222211()()121211211131.8851===224.0753n n n nn n n n n G s s s s G s s G s s G G s s s sωζωζωωωζωωωζωΦ==Φ==+++++++++,因为:所以:,4.3单位反馈控制系统开环传递函数为()(1)KG s s Ts =+,若116s =0.25s K T -=、,试求(1)动态性能指标%(0.05)s t σ∆=、.(2)欲使%=16%σ,当T 不变时,K 应取何值。

解:(1)对于单位反馈控制系统,已知开环传递函数可求出其闭环传递函数,并将其化为标准形式为:22()(1)()1()1(1)1180.25n n n K KG s K s Ts T s K s K G s Ts s K s s Ts T T T Tωζωζω+Φ====+++++++===⇒==即:;22 所以根据动态性能指标的计算公式将上述两参数带入后可得:0.2511ln( 1.5142%44.43%s n t s e e ζωσ--=-=-====(2)由于T=0.25s,所以可知: 1n n T ωζωζ==⇒=2将阻尼比带入超调量的计算公式中:%16%0.16ln 0.16ln 0.16 3.9388e e K σζπ--==⇒=⇒-==⇒=将阻尼比带入可得:4.4设控制系统如题4-4图所示,其中图(a)为无速度反馈系统,图(b)为带速度反馈系统,试确定系统阻尼比为0.5时K t 的值,并比较图(a)和图(b)系统阶跃响应的动态性能指标。

解:(1)根据系统结构图可求得两系统的闭环传递函数为:2210()101(1)(),101()()1021(1)10()10(1)(),101()()10101(1)(1)1+100.21622a n n b t tt n t n G s s s s G s H s s s s s G s s s s G s H s s s K s K s s s K K ωζωωζω+Φ=====++++++Φ===+++++++===所以所以因此(2)根据已经求得的两系统的阻尼比和无阻尼自然振荡角频率可分别计算两系统的动态性能指标:0.55351.0053%60.47%1ln( 6.0168ra pa a sa nt s t s e t sσζω-=======-=0.76481.1471%16.3%1ln( 1.9827rb pb b sb nt s t s e t sσζω-=======-=经对比可看出:采用速度反馈的b 系统虽然上升时间和峰值时间稍有延长,但超调量存在明显下降,系统振荡剧烈程度下降,另外调节时间也显著降低,即使说系统能够在较快的时间内达到稳定,系统动态性能得到了提高。

4.5某系统结构如题4-5图所示,试判断系统的稳定性。

解:根据系统结构图可利用梅森公式求解其传递函数,结构图中前向通道有一条,回路有两个,且两回路相关,因此有:2212111321010(+1)20101(2)1(1)(1)11011010()21101s s s s s s s sP s P s s s s s ∆=--⋅-=+++++=∆=∆+Φ==∆+++因此可得系统特征方程为:32211010s s s +++=列写其劳斯表为:3211021120902111s s s根据劳斯判据可知,劳斯表第一列系数符号均未发生变化,因此系统稳定。

4.6已知系统特征方程如下,用劳斯判据判断系统的稳定性,如不稳定求在s 右半平面的根数及虚根值。

另外用MATLAB 软件直接求其特征跟加以验证。

(1)5432543224312243248011232324484160124800480,148s=2js s s s s s s s s s +++++=+=±2劳斯表中出现全零行,根据上一行数据列写辅助方程为12s 对其求导得到全零行的新系数为:240,求解辅助方程可得通过分析劳斯表中第一列系数可知并没有符号变化,所以不存在位于s 右半平面的特征根,另外由于系统劳斯表中出现全零行,所以系统不稳定,存在对称于原点的根为s=2j ±用MATLAB 软件中函数roots 求特征方程根可得,系统特征方程根为如下所示,进一步验证了上述求解结果的正确性。

(2)5432543210312203525011235320251680(16)(80)033525000250,125s=s s s s s s s s s s +++++=+=±2同时乘以非正系数3可简化计算劳斯表中出现全零行,根据上一行数据列写辅助方程为5s 对其求导得到全零行的新系数为:100,求解辅助方程可得通过分析劳斯表中第一列系数可知并没有符号变化,所以不存在位于s 右半平面的特征根,另外由于系统劳斯表中出现全零行,所以系统不稳定,存在对称于原点的根为s=±用MATLAB 软件中函数roots 求特征方程根可得,系统特征方程根为如下所示,进一步验证了上述求解结果的正确性。

(3)6543265432(-20)(10)24447810014710448055100005100,510s=1290110s s s s s s s s s s s s s -+-+--+=------+=-±-4劳斯表中出现全零行,根据上一行数据列写辅助方程为-5s 对其求导得到全零行的新系数为:-20-10,求解辅助方程可得通过分析劳斯表中第一列系数符号变化两次,所以有两个位于s 右半平面的特征根,另外由于系统劳斯表中出现全零行,所以系统不稳定,存在对称于原点的根为s=1±±用MATLAB 软件中函数roots 求特征方程根可得,系统特征方程根为如下所示,进一步验证了上述求解结果的正确性。

(4)65432(4)(6)212741340134000340,34s=j,s=2250314s s s s s s s ----------=--±±--4劳斯表中出现全零行,根据上一行数据列写辅助方程为s 对其求导得到全零行的新系数为:4-6,求解辅助方程可得通过分析劳斯表中第一列系数符号变化1次,所以有1个位于s 右半平面的特征根,另外由于系统劳斯表中出现全零行,所以系统不稳定,存在对称于原点的根为s=j,s=2±±用MATLAB 软件中函数roots 求特征方程根可得,系统特征方程根为如下所示,进一步验证了上述求解结果的正确性。

4.7已知单位反馈系统的开环传递函数为2(0.51)()(+1)(0.51)K s G s s s s s +=++,试确定系统稳定时的K 值范围。

解:根据已知单位反馈系统开环传递函数可知系统特征方程为:4320.5 1.52(0.51)0s s s K s K +++++=,对系统列写劳斯表可得:4320.521.50.5100.51234.50.516(0.51)1s Ks K K s KK s K K K++-+--+欲使系统稳定则需满足:0.5120103K K +->⇒<4.50.510556(0.51)KK K K +->⇒--<<-+-+0K >汇总得:05 1.708K <<-+=4.8已知系统的特征方程为321340400s s s K +++=,试确定系统稳定时的K 值范围,若要求闭环系统极点均位于s=-1垂线之左,K 值该如何调整。

解:(1)对系统列写劳斯表为:321401340404013140s s K K s K-欲使系统稳定则需满足:40400;40013KK ->>,即:013K <<。

(2)将1s s =-带入系统原特征方程中得:323210174028011710402819.8414028s s s K s s K s K K +++-=---欲使系统稳定则需满足:19.840 4.95K K ->⇒<;402800.7K K ->⇒>;即0.7 4.95K <<4.9已知系统稳定,求2()1()2t r t t t =++的系统稳态误差。

相关主题