辅助分析与设计.ppt

第章自动控制系统计算机辅助分析 自动控制系统的计算机辅助分析是以理论分析为依据 在已经建立的自动控制系统数学模型的基础上 通过编程实现对系统稳定性 动态和稳态性能进行分析的一门应用技术 matlab以其方便灵活的编程 丰富的工具箱 以及强大的计算和绘图功能成为目前世界上最为流行的自动控制系统辅助分析软件 1自动控制系统的稳定性分析 matlab提供了求取特征方程根的函数roots 其调用格式为式中 p为特征多项式的系数向量 返回值v是特征根构成的列向量 matlab还提供一个可以直接求取矩阵特征值的函数eig 其调用格式为 其中d为矩阵a的特征值向量 调用该函数时 也可以给出两个返回值 其中v是由与特征值相对应的特征向量构成的变换矩阵 例2 某线性控制系统的状态方程为试求出系统特征多项式以及特征值 并且作线性变换要求变换后系统矩阵为对角阵 解题过程见教材第98页 1 2控制系统的能控性和能观性分析 在 现代控制理论 课程中 我们已经知道 线性定常系统 如果它的能控性矩阵为满秩 则该系统为状态完全能控 或称该系统是能控的 对于线性定常系统 如果它的能观性矩阵为满秩 则该系统为状态完全能观 或称该系统是能观的 matlab中有用于计算能控性矩阵的函数ctrb 其格式为 计算能观性矩阵的函数obsv matlab中还有计算矩阵秩的函数rank 这些函数可以帮助我们分析控制系统的能控性和能观性 例3 分析下面的线性系统是否能控 是否能观测 1 3利用传递函数的极点判别系统稳定性 控制系统的传递函数 或脉冲传递函数 以有理真分式形式给出时 matlab提供的函数tf2zp 和pzmap 可以用来求取系统的极点和零点 进而实现对系统稳定性的判断 例4 已知某控制系统如下图所示 试求出闭环系统的极点 并且判断闭环系统的稳定性 解输入命令 计算机显示 表示该系统的闭环传递函数为 再判断闭环极点 输入 计算机输出 显然 3个闭环极点全部位于左半复平面 因此 闭环系统稳定 1 4利用李亚普诺夫第二法判别系统稳定性 对于非线性系统 没有求lyapunov函数的一般方法 matlab也没有这个功能 只能判断齐次线性定常系统的稳定性 例5 齐次线性定常系统方程如下 试判断系统的稳定性 解编写matlab程序如下 计算机执行以后 输出 由于矩阵p的各阶主子式的行列式都为正 p为正定 因此本系统为大范围一致渐近稳定 2控制系统时域分析 2 1时域分析的一般方法 对于稳定的控制系统来说 其时域特性可以由暂态响应和稳态响应的性能指标来表示 最为常见的是用控制系统单位阶跃响应的特征来定义系统的动态时域性能指标 主要有 上升时间 峰值时间 超调量和调节时间等 需要指出 系统动态性能指标定义的前提是系统为稳定的 控制系统的稳态性能指标通常用系统的稳态误差来表示 2 2常用时域分析函数 在matlab中 常用的时域分析函数主要有以下几种 step 绘制连续系统的单位阶跃响应曲线 dstep 绘制离散系统的单位阶跃响应曲线 impulse 绘制连续系统的单位尖脉冲响应曲线 dimpulse 绘制系统的单位尖脉冲响应曲线 lsim 绘制连续系统的任意输入响应曲线 dlsim 绘制离散系统的任意输入响应曲线 例6 已知控制系统闭环传递函数如下 试用matlab绘制其单位阶跃响应曲线 解输入命令 计算机就绘制出该系统的单位阶跃响应曲线如图所示 再输入命令 计算机就绘制出该系统的单位脉冲响应曲线如下图所示 例7 已知二阶闭环控制系统如图所示 试在4个子图中绘出当无阻尼自然振荡频率 阻尼比 分别为0 2 0 5 1 0和2 5等不同值时 系统的单位阶跃响应曲线 2 3时域分析应用实例 摆杆长度为l 质量为m的单级倒立摆 摆杆的质心在杆的中心处 小车的质量为m 在水平方向施加控制力u 相对参考系产生位移为y 为了简化问题并且保其实质不变 忽略执行电机的惯性以及摆轴 轮轴 轮与接触面之间的摩擦力及风力 摆杆质心坐标为 在y轴方向上应用牛顿第二定律得以下方程 1 而 2 代入 1 式 化简为 在转动方向上 其转矩平衡方程为 3 4 或 不失一般性 不妨选取倒立摆的参数如下 选取状态变量 判断开环系统的稳定性 输入命令 计算结果 可见 有一个特征值位于右半复平面 开环系统不稳定 判断系统的能控性 计算机返回r 4 系统能控 则可以通过状态反馈配置系统极点 例如 我们希望通过状态反馈 将系统极点配置为 5 6 则使用命令place 可以求出状态反馈矩阵k 输入命令 计算机返回 建立该状态反馈控制系统的仿真模型 输出曲线 3控制系统频率域分析 稳定的线性定常系统 在正弦输入信号作用下 其输出的稳态分量是与输入同频率的正弦函数 进入稳态以后 输出正弦信号的振幅和输入正弦信号振幅之比称为幅频特性 而输出正弦信号的相位和输入正弦信号的相位之差称为相频特性 3 1频域分析的一般方法 在频率分析法中 判别闭环系统稳定性的最基本定理是nyquist判据 对于开环稳定的系统来说 开环传递函数的极点全部位于左半复平面以内 则闭环系统为稳定的充分必要条件为 开环频率特性的奈氏曲线不包围 1 j0 点 在半对数坐标纸上 分别绘制对数幅频特性和相频特性 就称为伯德图 在matlab中 为我们提供了nyquist bode 和margin 等命令 使我们非常方便地使用频率特性来分析系统 3 2频域分析应用实例 在matlab编程语言中 绘制奈氏曲线的命令是nyquist 其基本格式为 nyquist sys 或nyquist sys w 或nyquist sys wmin wmax 类似地 绘制伯德图的命令是bode 其基本格式为 bode sys 或bode sys w 或bode sys wmin wmax 例6 8 已知单位负反馈线性定常系统的开环传递函数为试绘制其奈氏曲线 并且判断闭环系统是否稳定 解首先判断开环系统是否稳定 输入命令 计算机返回 可见 开环系统稳定 再输入命令 计算机绘制出奈氏图 由于奈氏曲线不包围 1 j0 点 因此 闭环系统为稳定 输入命令 计算机绘制出bode图 并且计算出幅值裕度和相角裕度 显然 系统稳定 4根轨迹分析方法 控制系统的根轨迹分析方法就是利用系统的某个参数 通常是开环增益 从0变化到无穷大时 闭环系统特征根所留下的轨迹 即根轨迹 来分析系统性能以及参数变化对系统性能的影响 4 1幅值条件和相角条件 幅值条件方程 相角条件方程 复平面上满足相角条件的所有s点的集合就是系统的根轨迹 当被确定为某一数值时 根据幅值条件就可以确定闭环极点的位置 4 2绘制根轨迹的常用函数及其应用实例 在matlab编程语言中 有绘制根轨迹的命令rlocus 其基本格式为rlocus sys 和rlocus sys t 或者rlocus num den t 和rlocus num den t 执行该命令后 根轨迹图自动生成 如果给定参数t 则绘制当t从0变化到无穷大时的广义根轨迹 例6 12 已知系统开环传递函数如下 绘制系统根轨迹 并求出闭环系统临界稳定时的根轨迹增益值 解在matlab命令窗口键入命令 计算机绘制出系统根轨迹如图 再输入命令在图形窗口出现十字光标 因为闭环系统为临界稳定 所以选择闭环极点在虚轴上 点击鼠标左键 就确定了闭环极点 同时 在command窗口 计算机给出了相应的数值 5 基于计算机仿真的非线性定常控制系统新型稳定性判据 5 1问题的提出 对于一个控制系统来说 其最重要的属性就是稳定性 一个不稳定的系统是无法工作的 长期以来 对于非线性控制系统的稳定性分析 通常采用lyapunov第二法 但是 对于有些非线性系统构造合适的广义能量函数非常困难 到目前为止 仍然没有一个构造lyapunov函数的一般性的方法 近数十年来 计算机技术取得了突飞猛进的发展 而高配置的个人计算机以及像matlab这样优秀的计算与仿真软件越来越普及 使用计算机仿真来分析非线性控制系统的稳定性成为可能的解决方法之一 可见 计算机不仅在技术层面 而且在理论层面 都深刻地影响着控制理论与控制工程学科的发展 本书作者提出了一种新型的基于计算机仿真的非线性定常控制系统稳定性分析方法 5 2新型稳定性判据 首先 作为一个例子 我们考察以下非线性方程 选取正定的lyapunov函数如下 则其一阶导数为负定 并且当时 根据lyapunov稳定性理论 该系统状态空间原点为大范围一致渐近稳定的平衡状态 随时间变化曲线 随时间变化曲线 从以上观察得到启示 是否可以根据各状态分量的平方和函数是否收敛到零来判别非线性系统的稳定性 而不必构造出lyapunov函数 对于同一个非线性定常系统 如果选择各状态分量的平方和函数 非线性定常控制系统在平衡点的某个邻域内有以下4种运动形态 1 渐近稳定 2 发散 3 以极限环形式作自持振荡运动 4 在非平衡点的某些状态上驻留 定理 xiaojiangzhang scriterionforstability 对于非线性定常系统 其中为状态向量 设 1 该系统的平衡点为状态空间原点 如果平衡点不在原点 则通过变量代换坐标平移可以将平衡点平移至状态空间原点而不影响系统稳定性 2 当 其中为状态空间原点的半径为的邻域 时 为有界 则该系统为局部一致渐近稳定的充分必要条件是 通过仿真或数值计算 在充分长时间之后 各状态分量的平方和函数趋向于零 即 如果 则系统一致渐近稳定 如果不趋向零 则系统不是一致渐近稳定 如果在有界的范围内波动 则该非线性系统为自持震荡 需要指出 以上定理不适用于时变系统 因为时变系统的参数随时间变化而发生改变 有可能导致系统经过一段稳定状态或者驻留状态后又变成为发散 说明1 在lyapunov第二法中 条件 v为正定且为负定 可以确保得出结论 当时 则系统大范围渐近稳定 然而 本节提出的定理却不能得出大范围稳定性的这一结论 说明2 采用传统的lyapunov稳定性理论来判断系统稳定性时 其可信度取决于系统数学模型与真实系统的接近程度 而本节提出的基于计算机仿真的稳定性判据 其可信度取决于所建立的系统仿真模型与真实系统的接近程度 在matlab simulink环境下 系统仿真模型就是依据系统数学模型建立的 如果恰当地选择算法和步长 则两者并无显著差别 因此 两种判据的可信度是相当的 说明3 传统的lyapunov稳定性理论的实质是 对于某个控制系统 不去求它的解析解 对于许多非线性系统来说 无法得出解析解 转而去确定广义能量函数及其一阶导数的正定性 而xiaojiangzhang s稳定性判据的实质是 使用数字计算机解出控制系统的数值解 通过判断状态变量平方和函数是否收敛 就可以判断系统的稳定性 是一种基于数值解的稳定性判据 在matlab simulink环境下建立单级倒立摆模糊控制系统的仿真计算模型如图所示 随时间变化而振荡衰减收敛到零 根据本节的定理可以得知 对于实际的倒立摆模糊控制系统 参数设置和控制方法和仿真系统一致 该系统在状态空间原点为局域一致渐近稳定的平衡点 对该倒立摆模糊控制系统的实验也验证了该系统是稳定的 第章自动控制系统计算机辅助设计 1概述 使用matlab不仅可以解决控制系统的分析问题 还可以解决系统的设计问题 在掌握matlab以后 设计过程大大简化 设计效率大大提高 将人们从以往繁琐的计算绘图工作中彻底解放出来 自动控制系统设计变得方便 快捷 单输入单输出 siso 系统校正分为串联校正 并联校正和反馈校正等几种形式 在此我们仅以串联校正为例说明 2超前校正 滞后校正以及滞后 超前校正的bode图设计 在频率特性法中 由开环系统的bode图来分析闭环控制系统稳定性时 通常采用相角裕量和幅值裕量来描述闭环系统的相对稳定性 2 1超前校正器的bode图设计 设计超前校正器的步骤见教材第123页 在此举例说明 例7 1 某一个控制系统如图所示 设计超前校正器 使系统满足 1 在单位斜坡信号作用下 系统的稳态误差 0 001 2 校正后系统相角裕量的范围为 解 1 根据稳态误差要求 选取 则开环传递函数为 2 此时使用matlab中命令margin 来计算校正前系统的幅值裕量 相角裕量和穿越频率 输入命令 计算机绘制出该系统的bode图 并且计算出相应的幅值裕量和相角裕量 此时 幅值裕量 相角裕量 穿越频率 闭环系统将不稳定 需要进行超前校正 3 选取 建立一个m文件 不妨命名为fowrdgn m 如下 在command窗口键入该文件名fowrdgn并且回车 计算机就得出 即表示 于是 校正后系统的开环传递函数为 输入以下命令 计算机绘出校正以后系统的bode图 计算机同时计算出幅值增益裕度 相角裕度和穿越频率如下 显然 可以满足系统的性能指标要求 2 2滞后校正器的bode图设计 设计滞后校正器的步骤见教材第126页 在此举例说明 解 1 根据稳态误差要求 选取 则传函为 2 此时使用matlab中命令margin 来计算校正前系统的幅值裕量 相角裕量和穿越频率 输入命令 3 求滞后校正器的传递函数 根据设计要求 选取校正后的相角裕量 建立一个m文件 不妨命名为lagdgn m 如下 2 3滞后 超前校正器的bode图设计 当被校正的系统不稳定 并且要求校正后系统的响应速度 相角裕量和稳态精度较高时 以采用串联滞后 超前校正为宜 该方法是利用滞后 超前校正器的超前部分来增大系统的相角裕量 同时又利用滞后部分来改善系统的稳态性能 滞后 超前校正器的传递函数为 设计滞后 超前校正器的步骤见教材第128页 在此举例说明 解 1 根据稳态误差要求 因此 则 2 此时使用matlab中命令margin 来计算校正前系统的幅值裕量 相角裕量和穿越频率 输入命令 可知此时幅值裕量 相角裕量 穿越频率 系统不稳定 需要进行滞后 超前校正 3 求滞后校正器的传递函数 根据设计要求 选取校正后的相角裕量 建立一个m文件 不妨命名为laglead m 如下 得出滞后 超前校正器为 3pid控制器设计 pid 比例 积分 微分 控制器是目前在实际工程中应用最为广泛的一种控制策略 pid算法简单实用 不要求受控对象的精确数学模型 3 1pid控制器的传递函数 1 连续pid控制器的传递函数 连续系统pid控制器的表达式为 连续pid控制器的传递函数 为了避免纯微分运算 通常采用近似的pid控制器 其传递函数为 3 2pid控制器各参数对控制性能的影响 pid控制器的 和三个参数的大小决定了pid控制器的比例 积分和微分控制作用的强弱 例4 某直流电机速度控制系统如图所示 采用pid控制方案 使用期望特性法来确定 和这三个参数 建立该系统的simulink模型 观察其单位阶跃响应曲线 并且分析这三个参数分别对控制性能的影响 解使用期望特性法来设计pid控制器 假设pid控制器的传递函数为 系统闭环的传递函数为 不妨假设希望闭环极点为 和 则期望特征多项式为 对应系数相等 可求得 在commandwindow中输入这3个参数值 建立该系统的simulink模型如下 系统转速响应曲线如图所示 3 3使用ziegler nichols经验整定公式进行pid控制器设计 ziegler nichols经验整定公式是针对被控对象模型为带有延迟的一阶惯性传递函数提出的 例5 如图7 19所示的系统 被控对象为一个带有延迟的惯性环节 试用ziegler nichols经验整定公式 计算pid控制器的参数 并且绘制其仿真系统单位阶跃响应曲线 解由该系统传递函数可知 由ziegler nichols经验整定公式 可得 pid控制器的传递函数为 matlab simulink模型如图所示 仿真结果如图所示 可见系统可以稳定工作 4基于状态空间模型的控制器设计方法 状态空间表达式模型是最新型与最科学的描述方法 它能够全面地表达系统的全部状态信息 它不仅可以描述线性系统 而且可以描述非线性系统 状态空间模型既能够描述单输入单输出 siso 系统 也能够描述多输入多输出 mimo 系统 2 状态变量 动力学系统的状态变量是确定动力学系统状态的最小一组变量 3 状态向量 如果完全描述一个给定系统的动态行为需要n个状态变量 那么可以将这些状态变量看作是向量x t 的各个分量 即 则称为n维状态向量 4 状态空间 以各状态变量为坐标轴所组成的n维空间称为状态空间 在某一时刻的状态向量则可以用状态空间的某一个点来表示 5 状态空间表达式 描述系统输入 输出和状态变量之间关系的方程组称为系统的状态空间表达式 对于线性定常系统而言 具有以下形式 非齐次状态方程的解 例6 已知线性系统齐次状态方程如下 求系统状态方程的解 解用以下matlab程序计算齐次状态方程的解 程序执行后 表示 例7 已知系统状态方程如下 求系统状态方程的解 解用以下matlab程序计算状态方程的解 程序执行结果为 表示 4 2状态反馈极点配置控制器设计 线性系统是状态能控时 可以通过状态反馈来任意配置系统的极点 把极点配置到s左半平面所希望的位置上 则可以获得满意的控制特性 状态反馈的系统方程为 其中 p为一个行向量 其各分量为所希望配置的各极点 即 该命令计算出状态反馈阵k 使得 a bk 的特征值为向量p的各个分量 例8 线性控制系统的状态方程为 要求确定状态反馈矩阵 使状态反馈系统极点配置为 解首先判断系统的能控性 输入以下语句 执行结果为 说明系统能控性矩阵满秩 系统能控 可以应用状态反馈 任意配置极点 输入以下语句 语句执行结果为 用matlab simulink构造这一状态反馈控制系统模型如图所示 仿真