毕业论文---基于Matlab的系统的稳定性分析.doc
《毕业论文---基于Matlab的系统的稳定性分析.doc》由会员分享,可在线阅读,更多相关《毕业论文---基于Matlab的系统的稳定性分析.doc(21页珍藏版)》请在三一文库上搜索。
1、 编号编号 学学士士学学位位论论文文 基于基于 MatlabMatlab 的系统的稳定性分的系统的稳定性分 析析 学生姓名: 刘永辉 学 号: 20050104024 系 部: 物理系 专 业: 物理学 年 级: 2005 级 指导教师: 黄晓俊 完成日期: 年 月 日 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 摘要摘要 稳定性在系统的实际应用中非常的重要,本文介绍了系统的稳定性的概念, 论述了常用判定系统的稳定性的方法:奈奎斯特判据、根轨迹法、波特图法等, 也介绍了罗斯矩阵、朱里矩阵在稳定性分析中的作用。应用 MATLAB 编程来实 现奈奎斯特判据、根轨迹法、B
2、ode 图对稳定性的分析。 关键词:关键词:LTI 系统;稳定性;系统;稳定性;MATLAB Matlab-based analysis of system stability Abstract The stability of the systems practical application is very important, this paper introduces the concept of stability of the system, discusses the stability of the system used to determine the method: Ny
3、quist criterion, root locus method, Bode plots, such as law, Rose also introduced the matrix, where matrix Zhu at the role of stability analysis. Application of MATLAB programming to achieve the Nyquist criterion, root locus method, Bode diagram of the stability analysis. Key words:LTI system; stabi
4、lity; MATLAB 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 2 目 录 摘要 1 ABSTRACT1 引言 2 1.理论分析 2 1.1 概述2 1.1.1 MATLAB 语言介绍 2 1.1.2 LTI 系统的稳定性 3 1.2 LTI 连续时间系统的稳定性分析3 1.2.1 因果连续时间系统的稳定性准则 3 1.2.2 连续时间 LTI 反馈系统的奈奎斯特判据 6 1.3 LTI 离散时间系统的稳定性分析 .9 1.3.1 因果离散时间系统的稳定性准则 9 1.3.2 离散时间 LTI 反馈系统的奈奎斯特判据 11 2.基于 MTLAB 的稳定性分析
5、13 2.1 奈奎斯特图14 2.2 根轨迹15 2.3 波特图.16 3.结论 18 4.结语 18 致谢 19 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 2 引言引言 线性时不变系统通常被称为 LTI 系统,系统在不同的情况下有不同的函数 表达式。系统的稳定性对系统的输入输出行为至为重要。若系统稍微偏离其平 衡态,就可能会产生几种情况;若系统保持在平衡状态附近,则称系统是稳定 的;如果系统趋于返回平衡状态或一个极限状态,则称此系统为不稳定的。因 此,研究系统的稳定性的方法称为稳定性判据或稳定判据,如劳斯判据,胡尔 维茨(Hurwirz)稳定判据以及奈魁斯特稳定
6、判据等,在 MATLAB 未产生前,由 于系统的复杂性,判别计算量非常大,而用了 MATLAB 以后,稳定性分析将变的 很简单。 1.1.理论分析理论分析 线性系统的稳定性取决于系统本身的结构和参数,而与输出无关。线性系 统稳定的条件是其特征根均具有负实部。在实际工程系统中,为了避开对特征 方程的直接求解,就只好讨论特征根的分布,即看其是否全部具有负实部,并 以此来判别系统的稳定性。 1.11.1 概述概述 1.1.11.1.1 MATLABMATLAB 语言介绍语言介绍 MATLAB 是 MATrix LABoratory 的缩写,是 MathWork 公司于 1984 年推出的 一套面向工
7、程和科学运算的高性能软件。它具有强大的矩阵计算能力和良好的 图形可视化功能,为用户提供了非常直观和简洁的程序开发环境,因此被城为 第四代计算机语言。 MATLAB 发展至今,现已集成了许多工具箱,如控制系统集成箱(Control System Toolbox) 、信号处理工具箱(Single Processing Toolbox) 、模糊推理 系统工具箱(Fuzzy Logic Toolbox) 、Simulink 工具箱等。为此。MATLAB 语言 在控制工程领域已获得了广泛的应用。 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 3 1.1.21.1.2 LTILT
8、I 系统的稳定性系统的稳定性 LTI 系统的稳定性与其系统函数或有着密切的关系。)(sH)(zH 一个连续时间 LTI 系统,其冲激响应满足时,,而其系统函数0t0)(th 的 ROC 一定是 S 平面的右半部分: 。一个稳定连续时间 LTI 系)(sH 0 Res 统的充要条件是其单位冲激响应绝对可积。即:)(th (1) dtth )( 对应于系统函数则是其 ROC 包含轴。结合以上两种结果,可得稳)(sHj 定连续时间 LTI 系统,其系统函数的所以极点的实部都必须是负的。)(sH 离散时间 LTI 系统,也有类似的结果: (1)因果系统的充要条件是单位脉冲响应满足,其系统)(nh0,
9、0)(nnh 函数 H(z)的 ROC 为某内界圆的外部,即; 0 rz (2)稳定系统的充要条件是其单位脉冲响应绝对可和,或 k nh 系统函数 H(z)的 ROC 包含单位圆; (3)因果稳定离散时间系统 LTI 系统的系统函数 H(z)的所有极点必须落在 单位圆内部。因此,可以通过系统函数,很方便地了解系统的稳定性。不仅如 此,系统函数已经成为系统分析和综合的基本方法。 1.21.2 LTILTI 连续时间系统的稳定性分析连续时间系统的稳定性分析 1.2.11.2.1 因果连续时间系统的稳定性准则因果连续时间系统的稳定性准则 因果连续时间系统的系统函数 (2) )( )( )( sA s
10、B sH 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 4 式中 (3) 01 2 1 1 )(asasasasA n n n n 的极点就是的根,因此为判断系统是否稳定,亦即的极)(sH0)(sA)(sH 点是否都在左半开平面,只需判断的根,即特征根是否都在左半开平面,0)(sA 并不需要知道各特征根的确切位置。所有根均在左半开平面的多项式称为霍尔 维兹多项式。罗斯和霍尔维兹提出了判别多项式是否为霍尔维兹多项式的准则, 称为罗斯- 霍尔维兹准则。 对于特征根为实根和共轭复根,多项式可分解为许多一次aja )(sA 因子和二次因子的乘积。如果特征根都在左半开平面,则要
11、as 22 )( as 求各因子中,从而多项式的所有系数。也就是说,0a)(sA), 2 , 1 , 0(0niai 如果中任何一个或多个系数为零或负值,那么它就不是霍尔维兹多项式。)(sA 上述条件是必要条件,而不是充要条件。罗斯提出了一种列表的方法,常城罗 斯阵列。其方法如下表所示,将多项式的系数按下表的规律排列在 1,2 行)(sA 表 1 罗斯阵列 行 1n a 2n a 4n a 21n a 3n a 5n a 31n c 3n c 5n c 41n d 3n d 5n d n+1 罗斯阵列中第 3 行及以后的各行,按以上规则计算, , (4) 2 1 1 13 1 nn n n n
12、n aa c a aa 51 4 1 3 1 nn nn n n aa aa a c 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 5 , , (5) 31 31 1 1 1 nn nn n n cc aa c d 51 51 1 3 1 nn nn n n cc aa c d 依次类推,一直排列到第 n+1 行(以后各行为零) 。 罗斯准则指出:多项式是霍尔维兹多项式的充分和必要条件是罗斯阵)(sA 列中第一列元素的值均大于零,它保证了的根都在左半开平面。如果第0)(sA 一列元素的符号不完全相当,那么变号的次数就是在右半平面根的数目。 对于二阶系统 , (6) 2
13、 210 ( )A sa sa sa 若,根据上述稳定准则,可得为霍尔维兹多项式的冲要条件为 2 0a )(sA (7) 10 0,0aa 例 1:已知某系统的系统函数为的系数排列成罗斯阵列 32 ( )331H ssssk 13 31 8 3 1 k k k 如果系统是稳定的,根据罗斯准则,以上阵列中的第一列元素的值为正值, 即 和 8 0 3 k 10k 解得 18k 因此,当时,系统是稳定的。18k 1.2.21.2.2 连续时间连续时间 LTILTI 反馈系统的奈奎斯特判据反馈系统的奈奎斯特判据 根据反馈系统的稳定性,要求 (8))()(1sHsKG 学学 士士 学学 位位 论论 文文
14、 BACHELOR S THESIS 6 或等效函数 (9))()( 1 )(sHsG K sR 在 S 复平面的左半平面内没有零点。因此,可以考虑如图所示的一条半圆 围线。当 s 沿这条围线 C 顺时针旋转一周时,由的轨迹围线顺时针绕原点( )R s 的次数,可得出围线 C 内所包括的的零点个数和极点个数的差值。随着 M( )R s 增加至无穷大值时,对应的围线 C 就是沿 轴从到的半径为的半j 圆曲线,此时围线 C 包括了整个右半平面,且围线 C 变为整个虚轴。j o C M Im -jM jM Re 图1 包括右半平面的半圆闭合围线,当M 时,该围线包括了整个右半平面 为了保证随 M 增
15、加,围线 C 的半圆延伸至整个右半平面时,仍然是有( )R s 界的。该条件要求的极点数要大于等于它的零点数。这时( )R s 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 7 (10) n n n n n n n n n n ssM a b asasa bsbsb sRsR 0 1 1 0 1 1 lim)(lim)(lim 为常数。 当极点阶数大于零点阶数时,上述值也为零。因此,当 M 增大到无穷( )R s 大是。沿着这个围线 C 半圆部分的值不再变化,为一常数。)(sR 当时,图 1 所示的围线 C 与虚轴轴重合,对应的图就是当Mj)(sR 从变到时的图。如果
16、正向和反馈通路的系统函数是稳定的,那)(jR 么和分别是这两支路系统的频率响应函数。)(jG)(jH 注意到的围线只是复变函数的一个性质,不涉及 ROC 的问题。这样,)(sF 即使正向和反馈通路的系统不稳定,也可用上述方法,检查在)(jR 范围内的图,用于计算位于右半平面内的零点数和极点数之差。)(sR 再者,由式(9)可知,绕原点的次数,就是绕点-1/K 的次)(sR)()(sHsG 数,即绕原点的次数,就是围绕)(jR)()(jG)(jH)( 点-1/K 的次数。当从时,的图就称为奈奎斯特图。)(jG)(jH 注意到,的极点就是的极点,而的零点是闭环极点。因此,)(sR)(sG)(sH)
17、(sR 根据围线映射性质可得如下结论。 奈奎斯特图顺时钟绕-1/K 点的净次数等于右半平面内闭环极点数减去)(sG 在右半平面内的极点数。)(sH 由上述结果可得,如果反馈系统是稳定的,那么奈奎斯特图瞬时绕-1/K 点 的净次数等于在右半平面内的极点数,且是逆时针方向的。由此,就)(sG)(sH 学学 士士 学学 位位 论论 文文 BACHELOR S THESIS 8 可得出连续时间奈奎斯特稳定性判决。 例 2:设,试画出奈奎斯特图并确定使反馈系统 ) 15 . 0)(1( 1 )()( ss s sHsG 稳定的 K 的取值范围。 解:从的表示式,可知)()(jHjG (1)时, ;01)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 毕业论文 基于 Matlab 系统 稳定性 分析
链接地址:https://www.31doc.com/p-3941323.html