动力工程计算机控制matlab上机作业.docVIP

下载本文档

7
0
约2.63千字
约 4页
2017-04-24 发布于重庆
举报
版权申诉

动力工程计算机控制matlab上机作业.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文档。上传文档

动力工程计算机控制matlab上机作业

2-----12题已知系统如图所示。已知控制器为，对象特性为，采用零阶保持器，单位阶跃输入，采样周期1秒。 (1) 求出系统的闭环Z传递函数Gc(z)、误差Z传递函数E(z)、控制命令Z传递函数U(z)和输出的Z变换Y(z)； (2) 编程上机，算出y(kT)（k = 0, 1, ···, 20），画出响应曲线，求出超调量?p、过渡时间ts； Matlab程序如下： Ts=1;z=0;p=1;k=1; Rz=zpk(z,p,k,Ts) %按零点、极点、比例系数生成R(z)表达式（注意含采样周期Ts） z=[]; p=[0 -1]; k=10; Gs=zpk(z,p,k); %按零点、极点、比例系数生成G(s)表达式 Gd=c2d(Gs,Ts,zoh) %将G(s)按H0(s)离散化 z=[0.5 0.368]; p=[1 -0.718]; k=0.543; Dd=zpk(z,p,k,Ts) %按零点、极点、比例系数生成D(z)表达式（注意含采样周期Ts） G0d=Dd*Gd Gcd=feedback(G0d,1) %按单位负反馈生成闭环传递函数Gc(z)表达式 Yz=Rz*Gcd [yk t]=step(Gcd,0:1:20) %求Gc(z)对阶跃输入的响应，即y(kT) yk = 0 1.9976 1.0039 0.9978 1.0011 0.9992 1.0006 0.9996 1.0003 0.9998 1.0002 0.9999 1.0001 0.9999 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 t = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 pp=spline(t,yk); yy=ppval(pp,linspace(0,20)); plot(linspace(0,20),yy,b); hold on xlabel(t);ylabel(y); Ez=1/(1+G0d) [ek t]=step(Ez,0:1:20) %求E(z)对阶跃输入的响应 ek = 1.0000 -0.9976 -0.0039 0.0022 -0.0011 0.0008 -0.0006 0.0004 -0.0003 0.0002 -0.0002 0.0001 -0.0001 0.0001 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 t = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Uz=Dd/(1+G0d) [uk t]=step(Uz,0:1:20) %求U(z)对阶跃输入的响应 uk = 0.5430 -0.8599 0.7153 -0.5123 0.3671 -0.2631 0.1886 -0.1351 0.0968 -0.0694 0.0497 -0.0356 0.0255 -0.0183 0.0131 -0.0094 0.0067 -0.0048 0.0035 -0.0025 0.0018 t = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 pp=spline(t,uk); yy=ppval(pp,linspace(0,20)); plot(linspace(0,20),yy,b); hold on Ez=1/(1+G0d [ek t]=step(Ez,0:1:20) %求E(z)对阶跃输入的响应 ek = 1.0000 -0.9976 -0.0039 0.0022 -0.0011 0.0008