《线性系统时域响应分析.doc》由会员分享，可在线阅读，更多相关《线性系统时域响应分析.doc（9页珍藏版）》请在金锄头文库上搜索。
1、实验实验二二 线线性系性系统时统时域响域响应应分析分析 一 一 实验实验目的目的 1 熟练掌握 step 函数和 impulse 函数的使用方法 研究线性系统在 单位阶跃 单位脉冲及单位斜坡函数作用下的响应 2 通过响应曲线观测特征参量和对二阶系统性能的影响 n 二 基二 基础础知知识识及及 MATLAB 函数函数 一 基础知识 时域分析法直接在时间域中对系统进行分析 可以提供系统时间响应的全 部信息 具有直观 准确的特点 为了研究控制系统的时域特性 经常采用瞬 态响应 如阶跃响应 脉冲响应和斜坡响应 本次实验从分析系统的性能指标 出发 给出了在 MATLAB 环境下获取系统时域响应和分析系统
2、的动态性能和稳态 性能的方法 用 MATLAB 求系统的瞬态响应时 将传递函数的分子 分母多项式的系数分 别以 s 的降幂排列写为两个数组 num den 由于控制系统分子的阶次 m 一般小 于其分母的阶次 n 所以 num 中的数组元素与分子多项式系数之间自右向左逐 次对齐 不足部分用零补齐 缺项系数也用零补上 1 用 MATLAB 求控制系统的瞬态响应 1 阶跃响应 求系统阶跃响应的指令有 step num den 时间向量 t 的范围由软件自动设定 阶跃响应曲线 随即绘出 step num den t 时间向量 t 的范围可以由人工给定 例如 t 0 0 1 10 y x step nu
3、m den 返回变量 y 为输出向量 x 为状态向量 在 MATLAB 程序中 先定义 num den 数组 并调用上述指令 即可生成单位 阶跃输入信号下的阶跃响应曲线图 考虑下列系统 254 25 2 sssR sC 该系统可以表示为两个数组 每一个数组由相应的多项式系数组成 并且以 s 的降幂排列 则 MATLAB 的调用语句 num 0 0 25 定义分子多项式 den 1 4 25 定义分母多项式 step num den 调用阶跃响应函数求取单位阶跃响应曲 线 grid 画网格标度线 xlabel t s ylabel c t 给坐标轴加上说明 title Unit step Res
4、pinse of G s 25 s 2 4s 25 给图形加上标题名 则该单位阶跃响应曲线如图 2 1 所示 为了在图形屏幕上书写文本 可以用 text 命令在图上的任何位置加标注 例如 text 3 4 0 06 Y1 和 text 3 4 1 4 Y2 第一个语句告诉计算机 在坐标点 x 3 4 y 0 06 上书写出 Y1 类似地 第二个语句告诉计算机 在坐标点 x 3 4 y 1 4 上书写出 Y2 若要绘制系统 t 在指定时间 0 10s 内的响应曲线 则用以下语句 num 0 0 25 den 1 4 25 t 0 0 1 10 step num den t 即可得到系统的单位阶跃
5、响应曲线在 0 10s 间的部分 如图 2 2 所示 2 脉冲响应 求系统脉冲响应的指令有 impulse num den 时间向量 t 的范围由软件自动设定 阶跃响应曲 图 2 1 二阶系统的单位阶跃响应 图 2 2 定义时间范围的单位阶跃响应 线随即绘出 impulse num den t 时间向量 t 的范围可以由人工给定 例如 t 0 0 1 10 y x impulse num den 返回变量 y 为输出向量 x 为状态向量 y x t impulse num den t 向量 t 表示脉冲响应进行计算的时间 例 试求下列系统的单位脉冲响应 12 0 1 2 ss sG sR sC
6、 在 MATLAB 中可表示为 num 0 0 1 den 1 0 2 1 impulse num den grid title Unit impulse Response of G s 1 s 2 0 2s 1 由此得到的单位脉冲响应曲线如图 2 3 所示 求脉冲响应的另一种方法 应当指出 当初始条件为零时 G s 的单位脉冲响应与 sG s 的单位阶跃响 应相同 考虑在上例题中求系统的单位脉冲响应 因为对于单位脉冲输入量 R s 1 所以 sss s ss sGsC sR sC1 12 012 0 1 22 因此 可以将 G s 的单位脉冲响应变换成 sG s 的单位阶跃响应 图 2 3
7、二阶系统的单位脉冲响应 向 MATLAB 输入下列 num 和 den 给出阶跃响应命令 可以得到系统的 单位脉冲响应曲线如图 2 4 所示 num 0 1 0 den 1 0 2 1 step num den grid title Unit step Response of sG s s s 2 0 2s 1 3 斜坡响应 MATLAB 没有直接调用求系统斜坡响应的功能指令 在求取斜坡响应时 通常利用阶跃响应的指令 基于单位阶跃信号的拉氏变换为 1 s 而单位斜坡信 号的拉氏变换为 1 s2 因此 当求系统 G s 的单位斜坡响应时 可以先用 s 除 G s 再利用阶跃响应命令 就能求出系统
8、的斜坡响应 例如 试求下列闭环系统的单位斜坡响应 1 1 2 sssR sC 对于单位斜坡输入量 R s 1 s2 因此 sssssss sC 1 1 11 1 1 222 在 MATLAB 中输入以下命令 得到如图 2 5 所示的响应曲线 num 0 0 0 1 den 1 1 1 0 step num den title Unit Ramp Response Cuve for System G s 1 s 2 s 1 图 2 4 单位脉冲响应的另一种表示法 图 2 5 单位斜坡响应 2 特征参量和对二阶系统性能的影响 n 标准二阶系统的闭环传递函数为 22 2 2 nn n sssR sC
9、 二阶系统的单位阶跃响应在不同的特征参量下有不同的响应曲线 1 对二阶系统性能的影响 设定无阻尼自然振荡频率 考虑 5 种不同的值 1srad n 0 0 25 0 5 1 0 和 2 0 利用 MATLAB 对每一种求取单位阶跃响应曲线 分 析参数对系统的影响 为便于观测和比较 在一幅图上绘出 5 条响应曲线 采用 hold 命令实现 num 0 0 1 den1 1 0 1 den2 1 0 5 1 den3 1 1 1 den4 1 2 1 den5 1 4 1 t 0 0 1 10 step num den1 t grid text 4 1 7 Zeta 0 hold step num
10、 den2 t text 3 3 1 5 0 25 step num den3 t text 3 5 1 2 0 5 step num den4 t text 3 3 0 9 1 0 step num den5 t text 3 3 0 6 2 0 title Step Response Curves for G s 1 s 2 2 zeta s 1 由此得到的响应曲线如图 2 6 所示 2 对二阶系统性能的影响 n 同理 设定阻尼比时 当分别取 1 2 3 时 利用 MATLAB 求取25 0 n 单位阶跃响应曲线 分析参数对系统的影响 n num1 0 0 1 den1 1 0 5 1 t
11、 0 0 1 10 step num1 den1 t grid hold on text 3 1 1 4 wn 1 num2 0 0 4 den2 1 1 4 step num2 den2 t hold on text 1 7 1 4 wn 2 num3 0 0 9 den3 1 1 5 9 step num3 den3 t hold on text 0 5 1 4 wn 3 由此得到的响应曲线如图 2 7 所示 图 2 6 不同时系统的响应曲线 图 2 7 不同时系统的响应曲线 n 三 三 实验实验内容内容 1 观察函数 step 调用格式 假设系统的传递函数模型为 1464 73 234
12、2 ssss ss sG 绘制出系统的阶跃响应曲线 2 对典型二阶系统 2 2 2 2 nn n ss sG 1 分别绘出 分别取 0 0 25 0 5 1 0 和 2 0 时的单位阶跃响 2srad n 应曲线 分析参数对系统的影响 2 绘制出当 0 25 分别取 1 2 4 6 时单位阶跃响应曲线 分析参数 n 对系统的影响 n 四 四 实验报实验报告告 1 根据内容要求 写出调试好的 MATLAB 语言程序 及对应的 MATLAB 运算 结果 2 记录各种输出波形 根据实验结果分析参数变化对系统的影响 3 写出实验的心得与体会 五 五 预习预习要求要求 1 预习实验中基础知识 运行编制好的 MATLAB 语句 熟悉 MATLAB 指令及 step 和 impulse 函数 2 结合实验内容 提前编制相应的程序 3 思考特征参量和对二阶系统性能的影响 n
