原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
计算机仿真课程方案dongge
目录
一、课程设计的目的…………………………………………………………1
二、课程设计用的仪器和器件……………………………………………1
三、课程设计的具体内容………………………………………………… 1
1、离散系统………………………………………………………………… 1
①离散有波纹……………………………………………………………… 1
②离散无波纹……………………………………………………………… 2
2、连续系统 ………………………………………………………………… 4
四、课程设计的体会和发生的问题…………………………………… 6
五、致谢……………………………………………………………………… 6
六、参考文献………………………………………………………………… 7
一、课程设计的目的
本次计算机仿真课程设计的目的是熟悉并使用MATLAB在计算机中实现连续系统和离散系统的设计和仿真。
二、课程设计用的仪器和器件
硬件： PC计算机
软件： MATLAB7.0
三、课程设计的具体内容
1.离散系统：
离散控制系统建模、分析、设计和仿真[9号]
设被控对象的传递函数用零阶保持器离散化，采样周期取0.08秒，以MATLAB为工具，分别设计一单位速度信号输入时的最少拍有波纹控制器DY(z)和一阶跃信号输入时的最少拍无波纹控制器DW(z) 。要求通过仿真证实设计结果。
①离散有波纹
M文件：
sys=tf([999 8991 17982],[1 14 56 64 0 0])
gz=c2d(sys,0.08,zoh)
[a b k]=zpkdata(gz)
gz=zpk(a,b,k,0.08,variable,z^-1)
syms z h0 h1 g0 g1
gz=0.077519*z^-1*(1+3.398*z^-1)*(1-0.7866*z^-1)*(1-0.6188*z^-1)*(1+0.2411*z^-1)/(1-z^-1)^2/(1-0.8521*z^-1)/(1-0.7261*z^-1)/(1-0.5273*z^-1)
hz=z^-1*(1+3.398*z^-1)*(h0+h1*z^-1)
he=(1-z^-1)^2*(g0+g1*z^-1)
f1=subs(hz,z,1)-1
f2=subs(diff(hz),z,1)
[h0 h1]=solve(f1,f2)
h=[h0;h1]
hz=z^-1*(1+3.398*z^-1)*(h0+h1*z^-1)
f3=1-limit(he,z,Inf)
f4=1-subs(he,z,-3.398)
[g0 g1]=solve(f3,f4)
g=[g0;g1]
he=(1-z^-1)^2*(g0+g1*z^-1)
dz=hz/gz/he
sysy=zpk([-3.398 0.6394],[0 0 0],[0.6304],0.08)
sysyc=d2c(sysy,tustin)
hold on
td=0:0.08:2
ud=td
lsim(sysy,ud,td)
t=0:0.008:2
u=t
lsim(sysyc,u,t)
截图：
②离散无波纹：
M文件：
sys=tf([999 8991 17982],[1 14 56 64 0 0])
gz=c2d(sys,0.08,zoh)
[a b k]=zpkdata(gz)
gz=zpk(a,b,k,0.08,variable,z^-1)
syms z h0 h1 g0 g1
gz=0.077519*z^-1*(1+3.398*z^-1)*(1-0.7866*z^-1)*(1-0.6188*z^-1)*(1+0.2411*z^-1)/(1-z^-1)^2/(1-0.8521*z^-1)/(1-0.7261*z^-1)/(1-0.5273*z^-1)
hz=z^-1*(1+3.398*z^-1)*(1-0.7866*z^-1)*(1-0.6188*z^-1)*(1+0.2411*z^-1)*(h0+h1*z^-1)
syms g2 g3 g4
he=(1-z^-1)^2*(g0+g1*z^-1+g2*z^-2+g3*z^-3+g4*z^-4)
f1=subs(hz,z,1)-1
f2=subs(diff(hz),z,1)
[h0 h1]=solve(f1,f2)
h=[h0;h1]
hz=z^-1*(1+3.398*z^-1)*(1-0.7866*z^-1)*(1-0.6188*z^-1)*(1+0.2411*z^-1)*(h0+h1*z^-1)
f3=1-limit(he,z,Inf)
f4=1-subs(he,z,-3.398)
f5=1-subs(he,z,0.7866)
f6=1-subs(he,z,0.6188)
文档评论（0）