Matlab第十章作业.doc

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

Matlab第十章作业 10.4考虑一个单位负反馈控制系统,其前向通道传递函数为: 试应用Bode图法设计一个超前校正装置 使得校正后系统的相角裕度,幅值裕度,带宽。其中,。试问已校正系统的谐振峰值和谐振角频率的值各为多少? 解: 先建立超前校正函数fg_lead_pm(wc未知) 函数语句如下 function [ngc,dgc]=fg_lead_pm(ng0,dg0,Pm,w) [mu,pu]=bode(ng0,dg0,w); [gm,pm,wcg,wcp]=margin(mu,pu,w); alf=ceil(Pm-pm+5); phi=(alf)*pi/180; a=(1+sin(phi))/(1-sin(phi)) ; a1=1/a%求参数 dbmu=20*log10(mu); mm=-10*log10(a); wgc=spline(dbmu,w,mm); T=1/(wgc*sqrt(a)); ngc=[a*T,1]; dgc=[T,1]; 建立M文件l104其语句如下 ng0=[1];dg0=conv([1,0],conv([1,0],[1,5])); t=[0:0.01:5]; w=logspace(-3,2); g0=tf(ng0,dg0) b1=feedback(g0,1)%校正前系统闭环传函 [gm,pm,wcg,wcp]=margin(g0)%校正前参数 Pm=50; [ng1,dg1]=fg_lead_pm(ng0,dg0,Pm,w);%利用超前校正进行校正 g1=tf(ng1,dg1);%校正环节传递函数 g2=g0*g1;%校正后前向通道传函 [gm1,pm1,wcg1,wcp1]=margin(g2)%校正后系统参数 km1=20*log(gm1)%校正后系统幅值裕度工程表示 b2=feedback(g2,1);%校正后系统闭环传函 bode(b2) [mag,phase,w]=bode(b2);%对校正后系统闭环传函bode图进行离散化 figure, bode(g0,r--,g1,b--,g2,g,w), grid on ; Mr=max(mag)%求取Mr a=find(mag==Mr);%求Mr对应的脚标 wr=w(a)%求wr b=find(mag=0.707*mag(1));%求幅值小于0.707倍零频幅值的脚标所组成的数组 wb=w(b(1))%求wb 在命令窗口中得到以下结果 Transfer function: 1 ----------- s^3 + 5 s^2 Transfer function: 1 --------------- s^3 + 5 s^2 + 1 gm = 0 pm = -5.1009 wcg = 0 wcp = 0.4463 a1 = 0.0669 gm1 = 8.3933 pm1 = 51.0980 wcg1 = 3.8710 wcp1 = 0.8728 km1 = 42.5488 Mr = 1.2644 wr = 0.4695 wb = 1.7055 可得各项指标=0.0669,, , 10.5 考虑一个单位负反馈控制系统,其前向通道传递函数为: 试应用Bode图法设计一个校正装置,使得校正后系统的静态速度误差常数,相角裕度,幅值裕度。 解: 1 .先建立超前校正函数fg_lead_pm(wc未知) 函数语句如下 function [ngc,dgc]=fg_lead_pm(ng0,dg0,Pm,w) [mu,pu]=bode(ng0,dg0,w); [gm,pm,wcg,wcp]=margin(mu,pu,w); alf=ceil(Pm-pm+5); phi=(alf)*pi/180; a=(1+sin(phi))/(1-sin(phi)) ; a1=1/a%求参数 dbmu=20*log10(mu); mm=-10*log10(a); wgc=spline(dbmu,w,mm);

文档评论(0)

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

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

1亿VIP精品文档

相关文档