《计算机控制技术》课程设计设计题目:最少拍无纹波计算机控制系统设计及仿真实现学院:机械与电子工程专业:自动化班级:1121502学号:************姓名:***指导老师:***目录引言 (1)1、课题介绍 (1)1.1、课程设计目的 (1)1.2、课程设计内容 (2)2、课题论证 (2)2.1最少拍计算机控制系统 (3)2.2、最少拍有纹波控制系统 (3)2.3 最少拍无纹波控制系统 (4)3 系统设计 (5)3.1、单位阶跃信号输入时的最少拍无纹波设计 (5)3.2单位速度信号输入时的最少拍无纹波设计 (10)4 结果分析 (16)4.1有纹波系统和无纹波系统的对比 (16)4.2最少拍无纹波控制系统对典型输入的适应性问题 (16)4.3 K的变化对系统的影响 (18)d5 设计总结 (20)6参考文献 (20)引言《计算机控制技术》是一门理论性、实用性和实践性都很强的课程,课程设计环节就占有更加重要的地位。
本次课程设计通过对最少拍无纹波控制器的设计及仿真了解和掌握了最少拍无纹波设计和有纹波设计。
首先,通过学习和搜集资料了解和掌握了最少拍控制器的设计原理,设计了最少拍无纹波控制器,通过MATLAB程序验证了正确性。
最少拍系统设计是以采样点上误差为零或保持恒定为基础,采用Z变换方法进行设计并保证采样点之间的误差也为零或保持恒定值,因此在采样点之间可能存在波纹,即在采样点之间有误差存在,这就是有波纹设计。
而无波纹设计是指在典型输入信号的作用下,经过有限拍系统达到稳定,并且在采样点之间没有波纹,输入误差为零。
即要求采样点之间产生的波纹不能反映在采样点信号上,也就是对采样点之间的信号,不能形成闭环控制。
要得到无波纹系统设计,其闭环Z传递函数φ(z)必须包含被控制对象G(z)的所有零点。
设计的控制器D(Z)中消除了引起纹波振荡的所有极点,采样点之间的波纹也就消除了。
系统的闭环Z传递函数φ(z)中的Z-1的幂次增高,系统的调整时间ts就增长。
本文以实例来介绍最少拍无波纹控制的实现方法。
关键词:计算机控制技术最少拍无纹波 MATLAB单位阶跃信号单位速度信号1、课题介绍1.1、课程设计目的1)学习并掌握有纹波最少拍控制器的设计和Simulink 实现方法; 2)研究最少拍控制系统对典型输入的适应性及输出采样点间的纹波; 3)学习并掌握最少拍无纹波控制器的设计和Simulink 实现方法;4) 研究输出采样点间的纹波消除方法以及最少拍无纹波控制系统对典型输入的适应性。
5)编写算法MATLAB/simulink 仿真程序实现(模拟步进电机进给过程); 6)撰写设计说明书。
1.2、课程设计内容下面以一个具体实例介绍最少拍系统的设计和仿真。
如图1所示的采样-数字控制系统,被控对象: 120()(1)(2)dK G s T s T s =++ (1-1)零阶保持器: 1()Tse H s s--= (1-2)K d =2、T 1=1、T 2=1,选择采样周期T=1s , 试设计无纹波最少拍控制器,并分析仿真结果1. 分别在单位阶跃/单位速度输入下设计无纹波有限拍控制器2. 在Simulink 仿真环境画出仿真框图及得出仿真结果,画出数字控制器和系统输出波形。
3. 与有纹波系统进行对比分析(选用单位速度输入进行对比分析即可)4. 探讨最少拍无纹波控制系统对典型输入的适应性问题5.分析K d 的变化对系统的影响2、课题论证2.1最少拍计算机控制系统最少拍设计,是指系统在典型输入信号(如阶跃信号、速度信号、加速度信号等)作用下,经过最少拍(有限拍)使系统输出的系统稳态误差为零。
因此,最少拍控制系统也称最少拍无差系统或最少拍随动系统,它实质上是时间最优控制系统,系统的性能指标就是系统调节时间最短或尽可能短,即对闭环Z传递函数要求快速性和准确性。
图2-1 离散系统模型数字控制器模拟化设计方法是基于连续系统的设计,并在计算机上采用数字模拟方法来实现,选用的采样周期须足够小,且采样周期的变化对系统影响不大。
如图1 的数字离散控制系统中,GC(S)为被控对象,其中H(S)= (1-e-TS)/S 代表零阶保持器,D(Z)代表被设计的数字控制器,它是由单片机来实现的,D(Z)的输入输出均为离散信号。
上面图1示计算机控制系统框图中:Gc(s) ——被控对象的连续传递函数D(z) ——数字控制器的Z传递函数H(s) ——零阶保持器的传递函数,T ——采样周期广义对象的脉冲传递函数为:G(z)=Z[H(s)×G(s)] (2-1)2.2、最少拍有纹波控制系统数字控制器原理方框图如图所示图 2-2数字控制器原理方框图2.3 最少拍无纹波控制系统2.3.1产生纹波的原因:产生纹波的原因在于数字控制器u(k)经若干拍数后,不为常值或零,而是震荡收敛的。
它的作用在于保持器的输入端,保持器的输出必然也波动,使系统输出在采样点之间产生波纹。
2.3.2消除纹波的附加条件:最少拍无纹波系统设计,是在最少拍控制存在纹波时,对期望闭环响应φ(z)进行修正,以达到消除采样点之间波纹的目的。
无纹波系统的调整时间比有纹波系统的调整时间增加若干拍,增加的拍数等于G(z)在单位圆的零点数目。
2.3.3设计分析a)准确性要求。
系统对某种典型输入, 在采样点上无稳态误差,对特定的参考输入信号在到达稳态后系统输出在采样点的值准确跟踪输入信号即采样点上的输出不存在稳态误差。
b)快速性要求。
闭环系统过渡过程最短, 即最少采样点数内使采样点上稳态误差趋于零.即在各种使系统在有限拍内到达稳态的没计中系统准确跟踪输入量所需的采样周期数应为最少。
c)稳定性要求。
系统输出在采样点上不发散、不振荡, 且采样点之间也不能发散, 当广义对象G(Z) 含单位圆上或圆外零点或极点时, 前面两步设计出的φ(z),不能保证稳定性要求.数字控制器必须在物理上可实现且应该是稳定的闭环系统。
在采样点上的输出不存在稳态误差,但在采样点间的输出存在稳态误差的系统为有波纹最少拍控制系统。
3 系统设计3.1、单位阶跃信号输入时的最少拍无纹波设计⑴ 最少拍有纹波设计① 用matlab 求出该系统的广义脉冲传递函数G=tf([2],[1 3 2]) Transfer function: 2 ------------- s^2 + 3 s + 2 G1=c2d(G,1,'zoh') Transfer function: 0.3996 z + 0.147 ------------------------ z^2 - 0.5032 z + 0.04979 Sampling time: 1即: )1356.01)(3676.01()368.01(3996.0)(1111------+=z z z z z G (3-11)② z z e z 1)(1)(-=-=φφ (3-12))11()(--=z z e φ (3-13)③ )1368.01)(11(13996.0)11356.01)(13676.01(1)()()(1)(z z z z z z Z e Z G z e Z D -+--------=-=φφ (3-14) 即:2147.012526.03996.0204979.015032.01)(-----+--=ZZ Z Z Z D (3-15)④根据求出的结果进行仿真仿真原理图:图3-11 阶跃输入最少拍有纹波设计仿真原理图系统输出:图3-12 阶跃输入最少拍有纹波系统输出控制器输出:图3-13阶跃输入最少拍有纹波控制器输出误差输出:图3-14阶跃输入最少拍有纹波误差输出由图3-12、图3-13、图3-14可知:控制器输出不稳定,将使系统出现纹波⑵、最少拍无纹波设计① 求该系统的广义脉冲传递函数)1356.01)(3676.01()368.01(3996.0)(1111------+=z z z z z G (3-16)② 因G(Z)零点Z=-0.368,极点1356.0,3676.021==P P 。
闭环脉冲传递函数)(z Φ应选为包含G(Z)全部零点,所以:)1368.01(1)(1)(-+-=Φ-=Φz az z ez (3-17))(z eΦ由)(z G 的不稳定极点和)(z Φ的阶次来决定。
所以选择: )11)(11()(-+--=Φbz z z e(3-18) 因)(z eΦ=1-)(z Φ,将上述所得的)(z Φ和)(z e Φ的值代入后,可得)1368.01(11)11)(11(-+--=-+--z az bz z (3-19) 所以解得269.0,731.0==b a③ 求得控制器D(z)的传函为:)1269.01)(11()11356.01)(13676.01(83.1)()()(1)(-+------=ΦΦ-=z z z z Z e Z G z e Z D (3-20) 即:2269.01731.012091.0192.083.1)(-----+--=ZZ Z Z Z D (3-21)④根据求出的结果进行仿真仿真原理图:图3-15阶跃输入最少拍无纹波系统仿真仿真结果图:图3-16阶跃输入最少拍无纹波系统输出图3-17阶跃输入最少拍无纹波控制器输出图3-18阶跃输入最少拍无纹波误差输出由仿真结果可知:控制量稳定,将使系统输出无纹波。
无纹波系统的调整时间比有纹波系统的调整时间增加一拍,增加的拍数正好等于G(Z)在单位圆内的零点数目。
3.2单位速度信号输入时的最少拍无纹波设计⑴、最少拍有纹波设计① 求该系统的广义脉冲传递函数)1356.01)(3676.01()368.01(3996.0)(1111------+=z z z z z G (3-22)② Z Z z e z 212)(1)(---=-=φφ (3-23))11(2)(--=Z z eφ (3-24)③ )1z (12)z 10.368(1z 10.3996)11356.01)(13676.01)(212()()()(1)(---+--------=-=z z z z Z eZ G z e Z D φφ(3-25)即:Z Z Z Z Z Z Z D 3147.021056.016552.03996.0304979.026028.01006.22)(-+-+-----+--= (3-26)④ 根据求出的结果进行仿真仿真原理图:图3-19 速度输入时的有纹波系统仿真图仿真结果图:图3-20速度输入有纹波系统的输出图3-21速度输入有纹波系统的控制器输出图3-22速度输入有纹波系统的误差波形 由仿真结果可知:输出波形有纹波。
⑵、最少拍无纹波设计① 求该系统的广义脉冲传递函数)1356.01)(3676.01()368.01(3996.0)(1111------+=z z z z z G (3-27)② 因G(Z)有1-Z 零点Z=-0.368,极点1356.0,3676.021==P P 。