文档介绍：武汉理工大学《自动控制原理》课程设计说明书
学 号：
课 程 设 计
题 目
飞行器控制系统设计
学 院
自动化学院
专 业
自动化
班 级
自动化1002班
姓 名
2 校正前后系统的matlab仿真
武汉理工大学《自动控制原理》课程设计说明书
2
2。1。1 校正前系统bode图
校正前系统Bode图源程序如下:
〉〉 num=819000
〉〉 den=[1,，0］
〉> bode(num，den)
图1 校正前系统bode图
2。1。2 校正前系统奈奎斯特曲线
校正前系统奈奎斯特图源程序如下：
〉> num=819000
>〉 den=[1,,0]
>〉 nyquist(num,den）
图2 校正前系统奈奎斯特曲线
校正前系统单位阶跃响应曲线
校正前系统的闭环传递函数为：
校正前系统单位阶跃响应源程序如下：
〉〉 num=819000
〉〉 den=［1，,819000］
〉> step（num,den）
图3 校正前系统的单位阶跃响应曲线
由系统的响应曲线可知系统的调节时间为0。，系统的超调时间为0。00349s、，都比较大。
校正后系统bode图源程序如下：
>>G=tf（819000*［，1],conv(［1，]，[0。0001444,1,0]）,bode（G)
G =
武汉理工大学《自动控制原理》课程设计说明书
3
1818 s + 819000
—--——-—-——--—————---—-—-—----—-———-
0。0001444 s^3 + s^2 + s
图4 校正后系统bode图
由校正后系统bode图可以看出校正后系统相位裕度达到,满足系统设计要求。
校正后系统奈奎斯特曲线
校正后系统奈奎斯特图源程序如下：
>〉G=tf（819000＊［，1]，conv（［1，361。2],[，1,0］)),
nyquist（G)
G =
1818 s + 819000
--——————————--—---—-——-—-—----—-——-
s^3 + s^2 + 361。2 s
图5校正后系统奈奎斯特曲线
校正后系统闭环传递函数为：
校正后系统单位阶跃响应源程序如下:
〉〉num=［,819000]
〉>den=［0。000144,,,819000］
step（num,den)
由校正后系统单位阶跃响应曲线可知系统调节时间tss,小于给定值0。008s，能够满足系统动态性能指标要求。
图6 校正后系统单位阶跃响应曲线
3校正前后系统性能比较
3．1校正前后系统bode图比较
校正前后系统比较的bode图源程序如下，其中g1函数表示校正前系统,g2函数表示校正后系统：
武汉理工大学《自动控制原理》课程设计说明书
5
>> num=819000
〉〉den=［1，,0]
>>g1=tf(num，den)
g1 =
819000
-——--——---———
s^2 + 361。2 s
>〉 g2=tf(819000*［,1］，conv（［1，],［，1，0］))
g2 =
1818 s + 819000
———-——————————--—-———--————-—-—---—
0。0001444 s^3 + s^2 + 361。2 s
>> bode（g1，g2）
由校正前后系统bode图可以看出，校正后系统截止频率变大，系统的相相位裕度变大,满足系统设计要求。
图7 校正前后系统bode图比较
校正前后系统比较的奈奎斯特图源程序如下,其中g1函数表示校正前系统，g2函数表示校正后系统:
〉>