太阳观测控制系统设计课程设计.doc

学 号：
0121211350518

课 程 设 计

题 目
太阳观测控制系统设计
学 院
自 动 化
专 业
电气工程及其自动化
班 级
电气1205
姓 名
刘爽
指导教师
熊和金、李浩

2015
年
1
月
16
日

I
武汉理工大学《自动控制原理》课程设计说明书
课程设计任务书

学生姓名： 刘爽 专业班级： 电气1205班
指导教师： 熊和金、李浩 工作单位： 自动化学院
题 目 太阳观测控制系统设计
初始条件：
太阳观测控制系统的开环传递函数为：

要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）
1．分别用时域和频域方法设计该系统的控制器。
控制系统的时域性能指标为：
斜坡输入产生的稳态误差2%
阶跃响应的最大超调<5%
上升时间0.03秒
调节时间0.08秒
控制系统的频域性能指标为：
斜坡输入产生的稳态误差1%
相位裕量大于
2. 用Matlab对校正前后的系统进行仿真分析，画出阶跃响应曲线，计算其时域性能指标。

时间安排:
任务
时间（天）
审题、查阅相关资料
1
分析、计算
3
编写程序
1
撰写报告
2
论文答辩
0.5

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

I
武汉理工大学《自动控制原理》课程设计说明书
目录
摘要 I
1、以时域方法设计系统的控制器 2
1.1理论分析 2
1.1.1典型情况分析 2
1.1.2基本参数计算 3
1.1.3上升时间的计算 3
1.1.4超调量的计算 3
1.1.5调节时间的计算 3
1.1.6理论分析结论 4
1.2方案论证与选择 4
1.2.1方案1—比例-微分控制 4
1.2.2方案2—测速负反馈控制 5
1.2.3方案选择 6
1.3参数的确定 6
1.3.1稳态误差 6
1.3.2阶跃信号输入的超调量 7
1.3.3上升时间 7
1.3.4调节时间 7
1.4系统设计 7
2、以频域方法计系统的控制器 8
3、MATLAB仿真分析 11
3.1 系统的时域分析 11
3.1.1斜坡输入稳态误差 11
3.1.2阶跃输入性能 11
3.2系统的频域分析 13
结束语 14
参考文献 15
本科生课程设计成绩评定表 16

1
武汉理工大学《自动控制原理》课程设计
摘要
自动控制技术已广泛应用于制造业、农业、交通、航空及航天等众多产业部门，极大地提高了社会劳动生产率，改善了人们的劳动条件，丰富和提高了人民的生活水平。在今天的社会生活中，自动化装置无所不在，为人类文明进步做出了重要贡献。 
随着计算机技术的飞速发展，控制系统计算机辅助设计技术在工具、理论和算法上取得了巨大的进步，以前难于设计的控制系统现在可用新方法和新策略较容易地得到结果。在目前诸多控制系统设计方法中，处于主导地位的是频域设计方法和时域设计方法。 
计算机的不断普及和各工具软件的不断发展，Matlab已经作为一个很强大的软件而受到人们的广泛应用。结合Matlab来分析系统的动态性能，已经成为了未来发展的趋势。
关键词：太阳观测控制系统 时域分析 频域分析 Matlab仿真

I
武汉理工大学《自动控制原理》课程设计
太阳观测控制系统设计
1、 以时域方法设计系统的控制器
1.1理论分析
1.1.1典型情况分析
根据设计要求，若将系统设计为最基本的单位负反馈的二阶系统。单位负反馈的二阶系统的结构图如图1-1所示：

图1-1 单位二阶系统的结构图
其开环传递函数为：

其中，是二阶系统的阻尼比；是自然振荡频率或无阻尼振荡角频率，简称为自然振荡频率或无阻尼振荡频率。式②也可以写作：

其中，称为开环增益，为时间常数。
二阶系统的闭环传递函数为：

由于系统不存在闭环零点，上式所示的系统为无零点的二阶系统。显然二阶系统的响应取决于、两个参数，令式(4)等于零，得二阶系统特征方程为：
（5）
系统的闭环极点即特征方程式的根为：
（6）
从式(6)可以看出，当阻尼比，二阶系统的两个闭环极点为正实数或一对具有正实数部的共轭复数，位于右半S平面；当阻尼比，二阶系统的两个闭环极点为负实数或一对具有负实部的共轭复数，位于左半S平面。
1.1.2基本参数计算
由开环传递函数可得：
（7）
（8）
（9）
由式(7)、(8)、(9)可将、用含有K的式子表示：，，；阻尼角为：。
1.1.3上升时间的计算
欠阻尼二阶系统单位阶跃响应为：
（10）
上升时间是指第一次上升到稳态值所需要的时间。由于的稳态值为1，在计算上升时间时，令，可得到：
（11）
1.1.4超调量的计算
系统上升到最大值的峰值时间为：
（13）
将式（14）代入式（10），得到其峰值为：
（15）
所以可以计算其超调量为：
（16）
1.1.5调节时间的计算
从理论上来说要精确计算是十分困难的，因为曲线的收敛速度比比其包络线要快，所以可以用包络线代替响应曲线进行估算，所得结果略保守。设误差带为，包络线进入误差带的时间即可认为是调节时间设误差带为，则可得到：
（18）
对于或者 ，常采用以下公式计算调节时间：
（19）
（20）
本次设计采取误差为进行设计，所以调节时间用含有K的式子表示为：

1.1.6理论分析结论
按照最基本的单位负反馈二阶系统设计控制器，不能满足设计的要求，所以需要加入控制环节来增大系统的阻尼比，以达到控制调节的时间的目的。根据设计要求，添加比例-微分控制器或者测速负反馈控制器均能达到增大系统的阻尼比的目的。下面将对这两个控制方式进行对比及选择。
1.2方案论证与选择
1.2.1方案1—比例-微分控制
比例-微分控制是通过引入误差信号的倒数作为附加的控制信号，从而可以减小系统的超调量，如图1-2所示是比例-微分控制系统的结构图：

图1-2 比例-微分控制系统的结构图
通过调整的大小，可以改变附加信号的强弱，采用比例-微分控制时系统的开环传递函数为

闭环传递函数为

系统的阻尼比为

可见对于阻尼比为的欠阻尼二阶系统，采用比例-微分控制时其阻尼比增加为，从而使系统的超调量下降，调节时间变短。但是采用比例-微分控制后系统为有零点的二阶系统。
1.2.2方案2—测速负反馈控制
测速负反馈控制是通过引入输出量的倒数作为附加控制信号，来减小系统的超调量，如图1-3是测速负反馈系统的结构图：

图1-3测速负反馈系统的结构图
通过调整的大小，可以改变附加信号的强弱，采用测速负反馈时系统的开环传递函数为

闭环传递函数为

引入测速负反馈控制后，系统的阻尼比为

可见，测速负反馈控制增大了系统的阻尼比，从而有利于减小系统的超调量。采用测速负反馈控制后，由于闭环函数没有零点，性能指标的计算可以采用单位负反馈二阶系统的性能指标计算公式。
1.2.3方案选择
相比较比例-微分控制与测速负反馈两种控制方式，测速负反馈系统计算与典型无零点二阶系统性能指标的计算方式一样，比比例-微分控制分析更为简便，能满足本次设计的要求，所以采用测速负反馈控制。
1.3参数的确定
1.3.1稳态误差
当输入单位阶跃信号时，系统的静态误差为，得到，所以稳态误差。
当输入信号为单位斜坡函数时，系统输出量的拉氏变换式为

令，对上式取拉氏反变换可得单位斜坡响应

上式表明，欠阻尼二阶系统的单位斜坡响应由稳态分量和瞬态分组成。
控制系统的误差响应为:

当时间趋于无穷时，误差响应的稳态值称为稳态误差，以标志。对单位斜坡响应其稳态误差为:

将、和代入上式有

1.3.2阶跃信号输入的超调量
采用测速负反馈控制后，由于闭环函数没有零点，性能指标的计算可以采用单位负馈二阶系统的性能指标计算公式。所以此时超调量应为

1.3.3上升时间
欠阻尼二阶系统单位阶跃响应为：

上升时间是指第一次上升到稳态值所需要的时间。由于的稳态值为1，在计算上升时间时，令，可得到：

1.3.4调节时间
从理论上来说要精确计算是十分困难的，因为曲线的收敛速度比比其包络线要快，所以可以用包络线代替响应曲线进行估算，所得结果略保守。设误差带为，包络线进入误差带的时间即可认为是调节时间设误差带为，则可得到：

对于或者 ，常采用以下公式计算调节时间：

所以调节时间用含有K的式子表示为：

1.4系统设计
为满足设计要求，可分别得到

借助EXCEL计算，当和时能使太阳观测控制系统在时域下满足时域性能要求。此时的性能参数如下表所示：
表1-1 校正后系统的性能参数

60
23/6000
0.58%

0.014s
4.32%
0.005s
此时太阳观测控制系统开环传递函数为

2、 以频域方法设计系统的控制器
已知开环系统传递函数，其速度误差传递系数为

所以当输入为斜坡函数时，其稳态误差为

因为稳态误差，所以，取时，令，得到系统的频率特性为

所以其对数幅频特性为
设为系统的截止频率，可得

其相角裕度为

由上两式可求得,
因为系统要求，故系统需要使用串联超前校正来增大系统的相角裕度，从而满足系统要求。
采用串联超前校正时，串联超前校正的传递函数为

超前网络对频率在至之间的输入信号有明显的微分作用，在该频率范围内，输出信号相角比输入信号相角超前，在最大超前角频率处，具有最大超前角，且正好处于频率和的几何中心。
超前网络的相角为

将上式对求导并令其为零，得最大超前角频率

最大超前角

因为为控制系统要求的最小相角裕度，则有

为补偿角是用于补偿因超前校正装置的引入，考虑到校正后截止频率将比校正前截止频率略有增加，因而在计算公式中预先增加5°~10°。在此，取 。所以有

因为
，
由此可求得

设为校正后的截止频率，则有，即

求得

所以串联超前校正传递函数为

校正后的系统开环传递函数为

综上所述，可知当、、时太阳观测控制系统的开环传递函数为

此时系统的截止频率为，可计算相角裕度为。输入信号为单位斜坡信号时，其误差传递系数为

稳态误差为

可以使得控制系统在频域下满足以下两个条件：
（1）单位斜坡输入产生的稳态误差；
（2）控制系统的相角裕度
3、 MATLAB仿真分析
3.1 系统的时域分析
3.1.1斜坡输入稳态误差
Matlab中用函数dcgain（）计算系统的静态误差系数，通过对传递函数做适当处理，也可用dcgain（）来计算控制系统的静态位置误差系数，静态速度误差系数，静态加速度误差系数。
首先用Matlab求校正前的静态速度误差系数，代码与运行结果如下：
>> G=tf([2000],[1 30 0]);
>> sG=tf([2000 0],[1 30 0]);
>> kv=dcgain(sG)

kv1=
66.6667
>> ess1=1/kv1
ess1=
0.0150
可见系统在未校正之前，系统的稳态误差为1.5%。
用Matlab求校正后的静态速度误差系数，代码与运行结果如下：

>> G=tf([120000],[1 490 0]);
>> sG=tf([120000 0],[1 490 0]);
>> kv2=dcgain(sG)

kv2 =
244.8980
>> ess2=1/kv2
ess2 =
0.0041
可见在系统在校正之后，系统的稳态误差为0.41%<1%满足设计要求。
Matlab没有直接求系统斜坡响应的功能函数。在求取控制系统的斜坡响应时，通常利用阶跃响应功能函数。基于单位阶跃信号的拉氏变换为，而单位斜坡信号的拉氏变换为。所以在求取控制系统的单位斜坡响应时，可利用阶跃响应的功能函数step()求取传递函数为的系统的阶跃响应，则其结果就是原系统G(s)的斜坡响应。
3.1.2阶跃输入性能
首先用MATLAB分别建立校正前与校正后的simulink模型如图3-1和图3-2所示：

图3-1 系统校正前的simulink模型

图3-2 系统校正后的simulink模型
首先用MATLAB绘制未校正时系统的单位阶跃响应曲线，如图3-1所示，代码如下：
num=[2000];
den=[1,30,2000];
step(num,den)

图3-1 系统校正前的单位阶跃响应
由响应曲线图可以得出校正前的参数如下表所示：
表3-1 校正前系统的性能参数

0.0742s
0.344
33%
用MATLAB绘制未校正时系统的单位阶跃响应曲线如图3-2所示，代码如下：
num=[120000]
den=[1,490,120000]
step(num,den)

图3-2 系统校正后的阶跃响应
由响应曲线图可以得出校正前的参数如下表所示：
表3-1 校正前系统的性能参数

0.0113s
0.0206
4%
3.2系统的频域分析
首先用Matlab绘制校正前系统的伯德图并求幅值裕度、相角裕度如图3-3所示，其代码如下：
num=[3000];
den=[1,30,0];
g=tf(num,den);
margin(g);
grid；

图3-3 系统校正前的伯德图
用Matlab绘制校正后系统的伯德图并求幅值裕度、相角裕度如图3-4所示，其代码如下：
num=[3000];
den=[1,30,0];
g=tf(num,den);
margin(g);
grid；

图3-4 系统校正后的伯德图

结束语
自动控制技术广泛运用于生活中各行各业，极大地提高了生产率和改善了人们的生活水平，它已成为现代生活中不可缺少的重要组成部分。所谓自动控制，是指在没有人直接参与的情况下，利用外加的设备或装置，是机器、设备或生产过程的某个工作状态或参数自动地按照预定的规律运行。
回顾起此次自动控制课程设计，至今我仍感慨颇多，从审题到完成，从理论到实践，在整整两星期的日子里，熬了好多夜，经历了烦躁到平静，但是学到很多很多的的东西，不仅巩固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识，不仅是专业知识，更学会熟练的使用MATLAB这个非常实用的数学软件，它帮我解决了许多问题，但是它的编程也是一个让人头疼的问题，最终我还是攻克了这个难关。对于自动控制这门专业性非常强的课程，更需要我们打好扎实的基本功。在设计的过程中难免会遇到过各种各样的问题，比如有时候被一些细小的问题挡住了前进的步伐，有时一开始的计算错误会导致后面的无法进行，花费了大量时间这上面，还有的问题最后还要查阅其他的书籍才能找出解决的办法。
通过对课题的仔细分析，发现在用时域的方法设计太阳观测控制系统时，需要通过测速反馈控制来改善系统的性能从而使所设计的控制系统的时域性能指标满足要求。用频域地方法设计太阳观测控制系统时，通过改变原系统的开环增益可以使该系统在单位斜坡输入下产生的稳态误差满足要求，但不能增大系统的相角裕量，无法满足系统在频域下的设计要求。因此需对原系统增加串联超前校正来增大系统的相角裕量，而串联超前校正不改变系统在单位斜坡输入下产生的稳态误差，从而使所设计的控制系统满足频域下的性能指标。
通过这次自动控制课程设计，加深了对课本知识的理解，进一步熟练掌握了基础理论知识，对平时较难掌握的知识有了更深的了解。熟练运用了理论知识分析课题，并完成设计要求。这次课程设计终于顺利完成了，在设计中遇到了很多MATLAB编程问题通过和同学的讨论和参阅资料书，问题终于得以解决。在此，对给过我帮助的所有同学和指导老师表示忠心的感谢！

参考文献
[1]胡寿松.《自动控制原理》（第五版). 北京: 科学出版社，2001.
[2]王万良.《自动控制原理》[M]. 北京: 科学出版社，2008.
[3]张平.《MATLAB基础应用简明教程》[M]. 北京:北京航空航天大学，2001.
[4]林雪松.《MATLAB7.0应用集锦》[M]. 北京:机械工业出版社，2005.
[5]刘卫国.《MATLAB程序设计教程》[M]. 北京:水利水电出版社，2005.

15
武汉理工大学《自动控制原理》课程设计
本科生课程设计成绩评定表
姓 名
刘爽
性 别
男
专业、班级
电气1205班
课程设计题目：太阳观测控制系统设计
课程设计答辩或质疑记录：

成绩评定依据：
评 定 项 目
最高分限
评分成绩
1．选题合理、目的明确
10

2．设计方案正确、具有可行性、创新性
20

3．设计结果（例如：系统设计程序、仿真程序）
20

4．态度认真、学习刻苦、遵守纪律
15

5．设计报告的规范化、参考文献充分（不少于5篇）
10

6．答辩
25

总 分
100

最终评定成绩（以优、良、中、及格、不及格评定）

指导教师签字：

年 月 日目 录
第一章 总 论 1
1.1 项目概要 1
1.2 可行性研究报告编制依据 1
1.3 项目区简介及建设单位概况 2
1.4 可研报告研究内容 6
1.5 可研报告研究结论、问题及建议 7
第二章 项目背景及建设的必要性和可行性 9
2.1 建设背景 9
2.2 项目建设的必要性 13
2.3 项目建设的可行性 15
第三章 项目建设内容及规模 17
3.1 项目建设内容 17
3.2 建设规模 17
第四章 场址选择及建设条件 19
4.1 场址现状 19
4.2 建设条件 19
第五章 工程方案 24
5.1 方案设计指导思想及原则 24
5.2 建筑设计 24
5.3 结构设计 29
5.4 给排水设计 33
5.5 暖通设计 36
5.6 电气及弱电设计 39
5.7 消防设计 42
第六章 能源和资源节约措施 45
6.1 能源节约措施 45
6.2 建筑节能具体措施 46
6.3 给排水资源节约措施 46
6.4 电气节能具体措施 46
6.5 供热系统节能技术措施 47
第七章 环境影响评价 48
7.1 环境影响 48
7.2 保护措施 50
7.3 安全保护措施 51
7.4 环境影响评价结论 52
第八章 组织机构 53
8.1 管理机构 53
8.2 项目组织管理 53
第九章 工程管理及实施计划 54
9.1 项目建设管理原则 54
9.2 工程管理 54
9.3 项目实施步骤 55
9.4 项目实施进度计划 55
第十章 劳动安全、卫生与消防 58
10.1 劳动安全与卫生 58
10.2 消防安全 59
第十一章 投资估算及资金筹措 61
11.1 估算依据及内容 61
11.2 投资估算 63
11.3 资金筹措 64
12.1 招标依据 65
12.2 工程招标 65
第十三章 社会影响分析及评价 68
13.1 项目社会评价的依据 68
13.2 建设地区基本的社会环境情况 68
13.3 社会影响分析 68
13.4 项目在建设、运营中的社会风险分析 69
13.5 社会评价结论 70
第十四章 结论及建议 71
14.1 结论 71
14.2 项目建议 71
第十五章 附表、附图、附件 73
15.1 附表 73
15.2 附图 73
15.3 附件 73

18