最小拍系统matlab仿真(1)

1、10/11学年第一学期学院：专业：学生姓名：学号：课程设计题目：起迄日期：课程设计地点：指导教师：系 主任：下达任务书日期：2011年12月课程设计任务书1 .设计目的：1至,系第曩番薮菠羲方箱转编彝西簧开箱外而行藏羹市言东2 .研究最少拍控制系统在单位速度信号的作用下的适应性及输出采样点间的纹波2.设计内容和要求(包括原始数据、技术参数、条件、设计要求等):亘豆二不索获血禾国薪宗：: RO D (z) Gh(s) f G(s) /_其中O(z)离散控制器，Gx(s)是零阶保持器，受控对象为G(s)=也一式分别求出这个系统在 s(s + 1)单位速度信号作用下的数字控制器D(Z)和响应曲线，取

2、采样周期T=l。实验要求：1、所采用方法的基本思路和计算公式。2、将所得结果用SIMULINK进生仿真，并将输入，输出和误差三条曲线放置在同一图像内，比较它们的 跟随特性。3、仿真过程及仿真结果分析(1)系统的动态响应在几个采样周期内结束？(2)系统的稳态响应有无波纹，D(z)是什么控制器？3 .设计工作任务及工作量的要求(包括课程设计计算说明书(论文)、图纸、 实物样品等)：(1)查阅资料，确定系统设计方案(2)说明设计原理，构建系统,SIMULINK仿真(3)撰写课程设计说明书.4 .主要参考文献:1夏扬.计算机控制技术M.北京：机械工业出版社2顾德英，张健，马淑华.计算机控制技术M.北京

3、：北京邮电大学出版社3刘金琨.先进PID控制及其MATLAB仿真M.北京：电子工业出 版社4张艳兵，王忠庆，鲜浩，计算机控制技术M.北京：国防工业出 版社5黄一夫，微型计算机控制技术M,北京：机械工业出版社5 .设计成果形式及要求：一谏询吾祝丽吊T羽设计原理分析1份仿真结果1份.6 .工作计划及进度：2011年12月5日收集整理资料，确定控制方案12月6 012月7日 SIMULINK进行仿真.12月8日撰写课程设计说明书12月9日答辩系主任审查意见：签字：年 月 日一、基本原理最少拍设计，是指系统在典型输入信号（如阶跃信号、速度信 号、加速度信号等）作用下，经过最少拍（有限拍）使系统输出的系

4、 统稳态误差为零。因此，最少拍控制系统也称为最少拍无差系统或最 少拍随动系统，它实质上是时间最优控制系统，系统的性能指标就是 系统的调节时间最短或尽可能短，即对闭环Z传递函数要求快速性和 准确性。6b (s)被控对象的连续传递函数(z)数字控制器的Z传递函数(s)零阶保持器的传递函数，T采样周期广义对象的脉冲传递函数为：G(z)=Z H(S).GC(5)系统闭环脉冲传递函数为：6 (z)C= O(z)G(z) 而 l + D(z)G(z)系统误差脉冲传递函数为:6e(z)=名=1-小(z) =R1 + D(z)G(z)数字控制器脉冲传递函数:D(z)U二1 中 £(77 GU) l-

5、(b(z)若已知Gc(s),且可根据控制系统的性能指标要求构造(z), 则根据G(z)=(1-z Z跑和 D(z)二侬=5sE(z) G(z)l 6(z)J11、闭环Z传递函数(z)的确定；由图1：误差E(z)的Z传递函数为：(Mz)二华工1-6 (z)二1 R l + O(z)G(z)数字控制器脉冲传递函数:D(z)匹曲LE G(z) 1-奴 z)从上式看出，(力的求取主要取决于小(z),或者6e(z), 6(z) 的选择根据稳、准、快等指标设计。下面分析闭环传递函数6 (z)的 确定原则。1)由物理可实现性确定将(Z)写成分子分母关于Z-1有理多项式次暴相除的形式，即Z传递函数物理可实现的

6、条件是分子关于z"的基次低于分母关于一1 的嘉次，即rS。故(z)的展开式中不会出现z的正幕次项z d,其 展开式的最低次号为z若G(z)不含纯滞后环节。因为：6(z)二 2 QG'R l + O(z)G(z)所以6(z)应具有如下标准形式(Z)二iZ、2Z-2+中.'若G(z)包含纯滞后环节zd的因子。O (z)=(】z"+2Z'疝z-d2)由系统的准确性确定根据系统在采样点对稳态误差为0的要求确定6 (z),由得E(z) = R，=1 一R(z)可见误差E的大小与输入信号有关。典型输入信号有单位阶跃、单位速度、单位加速度。单位阶跃输入时r(t)=

7、l(t) ,R(z)=IT单位斜坡输入时r(t)=t ,(T)2单位加速度输入时r(t) = lt2 , R(t) Jfa + z：)22(15旅综合三种典型输入函数(单位阶跃、单位速度、单位加速度)R(z) =8(z)(T)gq=l、2、3, B(z)是不含(1 z ')因子的z1多项式, 阶次为q-lo根据Z变换的终值定理，系统的稳态误差为由于6(z)不含(1-z1)因子，因此要使稳态误差。(8)二0,必 须有(t)e(z)=(l-Z-，rFU)其中：尸= 1 + /1Z，+ f2Z-2 +. + fPz-p而由中，=1 一得=1-,所以：(t)e(z)=(i-r,rF(z)(z)

8、=iZ“+2Z，衣一' (z”的最高嘉次A'=p + q)3)由系统的快速性确定中(z)式 (Z)=l (1 - zT)"F(Z)= 1Z“+ 2Z+中 NZ ' 表明：系统闭环响应在采样点的值经N拍采样周期后可达到稳态误差 为0,要使误差尽快为0,应使表达式中拍数力最小。所以当叶0, 即尸(z)=l时，系统在采样点的输出可在最少拍(Afnin=g拍)内达 到稳态，即为最少拍控制。最少拍控制器设计时，闭环Z传递函数及误差传递函数为：4>e(z)=(i-r,r(z)=i-(i-rIr 最少拍控制器： D0)_u-1(t”e g(z)1-6 Ga)(i-z_

9、,r q=l、2、3,取决于输入信号的类型 2、最少拍控制器的可实现问题。前面根据最小拍定义确定的闭环误差脉冲传递函数只适用于对象不含有纯滞后环节的系统。实际中，很多对象都含 有纯滞后，为了使设计的控制器在物理上可实现，需对设计加以限制。 使闭环脉冲传递函数的零点包含纯滞后环节中(z)=z4 -(l一尸)” 3、最少拍控制的稳定性问题。前面根据最小拍定义确定的系统传递函数可使系统过渡过程实 现最小拍。但上述结论仅适用于被控对象G(z)是非常特殊的场合： G(z)的极点和零点都在单位圆内且不含纯滞后环节，即是说只有当被 控对象的零极点都在单位圆内且不含纯滞后环节时，前面的结论才正 确.实际中，被

10、控对象G(z) =此可能不满足上面条件， B也就是说它可能含有Z平面单位圆上和圆外的零极点，这时利用前面 结论来确定控制器就是错误的。为了保证闭环系统稳定，4)(Z)和小e(z)都不应含有单位圆上或 单位圆外的极点。从上式看出，对G(z)中位于单位圆外或圆上的极点，只能采用(力的零点来对消G(z)的不稳定极点,用6e(z的零点来对消G(z) 的不稳定极点。如果用(z)的零点来简单对消G(z)的不稳定极点，虽然从理论 上可以得到一个稳定的闭环系统，但是这种稳定是建立在零极点的完 全对消的基础上，当系统的参数产生漂移，或辨识参数有误差时，这 种绝对的对消不可能实现，从而引起系统的不稳定。因此，G(

11、z)的 不稳定极点不能用(z)的零点来对消。为了保证闭环系统稳定，只能采用第二种方法来将系统补偿成稳 定系统。因此，在选择6e(z)必须附加稳定约束条件。同理，当G(z)中含有单位圆外或圆上的零点时，由于从上式可以看出：G(z)位于单位圆上或圆外的零点可能成为控 制器(z)的不稳定极点，因此使得对象的输出不稳定。为了确保补 偿以后的系统稳定，用 的零点来对消G(z)圆上和圆外的零 点。因此，当G(z)含有单位圆外或圆上的零极点时，应对前面的 最小拍系统设计结论加两个限制条件。4、最少拍控制器的设计(稳定约束条件)假设Gc(z) = GcS)e-a, GJ(s)假设是不含滞后部分的传递函数；为纯

12、滞后时间，则G(z) = ZG'c (s)e-° = P-Z1-G'c (5) = z" (设 d =) 当 ssA T力0,被控对象不含纯滞后；当力0 ,被控对象含d个采样周期的纯 滞后。前面讨论最小拍控制器的稳定性问题知，只有当开环传递函数G口不含有单位圆外或圆上的零极点时才能利用前面根据最小拍含 义得出的结论。当含有单位圆外(上)的零极点时，为了使补偿 后的闭环系统稳定，应对闭环传递函数和闭环误差传递函数附加一些 约束条件：G(z)单位圆外(上)的全部极点由，的零点来对消；G® 单位圆外(上)的所有零点由(z)的零点来对消。约束条件设G(z)

13、含有：个在z平面的单位圆外或单位圆上的零点力、金、乩；/个在z平面的单位圆上和圆外的极点a、&、其中J个 极点在单位圆上；G，是G(z)中不含单位圆上或圆外的零极点部分。由分析知，选择闭环传递函数的约束条件如下：1)在，的零点中，必须包含G(z)在z平面单位圆外或单位圆上的 所有极点，即："(z)为小的多项式，且不包含G(z)中的不稳定极点a。K=1 +加一+/22"+一 +九广而根据上节最小拍控制的含义确定得我=(1 -综合最小拍控制含义和系统稳定性要求得若jWq文(z) = l-= 其声z"的财教力r-1v-j+q+m若J>q，=1 一(z) =

14、也(1一平"(1一/)/仁)其中的阶数为v+m2)在的零点中，必须包含G(z)在z平面单位圆外或单位圆上的所有零点，并包含滞后环节z -d即= zj(l4zT)F2(z)J=1F21公为力的多项式，且不包含G(z)中的不稳定零点b= f2iz'1 + f21z'2 +. + flnzn 中 Z" 的阶数为 d+u力73) Fl。)、咫(z)阶数的确定由于(z)和我的阶数相等，且有最低嘉次，因此：当m=d+un= v-j+q当 J> (/, HFU+dn=v4) Fl(z)和F2(z)中系数的确定根据：(z) = l-RC)比较系数得到其中参数，从而得到

15、最小拍 控制器如下。根据上述的约束条件设计的最少拍控制系统，保证了在最少的几 个采样周期后系统的响应在采样点时是稳态误差为零；但是，不能保 证任意两个采样点之间的稳态误差为零，也就是说在非采样点有波纹 存在，称此种控制器为最小拍有波纹控制器。二、参数计算该系统广义被控对象的脉冲传递函数G (z) =Z Gc (s) G (s) =Z L1 e 7 10S S (s 4 1)f 1 1 1二10(1-z i) z- + mi3二 10 (l-z 1)【(1 _ z 7) 2 一 1 - Z - 0.3679Z-43.679z 7(1+0. 718z - (1 - z - D (1 + 0. 71

16、8z- ') (1 - z- 1)闭环系统的输出为 12，C (z) =R (z)(z)=Z - 1 (2-Z - 1) N(1 - Z ) <1 - 0.3679z -)针对单位速度输入函数设计最少拍无差系统，m=l, v=l, u=0, q=2, q+v-l=2,因为G (z)包含一个单位圆上的极点z=l, w=l,所以可以 降一阶来处理。由(z)=UZ-nl n (1 - biz - 0 (0 + 逐一1 + + (i = 1)-q - v + w 4 1)q+v-Ll = l,待定系数有三个，选择闭环传递函数为(z) =z - 1 (o+z" 1)由中e(Z)二

17、(1 - Z 一 1) Pfl ：； ； (1 - AjZ 一 1) F (z),其中F (z) =l+f 1Z - 】+f2Z 一 2+fm + L q _ p _ 逐 一 m 一 U - Q P + 1p=2, n)+u+q-p-l=0,选择误差传递函数为e (z) = (l-z1) 2根据待定系数法解得：。=2,i二-1所以，(z) =z - 1 (2-Z - be(Z)= (2故数字控制器的传递函数为0. 5434 (1 - 0. 5z - 1) (1 - 0. 3679z '')系统偏差为2 - 1E (z) =R (z) e (z);777 (l-z-l) 2(1

18、- z )三、系统仿真本题用simulink来仿真最少拍无差系统，因此我们使用Mat lab 中的simulink工具画出其原理图，对输入、输出、误差曲线进行仿 真，观察它的曲线特性。其步骤如下：(1)打开simulink对话框，经历一个可以用来搭建的模型窗口；(2)分别在 Sources 模块库、Continus 模块库、Math Operations模块库、Discrete 模块库、Commonly Used Blocks 模块库、Sinks 模块库中找到step阶跃信号、Transfer fen、ADD、DiscreteTransfer Fen 和 Discrete Zero-Order Hold、Scope 模块;(3)对找到的模块进行参数设置;(4)按要求进行连线并保