自动控制理论课程设计自动控制系统建模、分析及校正.doc

课程设计（大作业）报告

课程名称： 自动控制理论
设计题目： 自动控制系统建模、分析及校正
院 系：自动控制与机械工程学院
班 级：电气工程及其自动化2014级 2班
设 计 者：
学 号：
指导教师：
设计时间：2016.12.12——2016.12.16

昆明学院课程设计（大作业）任务书
姓 名： 院（系）：自动控制与机械工程学院
专 业：电气工程及其自动化 学 号：201404170146
任务起止日期：第十五周（2016年12月12日——2016年12月16日）

课程设计题目：
自动控制系统建模、分析及校正

课程设计要求：
1．了解matlab软件的基本特点和功能,熟悉其界面、菜单和工具条；掌握线性系统模型的计算机表示方法、变换以及模型间的相互转换。了解控制系统工具箱的组成、特点及应用；掌握求线性定常连续系统输出响应的方法，运用连续系统时域响应函数（impulse,step,lsim），得到系统的时域响应曲线。
2．掌握使用MATLAB软件绘制开环系统的幅相特性曲线、对数频率特性曲线；观察控制系统的观察开环频率特性，对控制系统的开环频率特性进行分析；
3．掌握MATLAB软件中simulink工具箱的使用；熟悉simulink中的功能模块，学会使用simulink对系统进行建模；掌握simulink的仿真方法。

工作计划及安排：
在第15周内（2016.12.12——2016.12.16）内完成规定的题目。
1：2016.12.12
明确注意事项，完成题目1、2。

2：2016.12.13
完成题目3、4、5。

3：2016.12.14
完成题目6、7、8。

4：2016.12.15
完成题目9、10。

5：2016.12.16
完成题目11，打印并提交课程设计报告。

指导教师签字

2016年 12 月 8 日
课程设计（大作业）成绩
学号：201404170146 姓名：张雨 指导教师：杨祖元
课程设计题目：自动控制系统建模、分析及校正

总结：MATLAB 是一种直观、高效的计算机语言，同时也是一个科学计算平台。它的伴随工具Simulink 是用来对真实世界的动力学系统建模、模拟仿真和分析的软件。我们可将综合性和设计性实验项目通过MATLAB 在计算机上仿真，使系统的观察实验的动态过程。目前，MATLAB 已经成为我们当代大学生必须掌握的基本技能，在设计研究单位和工业部门，MATLAB 已经成为研究和解决各种具体工程问题的一种标准软件。在完成了验证性、综合性和设计性实验后，课程设计必不可少。课程设计是工科实践教学的一个重要的环节，目的是培养我们综合运用理论知识分析和解决实际问题的方法和能力,实现由知识向技能的初步化。所以课程设计是培养我们思维创造能力最有效的途径。

指导教师评语：

成绩：

填表时间： 指导教师签名：
课程设计（大作业）报告

1.用matlab语言编制程序，实现以下系统：
num =[4 56 288 72 720 288]
den =[1 6 14 21 24 17 5 0]
Transfer function:

Transfer function
1） 2）

2.两环节G1、G2串联，求等效的整体传递函数G(s)

Transfer function:

3.两环节G1、G2并联，求等效的整体传递函数G(s)

Transfer function:

4.已知系统结构如图，求闭环传递函数。其中的两环节G1、G2分别为

（2）负反馈：

Transfer function:

（1）正反馈：

Transfer function:

5.已知某闭环系统的传递函数为，求其单位阶跃响应曲线，单位脉冲响应曲线。
（1）单位阶跃响应曲线：

（2）单位脉冲响应曲线：

6.典型二阶系统的传递函数为 ， 为自然频率， 为阻尼比，试绘出当=0.5，分别取-2、0、2、4、6、8、10时该系统的单位阶跃响应曲线；分析阻尼比分别为–0.5、–1时系统的稳定性。

（1）当=0.5，分别取2、4、6、8、10时
w=0:2:10;
kosai=0.5;
figure(1)
hold on
for Wn=w
num=Wn^2;
den=[1,2*kosai*Wn,Wn^2];
step(num,den)
end
grid on;

w=-2

kosai=0.5;
figure(1)
Wn=w
num=Wn^2;
den=[1,2*kosai*Wn,Wn^2];
step(num,den)
grid on;
title('µ¥Î»½×Ô¾ÏìÓ¦')
xlabel('ʱ¼ä')
ylabel('Õñ·ù')
w =-2
Wn = -2

（2）分析阻尼比分别为–0.5、–1时系统的稳定性
1、当=-0.5时：
w=-2:2:10;
kosai=-0.5;
hold on
for Wn=w
den=[1,2*kosai*Wn,Wn^2];
v=roots(den)
end
v = -1.0000 + 1.7321i
-1.0000 - 1.7321i
v = 0
0
v = 1.0000 + 1.7321i
1.0000 - 1.7321i
v = 2.0000 + 3.4641i
2.0000 - 3.4641i
v = 3.0000 + 5.1962i
3.0000 - 5.1962i
v = 4.0000 + 6.9282i
4.0000 - 6.9282i
v = 5.0000 + 8.6603i
5.0000 - 8.6603i
由特征根可得，当w=-2时系统稳定，当w=0时系统临界稳定，当w=2、4、6、8、10时，系统不稳定。

2、当=-1时：
clear
w=-2:2:10;
kosai=-1;
hold on
for Wn=w
den=[1,2*kosai*Wn,Wn^2];
v=roots(den)
end
v = -2
-2
v = 0
0
v = 2
2
v = 4
4
v =6.0000 + 0.0000i
6.0000 - 0.0000i
v = 8
8
v = 10
10
由特征根可得，当w=-2时系统稳定，当w=0时系统临界稳定，当w=2、4、6、8、10时，系统不稳定。

7. 设有一高阶系统开环传递函数为，试绘制该系统的零极点图和闭环根轨迹图。
（1）系统的零极点：

num=[0.016 0.218 1.463 9.359];
den=[0.06 0.268 0.635 6.271];
[z,p,k]=tf2zp(num,den)
pzmap(num,den)
运行结果：
z = -10.2679
-1.6785 + 7.3587i
-1.6785 - 7.3587i
p = -5.7710
0.6522 + 4.2054i
0.6522 - 4.2054i
k = 0.2667
（2）系统的闭环根轨迹：
num=[0.016 0.218 1.463 9.359];
den=[0.06 0.268 0.635 6.271];
rlocus(num,den)

8.单位反馈系统前向通道的传递函数为： ，试绘制该系统的Bode图和Nyquist曲线，说明软件绘制曲线与手动绘制曲线的异同。
（1）bode图
num=[0,0,2,8,12,8,2];
den=[1,5,10,10,5,1,0];
bode(num,den)

（2）Nyquist曲线
num=[0,0,2,8,12,8,2];
den=[1,5,10,10,5,1,0];
nyquist(num,den)

9．已知某控制系统的开环传递函数，试绘制系统的开环频率特性曲线，并求出系统的幅值与相位裕量。
G=tf(1.5,conv(conv([1,0],[1,1]),[1,2]));
bode(G)
grid
[Gm,Pm,Wcg,Wcp]=margin(G)
Gm =4.0000
Pm =41.5340
Wcg =1.4142
Wcp =0.6118

10.在SIMULINK中建立系统，该系统阶跃输入时的连接示意图如下。k为学生学号后三位。绘制其单位阶跃响应曲线，分析其峰值时间tp、延迟时间td、上升时间tr、调节时间ts及超调量。

num1=[1];
den1=[1,0];
sys1=tf(num1,den1);
num2=[146];
den2=[1,9];
sys2=tf(num2,den2);
sys12=sys1*sys2;
G=feedback(sys12,1);
step(G)
grid

分析其峰值时间tp、延迟时间td、上升时间tr、调节时间ts及超调量,闭环系统的函数：

（1） 峰值时间：
G1=tf(146,[1,9,146]);
[y,t]=step(G1);
[Y,K]=max(y);
tp=t(K)
tp =0.2798
（2） 延迟时间
[y,t]=step(G1);
c=dcgain(G);
i=1;
while y(i)<0.5*c
i=i+1;
end
td=t(i)
td =0.1049
（3） 上升时间
G1=tf(146,[1,9,146]);
[y,t]=step(G1);
c=dcgain(G);
n=1;
while y(n)<c
n=n+1;
end;
tr=t(n)
tr =0.1749
（4） 调节时间
c=dcgain(G);
i=length(t);
while(y(i)>0.98*c)&(y(i)<1.02*c)
i=i-1;
end
ts=t(i)
ts = 0.8743
（5） 超调量
G1=tf(146,[1,9,146]);
[y,t]=step(G1);
c=dcgain(G1);
[Y,K]=max(y);
a=(Y-c)/c
a = 0.2835

*11. 给定系统如下图所示，试设计一个串联校正装置，使截止频率ωc=40rad/s、相位裕度γ≥45º

G=tf(100,[0.04,1,0]);
[Gw,Pw,Wcg,Wcp]=margin(G)
G1=tf(100,[0.04,1,0]);
G2=tf(100*[0.025,1],conv([0.04,1,0],[0.01,1]))
bode(G1)
hold
bode(G2,'r')
grid
figure
G1c=feedback(G1,1);
G2c=feedback(G2,1);
step(G1c)
hold
step(G2c,'r')
grid
Gm =4.0000
Pm =41.5340
Wcg =1.4142
Wcp =0.6118
Gw =Inf
Pw =28.0243
Wcg =Inf
Wcp =46.9701
Transfer function:

矫正前后的bode图

矫正前后的时域响应图目 录

第一章 总 论 1
一、项目提要 1
二、可行性研究报告编制依据 2
三、综合评价和论证结论 3
四、存在问题与建议 4
第二章 项目背景及必要性 5
一、项目建设背景 5
二、项目区农业产业化经营发展现状 11
三、项目建设的必要性及目的意义 12
第三章 建设条件 15
一、项目区概况 15
二、项目实施的有利条件 17
第四章 建设单位基本情况 19
一、建设单位概况 19
二、研发能力 20
三、财务状况 20
第五章 市场分析与销售方案 21
一、市场分析 21
二、产品生产及销售方案 22
三、销售策略及营销模式 22
四、销售队伍和销售网络建设 23
第六章 项目建设方案 24
一、建设任务和规模 24
二、项目规划和布局 24
三、生产技术方案与工艺流程 25
四、项目建设标准和具体建设内容 26
五、项目实施进度安排 27
第七章 投资估算和资金筹措 28
一、投资估算依据 28
二、项目建设投资估算 28
三、资金来源 29
四、年度投资与资金偿还计划 29
第八章 财务评价 30
一、财务评价的原则 30
二、主要参数的选择 30
三、财务估算 31
四、盈利能力分析 32
五、不确定性分析 33
六、财务评价结论 34
第九章 环境影响评价 35
一、环境影响 35
二、环境保护与治理措施 35
三、环保部门意见 36
第十章 农业产业化经营与农民增收效果评价 37
一、产业化经营 37
二、农民增收 38
三、其它社会影响 38
第十一章 项目组织与管理 40
一、组织机构与职能划分 40
二、项目经营管理模式 42
三、技术培训 42
四、劳动保护与安全卫生 43
第十二章 可行性研究结论与建议 46
一、可行性研究结论 46
二、建议 47