连续控制系统的仿真.ppt

第四章连续控制系统的仿真 4 1仿真模型的结构 1 面向闭环系统微分方程 传递函数 的数字仿真 方法 得到系统的闭环传递函数 再用Z变换法获得闭环系统的差分迭代格式 又称仿真模型 或者有微分方程组采用数值积分法 离散相似法离散成迭代格式后进行仿真 局限性 该方法要求用户上机前事先计算出闭环传递函数 在多回路系统中具有很大的局限性 开环系统往往有控制器 控制对象等多个环节组成 写出总的闭环传递函数后研究某个环节的参数变化对系统品质的影响不方便 特别是在多回路系统中研究小闭环中某环节参数变化对系统影响更麻烦 系统中含有非线性环节就写不出闭环传递函数 2 面向系统结构图的数字仿真将系统看作有许多典型环节组成 进行仿真时 将各个环节的参数以及各个环节的连接方式输入计算机让计算机求闭环系统的数值解 4 2典型环节的离散化系统及其差分方程 典型环节 比例环节 积分环节 比例 积分环节 惯性环节 比例 微分环节 滞后环节等 为了减少仿真系统的环节 一般引入二阶与三阶环节 1 三阶环节 展开成部分分式 令 得到状态方程为 其中 输出方程为 其中 1 确定状态变量的初始值 由 1 得 2 对 1 式两边取导数得 3 将式带入 3 式 取t 0得 对 3 式两边取导数得 取t 0得 联立三式得 状态方程标准化 将A B矩阵代入下式 其中 离散输出方程为 求初始值的方法2 5 状态变量的初始值 可以根据环节在初始状态的稳态输出和稳态输入直接由 5 所给定的参数确定 根据 得到 令S 0即t 根据 5 得状态变量得稳态初始值 2 滞后环节 传递函数 其中 滞后时间得整数部分 以采样周期T为基本单位 滞后时间得小数部分 取Z变换得 差分方程为 下面分两种情况讨论上式得解法及其程序实现 1 a为T的整数倍 仿真方法 在内存中设置一个区域 存放输入U过去的采样值 该区域占用 个单元 从左向右 其编号为0 1 2 每一次计算 首先将该时刻的输入存入单元 然后从第一单元取出 作为该时刻的输出 最后各单元向左平移一次 2 a不为T的整数倍 近似处理法 从1号单元和0号单元取出2个数值并加以差值 得到采用线性差值法 即假定1号和2号两个单元之间的数据是按线性规律变化 近似公式为 如 即表示 4 3控制系统结构图的数值表示方法 仿真方法 将系统看成是由若干环节构成 仿真时先将各典型环节的参数及各环节的连接关系输入计算机 由计算机程序构成系数矩阵 再用四阶龙格 库塔方法对构成的状态方程进行仿真计算 又称为连续系统结构图法仿真 优点 由于输入参数是各环节的参数 可以直接求出各环节的输出 方便地研究环节参数变化对系统地和环节本身的影响 避免了求整体传递函数的问题 1 仿真模型 1 典型环节的确定 一个控制系统的结构图通常由一些不同的环节组成 一般常见的环节如下 积分环节 比例 积分环节 惯性环节 比例 惯性环节 振荡环节 为了编程的方便 各典型环节可用一个典型环节来表示 选定的典型环节应具有如下特点 1 具有典型性 即能用它来描述较多的环节类型 保证所编制的环节有较大的通用性 2 环节的结构简单 用它编制的程序简单 易于计算 选定比例 惯性环节作为典型环节 比例 惯性环节的传递函数 通过改变系数可以转化为其它常用环节 振荡环节的特征方程因式分解后分为三种情况 1 实根时 用两个串联的惯性环节替代 2 共轭复根时 用两个串联的惯性环节加负反馈替代 3 一般复根时 采用积分环节和惯性环节串联加负反馈替代 以 为例 存在特征根 1 5 和 1 5 等效环节一 等效环节二 2 传递函数的矩阵表达式 用n个典型环节描述n阶控制系统的传递函数表达式为 将其转换成一阶微分方程组形式为 对每个环节列写出对应的传递函数为 改写成一阶微分方程组形式 矩阵方程形式 其中 C矩阵与D矩阵类似 A B C D都是n n维的系数矩阵 u和y分别为输入与输出的列向量 3 连接矩阵 连接矩阵 就是用矩阵形式表示各个环节之间的关系 该系统由5个环节构成 环节之间的连接关系 写成矩阵方程形式 向量形式 其中 W为系统的连接矩阵 反映了系统中各环节之间的连接关系 每个元素表示系统中第个环节的输出对第个环节的输入之间的连接关系 1表示直接连接 负号表示反馈 0表示没有直接连接关系 为外部连接矩阵 元素表示外加输入信号与第个环节输入端的连接关系 4 控制系统的状态方程 将 拉氏变换 代入 得到 展开合并得 拉氏反变换得到 令 则 当的逆矩阵存在时 当输入函数是阶跃函数时 在t 0瞬间出现一个无穷大的迈脉冲 方程右端第三项将趋于无穷大 系统求解不方便 为此 需要增加限制条件 当输入函数为脉冲函数时使式中的第三项为零 即 令 反映输入作用函数与被作用环节之间的连接关系 实际上只要限制输入作用函数作用的环节的 即该环节的传递函数不应具有导数项 则 令 标准型的状态方程 4 4非线性系统的仿真 非线性环节包括 饱和非线性 分段线性的非线性 失灵区非线性 间隙非线性 继电非线性 具有死区的继电非线性 具有滞环的继电非线性 具有死区和滞环的继电非线性 1 饱和 分段线性及失灵区非线性特性 表示第一段的斜率 表示第二段的斜率 为拐点的横坐标 分段线性的非线性特性 数学描述 当时 得到饱和非线性 当时 得到失灵区非线性 2 间隙非线性特性 间隙非线性特性 设为前一次输入 为前一次输出 间隙非线性数学描述 当 则 系统工作在右边的特性上 当 则 系统工作在左边的特性上 其他情况 系统工作在中间过渡的特性上 注意点 存储本次的输入与输出 3 继电非线性及具有死区的继电非线性 略 4 具有滞环的继电非线性 具有滞环的继电非线性特性 5 具有死区和滞环的继电非线性特性 具有死区和滞环的继电非线性特性 4 5仿真误差与数值稳定性 仿真误差和数值稳定性均与计算方法 机器精度 步长的选择有关 当计算方法和机器选定以后 则仅仅与步长的选择有关 所以在仿真中 步长是一个重要的参数 仿真误差有截断误差和舍入误差两种 计算步长愈小 截断误差也愈小 舍入误差愈大 在实际解决问题时 通常是根据系统响应的速度由经验确定 采样频率一般可取 其中为被仿真