MATLAB控制系统各种仿真例题(包括simulink解法)资料.doc

MATLAB控制系统各种仿真例题(包括simulink解法)资料.doc

  1. 1、本文档共21页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB控制系统各种仿真例题(包括simulink解法)资料.doc

控制系统的模型与转换 请将下面的传递函数模型输入到matlab环境。 ,T=0.1s s=tf(s); G=(s^3+4*s+2)/(s^3*(s^2+2)*((s^2+1)^3+2*s+5)); G Transfer function: s^3 + 4 s + 2 ------------------------------------------------------ s^11 + 5 s^9 + 9 s^7 + 2 s^6 + 12 s^5 + 4 s^4 + 12 s^3 num=[1 0 0.56]; den=conv([1 -1],[1 -0.2 0.99]); H=tf(num,den,Ts,0.1) Transfer function: z^2 + 0.56 ----------------------------- z^3 - 1.2 z^2 + 1.19 z - 0.99 请将下面的零极点模型输入到matlab环境。请求出上述模型的零极点,并绘制其位置。 ,T=0.05s z=[-1-j -1+j]; p=[0 0 -5 -6 -j j]; G=zpk(z,p,8) Zero/pole/gain: 8 (s^2 + 2s + 2) -------------------------- s^2 (s+5) (s+6) (s^2 + 1) pzmap(G) z=[0 0 0 0 0 -1/3.2 -1/2.6]; p=[1/8.2]; H=zpk(z,p,1,Ts,0.05) Zero/pole/gain: z^5 (z+0.3125) (z+0.3846) ------------------------- (z-0.122) Sampling time: 0.05 pzmap(H) num=[0,7.1570,-6.4875 ]; den=[1,-2.2326,1.7641,-0.4966]; sysd=tf(num,den,0.05,variable,z^-1) Transfer function: 7.157 z^-1 - 6.487 z^-2 ----------------------------------------- 1 - 2.233 z^-1 + 1.764 z^-2 - 0.4966 z^-3 Sampling time: 0.05 线性系统分析 请分析下面传递函数模型的稳定性。 num=[1]; den=[1 2 1 2]; G=tf(num,den); eig(G) ans = -2.0000 0.0000 - 1.0000i 0.0000 + 1.0000i 可见,系统有两个特征根在虚轴上,一个特征根在虚轴左侧,所以系统是临界稳定的。 num=[3 1]; den=[300 600 50 3 1]; G=tf(num,den); eig(G) ans = -1.9152 -0.1414 0.0283 - 0.1073i 0.0283 + 0.1073i 可见,有两个特征根在虚轴右侧,所以系统是不稳定的。 请判定下面离散系统的稳定性。 num=[-3 2]; den=[1 -0.2 -0.25 0.05]; H=tf(num,den,Ts,0.1); [eig(H) abs(eig(H))] ans = -0.5000 0.5000 0.5000 0.5000 0.2000 0.2000 可以看出,由于各个特征根的模均小于1,所以可以判定闭环系统是稳定的。 z=tf(z,0.1); H=(2.12*z^-2+11.76*z^-1+15.91)/…; (z^-5-7.368*z^-4-20.15*z^-3+102.4*z^-2+80.39*z-1-340); [eig(H) abs(eig(H))] ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0

文档评论(0)

四娘 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档