《PX单位负反馈系统的校正装置设计》由会员分享,可在线阅读,更多相关《PX单位负反馈系统的校正装置设计(20页珍藏版)》请在金锄头文库上搜索。
1、武汉理工大学自动控制原理课程设计说明书附件 1:学 号: 0121111360126课 程 设 计题 目 单位负反馈系统的校正装置设计学 院 自动化学院专 业 自动化班 级 自动化 1001 班姓 名 彭玄指导教师 谭思云2013 年 12 月 24 日武汉理工大学自动控制原理课程设计说明书课程设计任务书学生姓名: 彭玄 专业班级: 自动化 11 指导教师: 谭思云 工作单位: 自动化学院 题 目: 单位负反馈系统的校正装置设计 初始条件:已知某控制系统结构如图所示,要求设计校正环节 Gc(s ) ,使系统对于阶跃输入的稳态误差为 0,使系统校正后的相角裕量 ,幅值裕量 .045dBh10R(
2、t)-)1.0)(ss要求完成的主要任务: (包括课程设计工作量及其技术要求,以及说明书撰写等具体要求)(1) 用 MATLAB 作出原系统的系统伯德图,计算系统的幅值裕度和相度。(2) 在系统前向通路中插入一校正装置,确定校正网络的传递函数,用MATLAB 进行验证。(3) 用 MATLAB 画出未校正和已校正系统的根轨迹。(4) 对上述任务写出完整的课程设计说明书,说明书中必须进行原理析,写清楚分析计算的过程及其比较分析的结果,并包含 Matlab 源程或 Simulink 仿真模型,说明书的格式按照教务处标准书写。时间安排:1、 课程设计任务书的布置,讲解 (半天)2、 根据任务书的要求
3、进行设计构思。 (半天)3、 熟悉 MATLAB 中的相关工具(一天)4、 系统设计与仿真分析。 (三天)5、 撰写说明书。 (二天)6、 课程设计答辩(半天)指导教师签名: 年 月 日系主任(或责任教师)签名: 年 月 武汉理工大学自动控制原理课程设计说明书日武汉理工大学自动控制原理课程设计说明书目录摘要 .11设计题目 .22设计要求 .23设计的总体思路 .34用 MATLAB 作出原系统的系统伯德图和根轨迹 .35 校正过程 .65.1 超前校正过程 .65.2滞后校正过程 .86simulink 仿真 .137、设计总结 .15参考文献 .16武汉理工大学自动控制原理课程设计说明书0
4、摘要在科技高速发展的今天,自动控制技术已广泛运用于制造,农业,交通,航空航天等众多部门,极大的提高了社会劳动生产率,改善了人们的劳动环境,丰富和提高了人民的生活水平,在今天的社会生活中,自动化装置已经无所不在,为人类文明的进步作出了重要贡献,自动控制系统的课程设计就是检验我们学过知识扎实程度的好机会,也让我们的知识体系更加系统,更加完善。在不断学习新知识的基础上得到动手能力的训练,启发新思维及独立解决问题的能力,提高设计、装配、调试能力。如果系统设计要求满足性能指标属频域特征量,则通过采用频域校正方法。在开环系统对数频率特性基础上,满足稳态误差、开环系统截止频率和相角裕度等要求出发,进行串联校
5、正的方法。在 bode 图上虽然不能严格定量的给出系统的动态性能,但是却能方便的根据频域指标确定校正装置的形式和参数,特别是对校正系统的高频特性有要求时,采用频率校正较其他方法更方便。串联滞后校正-超前校正,其基本原理是利用滞后超前网络的超前部分来郑大系统的相角裕度,同时利用滞后部分来改善系统的稳定性能。关键词:控制系统,校正,传递函数,频率,幅值相武汉理工大学自动控制原理课程设计说明书1单位负反馈系统的校正装置设计1设计题目已知某控制系统结构如图所示,要求设计校正环节 Gc(s) ,使系统对于阶跃输入的稳态误差为 0,使系统校正后的相角裕量 ,幅值裕量 。045dBh10R(t)-)1.0)
6、(ss2设计要求1.用 MATLAB 作出原系统的系统伯德图,计算系统的幅值裕度和相位裕度。2.在系统前向通路中插入一校正装置,确定校正网络的传递函数,并用MATLAB 进行验证。3.用 MATLAB 画出未校正和已校正系统的根轨迹。4.对上述任务写出完整的课程设计说明书,说明书中必须进行原理分析,写清楚分析计算的过程及其比较分析的结果,并包含 Matlab 源程序或Simulink 仿真模型。武汉理工大学自动控制原理课程设计说明书23设计的总体思路(1)根据开环传递函数得 E(s )= ,ess = )1(GsH)(lim0sEs得 ess=0。(2)根据开环传递函数画出 bode 图和根轨
7、迹图,求出幅值裕度和相位裕度。从图上看出校正前系统的相位裕度 和剪切频率c(3)根据相位裕度 的要求,计算出滞后校正装置的参数 a 和 T。即得校正装置的传递函数,然后得到校正后系统的开环传递函数。(4)验证已校正系统的相位裕度 和幅值裕度 h。4用 MATLAB作出原系统的系统伯德图和根轨迹根据开环传递函数画出 bode 图以及根轨迹,得出幅值裕度和相位裕度。=)(sG)1.0)(1s校正前的 bode 图,如图 1 所示源程序:num=100;den=1,11,10,0;margin(num,den);grid; 图 1 校正前的 bode 图武汉理工大学自动控制原理课程设计说明书3由图
8、1 得出:幅值裕度为 0.828db相位裕度为 1.58deg幅值截止频率 Wcg=3.16rad/s相位截止频率 Wcp=3.01rad/s校正前的根轨迹图如图 2 所示。源程序:num = 100;den = 1 11 10 0;rlocus(num,den) -35 -30 -25 -20 -15 -10 -5 0 5 10 15-25-20-15-10-50510152025 Root LocusReal Axis (seconds-1)Imaginary Axis (seconds-1)图 2 校正前的根轨迹图单位阶跃响应分析如图 3 所示。源程序:G=tf(100,1 11 10
9、0 );G1=feedback(G,1);武汉理工大学自动控制原理课程设计说明书4t=0:0.01:100;step(G1,t);gridxlabel(t );ylabel(c(t)0 10 20 30 40 50 60 70 80 90 10000.20.40.60.811.21.41.61.82 Step ResponseTime (seconds)Amplitude图 3 校正前单位阶跃响应图武汉理工大学自动控制原理课程设计说明书55校正过程5.1 超前校正过程(1). 得 =18090actgWc=17.5611m(2).则有 =4517.56+10=37.33得出 a=0.245(3
10、).然后未校正系统的 特性曲线上查出其幅值等于 10lg(1/a)对应0()L的频率 。如图 4。 m源程序:num=100;den=1,11,10,0;margin(num,den);10-2 10-1 100 101 102 103-270-225-180-135-90Phase (deg)Bode DiagramGm = 0.828 dB (at 3.16 rad/s) , Pm = 1.58 deg (at 3.01 rad/s)Frequency (rad/s)-150-100-50050100System: untitled1Frequency (rad/s): 2.09Magni
11、tude (dB): 6.11Magnitude (dB)图 4 超前校正的 bode 图(4).这就是校正后系统的截止频率 =2.09,且 。cmc(5).确定校正网络的传递函数。根据步骤 3 所求得的 和 a 两值,可求时间常数 T。得到 T=0.9671mTa武汉理工大学自动控制原理课程设计说明书6(6).校正系统的传递函数为 校正系统的 bode 图,如图 5源程序:G0=tf(0.967 1, 0.237 1);margin(G0); 051015Magnitude (dB)10-2 10-1 100 101 102010203040Phase (deg)Bode DiagramGm
12、 = Inf , Pm = -180 deg (at 0 rad/s)Frequency (rad/s)图 5 校正系统的 bode 图校正后系统 bode 如图 6源程序:n1=100;d1=1 11 10 0s1=tf(n1,d1);s2=tf(0.967 1, 0.237 1);s=s1*s2;Gm,Pm,Wcm,Wcp=margin(s)margin(s)武汉理工大学自动控制原理课程设计说明书7-150-100-50050Magnitude (dB)10-1 100 101 102 103-270-225-180-135-90Phase (deg)Bode DiagramGm = 3.
13、24 dB (at 6.46 rad/s) , Pm = 10.1 deg (at 5.31 rad/s)Frequency (rad/s)图 6 校正后系统的 bode 图由图得到的相位裕度跟幅值裕度都不满足要求,超前校正无法满足要求。5.2滞后校正过程(1) 根据相角裕量 45 0的要求,再考虑到串接滞后校正装置的相角滞后,从未校正系统的频率特性曲线图 1 上,找出对应相角180 0+(45 0+100)=125 0处的频率 wc0.612rad/s。w c将作为校正后系统的增益交界频率。(2) 确定滞后装置的传递函数为 : G c=(1+aTs)/(1+Ts)根据滞后校正装置的最大幅值和原系统在 wc上的幅值相等条件,求出 a 值。在 wc=wc处,从未校正的对数幅频特性曲线上得 :20lg|G 0(jwc)|=22.9dB再由 20lg 1/a=22.9dB 计算出 a=0.0716由 1/(aT)=1/10w c所以当 wc 0.612rad/s a=0.0716 时,可求得:T=228.21s将所求的 a 值和 T 值代入式得校