自动控制原理校正课程设计4.docx
《自动控制原理校正课程设计4.docx》由会员分享,可在线阅读,更多相关《自动控制原理校正课程设计4.docx(21页珍藏版)》请在冰点文库上搜索。
自动控制原理校正课程设计4
自动控制原理校正课程设计4
课程设计报告
课程名称自动控制原理
院部名称机电工程学院
专业电气工程及其自动化
班级
指导教师
金陵科技学院教务处制
设计任务...............................................3
设计要求...............................................3
设计步骤..............................................3
未校正前系统的性能分析....................................3
1.1开环增益
...........................................3
1.2校正前系统的各种波形图...............................4
1.3由图可知校正前系统的频域性能指标....................7
1.4特征根..............................................8
1.5判断系统稳定性......................................8
1.6分析三种曲线的关系.................................8
1.7求出系统校正前动态性能指标及稳态误差................8
1.8绘制系统校正前的根轨迹图...........................9
1.9绘制系统校正前的Nyquist图.........................9
校正后的系统的性能分析.........................................10
2.1滞后超前校正.......................................10
2.2校正后系统的各种波形图............................11
2.3由图可知校正后系统的频域性能指标..................14
2.4特征根..............................................14
2.5判断系统稳定性.....................................14
2.6分析三种曲线的关系................................15
2.7求出系统校正后动态性能指标及稳态误差..............15
2.8绘制系统校正后的根轨迹图和Nyquist图..............15
心得体会..............................................18
主要参考文献...........................................19
一、设计任务:
已知单位负反馈系统的开环传递函数
,试用频率法设计串联滞后——超前校正装置。
(1)使系统的相位裕度
(2)静态速度误差系数
(3)幅值穿越频率
二、设计要求:
(1)首先,根据给定的性能指标选择合适的校正方式对原系统进行校正,使其满足工作要求。
要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数,校正装置的参数T,
等的值。
(2)利用MATLAB函数求出校正前与校正后系统的特征根,并判断其系统是否稳定,为什么?
(3)利用MATLAB作出系统校正前与校正后的单位脉冲响应曲线,单位阶跃响应曲线,单位斜坡响应曲线,分析这三种曲线的关系?
求出系统校正前与校正后的动态性能指标σ%、tr、tp、ts以及稳态误差的值,并分析其有何变化?
(4)绘制系统校正前与校正后的根轨迹图,并求其分离点、汇合点及与虚轴交点的坐标和相应点的增益
值,得出系统稳定时增益
的变化范围。
绘制系统校正前与校正后的Nyquist图,判断系统的稳定性,并说明理由?
(5)绘制系统校正前与校正后的Bode图,计算系统的幅值裕量,相位裕量,幅值穿越频率和相位穿越频率。
判断系统的稳定性,并说明理由?
三、设计步骤:
开环传递函数
1、未校正前系统的性能分析
1.1开环增益
已知系统中只有一个积分环节,所以属于I型系统
由静态速度误差系数
可选取
开环传递函数为
1.2通过matlab求未校正装置的传递函数:
程序代码如下:
>>num=260;
den=conv([10],conv([0.11],[0.011]));
G=tf(num,den);
[gm1,pm1,wcg1,wcp1]=margin(G)
margin(G)
Warning:
Theclosed-loopsystemisunstable.
>Inlti.marginat89
gm1=
0.4231
pm1=-13.8089
wcg1=31.6228
wcp1=47.9095
>>
1.3输出校正前系统的幅值裕度、相角裕度、相角穿越频率、幅值穿越频率:
Kg=0.4231、r=-13.8089、Wg=31.6228;Wc=47.9095.
k=260;n1=1;d1=conv(conv([10],[0.11]),[0.011]);s1=tf(k*n1,d1);
figure
(1);sys=feedback(s1,1);step(sys)
figure
(2);margin(s1);holdon
>>figure(3);sys=feedback(s1,1);impulse(sys)
>>figure(4);step(k*n1,[d1,0])
图1-1校正前系统的bode图
1-2校正前系统的单位阶跃响应
图1-3校正前系统的单位脉冲响应
图1-4校正前系统的单位斜坡响应
1.4特征根:
>>s=solve('(s*(0.1*s+1)*(0.01*s+1)+260=0)')
s=
-119.77456778288007033207608505346
4.8872838914400351660380425267284-46.334210769402293735695159755370*i
4.8872838914400351660380425267284+46.334210769402293735695159755370*i
1.5判断系统稳定性:
(1)由图1可以看出,
之前对数幅频渐近特性曲线所对应的相频特性曲线穿越了
由程序输出结果可知系统校正前有三个特征根,且有两个特征根的实部为正值。
系统稳定的充分必要条件为系统特征方程的所有根都具有负实部,或者说都位于s平面的虚轴之左。
所有校正前的系统是不稳定的。
1.6分析校正前单位脉冲响应曲线,单位阶跃响应曲线,单位斜坡响应曲线,这三种曲线的关系
单位斜坡响应的一次导数是阶跃响应曲线,阶跃响应的一次导数是冲击响应。
1.7求出系统校正前动态性能指标σ%、tr、tp、ts以及稳态误差的值:
ess=0;tr=0.0188;tp=0.4806;ts=0.4429.1.8绘制系统校正前的根轨迹图,
k=260;n1=1;d1=conv(conv([10],[0.11]),[0.011]);s1=tf(k*n1,d1);
figure
(1);sys=feedback(s1,1);step(sys)
c=dcgain(sys);[y,t]=step(sys);[max_y,k]=max(y);tp=t(k)
max_overshoot=100*(max_y-c)/c
r1=1;
while(y(r1)<0.1*c)
r1=r1+1;
end
r2=1;
while(y(r2)<0.9*c)
r2=r2+1;
end
tr=t(r2)-t(r1)
s=length(t);
whiley(s)>0.98c&&y(s)<1.02*c
s=s-1;
end
ts=t(s)
figure
(2);margin(s1);holdon
>>figure(3);sys=feedback(s1,1);impulse(sys)
>>figure(4);step(k*n1,[d1,0])
>>ess=1-dcgain(sys)
tp=0.4806
max_overshoot=947.9917
tr=0.0188
ans=0
ts=0.4429
k=260;n1=1;d1=conv(conv([10],[0.11]),[0.011]);s1=tf(k*n1,d1);
k=0:
0.05:
200;
figure
(1);rlocus(s1,k)
figure
(2);nyquist(s1)
[k,poles]=rlocfind(s1)
Selectapointinthegraphicswindow
selected_point=
-7.7133-1.5528i
k=
0.0077
poles=
-100.2220
-6.8571
-2.9209
图1-5校正前系统的根轨迹
分离点-10与虚轴的交点为-7.4289
1.5528i,当取-7.4289-1.5528i点时,k=(00.0081)1.9绘制系统校正前的Nyquist图,判断系统的稳定性,并说明理由。
图1-6校正前系统的耐奎斯特曲线
因为系统的耐奎斯特曲线的分离点为(-7.43,j0)顺时针包围(-1,j0)点1圈,所以R=-1,没有实部为正的极点所以P=0,Z=P-R=1,闭环系统不稳定。
2、校正后的系统的性能分析2.1滞后超前校正
题目中要求
取
=34rad/s,过
处作一斜率为-20dB/dec的直线作为期望特性的中频段。
为使校正后系统的开环增益不低于250rad/s,期望特性的低频段应与未校正系统特性一致。
而未校正系统的低频段斜率与期望特性的中频段斜率同为-20dB/dec,即两线段平行,为此,需在期望特性的中频段与低频段之间用一斜率为-40dB/dec的直线作连接线。
连接线与中频段特性相交的转折频率取w=1/10*wc=3.4rad/s.
图2-1
为使校正装置不过于复杂,期望特性的高频段与未校正系统特性一致。
由于未校正系统高频段特性的斜率是-60dB/dec,故期望特性中频段与高频段之间也应有斜率为-40dB/dec的直线作为连接线。
其传递函数为:
根据图2-1知.
因此,串联滞后-超前校正装置的传递函数为:
校正后系统的开环传递函数为:
2.2通过MATLAB绘制出校正后系统的bode图和校正后系统的单位阶跃响应图分别如:
MATALAB程序为:
图2-2校正后系统的bode图
图2-3校正后系统的单位阶跃响应图
图2-4校正后系统的单位脉冲响应图
图2-5校正后系统的单位斜坡响应图
2.3由图可知校正前系统的频域性能指标如下:
Kg=13.6、r=45.6、Wg=84.3;Wc=30.6.
满足相位裕度
,幅值穿越频率
两个条件,完成系统校正。
2.4特征根:
>>s=solve('(s*(2.227*s+1)*(0.01*s+1)*(0.013*s+1)+260*(0.294*s+1)=0)','s')
s=
-133.86874850086305599109243800163
-19.887924684734516793686156938933-37.472882756767467263003301646349*i
-19.887924684734516793686156938933+37.472882756767467263003301646349*i-3.7275136284071594975602558922593
2.5判断系统稳定性
由图可以看出,
之前对数幅频渐近特性曲线所对应的相频特性曲线没有穿越
,或者看特征根,可以看出所有的根都在左半平面,因此系统稳定。
2.6分析校正后单位脉冲响应曲线,单位阶跃响应曲线,单位斜坡响应曲线,这三种曲线的关系,并分析其与校正前相比有何变化
单位斜坡响应的一次导数是阶跃响应曲线,阶跃响应的一次导数是冲击响应。
与校正前作比较,可以发现,校正后阶跃响应波形上升时间变长,峰值时间变短,
调节时间变长,超调量大大减小,稳态误差保持不变。
2.7求出系统校正后动态性能指标σ%、tr、tp、ts以及稳态误差的值
n1=260;d1=conv(conv([10],[0.11]),[0.011]);
s1=tf(n1,d1);
s2=tf([1/3.41],[2.2271]);
s3=tf([0.11],[0.0131]);
sope=s1*s2*s3;
figure
(1);margin(sope);holdon
figure
(2);sys=feedback(sope,1);step(sys)
[y,t]=step(sys);
c=dcgain(sys);
[max_y,k]=max(y);
tp=t(k)
max_overshoot=100*(max_y-c)/c
r1=1;
while(y(r1)<0.1*c)
r1=r1+1;
end
r2=1;
while(y(r2)<0.9*c)
r2=r2+1;
end
tr=t(r2)-t(r1)
s=length(t);
whiley(s)>0.98*c&&y(s)<1.02*c
s=s-1;
end
ts=t(s)
figure(3);sys=feedback(s1,1);impulse(sys)
figure(4);step(k*n1,[d1,0])
ess=1-dcgain(sys)
tp=0.0916
max_overshoot=27.2139
tr=0.0363
ts=0.4503
ess=0
2.8绘制系统校正后的根轨迹图,并求其分离点、汇合点及与虚轴交点的坐标和相应点的增益
值,得出系统稳定时增益
的变化范围。
绘制系统校正后的Nyquist图,判断系统的稳定性,并说明理由。
分离点-32-2.33i与虚轴的交点为-0.7109
31.0559i,当取-0.7109-31.0559i点时,k=(01.0285)
n1=260;d1=conv(conv([10],[0.11]),[0.011]);
s1=tf(n1,d1);
s2=tf([1/3.41],[2.2271]);
s3=tf([0.11],[0.0131]);
sope=s1*s2*s3;
k=0:
0.05:
200;
figure
(1);rlocus(sope,k)
figure
(2);nyquist(sope)
[k,poles]=rlocfind(sope)
Selectapointinthegraphicswindow
selected_point=-76.0664-31.0559i
k=0.3920
poles=
-119.7115
-35.5265
-17.3689
-10.0000
-4.7652
图2-6校正后系统的根轨迹
图2-7校正后系统的耐奎斯特曲线
因为系统的耐奎斯特曲线顺时针包围(-1,j0)点0圈,所以R=0,没有实部为正的极点所以P=0,Z=P-R=0,闭环系统稳定。
4、心得体会
通过此次的课程设计,我深刻体会到对知识理解和应用的重要性,学习书本知识固然重要,但是遇到课程设计这种题目的时候还是不够,它需要全方面的思考和验证,在选取一个数值时,要前后分析计算出理论值,再用MATLAB进行验证。
在这次自控的课程设计中我受益匪浅。
这次课程设计让我学会了使用MATLAB软件作图,对我以后的实习和工作都有很大的帮助。
我还了解了自动控制中校正设计的一般步骤与要点,以后遇到这种问题,应该也比较容易解决了。
而且我发现哪怕看起来多么困难、觉得无从下手的题目,只要静下心来,多去查阅资料,多去动手操作尝试,而不是单纯的盯着书上那些我们没学过的程序语句,在那深思那语句到底什么意思,那么我们一定能从实际操作中领悟出很多真谛。
总而言之,实践是检验一切的真理,只要真正用心,认真地去对待了课程设计,最终我们一定会有不错的收获的。
五、主要参考文献
【1】程鹏.自动控制原理[M].北京:
高等教育出版社,2009
【2】徐薇莉.自动控制理论与设计[M]上海:
上海交通大学出版社,2001
【3】欧阳黎明.MATLAB控制系统设计[M].北京:
国防工业出版社,2001
【4】黄忠霖编著.自动控制原理的MATLAB实现.北京:
国防教育出版社,2007
【5】杨庚辰主编.自动控制原理.西安电子科技大学出版社,2009