资源描述:
《控制系统的MATLAB仿真与设计课后习题答案(王海英,袁丽英,吴勃编著)高等教育出版社》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、第八章8.1创建连续二阶系统和离散系统的传递函数模型。ˉ2′(3)C(z)=5z(1)C(s)=挠:s2十)〓(2)Cls·2s+2εz2-1.5z+0.5程序如下:C△ear`c1c`G1=t£(5`[1`2`2])G2=tf(5`[1`2`2]`li∶j∵;`2)∶o∷
2、0△∶G3=tf([0.5`0]`[1`-1.5`0.5]`0.1)2(s-+0.5)8.2已知系统的传递函数为G$)=建立系统的传递函数模型,并转换为零极点模型和状(s+2)2+1态空间模型。程序如下:C△ear氵clc氵num=2*[1`0.5]氵den=[1`0.2`
3、1.01]氵sys1=tf(num`den)[z`p`k]=tf2zp(num`den)氵sys2=zpk(z'p`k)[A、B`C`D]=t£2sS(num`den)氵sys3=ss(A`B`C`D)ⅢElkg系统的方框图如图所示。氏=L凡=⒉G=⒊G=仉计算系统的Φ⑴=号=r。黯←(1)结构图如下:TmnζerFcn2TranζerFcn3(2)程序如下:Clear`clc`R1=1`R2=2氵C1=3`C2=4氵l)氵[A`B`C'D]≡linm°d(l^∵:∶[num`den]=SS2tf(A`B`C`D)`sys=tf(num、den
4、)·Ι110nυ0l000一00^υB〓〓8.5己知系统的系数矩阵分别为/0乙,试判断系统的能控性。2程序如下:C△ear'c△c`A=[1、1、0氵0'1'0'0、0、2]'B=[0'0`l、0`O'-2]丿n=size(A)Tc=ctrb(A`B)`i£n==rank(TC)d⊥sp(帚乐-∶∵已-I茫i}∮)氵elsedisp(I系铰∷彳Ⅱt个走扌}’)'end第九章9.1系统的传递函数为Gls)=午’绘制出其根轨迹、伯德图和奈奎斯特图。毛毛⒈足÷程序如下:Clear`clc`num=[2`5`1]氵den=[1`2`3]'sys=t£(
5、num`den)氵subplot(311)rlocus(sys)subplot(312)b°de(Sys)subp1ot(313)nyquist(syS)9.2系统的传递函数为C$)=,绘制系统的单位阶跃响应和单位冲激响应曲线。程序如下:C△ear`c△c氵num=5*[1`5`6];den=[1`6`1O`8]氵step(num`den)`grid●∴氵figure`impulse(num`den)氵gridi=∴'9.3(改)典型二阶系统的传递函数为C(s)=,试绘制当c=⒍7,ω刀=6时系统的单位冲激响应和阶跃响应曲线。另求系统mp,屺“
6、,tp程序如下:C△ear氵c△c冫zeta=0.7氵wn=6氵num=wn^2`den=[1`2*Zeta*wn`wn^2]'sys=t£(num'den)`[mp`tr`ts`tp]=stepchar(sys`0.02)step(Sys)氵grid=Ⅱ氵figure氵impu△se(sys)氵grid=∷`9.4系统的特征方程为s6+2s5+Bs4+12s3+20s2+16s十16=0,试应用Routh稳定判据判定该系统的稳定性。程序如下:Clear氵clC`den=[1`2`8`△2`20`△6`16〗氵[rtab`msg]=routh(
7、den)a=rtab(:`1)ifa11(a)0)disp(’系坑廴称t的l)氵elsel)`disp(、系坨走下稳t呐end7(s+5)9.5系统的开环传递函数为G(s)=计算系统的幅值裕度和相角裕度。s2(s+10)(s+1)程序如下:C△ear`c△c`num≡7*[1、5]氵den=conv([1'0`0]`conv([1`10]`[1`1]))`sys=tf(num`den)氵[gm`pm`wg`wC]=margin(sys)例9.1.5二阶系统的传递函数为Gls·)=,当系统的输入信号是幅值为1,周期为8s的方旒波时,绘制系统的输出
8、响应曲线。程序如下:C△ear`clc`sys=tf([16]'[1'3'16])`∮[u't]=gens土g(lFqIⅠ91△`8`32`0.1)氵1Sim(sys'u't)`grid●9Κ例9.2.1负反馈系统的开环传递函数为G←)=绘制系统的根轨迹图,并求使系统稳定s(s2+4s+5)的Κ值范围和使系统无超调的K值范围。程序如下:C△ear`c△c氵sys=t£(1`[1`恐`5`0])r1ocus(sys)[x`y]=ginput(3)`p=x+i*yK=rlocfind(sys`p)1000例⒐3.1系统的开环传递函数为C$)〓,绘制
9、系统的Nyqu巛图,并讨论其稳定性。(s2+3s+2)(s+5)程序如下:C△ear`c1c`G=tf(1000'Conv([1`3、2]'[1`5]