1、1,5.1 连续域离散化设计 5.2 数字PID控制器设计 5.3 控制系统z平面设计性能指标要求 5.4 z平面根轨迹设计 5.5 w变换及频率域设计,2,5.1.1 设计原理和步骤,实质是将数字控制器部分看成是一个整体，其输入和输出都是模拟量，因而可等效为连续传递函数De(s)。,A/D输出与输入关系：,系统低通且 采样频率较高,D/A的频率特性：,等效连续 传递函数,计算机实现算法D(z)的计算表示：,设计时常近似为：,图5-1计算机控制系统典型组成,3,连续域-离散化设计的步骤如下：,第1步：根据系统的性能，选择采样频率，并设计抗混叠前置滤波器。 第2步：考虑ZOH的相位滞后，根据系统
2、的性能指标和连续域设计方法，设计数字控制算法等效传递函数Ddc(s)。 第3步：选择合适的离散化方法，将Ddc(s)离散化，获得脉冲传递函数D(z)，使两者性能尽量等效。 第4步：检验计算机控制系统闭环性能。若满足指标要求，进行下一步；否则，重新进行设计。 改进设计的途径有： 选择更合适的离散化方法。 提高采样频率。 修正连续域设计，如增加稳定裕度指标等。 第5步：将D(z)变为数字算法，在计算机上编程实现。,4,5.1.2 各种离散化方法,最常用的表征控制器特性的主要指标： 零极点个数； 系统的频带； 稳态增益； 相位及增益裕度； 阶跃响应或脉冲响应形状； 频率响应特性。,等效离散,D(z)
3、,D(s),数值积分法 一阶向后差法 一阶向前差法 双线性变换法及修正双线性变换法 零极点匹配法 保持器等价法 z变换法(脉冲响应不变法),离散化方法,5,1. 一阶向后差分法,(1)离散化公式,实质：将连续域中的微分用一阶向后差分替换,做z变换，得,s与z之间的变换关系,比较,图5-3向后差分(矩形积分)法,6,1. 一阶向后差分法,(2)主要特性 s平面与z平面映射关系 当=0 （s平面虚轴），s平面虚轴映射到z平面为该小圆的圆周。 当 0（s右半平面），映射到z平面为上述小圆的外部。 当 0（s左半平面），映射到z平面为上述小圆的内部。 若D(s)稳定，则D(z)一定稳定 变换前后，稳态
4、增益不变。 离散后控制器的时间响应与频率响应，与连续控制器相比有相当大的畸变。,图5-4 向后差分法的映射关系,(3) 应用由于这种变换的映射关系畸变严重，变换精度较低。所以，工程应用受到限制，用得较少。,7,2. 一阶向前差分法,(1)离散化公式,实质：将连续域中的微分用一阶向前差分替换,做z变换，得,s与z之间的变换关系,比较,图5-7 向前差分矩形积分法,8,2.一阶向前差分法,(2)主要特性 s平面与z平面映射关系若D(s)稳定，则D(z)不一定稳定改进方法是适当减少采样周期T 。,向前差分法的映射关系图,(3) 应用由于这种变换的映射关系畸变严重，不能保证D(z)一定稳定，或者如要保
5、证稳定，要求采样周期较小，所以应用较少。,平移放大关系,9,3.双线性变换法（突斯汀-Tustin变换法）,(1)离散化公式,实质：将梯形面积近似代替积分,进行z变换，得,s与z之间的变换关系,比较,图5-9 梯形积分法,10,3.双线性变换法,(2)主要特性 s平面与z平面映射关系 当=0（s平面虚轴)映射为z平面的单位圆周。 当 0（s右半平面），映射到z平面单位圆外 。 当 0（s左半平面），映射到z平面单位圆内 。 若D(s)稳定，则D(z)一定稳定,图5-10 双线性变换映射关系,s域角频率,z域角频率为D,11,3.双线性变换法,(2)主要特性 频率畸变：双线性变换的一对一映射，保
6、证了离散频率特性不产生频率混叠现象，但产生了频率畸变。,图5-11双线性变换的频率关系,图5-12双线性变换的频率关系,当采样频率较高,足够小,12,3. 双线性变换法,(2)主要特性 变换前后，稳态增益不变。 双线性变换后D(z)的阶次不变，且分子、分母具有相同的阶次。并有下式成立：,(3) 应用 这种方法使用方便，且有一定的精度和前述一些好的特性，工程上应用较为普遍。 这种方法的主要缺点是高频特性失真严重，主要用于低通环节的离散化，不宜用于高通环节的离散化。,13,4. 修正双线性变换,(1)离散化方法,解决“双线性变换产生频率轴非线性畸变 ”问题的方法。,1是设计者选定的特征角频率,依据
7、连续域与双线性变换后频率的非线性关系，首先修正原连续域传递函数，然后再进行双线性变换的结果。,(2) 主要特性 该方法本质上仍为双线性变换法，因此具有双线性变换法的各种特性。但由于采用了频率预修正，故可以保证在处连续频率特性与离散后频率特性相等，即满足(3) 应用 由于该方法的上述特性，所以主要用于原连续控制器在某些特征频率处要求离散后频率特性保持不变的场合。,14,5. 零极点匹配法,(1)离散化方法,特点： 零、极点分别按 匹配 若分子阶次m小于分母阶次n，离散变换时，在D(z)分子上加(z+1)n-m因子 确定D(z)的增益k1的方法： 按右式来匹配 若D(s)分子有s因子，可依高频段增
8、益相等原则确定增益,即也可选择某关键频率处的幅频相等，即,15,5. 零极点匹配法,(2)主要特性 零极点匹配法要求对D(s)分解为极零点形式，且需要进行稳态增益匹配，因此工程上应用不够方便。 由于该变换是基于z变换进行的，所以可以保证D(s)稳定，D(z)一定稳定。 当D(s)分子阶次比分母低时，在D(z)分子上匹配有(z+1)因子，可获得双线性变换的效果，即可防止频率混叠。,16,6. 其他离散方法,(1) z变换法（脉冲响应不变法）,(2) 带保持器的z变换 带零阶保持器z变换法（阶跃响应不变法） 一阶保持器z变换法（斜坡响应不变法）,这种方法可以保证连续与离散环节脉冲响应相同，但由于z
9、变换比较麻烦，多个环节串联时无法单独变换以及产生频率混叠和其他特性变化较大，所以应用较少。,注意，这里的零阶保持器是假想的，并没有物理的零阶保持器。这种方法可以保证连续与离散环节阶跃响应相同，但要进行z变换，同样具有z变换法的一系列缺点，所以应用亦较少。,由于和零阶保持器z变换法类似的原因，这种方法应用得较少。,17,5.1 连续域离散化设计 5.2 数字PID控制器设计 5.3 控制系统z平面设计性能指标要求 5.4 z平面根轨迹设计 5.5 w变换及频率域设计,18,5.2 数字PID控制器设计,根据偏差的比例(P)、积分(I)、微分(D)进行控制(简称PID控制)，是控制系统中应用最为广
10、泛的一种控制规律。 优点： 原理简单 通用性强,19,5.2.1 数字PID基本算法,1模拟PID控制算法的离散化,模拟PID控制器的基本规律：,离散化,kT均用k简化表示,20,2PID的增量式算法,仅对应执行机构（如阀门）位置的改变量,算法优点： （1）较为安全。因为一旦计算机出现故障，输出控制指令为零时，执行机构的位置（如阀门的开度）仍可保持前一步的位置，不会给被控对象带来较大的扰动。 （2）计算时不需进行累加，仅需最近几次误差的采样值。,主要问题：执行机构的实际位置也就是控制指令全量的累加需要用计算机外的其他的硬件(如步进电机)实现。,21,图5-21 PID计算机控制系统,22,5.
11、2.2 数字PID控制算法改进,1. 抗积分饱和算法 (1)积分饱和的原因及影响 因长时间出现偏差或偏差较大，计算出的控制量有很大，超出D/A转换器所能表示的数值范围。这时的执行机构已到极限位置，仍不能消除偏差，且由于积分作用，尽管计算PID差分方程式所得的运算结果继续增大或减小，但执行机构已无相应的动作，这就称为积分饱和。当控制量达到饱和后，闭环控制系统相当于被断开。,23,小信号控制下，积分器没有饱和的响应曲线。,控制饱和值不变，但系统给定值加大，使控制作用出现饱和时的仿真曲线,在同样给定值时，控制作用没有饱和限制时的仿真曲线。,24,1.抗积分饱和算法 (2) 积分饱和抑制,积分分离法：
12、 系统加入积分控制的主要作用提高稳态精度，减少或消除误差。 积分分离法的基本控制思想：,某个规定的门限值 ; 当e(k)，=0（即取消积分）； 当误差e(k)=，=1（即引入积分）,无积分分离的响应曲线,有积分分离的响应曲线,图5-23 积分分离法,25,(2) 积分饱和抑制,遇限削弱积分法 ： 基本思想： 当控制量进入饱和区后，只执行削弱积分项的累加，不进行增大积分项的累加。即系统在计算u(k)时，先判断u(k-1)是否超过门限值。若超过某个方向门限值时，积分只累加反方向的e(k)值。,具体算式为：,若,且,不进行积分累加；,进行积分累加。,若,且,不进行积分累加；,若,进行积分累加。,若,
13、26,(2) 积分饱和抑制,饱和停止积分法 ： 基本思想： 当控制作用达到饱和时，停止积分器积分，而控制器输出未饱和时，积分器仍正常积分。 特点： 简单易行，但不如上一种方法容易使系统退出饱和,具体算式为：,若,不进行积分运算；,进行积分运算。,若,27,(2) 积分饱和抑制,反馈抑制积分饱和法： 基本思想： 测量执行机构的输入与输出，并形成误差es，将该信号经过增益1/Tt反馈至积分器输入端，降低积分器输出。 当执行机构未饱和时，es=0； 当执行机构饱和时，附加反馈通道使误差信号es趋于零，使控制器输出处于饱和极限。,图5-24反馈抑制积分饱和法,方案要求： 系统可以测量执行机构的输出。若
14、无法测量执行机构的输出，可以在执行机构之前加入执行机构带饱和限幅的静态数学模型，利用该模型形成误差es，并构成附加反馈通道。,28,2. 防积分整量化误差的方法,当采样周期较小而积分时间常数较大时，积分项的数值很小，有可能使微型机二进制数字最低位无法表示，产生整量化误差，发生积分项丢失的现象。 为了防止积分项由于数的整量化误差所导致的丢失现象，在控制算法及编程方面应采取一定的改进措施： 在积分项运算法时，可以将其结果用双字长单元存贮，若积分项小于单字长时，其积分结果存放在低字节单元中，经过若干次累加后，当其值超过低字节表示时，则在高字节最低位加1，从而消除了有限字长造成的量化截尾误差。,图5-
15、25双字节积分累加,积分项,29,3. 微分算法的改进,(1)不完全微分的PID算式（采用带惯性环节的实际微分器）,引入微分改善了系统的动态特性，但由于微分放大噪声的作用也极易引进高频干扰。,不完全微分PID位置算法,30,不完全微分PID 与基本PID控制作用比较,在e(k)发生阶跃突变时， 完全微分作用仅在控制作用发生的一个周期内起作用； 不完全微分作用则是按指数规律逐渐衰减到零，可以延续几个周期，且第一个周期的微分作用减弱。,图5-26 不完全微分的阶跃响应,31,3. 微分算法的改进,(2)微分先行PID,图5-27 微分先行结构图,适用于给定值频繁升降的场合，可以避免因输入变动而在输
16、出上产生跃变,32,4. 带非灵敏区的PID控制,控制算法 非灵敏区设置值,若,若,，则,，则,图5-28带非灵敏区的PID控制,33,5. 自动与手动无扰转换的PI算法,工业上通过PID控制的被控对象常常有手动与自动两种控制方式，转换时要求实现无扰转换。,图5-29 自动与手动无扰转换,系统处于自动状态时,PI控制,开关处于手动位置时,实现从自动到手动的无扰切换,34,5.2.3 PID调节参数的整定,1扩充临界比例度法 (1)选择一个足够短的采样周期T，通常可选择采样周期为被控对象纯滞后时间的1/10。 (2)用选定的T使系统工作。这时，去掉数字控制器的积分作用和微分作用，只保留比例作用。
17、然后逐渐减小比例度(=1/KP)，直到系统发生持续等幅振荡。记下此时的临界比例度k及系统的临界振荡周期Tk (即振荡波形的两个波峰之间的时间) 。 ,控制度=,(3)选择控制度,(4)根据选定的控制度，查表5-1 ，求得T、KP、TI、TD的值。,(5)按计算所得参数投入在线运行，观察效果，如果性能不满意，可根据经验和对P、I、D各控制项作用的理解，进一步调节参数，直到满意为止。,图5-30等幅振荡曲线,35,表5-1 扩充临界比例度法整定参数,36,2扩充阶跃响应曲线法,用扩充响应曲线法整定T和KP、TI、TD的步骤如下: (1)数字控制器不接入系统，将被控对象的被控制量调到给定值附近，并使
18、其稳定下来，然后测出对象的单位阶跃响应曲线 。 (2)在对象响应曲线的拐点处作一切线，求出纯滞后时间和时间常数Tm以及它们的比值Tm/ 。 (3)选择控制度 (4)查表5-2，即可求得数字控制器的KP、TI、TD及采样周期T。,图5-31 对象的响应曲线,37,表5-2 扩充阶跃响应曲线法PID参数,38,3试凑法确定PID参数,整定步骤： (1)首先只整定比例部分。比例系数KP由小变大，观察相应的系统响应，直到得到反应快，超调小的响应曲线。系统若无静差或静差已小到允许范围内，并且响应效果良好，那么只须用比例调节器即可。 (2)若稳态误差不能满足设计要求，则需加入积分控制。整定时先置积分时间T
19、I为一较大值，并将经第1步整定得到的KP减小些，然后减小TI ，并使系统在保持良好动态响应的情况下，消除稳态误差。这种调整可根据响应曲线的状态，反复改变KP及TI ，以期得到满意的控制过程。 (3)若使用PI调节器消除了稳态误差，但动态过程仍不能满意，则可加入微分环节。在第2步整定的基础上，逐步增大TD，同时相应地改变KP和TI，逐步试凑以获得满意的调节效果。,39,5.1 连续域离散化设计 5.2 数字PID控制器设计 5.3 控制系统z平面设计性能指标要求 5.4 z平面根轨迹设计 5.5 w变换及频率域设计,40,5.3.1 时域性能指标要求,(1) 稳定性要求 (2) 系统稳态特性的要
20、求： 主要以系统在一定指令信号及干扰信号作用下稳态误差的大小来衡量。 影响稳态误差的主要因素是系统的类型及开环放大系数。 (3) 系统动态特性要求： 主要以系统单位阶跃响应的升起时间、峰值时间、超调量和调节时间来表示。 任意高阶系统动态指标是由系统的零极点分布决定的，并且很难计算。但在很多情况下，高阶系统中都有一对主导极点，这时可把高阶系统近似看作二阶系统来研究。,由于多数计算机控制系统的被控对象是连续的，设计时所给定的性能指标要求，基本上与连续系统设计时相同。因此，若在z平面上直接进行离散系统设计，需要考虑如何将连续系统的性能指标转换为z平面的描述。,41,动态指标的求取,动态指标如下： 超
21、调量上升时间峰值时间调节时间（5%误差带）,特征根为,实部和虚部的绝对值,单位阶跃响应,42,5.3.2 频域性能指标要求,开环频率特性低频段的形状： 低频段的形状及幅值大小充分反映了系统的稳态特性，其结论与连续系统类似。2. 开环频率特性中频段的形状： 主要反映系统动态特性要求。通常以开环系统的截止频率、相位稳定裕度、增益稳定裕度h以及在附近幅频特性的斜率要求来描述。3. 开环频率特性高频段的形状： 主要反映系统抑制高频噪声的能力，通常要求开环频率特性高频段幅值衰减要多、要快。,鉴于离散系统频率特性的特点(即，频率特性不是的有理函数) ，因此，频率域设计时，并不直接利用z平面的频率特性，而是
22、将其变换到其他更有利的平面上进行，这时相关的性能要求也应进行变换。,43,5.1 连续域离散化设计 5.2 数字PID控制器设计 5.3 控制系统z平面设计性能指标要求 5.4 z平面根轨迹设计 5.5 w变换及频率域设计,44,5.4.1 z平面根轨迹,系统闭环脉冲传函,图5-33 离散控制系统,D(z)为数字控制器,G(z)为广义被控对象,闭环系统特征方程,连续系统闭环特征方程,结论：离散系统与连续系统的闭环特征方程形式完全一样。连续系统中根轨迹的定义及绘制法则，在z域完全适用.,Z平面根轨迹应相对于单位圆来分析,45,离散系统中根轨迹的绘制法则,z平面根轨迹的特殊性： (1) z平面极点
23、的密集度很高，在用根轨迹分析系统性能时，要求根轨迹的计算精度较高。 (2) z平面的临界放大系数由根轨迹与单位圆的交点求得。 (3) 离散系统脉冲传递函数的零点多于相应的连续系统，只考虑闭环极点位置对系统动态性能的影响是不够的，还需考虑零点对动态响应的影响。,根轨迹方程,开环传递函数,46,5.4.2 z平面根轨迹设计方法,1. 设计步骤 (1)根据给定的时域指标，在z平面画出期望极点的允许范围。 (2)设计数字控制器D(z)。 先求出广义对象脉冲传递函数 然后确定控制器D(z)的结构形式 若要求数字控制器不影响系统的稳态性能，则要求： (3)进行数字仿真研究，检验闭环系统的动态响应。 (4)
24、在计算机上编程实现D(z)算法。,根轨迹法实质上是一种闭环极点的配置技术，即通过反复试凑，设计控制器的结构和参数，使整个闭环系统的主导极点配置在期望的位置上。,常用控制器有一阶相位超前及相位滞后环节：,47,采用一阶控制器的设计流程：,48,2. 设计举例 天线伺服系统离散域根轨迹设计。,系统设计指标（设采样周期T=0.1s ） 超调量 % =15%； 升起时间 tr 0.55s；,(1)设计指标与z平面期望极点位置,根据设计指标计算得：,z域同心圆半径r 0.5,z域射线,图5-36 理想的 z平面极点范围,调节时间 ts1s； 静态速度误差 Kv5。,49,(2) 设计数字控制器D(z),
25、被控对象的脉冲传递函数,Matlab指令,num=20; den=1 10 0; n,d=c2dm(num,den,0.1,zoh),n= 0 0.0736 0.0528 d=1.0000 -1.3679 0.3679,先可取控制器为纯比例环节,绘制系统的根轨迹,结论:根轨迹没有进入期望极点范围,50,改进控制器D(z)的设计,采用零极对消法，选用,得到：,系统的根轨迹,利用Matlab指令K,pole=rlocfind(num, den)，可在选定极点位置后自动计算得：,满足性能指标要求,希望极点：,根轨迹增益,控制器增益,控制器传函,系统静态速度误差系数,51,(3)系统时域仿真,结论：
26、该系统较好地满足了给定的时域动态性能要求。,图5-39 单位阶跃响应曲线,52,5.1 连续域离散化设计 5.2 数字PID控制器设计 5.3 控制系统z平面设计性能指标要求 5.4 z平面根轨迹设计 5.5 w变换及频率域设计,53,5.5.1 W变换,1. 方法,2. W 变换主要特性 (1) 映射关系,图5-40 s、z、w域之间的映射关系,54,(2) s域和w域频率对应关系,s域和z域的频率都用来表示，是系统的真实频率，变换至w域后得到的频率为虚拟频率，以表示。,图5-41 s域和w域的频率变换关系,55,(3) w域传递函数与z传递函数的关系,结论： 如果(n+k)m，则变换后分子
27、填加有新的零点 w传递函数是w的有理分式函数，故G(jv)是虚拟频率v的有理分式函数。 分子分母一般是同阶的。 变换前后稳态增益不变。,56,(4) s域和w域传递函数的关系,当采样周期T减小时，复变量w近似等于复变量s； 传递函数G(s)与G(w )的相似性；,G(s)与G(w )稳态增益维持不变。,若a=5，T=0.1s，则有,因为带ZOH的z变换与双线性变换都能维持稳态增益的不变。,57,(5) w变换与突斯汀变换,结论： 如果(n+k)m，则变换后分子填加有新的零点 w传递函数是w的有理分式函数，故G(jv)是虚拟频率v的有理分式函数。 分子分母一般是同阶的。 变换前后稳态增益不变。,
28、58,5.5.2 w域设计法（步骤）,(1) 给定连续被控对象G(s)，求出z域的广义对象的脉冲传递函数G(z)(2) 将G(z)变换到w平面上 (3) 在w平面设计控制器D(w) 由于平面和s平面的相似性，s平面上的设计技术，如频率法、根轨迹法等均可应用到w平面。 (4) 进行w反变换，求得z域控制器 D(z)(5) 检验z域闭环系统的品质 (6) D(z)控制器在计算机上编程实现。,59,5.5.3 设计举例天线转角 计算机伺服控制系统w域设计,系统设计指标（设采样周期T=0.1s ） 超调量 % =15%； 相稳定裕度m 50o，增益稳定裕度Lh6dB,- 调节时间 ts1s； - 静态
29、速度误差 Kv5,Matlab指令,numz=0.0736 0.0528; denz=1 -1.3674 0.3674; nw dw=d2cm(umz,denz, 0.1,tustin),nw=-.0076 -0.7723 18.4876 dw=1.000 9.2526 0.0000,(1) 求被控对象传函,60,(2) 在w域设计数字控制器,系统开环放大系数设计,数字控制器D(w)设计,w平面的开环传递函数,在w域检查开环稳定裕度要求,Matlab指令,nw=-0.0189 -1.9318 46.2117; dw=1.0000 9.2423 0; figure(1);margin(nw,dw
30、); grid,先取,满足指标要求，但截止频率较低。,61,数字控制器D(w)设计,时域响应特性检查 Simulink仿真结果 超调量大于给定要求，调节时间虽满足要求，但余量不大。,闭环单位阶跃响应,结论：不能满足要求，故需进一步设计动态控制器，其目的是在保证稳定裕度的条件下，进一步增大开环截止频率。,62,数字控制器D(w)设计,为了实现提高截止频率的目的，在正向通道引入超前-迟后环节是合适的。利用连续系统控制理论方法，依据开环频率响应的特点，通过试凑，可以确定超前-迟后环节的分子及分母的时间常数和增益。通过23次修正，最后取,dn=0.2 2; dd=0.02 1; nw=-0.0189
31、-1.9318 46.2117; dw=1.0000 9.2423 0;dgn,dgd=series(nw,dw,dn,dd),Matlab环节串联指令,dgn=-0.0038 -0.4242 5.3787 92.4234; dgd=0.0200 1.1848 9.2423 0;,在w域检查开环稳定裕度要求,满足要求,注意上述两个频率均为虚拟频率。,63,(3) 获取z平面的控制器D(z),进行w反变换,wdd=0.02 1; wdn=0.2 2; zdn,zdd=c2dm(wdn,wdd,0.1,tustin),Matlab指令,zdn=4.2857 -1.4286; zdd=1.0000 0.4286;,控制器的稳态增益,静态设计时要求,故最终z平面的控制器D(z) 应增大稳态增益1.25倍：,64,(4) 闭环系统仿真,图5-46系统单位阶跃响应,图5-47 系统单位斜坡响应,系统无超调， 调节时间小于0.6s。,稳态误差：,均满足要求。,65,第5章 内容结束！,
