自动控制原理MATLAB课程设计--滞后-超前校正
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
滞后-超前校正
——课程设计
一、设计目的:
1. 了解控制系统设计的一般方法、步骤。
2. 掌握对系统进行稳定性的分析、稳态误差分析以及动态特性分析的方
法。
3. 掌握利用MATLAB 对控制理论内容进行分析和研究的技能。
4. 提高分析问题解决问题的能力。
二、设计内容与要求:
设计内容:
1. 阅读有关资料。
2. 对系统进行稳定性分析、稳态误差分析以及动态特性分析。
3. 绘制根轨迹图、Bode 图、Nyquist 图。
4. 设计校正系统,满足工作要求。
设计条件:
1、被控制对象的传递函数是
m m 1m 2012m
n sn 1n 2012n
b s b s b s b ()a s a a s a G S ----+++⋯+=+++⋯+(n≥m)
2、参数a0,a1,a2,...an和b0,b1,b2,...bm因小组而异。
设计要求:
1. 能用MATLAB 解复杂的自动控制理论题目。
2. 能用MATLAB 设计控制系统以满足具体的性能指标。
3. 能灵活应用MATLAB 的CONTROL SYSTEM 工具箱和SIMULINK 仿真软件,
分析系统的性能。
三、设计步骤:
1、自学MATLAB软件的基本知识,包括MATLAB的基本操作命令。控制系统工具箱的用法等,并上机实验。
2、基于MALAB用频率法对系统进行串联校正设计,使其满足给定的领域
性能指标。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数,校正装置的参数T,α等的值。
已知开环传递函数为G(S)= 0
(2)(40)k s s s ++,使用频率法设计串联滞后—超前校
正装置,使系统的相角裕度大于等于40°,静态速度误差系数等于20。 校正前
根据上式可化简G(S)= 0
0.0125(0.51)(0.0251)
k s s s ++,所以公式G(S)=
20
(0.51)(0.0251)s s s ++,
所以=1,则c w = 6.1310,
相角裕度γ为9.3528。 校正后
串联校正滞后—超前校正装置2121(1)(1)
()(1)(1)
C bT S aT S G S T S T S ++=
++, 由
111(~)510
C W bT =,取
2211
*60.8335
bT bT =⇒=,又由20lg 2020lg520lg 4
a a -=⇒=
,
121aT =,10.125T =;10.25ab b =⇒=;1 1.22 3.332
T bT =⇒= 所以校正后的函数为 (10.833)(10.5)20
()(10.125)(1 3.33)(0.51)(0.0251)C S S G S S S s s s ++=
++++
3、利用MATLAB函数求出校正前与校正后系统的特征根,并判断其系统是否稳定,为什么? 校正前
>> % MATLAB PROGRAM j005.m >> %
>> num=[20];
den= conv([1 0],conv([0.5 1],[0.025 1])); g=tf(num,den);
sys=feedback(g,1); pzmap(g);
den= conv([1 0],conv([0.5 1],[0.025 1])); t=tf(num,den); pzmap(t);
[p,z]=pzmap(g); den=sys.den{1}; r=roots(den);
disp(r)
-41.0006
-0.4997 + 6.2269i
-0.4997 - 6.2269i
系统没有零极点在右边,所以系统开环稳定。
校正后
>> % MATLAB PROGRAM j005.m
%
num=conv([0.833 1],[6.6 20]);
d3=conv([0.5 1],conv([0.025 1],[1 0]));d2=conv([0.125 1],[3.33 1]);den=conv([0 d3],[0 d2]);
g=tf(num,den);
sys=feedback(g,1);
pzmap(g);
conv([0.5 1],conv([0.025 1],[1 0]));d2=conv([0.125 1],[4 1]);den=conv([0 d3],[0 d2]);
t=tf(num,den);
pzmap(t);
[p,z]=pzmap(g);
den=sys.den{1}; r=roots(den); disp(r)
-40.7534 -4.8304 -1.6921 + 3.4339i -1.6921 - 3.4339i -1.3324
极点和特征根都在左半平面,故系统稳定
4、利用MATLAB作出系统校正前与校正后的单位脉冲响应曲线,单位阶跃响应曲线,单位斜坡响应曲线,分析这三种曲线的关系?求出系统校正前与校正后的动态性能指标Ϭ%,tr,tp,ts,e ss ,并分析其有何变化?
脉冲响应
校正前
>> % MATLAB PROGRAM j005.m >> %
>> k=20;n1=1;d1=conv([1 0],conv([0.5 1],[0.025 1])); >> s1=tf(k*n1,d1);