控制系统计算机仿真上机报告.docVIP

下载本文档

52
0
约8.95千字
约 16页
2019-08-14 发布于江西
举报
版权申诉

控制系统计算机仿真上机报告.doc

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语言求下列系统的状态方程、传递函数、零极点增益和部分分式形式的模型参数，并分别写出其相应的数学模型表达式：（1）（2）解：（1）源程序： num=[1 7 24 24]; den=[1 10 35 50 24]; G=tf(num,den) [Z,P,K]=tf2zp(num,den) [A,B,C,D]=tf2ss(num,den) [R,P,H]=residue(num,den) 结果： Transfer function: s^3 + 7 s^2 + 24 s + 24 --------------------------------- s^4 + 10 s^3 + 35 s^2 + 50 s + 24 Z = -2.7306 + 2.8531i -2.7306 - 2.8531i -1.5388 P = -4.0000 -3.0000 -2.0000 -1.0000 K = 1 A = -10 -35 -50 -24 1 0 0 0 0 1 0 0 0 0 1 0 B = 1 0 0 0 C = 1 7 24 24 D = 0 R = 4.0000 -6.0000 2.0000 1.0000 P = -4.0000 -3.0000 -2.0000 -1.0000 H = [] （2）源程序： 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 2] D=[0] [num,den]=ss2tf(A,B,C,D) [Z,P,K]=tf2zp(num,den) [R,P,H]=residue(num,den) 结果： A = 2.2500 -5.0000 -1.2500 -0.5000 2.2500 -4.2500 -1.2500 -0.2500 0.2500 -0.5000 -1.2500 -1.0000 1.2500 -1.7500 -0.2500 -0.7500 B = 4 2 2 0 C = 0 2 0 2 D = 0 num = 0 4.0000 14.0000 22.0000 15.0000 den = 1.0000 4.0000 6.2500 5.2500 2.2500 Z = -1.0000 + 1.2247i -1.0000 - 1.2247i -1.5000 P = -1.5000 -1.5000 -0.5000 + 0.8660i -0.5000 - 0.8660i K = 4.0000 R = 4.0000 0.0000 0.0000 - 2.3094i 0.0000 + 2.3094i P = -1.5000 -1.5000 -0.5000 + 0.8660i -0.5000 - 0.8660i H = [] 2-3 用殴拉法求下列系统的输出响应在上，时的数值解。，要求保留4位小数，并将结果以图形的方式与真解比较。解：源程序： t=0:0.1:1; h=0.1; y(1)=1; for i=1:1:10; y(i+1)=y(i)+h*(-2*y(i));