matlab最少拍有波纹控制器控制系统建模分析设计和仿真.docx
《matlab最少拍有波纹控制器控制系统建模分析设计和仿真.docx》由会员分享，可在线阅读，更多相关《matlab最少拍有波纹控制器控制系统建模分析设计和仿真.docx（20页珍藏版）》请在冰豆网上搜索。
matlab最少拍有波纹控制器控制系统建模分析设计和仿真
北京理工大学珠海学院
课程设计任务书
2010～2011学年第2学期
学生姓名：
专业班级：
08自动化1班
指导教师：
钟秋海工作部门：
信息学院
一、课程设计题目
《控制系统建模、分析、设计和仿真》
本课程设计共列出10个同等难度的设计题目，编号为：
[0号题]、[1号题]、[2号题]、[3号题]、[4号题]、[5号题]、[6号题]、[7号题]、[8号题]、[9号题]。
学生必须选择与学号尾数相同的题目完成课程设计。
例如，学号为82的学生必须选做[2号题]。
二、课程设计内容
（一）《控制系统建模、分析、设计和仿真》课题设计内容
[2号题]控制系统建模、分析、设计和仿真
设连续被控对象的实测传递函数为：
用零阶保持器离散化，采样周期取0.1秒，分别设计一单位加速度信号输入时的最少拍有波纹控制器（z）和一单位速度信号输入时的最少拍无波纹控制器（z）。
具体要求见
（二）。
（二）《控制系统建模、分析、设计和仿真》课题设计要求与评分标准【共100分】
1、求被控对象传递函数G（s）的描述。
（2分）
2、求被控对象脉冲传递函数G（z）。
（4分）
3、转换G（z）为零极点增益模型并按1形式排列。
（2分）
4、确定误差脉冲传递函数（z）形式,满足单位加速度信号输入时闭环稳态误差为零和实际闭环系统稳定的要求。
（6分）
5、确定闭环脉冲传递函数（z）形式，满足控制器（z）可实现、最少拍和实际闭环系统稳定的要求。
（8分）
6、根据4、5、列写方程组，求解（z）和（z）中的待定系数并最终求解（z）和（z）。
（12分）
7、求针对单位加速度信号输入的最少拍有波纹控制器（z）并说明（z）的可实现性。
（3分）
8、用程序仿真方法分析加速度信号输入时闭环系统动态性能和稳态性能。
（7分）
9、用图形仿真方法（）分析单位加速度信号输入时闭环系统动态性能和稳态性能。
（8分）
10、确定误差脉冲传递函数（z）形式,满足单位速度信号输入时闭环稳态误差为零和实际闭环系统稳定的要求。
（6分）
11、确定闭环脉冲传递函数（z）形式，满足控制器（z）可实现、无波纹、最少拍和实际闭环系统稳定的要求。
（8分）
12、根据10、11、列写方程组，求解（z）和（z）中的待定系数并最终求解（z）和（z）。
（12分）
13、求针对单位速度信号输入的最少拍无波纹控制器（z）并说明（z）的可实现性。
（3分）
14、用程序仿真方法分析单位速度信号输入时闭环系统动态性能和稳态性能。
（7分）
15、用图形仿真方法（）分析单位速度信号输入时闭环系统动态性能和稳态性能。
（8分）
16、根据8、9、14、15、的分析，说明有波纹和无波纹的差别和物理意义。
（4分）
三、进度安排
6月13至6月14：
下达课程设计任务书；复习控制理论和计算机仿真知识，收集资料、熟悉仿真工具；确定设计方案和步骤。
6月14至6月16：
编程练习，程序设计；仿真调试，图形仿真参数整定；总结整理设计、仿真结果，撰写课程设计说明书。
6月16至6月17：
完成程序仿真调试和图形仿真调试；完成课程设计说明书；课程设计答辩总结。
四、基本要求
1．学生应按照课程设计任务书的要求独立分析、解决问题，按计划完成课程设计任务；
2．不得抄袭或找人代做，否则按考试作弊处理；
3.学生在完成课程设计时须提交不少于3000字课程设计说明书；说明书结构为：
（1）封面，
（2）任务书，（3）摘要，（4）关键词，（5）目录，（6）正文，（7）参考文献；
教研室主任签名：
年月日
摘要
本次课程的目的是让自动化的同学们学习并熟悉使用计算机软件6.5（部分同学是7.0的版本）去建模、分析、设计和仿真一个控制系统。
控制系统中有的需要使用零阶保持器去将连续输入函数离散化，有的则要求使用一阶保持器，从而设计一单位加速度信号输入时的最少拍有波纹控制器（z）和一单位速度信号输入时的最少拍无波纹控制器（z）。
在这个研究过程当中，同学们要先重温《自动控制理论》和《计算机仿真》两门关键学科，才能在使用时减少许多不解与疑问。
设计过程当中，需要将G（S）进去Z变换，按要求实现1形式的展开，以与设定能满足于最少拍数字无波纹（有波纹）控制器的（Z）和（Z），列出方程组，解出待定系数，再求得所需的数字控制器。
得出的误差脉冲传递函数和闭环脉冲传递函数在后边工作中还需要进去多次调整，从而获得最佳表达形式。
最后分别使用程序仿真方法和去分析系统在加速度和速度两种输入信号下的动态性能和稳态特性，比较有波纹和无波纹的差别并理解物理意义。
在输出的波形，理解曲线进入稳定所需的时间，即最小周期，一般是小于等于3T。
程序仿真中注意对和两个指令的使用，否则对同时展示多个参数曲线有所影响。
中选择正确适当的程序框图拖入到新建的模型文件，执行对各部件的编辑，纠正程序框图中存在的错误，输出系统所求波形。
关键词：
仿真建模离散化零阶保持器波纹最小拍Z变换待定系数模型
一、设计任务书1
二、摘要4
三、关键词4
四、关于与控制系统
1.该仿真软件的特点6
2.零阶保持器6
3.最小拍数字控制器6
五、最小拍系统设计
1.单位加速度输入有波纹6
2.单位速度输入无波纹7
六、题目详解与程序指令8
七、原程序指令的运行结果14
八、课程设计的体会和发生的问题16
九、参考文献17
7.0软件的特点：
1）高效的数值计算与符号计算功能,能使用户从繁杂的数学运算分析中解脱出来
2）具有完备的图形处理功能,实现计算结果和编程的可视化
3）友好的用户界面与接近数学表达式的自然化语言,使学者易于学习和掌握
4）功能丰富的应用工具箱（如信号处理工具箱、通信工具箱等）,为用户提供了大量方便实用的处理工具
零阶保持器：
（），它的作用是在信号传递过程中，把第时刻的采样信号值一直保持到第
（1）T时刻的前一瞬时，把第
（1）T时刻的采样值一直保持到
（2）T时刻，依次类推，从而把一个脉冲序列变成一个连续的阶梯信号。
因为在每一个采样区间内连续的阶梯信号的值均为常值，亦即其一阶导数为零，故称为零阶保持器。
最少拍数字控制器：
最少拍，即让控制器的输出在最少的周期内达到稳定。
数字控制器，电子控制器的一类，计算机控制系统的核心部分，与系统中反馈部分的元件、设备相连，该系统中的其他部分可能是数字的也可能是模拟的。
数字控制器通常是利用计算机软件编程，如本次课程设计中使用的7.0，从而完成特定的控制算法。
最小拍系统设计
1.单位加速度有波纹
M文件指令：
668*（[12],[16]）;
（[100]（[11]（[15],[18]）））;
（）%求得
2d（,0.1,''）%用零阶保持器离散化，求得
[abk]（）%零极点增益模型
（,0.1,'','z^-1'）%按1展开
za0a1a2b0
0.096757*z^-1*（1+3.252*z^-1）*（1-0.8187*z^-1）*（1-0.5488*z^-1）*（1+0.2281*z^-1）/（1^-1）^2/（1-0.9048*z^-1）/（1-0.6065*z^-1）/（1-0.4493*z^-1）
^-1*（1+3.252*z^-1）*（a01*z^-12*z^-2）
（1^-1）^3*（10*z^-1）%初步定出形式
f1（,1）-1
f2（（,1）,1）
f3（（,2）,1）
[a0ja1ja2j]（f123）
（[a0ja1ja2j]）
f4（3.252）-1
b0（f4）
（b0j）%解出待定系数
（0）
（,[a0a1a2]）%求得和
%求得
[]（（））;
2（N）
2（D）
[]（（））;
2（N）
2（D）
0:
0.1:
1
.*2
（）
（）
[]（（））;
2（N）
2（D）
2.单位加速度有波纹
M文件指令：
c0c1c2c3c4d0d1
（1^-1）^2*（c01*z^-12*z^-23*z^-34*z^-4）%初定和
^-1*（1+3.252*z^-1）*（1-0.8187*z^-1）*（1-0.6065*z^-1）*（1-0.4493*z^-1）*（d01*z^-1）
f6（,1）-1
f7（（,1）,1）
[d0jd1j]（f67）
（,[d01],[d01j]）%求得
f8（,2）（,2）-1
f9（3.773）-1
f10（,0.6703）-1
f11（,0.5488）-1
f12（,0.3679）-1
[c01234j]（f89101112）
（,[c01234],[c01234j]）%求得
0.096757*z^-1*（1+3.252*z^-1）*（1-0.8187*z^-1）*（1-0.5488*z^-1）*（1+0.2281*z^-1）/（1^-1）^2/（1-0.9048*z^-1）/（1-0.6065*z^-1）/（1-0.4493*z^-1）
%求得
[]（（））;
2（N）
2（D）
[]（（））;
2（N）
2（D）
0:
0.1:
2
（）
（）
[]（（））;
2（N）
2（D）
题目详细解答如下：
1、求被控对象传递函数G（s）的描述。
（2分）
668*（[12],[16]）;
（[100]（[11]（[15],[18]）））;
（）
2、求被控对象脉冲传递函数G（z）。
（4分）
2d（,0.1,''）
3、转换G（z）为零极点增益模型并按1形式排列。
（2分）
[abk]（）
（,0.1,'','z^-1'）
4、确定误差脉冲传递函数（z）形式,满足单位加速度信号输入时闭环稳态误差为零和实际闭环系统稳
定的要求。
（6分）
za0a1a2b0
0.096757*z^-1*（1+3.252*z^-1）*（1-0.8187*z^-1）*（1-0.5488*z^-1）*（1+0.2281*z^-1）/（1^-1）^2/（1-0.9048*z^-1）/（1-0.6065*z^-1）/（1-0.4493*z^-1）
（1^-1）^3*（10*z^-1）
5、确定闭环脉冲传递函数（z）形式，满足控制器（z）可实现、最少拍和实际闭环系统稳定的要求。
（8分）
^-1*（1+3.252*z^-1）*（a01*z^-12*z^-2）
6、根据4、5、列写方程组，求解（z）和（z）中的待定系数并最终求解（z）和（z）。
（12分）
f1（,1）-1
f2（（,1）,1）
f3（（,2）,1）
[a0ja1ja2j]（f123）
（[a0ja1ja2j]）
f4（3.252）-1
b0（f4）
（b0j）
（0）
（,[a0a1a2]）
7、求针对单位加速度信号输入的最少拍有波纹控制器（z）并说明（z）的可实现性。
（3分）
8、用程序仿真方法分析加速度信号输入时闭环系统动态性能和稳态性能。
（7分）
[]（（））;
2（N）
2（D）
[]（（））;
2（N）
2（D）
0:
0.1:
1
.*2
（）
（）
图为程序仿真下有波纹控制器的波形
9、用图形仿真方法（）分析单位加速度信号输入时闭环系统动态性能和稳态性能。
（8分）
[]（（））;
2（N）
2（D）
结构图如下：
图为以上程序框图在仿真下的波形：
10、确定误差脉冲传递函数（z）形式,满足单位速度信号输入时闭环稳态误差为零和实际闭环系统稳定的要求。
（6分）
c0c1c2c3c4d0d1
（1^-1）^2*（c01*z^-12*z^-23*z^-34*z^-4）
11、确定闭环脉冲传递函数（z）形式，满足控制器（z）可实现、无波纹、最少拍和实际闭环系统稳定的要求。
（8分）
^-1*（1+3.252*z^-1）*（1-0.8187*z^-1）*（1-0.6065*z^-1）*（1-0.4493*z^-1）*（d01*z^-1）
12、根据10、11、列写方程组，求解（z）和（z）中的待定系数并最终求解（z）和（z）。
（12分）
f6（,1）-1
f7（（,1）,1）
[d0jd1j]（f67）
（,[d01],[d01j]）
f8（,2）（,2）-1
f9（3.773）-1
f10（,0.6703）-1
f11（,0.5488）-1
f12（,0.3679）-1
[c01234j]（f89101112）
（,[c01234],[c01234j]）
13、求针对单位速度信号输入的最少拍无波纹控制器（z）并说明（z）的可实现性。
（3分）
0.096757*z^-1*（1+3.252*z^-1）*（1-0.8187*z^-1）*（1-0.5488*z^-1）*（1+0.2281*z^-1）/（1^-1）^2/（1-0.9048*z^-1）/（1-0.6065*z^-1）/（1-0.4493*z^-1）
14、用程序仿真方法分析单位速度信号输入时闭环系统动态性能和稳态性能。
（7分）
[]（（））;
2（N）
2（D）
[]（（））;
2（N）
2（D）
0:
0.1:
2
（）
（）
图为程序仿真下无波纹控制器的波形
15、用图形仿真方法（）分析单位速度信号输入时闭环系统动态性能和稳态性能。
（8分）
[]（（））;
2（N）
2（D）
结构图如下：
图为以上程序框图在仿真下的波形
16、根据8、9、14、15、的分析，说明有波纹和无波纹的差别和物理意义。
（4分）
有波纹和无波纹的差别在于有波纹控制器经过一定的采样周期之后，其输入跟踪输出，但其数字控制器的曲线会出现小的波动，该波动为波纹，而无波纹数字控制器D（z）的输出曲线经过一定采样周期之后会为一常数，不会产生波纹。
最少拍控制系统的系统结构简单，设计完成后容易在计算机上实现。
但它对系统的适应性差，只能保证在采样点上的输出可以跟踪输入，有可能在采样点之间呈现波纹，会引起系统的振荡。
产生波纹的原因是U（z）不能在优先个采样周期内变为0，即u（）不等于零，使系统的输出y（t）产生波动。
实现最少拍无波纹的控制，必须要系统在典型信号的作用下，经有限个采样周期后，系统的稳态误差保持恒值或为0，系统的数字控制器D（z）的输出u（）也必须保持为0。
原程序指令运行结果
:
668s^2+5344s+8016
s^5+14s^4+53s^3+40s^2
:
0.09676z^4+0.2044z^3-0.3452z^2+0.05314z+0.03225
z^5-3.961z^4+6.149z^3-4.663z^2+1.721z-0.2466
:
0.1
a=[4x1]
b=[5x1]
k=0.0968
:
0.096757z^-1（1+3.252z^-1）（1-0.8187z^-1）（1-0.5488z^-1）（1+0.2281z^-1）
（1^-1）^2（1-0.9048z^-1）（1-0.6065z^-1）（1-0.4493z^-1）
:
0.1
69723/727936*（1+813/250）*（1-8187/10000）*（1-343/625）*（1+2281/10000）/（1-1）^2/（1-1131/1250）/（1-1213/2000）/（1-4493/10000）
=1*（1+813/250）*（a012^2）
=（1-1）^3*（10）
f1=1063/250*a0+1063/250*a1+1063/250*a2-1
f2938/125*a0-2939/250*a1-2001/125*a2
f3=2689/125*a0+5628/125*a1+1926/25*a2
a0j=1444828500/1201157047
a1j1826125500/1201157047
a2j=663789250/1201157047
A=1.2029-1.52030.5526
f4=663789250/5373661750/436880018961*b0
b0j=2158642641/1201157047
B=1.7971
=（1-1）^3*（1+85687/45496）
1*（1+813/250）*（5417217599035801/4564948851/4549649776/90992^2）
727936/69723*（5417217599035801/4564948851/4549649776/90992^2）/（1-8187/10000）/（1-343/625）/（1+2281/10000）/（1-1）*（1-1131/1250）*（1-1213/2000）*（1-4493/10000）/（1+85687/45496）
=1.0018*
2.70865.3850-9.88864.0468
=1.0018*
2.25180000
=1.0028*
0.5417-2.83026.2917-7.72175.6551-2.47230.5971-0.0614
=1.0026*
4.3575-4.96500.59860.44660000
=1.0044*
0.2440-0.78671.0162-0.65850.2136-0.0276
=1.0042*
1.9625-0.6717-5.03964.4345-0.3242-0.3614
=（1-1）^2*（c012^23^34^4）
1*（1+813/250）*（1-8187/10000）*（1-1213/2000）*（1-4493/10000）*（d01）
f6=8352584490971/500*d0+8352584490971/500*d1-1
f7=8533117278079/100*d0+17/00*d1
d0j172000/697656676788841
d1j=297500/697656676788841
f8=1/4*c0+1/8*c1+1/16*c2+1/32*c3+1/64*c4-2757/562949953421312
f9=22781529/14235529*c0-22781529000/53710650917*c1+22781529000000/2841*c2-227800/764599528737830093*c3+227800000/28848340889*c4-1
f10=10870209/44930209*c0+1/3*c1+10000/2681*c2+100000000/293013743*c3+10/960071119329*c4-1
f11=79524/117649*c0+49702500/40353607*c1+31064062500/*c2+194/4747561509943*c3+122500/16284*c4-1
f12=39955041/13535041*c0+399550410000/49795415839*c1+39955/1681*c2+399550000/673982994992914399*c3+3995500000000/2479583438578932073921*c4-1
c39985825724072291/1416
c116118126287252682/47253242735819440591601664
c227487343919321963360445596793/3366661984722944
c325629/752712953061376
c46674674584926934716646131/14394968429824
=1.0018*
2.70865.3850-9.88864.0468
=1.0018*
2.25180000
=1.0028*
0.5417-2.83026.2917-7.72175.6551-2.47230.5971-0.0614
=1.0026*
4.3575-4.96500.59860.44660000
=1.0095*
-0.56872.4226-4.16353.7060-1.80720.4589-0.0475
=1.0094*
0.32730.6583-1.68011.0274-0.1459-0.05130.0130
我的体会：
设计最小拍数字控制系统是之前的《计算机控制系统》中学习过的知识，我基本能够熟练地掌握，但由于我之前从未学习过软件，一开始连的基本的功能和函数都根本不清楚，要用设计控制系统并进行仿真，对我是一个非常大的挑战。
但我相信，只要认真查阅相关书籍学习，虚心向老师和同学请教，一定能很快地学会使用软件的基本功能。
我借来一本《仿真技术》，从最基础的了解菜单界面窗口、建立M文件学起，到学习用它解方程，表达控制系统的传递函数，再到使用一些常用的函数命令，比如离散化、建立零极点增益模型等。
一步一步地亲身实践尝试，慢慢就掌握了这些基本的软件操作能力。
在这一基础上，我总于可以用软件设计最小拍控制系统。
通过利用待定系数，到实现程序仿真，到最后形象的图形仿真，我意识到软件的强大，它可以大大方便我们自动化专业的学生在控制系统设计中的计算和设计。
此外课程设计使我基本掌握了一个很有用的软件，也使我更加意识到学习知识要有信心，从基础打起，一点一点不断累积和实践尝试。
多和同学探讨在学习中遇到的错误和经验，也能很好提高学习效率，少犯错误。
这是一次很有意义的实践，我会把收获用在以后的学习中,不断努力学习新的知识！
发生的问题与解决方法：
1、系统传递函数的表达式不知道用什么命令表达。
解决方法：
翻阅资料，掌握了常用的几种表达方法。
2、在待定和的参数时，列出的方程有错误而导致无解。
解决方法：
利用的错误提示功能，能方便地找到错误并改正。
3、程序仿真时，图像曲线只看到整体曲线，局部的微小区间部分不够清晰、具体。
解决方法：
改变设定中终点步，可改变整体的曲线分辨率。
如将0:
0.1:
10改为0:
0.1:
1，输入，则局部部分非常具体清晰。
4、在图形仿真时，D（Y）或D（W）的参数不知道如何设定。
解决方法：
翻阅资料，用（）化简，用（）函数和2（）函数提取系数，可方便地设定D（Y）或D（W）的参数.。
5、图形仿真时，响应曲线出现震荡。
解决方法：
经过细心检查，发现是太粗心设置了错误的传递函数参数导致，重新修正后系统响应曲线满足要求。
参考文献
【1】魏克新编《语言与自动控制系统设计》机械工业出版社2004年
【2】周品编《数学建模与仿真》国防工业出版社 2009年
升级会员 