自动控制原理课程设计 滞后校正

课程设计报告

题 目 控制系统的设计与校正

课 程 名 称 自动控制原理 院 部 名 称 专 业 自动化 班 级 学 生 姓 名 姜文华 学 号 课程设计地点 C214 课程设计学时 一周 指 导 教 师 陈丽换

成绩

金陵科技学院教务处制

一、 设计目的：

1、 了解控制系统设计的一般方法、步骤。

2、 掌握对系统进行稳定性分析、稳态误差分析以及动态特性分析的方法。

3、 掌握利用MATLAB 对控制理论内容进行分析和研究的技能。

4、 提高分析问题解决问题的能力。

二、 设计内容与要求：

设计内容：

1、阅读有关资料。

2、对系统进行稳定性分析、稳态误差分析以及动态特性分析。

3、绘制根轨迹图、Bode 图、Nyquist 图。

4、设计校正系统，满足工作要求。

设计条件：

1、已知单位负反馈系统被控制对象的传递函数为

m m 1m 2012m n n 1n 2012n b b b b ()s s s G s a s a s a s a ----++++=+++

+ （n m ≥）。 2、参数,,,012a a a a n 和b ,b ,b ,b 012m 因小组而异。

设计要求：

1、能用MATLAB 解复杂的自动控制理论题目。

2、能用MATLAB 设计控制系统以满足具体的性能指标。

3、能灵活应用MATLAB 的CONTROL SYSTEM 工具箱和SIMULINK 仿

真软件，分析系统的性能。

三、 设计方法，步骤，时间分配

1、自学MATLAB 软件的基本知识。包括MATLAB 的基本操作命令、控制系统

工具箱的用法等，并上机实验。 1天

2、基于MATLAB 用频率法对系统进行串联校正设计，使其满足给定的频域

性能指标。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函

数，校正装置的参数T ，α等的值。 1天

3、利用MATLAB 函数求出校正前与校正后系统的特征根，并判断其系统是否

稳定，为什么? 0.5天

4、利用MATLAB 作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响

应曲线，单位斜坡响应曲线，分析这三种曲线的关系？求出系统校正前与校正后

的动态性能指标%σ， r t ，p t ，s t 以及稳态误差的值，并分析其有何变化？

0.5天

5、绘制系统校正前与校正后的根轨迹图，并求其分离点、汇合点及与虚轴

交点的坐标和相应点的增益K *值，得出系统稳定时增益K *的变化范围。绘制

系统校正前与校正后的Nyquist 图，判断系统的稳定性，并说明理由？ 0.5

天

6、绘制系统校正前与校正后的Bode 图，计算系统的幅值裕量，相位裕量，

幅值穿越频率和相位穿越频率。判断系统的稳定性，并说明理由？ 0.5天

7、整理资料，撰写设计报告，准备答辩。 0.5天

8、课程设计答辩 。 0.5天 四、 成绩考核方法及标准

1、评判设计报告 50%

2、答辩成绩 30%

3、设计期间表现 20%

4、评分标准按本校教务处有关规定执行。

五、 设计主要参考资料

1、《自动控制原理》教材 程鹏 主编 机械工业出版社

2、《机电控制工程》 王积伟 主编 机械工业出版社

3、《MATLAB 控制系统设计》 欧阳黎明 主编 国防工业出版社

一、设计条件：

已知单位负反馈系统的开环传递函数0K G(S)S(0.0625S 1)(0.2S 1)

=++， 试用频率法设计串联滞后校正装置，使系统的相位裕度050γ=，静态速度误差系数1v K 40s -=，增益欲度30—40dB 。

二、设计方法和步骤

1、 基于MATLAB 用频率法对系统进行串联校正设计，使其满足给定的频域性能指标。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数，校正装置的参数T ，α等的值。如下图所示。

根据系统校正的一般步骤，先确定K0=40，然后输入如下程序得到校

正前的波特图。

>> clear

>> k0=40;

>> n1=1;

>> d1=conv(conv([1 0],[0.0625 1]),[0.2 1]);

>> s1=tf(k0*n1,d1);

>> margin(s1)

根据相角裕量γ的要求，γ=γ0（Wc1 ）+φc （Wc1 ），φc （Wc1 ）取5°，求出γ0=55 180+x=55 x=-125 由上图可知Wc1=2.46 L0（wc1）=23.2 根

据公式20lgb+ L0（wc1）=0 求得b= 0.069 bT 1

=0.1Wc1 求得T=59.1

所以校正装置的传递函数为

159.1s 1

4.081s )(++=s Gc

输入如下程序

>> clear

>> n1=conv([0 40],[4.081 1]);

>> d1=conv(conv(conv([1 0],[0.0625 1]),[0.2 1]),[59.1 1]);

>> Gh=tf(n1,d1);

>> margin(Gh)

由上图可知Wc2=0.211 L0（wc2）=26.1 根据公式20lgb+ L0（wc2）=0 求

得b= 0.828 bT 1

=0.1Wc2 求得T=5.75

所以校正装置的传递函数为

15.75s 1

4.76s )(++=s Gc

输入如下程序

clear

>> n1=conv([0 40],[4.081 1]);

>> d1=conv(conv(conv([1 0],[0.0625 1]),[0.2 1]),[59.1 1]);

D1=[4.76 1];

D2=[5.75 1];

Den1=conv(D1,n1);1

Den2=conv(D2,d1);

>> Gh=tf(Den1,Den2);

>> margin(Gh)