自动控制原理课程设计报告开环传递函数docWord下载.docx
《自动控制原理课程设计报告开环传递函数docWord下载.docx》由会员分享，可在线阅读，更多相关《自动控制原理课程设计报告开环传递函数docWord下载.docx（27页珍藏版）》请在冰豆网上搜索。
（1）根据要求的稳态品质指标，求系统的开环增益值；
（2）根据求得的值，画出校正前系统的Bode图，并计算出幅值穿越频率、相位裕量（要求利用MATLAB软件编程进行辅助设计）,以检验性能指标是否满足要求。
若不满足要求，则执行下一步；
（3）画出串联校正结构图，分析并选择串联校正的类型（超前、滞后和滞后-超前校正）。
（若可以采用多种方法，分别设计并进行比较）
（4）确定校正装置传递函数的参数；
（5）画出校正后的系统的Bode图，并校验系统性能指标（要求利用MATLAB软件编程进行辅助设计）。
若不满足，跳到第（4）步。
否则进行下一步。
（6）提出校正的实现方式及其参数。
（要求实验实现校正前、后系统并得到的校正前后系统的阶跃响应）
（7）若采用不同串联校正方法进行设计，比较不同串联校正方法的特点。
（稳定性、稳态性能、动态性能和实现的方便性的比较）
二、设计正文
、第一问
1-1题目的处理
设
K=10时
K*=100
GK=100/[s（s+1）]
G=100/[s^2+10s+100]
1-2题目的代码
在MATLAB中运行以下代码：
GK=tf（[100],[1100]）;
G=tf（[100],[110100]）;
step（G）
Wn=10;
Kexi=10/（2*Wn）;
Beta=acos（Kexi）;
Wd=Wn*sqrt（1-Kexi^2）;
Tr=（pi-Beta）/Wd;
Tp=pi/Wd;
Ts=4.4/（Kexi*Wn）;
Ct=exp（-pi*Kexi/sqrt（1-Kexi^2））*100;
1-3运行的的结果
Transferfunction——G:
100
-----------------------
s^2+10s+100
Wn=
10
Kexi=
0.5000
Beta=
1.0472
Wd=
8.6603
Tr=
0.2418
Tp=
0.3628
Ts=
0.8800
Ct=
16.3034
1-4比较
以上数据与图中数据相吻合
、第二问
2-1题目的处理
GK=
K*
-------------
s^2+10s
分K*>
0和K*<
0两项
2-2题目的代码
当K*>
0时；
num=[1];
den=[1100];
rlocus（num,den）;
[Kp]=rlocfind（num,den）
title（‘ROOTLOCUS’）;
K*<
0时
num=[-1];
[K*p]=rlocfind（num,den）；
2-3运行的结果
K*=72.8796即K=7.28396
p=-14.89344.8934
K>
K<
2-4结论
由图可看出K*>
0即K>
0可满足系统稳定；
由图可看出K*<
0即K<
0时，系统不稳定
、第三问
3-1题目的处理
由以上可知K>
0时,系统都稳定
故设K=10,K*=100;
则GK=100/[s（s+10）];
3-2题目的代码
在MATLAB中代码为：
num=[100];
[magphasew]=bode（num,den）
margin（mag,phase,w）
nyquist（num,den）
grid
3-3运行的结果
3-4结论
根据nyquist稳定判据，R=0，P=0，Z=0
可知系统稳定
、第四问
第四问中的第一小题
4-1-1题目的分析
由单位斜坡信号作用下
<
=0.01得出K=100
GK=100/[s（0.1s+1）];
先做出未校正前的相关图形；
4-1-2题目的代码
Gk=tf（[100],[0.110]）;
G=feedback（Gk,1）;
num1=[100];
den1=[0.110];
figure
（1）;
[mag1,phase1,w]=bode（num1,den1）;
margin（mag1,phase1,w）
4-1-3运行结果
看图可知不满足相位裕度>
45
（图见附录）
4-1-4系统的校正
4-1-4-1超前校正
4-1-4-1-1问题的计算
由计算可得T=0.04sa=3.36;
Gc=（0.s+1）/（0.04s+1）;
结构图
原理图
4-1-4-1-2题目的代码
代码；
Gk=tf（num0,den0）;
代码；
num2=[0.0311];
den2=[0.0621];
num0=conv（num2,num1）;
den0=conv（den2,den1）;
[mag0,phase0,w]=bode（num0,den0）;
figure;
margin（mag0,phase0,w）
4-1-4-1-3运行的结果
4-1-4-2滞后校正
4-1-4-2-1问题的计算
计算可得T=11.87b=0.104;
Gc=（1+1.23s）/（1+11.87s）
4-1-4-2-2问题的代码
figure
num3=[1.231];
den3=[11.871];
num4=conv（num3,num1）;
den4=conv（den3,den1）;
[mag4,phase4,w]=bode（num4,den4）;
margin（mag4,phase4,w）;
由代码
Gk=tf（num4,den4）;
4-1-4-2-3运行的结果
第四问中的第二小题
4-2-1题目的分析
由题在单位斜坡信号作用下，系统的稳态误差<
=0.005得出K=200;
GK=200/[s（0.1s+1）];
4-2-2题目的代码
Gk=tf（[200],[0.110]）;
num1=[200];
margin（mag1,phase1,w）
4-2-3运行的结果
4-2-4系统校正
4-2-4-1超前校正
4-2-4-1-1题目的计算
由计算可得T=0.016sa=1.63;
Gc=（1+0.026s）/（1+0.016s）;
4-2-4-1-2题目的代码
由代码
num2=[0.0261];
den2=[0.0161];
margin（mag0,phase0,w）;
4-2-4-1-3运行的结果
（见附录）
4-2-4-2滞后校正
4-2-4-2-1问题的分析
计算可得T=23.52b=0.052;
Gc=（1+1.22s）/（1+23.52s）
4-2-4-2-2问题的代码
Figure
num3=[1.221];
den3=[23.521];
4-2-4-2-3运行的结果
4-2-4-3滞后—超前校正
4-2-4-3-1问题的分析
计算可得Tb=0.1
=50a=40
=0.51rad/s
Gc=[（1+1.28s）（1+0.1s）]/[（1+5.02s）（1+0.03s）];
4-2-4-3-2问题的代码
num5=conv（[201],[0.11]）;
den5=conv（[78.41],[0.00251]）;
num6=conv（num5,num1）;
den6=conv（den5,den1）;
[mag6,phase6,w]=bode（num6,den6）;
margin（mag6,phase6,w）;
Gk=tf（num6,den6）;
4-2-4-3-3运行的结果
三、附录
第一问
第二问
第三问
第四问
第一题
未校正的单位阶跃响应图
未校正前的BODE图
超前校正
可得校正后的BODE图
可得超前校正后的单位阶跃响应图
滞后校正
可得校正后的BODE图
可得滞后校正后的单位阶跃响应图
第二题
可得未校正前的BODE图
滞后—超前校正
可得超前—滞后校正后的单位阶跃响应图
四、总结
此次课程设计使我对自动控制原理这门学科有了更深层次的理解与认识，在学习和操作过程中收获颇多，知道了自动控制原理的实际用途与用法，为我以后工作中使用所学过的自动控制原理知识奠定了良好的基础，使我了解并掌握了一定的软件运用，对一些数据的更能深入的分析，通过画各种图，让我对于数据的变化更加清晰，对这个课题的掌握更加全面，此次课程设计，为我以后的实际应用更打下了坚实的基础。
最后，对此次帮助我的同学表示感谢，特别要致谢老师这个引路人，谢谢大家对我的帮助和教导。
五、参考文献
1、程鹏.自动控制原理[M].北京：
高等教育出版社，2009
2、徐薇莉.自动控制理论与设计[M]上海：
上海交通大学出版社,2001
3、欧阳黎明.MATLAB控制系统设计[M].北京：
国防工业出版社,2001
升级会员 