《控制系统数字仿真与CAD》张晓华版课后答案.docxVIP

下载本文档

323
0
约4.65千字
约 6页
2019-02-26 发布于广东
举报
版权申诉

《控制系统数字仿真与CAD》张晓华版课后答案.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

2?2.用MATLAB语言求下列系统的状态方程状态方程模型参数:编写matlab程序如下 ?num=[l 7 24 24]; ?den=[l 10 35 50 24]; ?[AB C D]=tf2ss(num,den) 零极点增益:编写程序 ? num=[l 7 24 24]; ? den=[l 10 35 50 24]; ? fZPK]=tf2zp(num,den) 得到结果 Z=-2.7306 + 2.8531 ,-2.7306 -2.853li ,-1.5388 P=-4, -3 ,?2 ,-l K=1 ⑶部分分式形式:编写程序? num=[l 7 24 24]; ? den=[l 10 35 50 24]; ? [R PH]=residue(num,den) 得到结果 R= 4.0000 6.0000, 2.0000, 1.0000 P= -4.0000, -3.0000 ,-2.0000 ,-1.0000 H=[] (2)解：(1)传递函数模型参数：编写程序?A=[2.25 -5 -1.25 -0.5;2.25 -4.25 -1.25 -0.25;0.25 -0.5 -1.25 ?1;1.25 -1.75 -0.25 -0.75]； ? B=[4 2 2 0]; ? C=[0 2 0 2J; ? D=[0]; ? fnum den]=ss2tf(A,B,C,D) 得到结果 num = 0 4.0000 14.0000 22.0000 15.0000 den= 1.0000 4.0000 6.2500 5.2500 2.2500 (2)零极点增益模型参数： A=[2.25 -5 -1.25 -0.5: 2.25 ?4.25 -1.25 ?0.25 :0.25 -0.5 -1.25 1.25 ? 1.75 ?0.25 -0.75]； B=[4 2 2 0]; C=[0 2 0 2]; D=IOJ; [Z,P,K]=ss2zp(A,B,C,D) 部分分式形式的模型参数： A=[2.25 -5 -1.25 -0.5 ?4.25 -1.25 ?0.25 0.25 ?0.5-1.25 ?1 -1.75 -0.25 -0.75]； B=f4 2 2 0]： C=[0 2 0 2); D二[0]; [num denJ=ss2tf(A,B,C,D) [R,P,H]=residue(num,den) 2-6 (1)解：m文件程序为h=0.1; dispCa数的数值解为);％显示中间的文字％ dispCy=); %同上％ y 二 1; fort=0:h:l m=y; disp(y); %显示y的当前值％ y=m-m*h; end 保存文件da2.m 在matalb命令行屮键入 da 得到结果函数的数值解为 y= 1 0.9000 0.8100 0.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.3487 (2)另建一个m文件求解y二e在点[0,1]的数值程序为h=0.1; disp(屈数的离散时刻解为)； disp(yM); fort=0:h:l y=exp(-t); disp(y); end保存文件da3.m 在matalb命令行中键入? da3 函数的离散时刻解为 y= 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.4066 0.3679 2-5. 解： ⑴编辑m文件程序h=0.1; dispC四阶龙格■库塔方法求解两数数值解为 dispO*); y=i; fort=0:h:l disp(y); kl=-y; k2=-(y+kl *h/2); k3=-(y+k2*h/2); k4=-(y+k3*h); y=y+(k 1 +2*k2+2*k3+k4)*h/6; end保存文件q5.m 在matlab命令行里键入q5 得到结果四阶龙格■库塔方法求解函数数值解为 y= 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.4066 0.3679 2-7 ?a=ll 0]; ? b二[14.6]; ? c=[l 3.4 16.35]; ? d=conv(a,b)； ? e=conv(d,c) e= 1.0000 8.0000 31.9900 75.2100 0 ? f=[oo 0 5 100]; ? g=e+f g= 1.0000 8.0000 31.9900 80.2100 100.0000 %以上是计算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式％ ? p=roo