零极点匹配法.ppt

1 5 1连续域 离散化设计5 2数字PID控制器设计5 3控制系统z平面设计性能指标要求5 4z平面根轨迹设计5 5w 变换及频率域设计 2 5 1 1设计原理和步骤 实质是将数字控制器部分看成是一个整体 其输入和输出都是模拟量 因而可等效为连续传递函数De s A D输出与输入关系 系统低通且采样频率较高 D A的频率特性 等效连续传递函数 计算机实现算法D z 的计算表示 设计时常近似为 图5 1计算机控制系统典型组成 3 连续域 离散化设计的步骤如下 第1步 根据系统的性能 选择采样频率 并设计抗混叠前置滤波器 第2步 考虑ZOH的相位滞后 根据系统的性能指标和连续域设计方法 设计数字控制算法等效传递函数Ddc s 第3步 选择合适的离散化方法 将Ddc s 离散化 获得脉冲传递函数D z 使两者性能尽量等效 第4步 检验计算机控制系统闭环性能 若满足指标要求 进行下一步 否则 重新进行设计 改进设计的途径有 选择更合适的离散化方法 提高采样频率 修正连续域设计 如增加稳定裕度指标等 第5步 将D z 变为数字算法 在计算机上编程实现 4 5 1 2各种离散化方法 最常用的表征控制器特性的主要指标 零极点个数 系统的频带 稳态增益 相位及增益裕度 阶跃响应或脉冲响应形状 频率响应特性 等效离散 D z 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 一定稳定 变换前后 稳态增益不变 离散后控制器的时间响应与频率响应 与连续控制器相比有相当大的畸变 图5 4向后差分法的映射关系 3 应用由于这种变换的映射关系畸变严重 变换精度较低 所以 工程应用受到限制 用得较少 7 2 一阶向前差分法 1 离散化公式 实质 将连续域中的微分用一阶向前差分替换 做z变换 得 s与z之间的变换关系 比较 图5 7向前差分矩形积分法 8 2 一阶向前差分法 2 主要特性 s平面与z平面映射关系 若D s 稳定 则D z 不一定稳定 改进方法是适当减少采样周期T 向前差分法的映射关系图 3 应用由于这种变换的映射关系畸变严重 不能保证D z 一定稳定 或者如要保证稳定 要求采样周期较小 所以应用较少 平移放大关系 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 主要特性 频率畸变 双线性变换的一对一映射 保证了离散频率特性不产生频率混叠现象 但产生了频率畸变 图5 11双线性变换的频率关系 图5 12双线性变换的频率关系 当采样频率较高 足够小 12 3 双线性变换法 2 主要特性 变换前后 稳态增益不变 双线性变换后D z 的阶次不变 且分子 分母具有相同的阶次 并有下式成立 3 应用 这种方法使用方便 且有一定的精度和前述一些好的特性 工程上应用较为普遍 这种方法的主要缺点是高频特性失真严重 主要用于低通环节的离散化 不宜用于高通环节的离散化 13 4 修正双线性变换 1 离散化方法 解决 双线性变换产生频率轴非线性畸变 问题的方法 1是设计者选定的特征角频率 依据连续域与双线性变换后频率的非线性关系 首先修正原连续域传递函数 然后再进行双线性变换的结果 2 主要特性该方法本质上仍为双线性变换法 因此具有双线性变换法的各种特性 但由于采用了频率预修正 故可以保证在处连续频率特性与离散后频率特性相等 即满足 3 应用由于该方法的上述特性 所以主要用于原连续控制器在某些特征频率处要求离散后频率特性保持不变的场合 14 5 零极点匹配法 1 离散化方法 特点 零 极点分别按匹配若分子阶次m小于分母阶次n 离散变换时 在D z 分子上加 z 1 n m因子确定D z 的增益k1的方法 按右式来匹配若D s 分子有s因子 可依高频段增益相等原则确定增益 即也可选择某关键频率处的幅频相等 即 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变换比较麻烦 多个环节串联时无法单独变换以及产生频率混叠和其他特性变化较大 所以应用较少 注意 这里的零阶保持器是假想的 并没有物理的零阶保持器 这种方法可以保证连续与离散环节阶跃响应相同 但要进行z变换 同样具有z变换法的一系列缺点 所以应用亦较少 由于和零阶保持器z变换法类似的原因 这种方法应用得较少 17 5 1连续域 离散化设计5 2数字PID控制器设计5 3控制系统z平面设计性能指标要求5 4z平面根轨迹设计5 5w 变换及频率域设计 18 5 2数字PID控制器设计 根据偏差的比例 P 积分 I 微分 D 进行控制 简称PID控制 是控制系统中应用最为广泛的一种控制规律 优点 原理简单通用性强 19 5 2 1数字PID基本算法 1 模拟PID控制算法的离散化 模拟PID控制器的基本规律 离散化 kT均用k简化表示 20 2 PID的增量式算法 仅对应执行机构 如阀门 位置的改变量 算法优点 1 较为安全 因为一旦计算机出现故障 输出控制指令为零时 执行机构的位置 如阀门的开度 仍可保持前一步的位置 不会给被控对象带来较大的扰动 2 计算时不需进行累加 仅需最近几次误差的采样值 主要问题 执行机构的实际位置也就是控制指令全量的累加需要用计算机外的其他的硬件 如步进电机 实现 21 图5 21PID计算机控制系统 22 5 2 2数字PID控制算法改进 1 抗积分饱和算法 1 积分饱和的原因及影响因长时间出现偏差或偏差较大 计算出的控制量有很大 超出D A转换器所能表示的数值范围 这时的执行机构已到极限位置 仍不能消除偏差 且由于积分作用 尽管计算PID差分方程式所得的运算结果继续增大或减小 但执行机构已无相应的动作 这就称为积分饱和 当控制量达到饱和后 闭环控制系统相当于被断开 23 小信号控制下 积分器没有饱和的响应曲线 控制饱和值不变 但系统给定值加大 使控制作用出现饱和时的仿真曲线 在同样给定值时 控制作用没有饱和限制时的仿真曲线 24 1 抗积分饱和算法 2 积分饱和抑制 积分分离法 系统加入积分控制的主要作用 提高稳态精度 减少或消除误差 积分分离法的基本控制思想 某个规定的门限值 当e k 0 即取消积分 当误差e k 1 即引入积分 无积分分离的响应曲线 有积分分离的响应曲线 图5 23积分分离法 25 2 积分饱和抑制 遇限削弱积分法 基本思想 当控制量进入饱和区后 只执行削弱积分项的累加 不进行增大积分项的累加 即系统在计算u k 时 先判断u k 1 是否超过门限值 若超过某个方向门限值时 积分只累加反方向的e k 值 具体算式为 若 且 不进行积分累加 进行积分累加 若 且 不进行积分累加 若 进行积分累加 若 26 2 积分饱和抑制 饱和停止积分法 基本思想 当控制作用达到饱和时 停止积分器积分 而控制器输出未饱和时 积分器仍正常积分 特点 简单易行 但不如上一种方法容易使系统退出饱和 具体算式为 若 不进行积分运算 进行积分运算 若 27 2 积分饱和抑制 反馈抑制积分饱和法 基本思想 测量执行机构的输入与输出 并形成误差es 将该信号经过增益1 Tt反馈至积分器输入端 降低积分器输出 当执行机构未饱和时 es 0 当执行机构饱和时 附加反馈通道使误差信号es趋于零 使控制器输出处于饱和极限 图5 24反馈抑制积分饱和法 方案要求 系统可以测量执行机构的输出 若无法测量执行机构的输出 可以在执行机构之前加入执行机构带饱和限幅的静态数学模型 利用该模型形成误差es 并构成附加反馈通道 28 2 防积分整量化误差的方法 当采样周期较小而积分时间常数较大时 积分项的数值很小 有可能使微型机二进制数字最低位无法表示 产生整量化误差 发生积分项丢失的现象 为了防止积分项由于数的整量化误差所导致的丢失现象 在控制算法及编程方面应采取一定的改进措施 在积分项运算法时 可以将其结果用双字长单元存贮 若积分项小于单字长时 其积分结果存放在低字节单元中 经过若干次累加后 当其值超过低字节表示时 则在高字节最低位加1 从而消除了有限字长造成的量化截尾误差 图5 25双字节积分累加 积分项 29 3 微分算法的改进 1 不完全微分的PID算式 采用带惯性环节的实际微分器 引入微分改善了系统的动态特性 但由于微分放大噪声的作用也极易引进高频干扰 不完全微分PID位置算法 30 不完全微分PID与基本PID控制作用比较 在e k 发生阶跃突变时 完全微分作用仅在控制作用发生的一个周期内起作用 不完全微分作用则是按指数规律逐渐衰减到零 可以延续几个周期 且第一个周期的微分作用减弱 图5 26不完全微分的阶跃响应 31 3 微分算法的改进 2 微分先行PID 图5 27微分先行结构图 适用于给定值频繁升降的场合 可以避免因输入变动而在输出上产生跃变 32 4 带非灵敏区的PID控制 控制算法非灵敏区设置值 若 若 则 则 图5 28带非灵敏区的PID控制 33 5 自动与手动无扰转换的PI算法 工业上通过PID控制的被控对象常常有手动与自动两种控制方式 转换时要求实现无扰转换 图5 29自动与手动无扰转换 系统处于自动状态时 PI控制 开关处于手动位置时 实现从自动到手动的无扰切换 34 5 2 3PID调节参数的整定 1 扩充临界比例度法 1 选择一个足够短的采样周期T 通常可选择采样周期为被控对象纯滞后时间的1 10 2 用选定的T使系统工作 这时 去掉数字控制器的积分作用和微分作用 只保留比例作用 然后逐渐减小比例度 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 数字控制器不接入系统 将被控对象的被控制量调到给定值附近 并使其稳定下来 然后测出对象的单位阶跃响应曲线 2 在对象响应曲线的拐点处作一切线 求出纯滞后时间 和时间常数Tm以及它们的比值Tm 3 选择控制度 4 查表5 2 即可求得数字控制器的KP TI TD及采样周期T 图5 31对象的响应曲线 37 表5 2扩充阶跃响应曲线法PID参数 38 3 试凑法确定PID参数 整定步骤 1 首先只整定比例部分 比例系数KP由小变大 观察相应的系统响应 直到得到反应快 超调小的响应曲线 系统若无静差或静差已小到允许范围内 并且响应效果良好 那么只须用比例调节器即可 2 若稳态误差不能满足设计要求 则需加入积分控制 整定时先置积分时间TI为一较大值 并将经第1步整定得到的KP减小些 然后减小TI 并使系统在保持良好动态响应的情况下 消除稳态误差 这种调整可根据响应曲线的状态 反复改变KP及TI 以期得到满意的控制过程 3 若使用PI调节器消除了稳态误差 但动态过程仍不能满意 则可加入微分环节 在第2步整定的基础上 逐步增大TD 同时相应地改变KP和TI 逐步试凑以获得满意的调节效果 39 5 1连续域 离散化设计5 2数字PID控制器设计5 3控制系统z平面设计性能指标要求5 4z平面根轨迹设计5 5w 变换及频率域设计 40 5 3 1时域性能指标要求 1 稳定性要求 2 系统稳态特性的要求 主要以系统在一定指令信号及干扰信号作用下稳态误差的大小来衡量 影响稳态误差的主要因素是系统的类型及开环放大系数 3 系统动态特性要求 主要以系统单位阶跃响应的升起时间 峰值时间 超调量和调节时间来表示 任意高阶系统动态指标是由系统的零极点分布决定的 并且很难计算 但在很多情况下 高阶系统中都有一对主导极点 这时可把高阶系统近似看作二阶系统来研究 由于多数计算机控制系统的被控对象是连续的 设计时所给定的性能指标要求 基本上与连续系统设计时相同 因此 若在z平面上直接进行离散系统设计 需要考虑如何将连续系统的性能指标转换为z平面的描述 41 动态指标的求取 动态指标如下 超调量上升时间峰值时间调节时间 5 误差带 特征根为 实部和虚部的绝对值 单位阶跃响应 42 5 3 2频域性能指标要求 开环频率特性低频段的形状 低频段的形状及幅值大小充分反映了系统的稳态特性 其结论与连续系统类似 2 开环频率特性中频段的形状 主要反映系统动态特性要求 通常以开环系统的截止频率 相位稳定裕度 增益稳定裕度h以及在附近幅频特性的斜率要求来描述 3 开环频率特性高频段的形状 主要反映系统抑制高频噪声的能力 通常要求开环频率特性高频段幅值衰减要多 要快 鉴于离散系统频率特性的特点 即 频率特性不是的有理函数 因此 频率域设计时 并不直接利用z平面的频率特性 而是将其变换到其他更有利的平面上进行 这时相关的性能要求也应进行变换 43 5 1连续域 离散化设计5 2数字PID控制器设计5 3控制系统z平面设计性能指标要求5 4z平面根轨迹设计5 5w 变换及频率域设计 44 5 4 1z平面根轨迹 系统闭环脉冲传函 图5 33离散控制系统 D z 为数字控制器 G z 为广义被控对象 闭环系统特征方程 连续系统闭环特征方程 结论 离散系统与连续系统的闭环特征方程形式完全一样 连续系统中根轨迹的定义及绘制法则 在z域完全适用 Z平面根轨迹应相对于单位圆来分析 45 离散系统中根轨迹的绘制法则 z平面根轨迹的特殊性 1 z平面极点的密集度很高 在用根轨迹分析系统性能时 要求根轨迹的计算精度较高 2 z平面的临界放大系数由根轨迹与单位圆的交点求得 3 离散系统脉冲传递函数的零点多于相应的连续系统 只考虑闭环极点位置对系统动态性能的影响是不够的 还需考虑零点对动态响应的影响 根轨迹方程 开环传递函数 46 5 4 2z平面根轨迹设计方法 1 设计步骤 1 根据给定的时域指标 在z平面画出期望极点的允许范围 2 设计数字控制器D z 先求出广义对象脉冲传递函数然后确定控制器D z 的结构形式若要求数字控制器不影响系统的稳态性能 则要求 3 进行数字仿真研究 检验闭环系统的动态响应 4 在计算机上编程实现D z 算法 根轨迹法实质上是一种闭环极点的配置技术 即通过反复试凑 设计控制器的结构和参数 使整个闭环系统的主导极点配置在期望的位置上 常用控制器有一阶相位超前及相位滞后环节 47 采用一阶控制器的设计流程 48 2 设计举例 天线伺服系统离散域根轨迹设计 系统设计指标 设采样周期T 0 1s 超调量 15 升起时间tr 0 55s 1 设计指标与z平面期望极点位置 根据设计指标计算得 z域同心圆半径r 0 5 z域射线 图5 36理想的z平面极点范围 调节时间ts 1s 静态速度误差Kv 5 49 2 设计数字控制器D z 被控对象的脉冲传递函数 Matlab指令 num 20 den 1100 n d c2dm num den 0 1 zoh n 00 07360 0528 d 1 0000 1 36790 3679 先可取控制器为纯比例环节 绘制系统的根轨迹 结论 根轨迹没有进入期望极点范围 50 改进控制器D z 的设计 采用零极对消法 选用 得到 系统的根轨迹 利用Matlab指令 K pole rlocfind num den 可在选定极点位置后自动计算得 满足性能指标要求 希望极点 根轨迹增益 控制器增益 控制器传函 系统静态速度误差系数 51 3 系统时域仿真 结论 该系统较好地满足了给定的时域动态性能要求 图5 39单位阶跃响应曲线 52 5 1连续域 离散化设计5 2数字PID控制器设计5 3控制系统z平面设计性能指标要求5 4z平面根轨迹设计5 5w 变换及频率域设计 53 5 5 1W 变换 1 方法 2 W 变换主要特性 1 映射关系 图5 40s z w 域之间的映射关系 54 2 s域和w 域频率对应关系 s域和z域的频率都用 来表示 是系统的真实频率 变换至w 域后得到的频率为虚拟频率 以 表示 图5 41s域和w 域的频率变换关系 55 3 w 域传递函数与z传递函数的关系 结论 如果 n k m 则变换后分子填加有新的零点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的有理分式函数 分子分母一般是同阶的 变换前后稳态增益不变 58 5 5 2w 域设计法 步骤 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 增益稳定裕度Lh 6dB 调节时间ts 1s 静态速度误差Kv 5 Matlab指令 numz 0 07360 0528 denz 1 1 36740 3674 nwdw d2cm umz denz