控制系统校正课程设计.doc

  1. 1、本文档共18页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
控制系统校正课程设计 PAGE 1 - 一、设计目的: 、了解控制系统设计的一般方法、步骤。 2、掌握对系统进行稳定性分析、稳态误差分析以及动态特性分析的方法。 3、掌握利用MATLAB对控制理论内容进行分析和研究的技能。 4、提高分析问题解决问题的能力。 二、设计内容与要求: 设计内容: 1、阅读有关资料。 2、对系统进行稳定性分析、稳态误差分析以及动态特性分析。 3、绘制根轨迹图、Bode图、Nyquist图。 4、设计校正系统,满足工作要求。 设计条件: 已知单位负反馈系统的开环传递函数: , 试用频率法设计串联滞后校正装置,使系统的相位裕度为,增益裕度不低于10dB,静态速度误差系数,剪切频率不低于1rad/s。 设计要求: 能用MATLAB解复杂的自动控制理论题目。 能用MATLAB设计控制系统以满足具体的性能指标。 能灵活应用MATLAB的CONTROL SYSTEM 工具箱和SIMULINK仿真软件,分析系统的性能。 三、设计步骤: 1、自学MATLAB软件的基本知识。包括MATLAB的基本操作命令、控制系统工具箱的用法等,并上机实验。 2、基于MATLAB用频率法对系统进行串联校正设计,使其满足给定的频域性能指标。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数,校正装置的参数,等的值。 解: 求满足稳态误差要求的系统开环增益. 即被控对象的传递函数为: = = ⑵.利用已经确定的开环增益k,画出未校正的系统的Bode图,求出 相位裕度和幅值裕度。 在MATLAB中输入: G=tf(7,[0.083 0.667 1 0]);[kg,r]=margin(G) Margin(G) 该程序执行后,得系统的Bode图及性能指标 kg = 1.1480 r = 3.4614 即模稳定裕度:= 穿越频率:= 相稳定裕度: 剪切频率:= ⑶.根据题意,选择校正装置:采用串联滞后校正。 截止频率 相角滞后通常取(-6~-14) 取-6 =180+= =- 得=1.6 rad/s ⑷.滞后校正装置的传递函数为:G(S)= , =(1/5~1/10)* -20lgb=20lg7/1.92 得b=0.23 取T=27 所以滞后校正传递函数为 校正后开环传递函数为* ⑸.检验: 在MATLAB中输入: G=tf(7*[6.21 1],conv([0.083,0.667,1,0],[27 1])) [Gm,Pm,Wcp,Wcg]=margin(G); margin(G) Transfer function: 43.47 s + 7 2.241 s^4 + 18.09 s^3 + 27.67 s^2 + s Gm =4.5769 Pm =38.811 Wcp =3.3245 Wcg = 1.3206 即,模稳定裕度:= 穿越频率: 相稳定裕度: 剪切频率: = 由以上数据可以看出,在串联了一个滞后校正装置后系统稳定, ,满足,增益裕度不低于10dB。满足要求。 3、利用MATLAB函数求系统的特征根,并判断系统的稳定性,为什么? 校正前: 程序如下: num=[7]; den=[0.083 0.667 1 0]; g=tf(num,den); sys=feedback(g,1); pzmap(g); den=[0.083 0.667 1 7]; t=tf(num,den); pzmap(t); [p,z]=pzmap(g); den=sys.den{1}; r=roots(den); disp(r) -7.8673 -0.0844 + 3.2730i -0.0844 - 3.2730i 开环零极点分布图 征根中无实部为正的根,所以闭环系统是稳定的。 校正后: 程序如下所示: num=7; den=conv([0.5 1 0],[0.167 1]); Gs=tf(num,den); n1=[6.21 1]; d1=[27 1]; Gc=tf(n1,d1); sope=Gs*Gc; Go=feedback(sope,1,-1) eig(Go) Transfer function: 43.47 s + 7 2.255 s^4 + 18.09 s^3 + 27.67 s^2 + 44.47 s + 7 ans = -6.6091 -0.6209 + 1.5204i -0.6209 - 1.5204i -0

文档评论(0)

153****9595 + 关注
实名认证
内容提供者

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

1亿VIP精品文档

相关文档