用MATLAB进行控制系统滞后校正设计

学号：0121011360128

课程设计

用MATLAB进行控制系统的滞题目

后校正设计

学院自动化

专业自动化

班级

姓名

指导教师肖纯

2012 年12 月23 日

课程设计任务书

学生姓名： 专业班级： 自动化1001班 指导教师： 肖 纯 工作单位： 自动化学院

题 目: 用MATLAB 进行控制系统的滞后校正设计。 初始条件：已知一单位反馈系统的开环传递函数是

)

1)(01(10)(s s s K

s G ++=

要求系统的静态速度误差系数1100-=s K v ，相角裕度 40≥γ，并且幅值裕度不小于10分贝。

要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写

等具体要求）

（1） 用MATLAB 作出满足初始条件的K 值的系统伯德图，计算系统的幅值

裕度和相位裕度。

（2） 系统前向通路中插入一相位滞后校正，确定校正网络的传递函数，

并用MATLAB 进行验证。

（3） 用MATLAB 画出未校正和已校正系统的根轨迹。

（4） 课程设计说明书中要求写清楚计算分析的过程，列出MATLAB 程序和

MATLAB 输出。说明书的格式按照教务处标准书写。

时间安排：

指导教师签名： 年 月 日 系主任（或责任教师）签名： 年 月 日

目录

1．正特性及校正方法 (1)

1.1滞后校正特性 (1)

1.2滞后校正设计的一般步骤与方法 (1)

2．未校正时系统分析 (3)

2.1伯德图绘制 (3)

2.2未校正系统的相位裕度和幅值裕度 (4)

3.确定滞后校正传递函数 (5)

4.系统校正前后根轨迹图 (7)

4.1未校正系统根轨迹图 (7)

4.2校正后系统根轨迹图 (8)

5.心得体会 (10)

参考文献 (11)

1．正特性及校正方法

1.1滞后校正特性

滞后校正就是在前向通道中串联传递函数为)(s G c 的校正装置来校正控制系统，)(s G c 的表达式如下所示。

1

,11)(<++=

a Ts aTs

s G c （1-1）

其中，参数a 、T 可调。滞后校正的高频段是负增益，因此，滞后校正对系统中高频噪声有削弱作用，增强了抗干扰能力。可以利用滞后校正的这一低通滤波所造成的高频衰减特性，降低系统的截止频率，提高系统的相位裕度，以改善系统的暂态性能。

滞后校正的基本原理是利用滞后网络的高频幅值衰减特性使系统截止频率下降，从而使系统获得足够的相位裕度。或者，是利用滞后网络的低通滤波特性，使低频信号有较高的增益，从而提高了系统的稳态精度。

可以说，滞后校正在保持暂态性能不变的基础上，提高开环增益。也可以等价地说滞后校正可以补偿因开环增益提高而发生的暂态性能的变化。

1.2滞后校正设计的一般步骤与方法

（1）按稳态性能指标要求的开环放大系数绘制未校正系统的伯德图。

如果未校正系统需要补偿的相角较大，或者在截止频率附近相角变化大，具有这样特性的系统一般可以考虑用滞后校正。

（2）在未校正系统的伯德图上找出相角为）（εγ--︒180-的频率作为校正后系

统的截止频率'c w ，其中γ为要求的相位裕度，ε为补偿滞后校正在'

c w 产生的相

位滞后，一般取︒︒10~5。

ε的选取：ε是为了补偿滞后校正的相位滞后的，一般限制滞后校正的滞后

相角小于︒10，所以可以取小于︒10的值。ε应取一个尽量小，但又能补偿滞后校

正在'c w 处的滞后相角的值。一般，若'c w 较大，ε可取小一些。反之，若'

c w 小，

则ε取大一些。

（3）在未校正系统的伯德图上量取量取）（'0c w L （或由|)(|lg 20'0c jw G 求取）的分贝值，并令a

jw G c 1

lg

20|)(|lg 20'0=，由此确定参数a(a<1)。 这一步的意思是，在'

c w w =处，设计滞后校正的幅值与原系统的幅值反向相

等才能相互抵消，使校正后系统的截止频率为'

c w 。

（4）取'

)10

1~51(1c

w aT =，并由a 求参数T 。 选

aT

1

的原则是使)(c c w ϕ不超过（2）中所选的ε值，即εϕ≤-=T T 'c 'c 'c c arctanaw arctanw |)(w |。校正时只需要近似值，aT

1

可近似由下式

选取。

ε

εtan )90cot(1'

'c c w w aT =-︒= (1-2)

（5）绘制校正后系统的伯德图，校验各项性能指标，若不满足，可重新选

择'

c w 或aT 1的值

2．未校正时系统分析

2.1伯德图绘制

(1)根据稳态误差的要求，确定系统的开环增益K

100

)1)(110(10lim 0

=++=→s s s K

s

K s v （2-1）

则解得100k =

s s s s s s K S G 1021131000

)1)(10(10）(+∧+∧=++=

（2-2）

（2）由于开环增益K=100,在MATLAB 中输入以下命令：

G=tf(1000,[1,11,10,0]); [kg,r]=margin(G);

[mag,phase,w]=bode(G); %绘制G 的幅频特性和相频特性曲线 margin(G)

matlab 得出如下结果：

Warning: The closed-loop system is unstable. > In warning at 26

In DynamicSystem.margin at 63

kg = 0.1100

r = -34.3247