《自动控制原理实验用Matlab软件编制劳斯判据程序并解题(《学习辅导》例4.3.5)》由会员分享，可在线阅读，更多相关《自动控制原理实验用Matlab软件编制劳斯判据程序并解题(《学习辅导》例4.3.5)(8页珍藏版)》请在人人文库网上搜索。

1、上海电力学院实验报告自动控制原理实验课程主题：用Matlab软件编制劳斯的标准程序，解决问题(学习辅导例4.3.5 )。班级：名字：学号：时间： 2012年11月4日自动工程学院自动化(电站自动化)专业实验报告目录。一、问题的说明二、理论方法分析三、实验设计和实现3四、实验结果和分析5五、结论和讨论六、实验须知6七、参考文献7八附录七一、问题的说明使用MATLAB建立劳斯标准，给出其劳斯矩阵，确定相应系统的稳定性二、理论方法的分析使用m文件实现Matlab编程。1)M文件的制作和调用从Matlab桌面上的“文件”菜单中选择“New”菜单项，然后选择“M-file”命令，将显示Matlab文本编。

2、辑器窗口。在Matlab命令窗口的“文件”菜单中选择“打开”命令后，将显示“打开”对话框，您可以在“文件名”对话框中选择要打开的m文件名。2)M文件的调试文件编辑器窗口的菜单栏和工具栏下方有三个区域，右侧的大区域是程序窗口。最左边的区域显示行号，每行包含数字，包含空白行，行号自动出现。 在行号的增加而增加的行号和程序窗口之间的区域中，只在可执行行中有细横线，在非执行行(如空白行、注释行、函数定义行)之前什么也没有。 在调试程序时，可以直接在这些程序上单击鼠标来设置或删除断点。三、实验设计和实现(1)程序%RouthMatrix*确定劳斯矩阵(带参数的特征多项式)和对应系统的稳定性*清除；sym。

3、s k z q %定义变量k z qp=input (特征多项式的参数=)； 系统提示输入%参数得到n=length(p) %p的长度for i=0:ceil(n/2)-1 %多项式按劳斯矩阵排序a(1，i 1)=p(2*i 1)PK2* (P1 ) na(2，i 1)=0；break结束a(2，i 1)=p(2*(i 1) )；结束for k=3:n %计算从第三行开始的劳斯矩阵的内容for j=1:ceil(n-k 1)/2 )判定ifsa (k-1，1 )=0%共轭虚根的有无disp (系统中有共轭虚根)中断符号=1；break结束(k，j )=(a (k-1，1 ) * a (k-2，。

4、j 1)-a(k-1，J1 ) * a (k-2，1 )/a (k-1，1 )；结束结束支持disp (劳斯矩阵) %输出的劳斯矩阵disp (双倍(a ) )for i=3:k %是通过劳斯的判断标准判断系统的稳定性的IFA (I-1，1 )=0%判断第一列的元素是否小于或等于0q=1；break结束结束PS=1disp (系统不稳定)else输出disp (系统稳定) %系统稳定性判定结果结束(2)设计劳斯表的内容是根据书的原理设计的，把多项式排列在劳斯矩阵上，可以得到其前两行的内容，剩下的几行必须根据书的原理写。 是否其稳定的判断取决于第1列元素是否有符号的变化，如果从第2行第1列的元素。

5、大于0则判断endfor循环、i 1、第3行第1列的元素，以下，也将q分配1，直到某行第1列的元素成为零以下，并跳过循环四、实验结果和分析如图所示的m文件代码：创建m文件后，可以保存并运行m文件，然后通过MATLAB接口进行调试。自由选择特征方程式的系数，执行后的结果如下所示五、结论和讨论上述实验过程表明，劳斯的判断标准的前两行是根据特征方程的参数写的，接下来的几行是根据前两行计算的，可以用劳斯的判断标准来判断稳定性，有零行则有共轭虚根，第一列的元素符号有变化，系统就可以如果有未知的参数，也可以判断该值的范围。六、实验体会实践是检验真理的唯一标准，死板的书交给我们的是乍一看的计算方法，实践不仅。

6、反馈理论的实践者，还反馈理论学习的本质。 不是第一次听MATLAB，而是在自习原理的课上老师也使用它展示了几个程序。 那时，MATLAB留下了非常强大，操作非常简单的好印象。 之后的实验课自己操作，觉得那只是按照老师说的顺序一步一步地操作。 老实说，看到这个主题的瞬间有点茫然。 幸而主题可以作为参考，解读了书上写的程序后，自己摸索练习，终于用MATLAB做了完整的主题。 实际感受到了这个软件的实用性。 这次实验接触了只有功能模块的冰山一角，只学习了基本操作，只是初步理解了这个软件，但对更高级的东西还不熟悉。 虽然并不是所有用于报告的编程都学到了，但是我已经切实感受到了MATLAB的魅力。 我必。

7、须佩服好的程序软件有助于学习自动控制的原理。 好的开始是成功的一半，今后的学习将理论和实践结合起来，为MATLAB提供学习自主原理的服务。七、参考文献自动控制原理学习辅导杨平，如翁所示，编着王志萍。中国电力出版社如杨平、馀洁、冯照坤、翁所示编纂。中国电力出版社八、附录%RouthMatrix*确定劳斯矩阵(带参数的特征多项式)和对应系统的稳定性*清除；syms k z q %定义变量k z qp=input (特征多项式的参数=)； 系统提示输入%参数得到n=length(p) %p的长度for i=0:ceil(n/2)-1 %多项式按劳斯矩阵排序a(1，i 1)=p(2*i 1)PK2* 。

8、(P1 ) na(2，i 1)=0；break结束a(2，i 1)=p(2*(i 1) )；结束for k=3:n %计算从第三行开始的劳斯矩阵的内容for j=1:ceil(n-k 1)/2 )判定ifsa (k-1，1 )=0%共轭虚根的有无disp (系统中有共轭虚根)中断符号=1；break结束(k，j )=(a (k-1，1 ) * a (k-2，j 1)-a(k-1，J1 ) * a (k-2，1 )/a (k-1，1 )；结束结束支持disp (劳斯矩阵) %输出的劳斯矩阵disp (双倍(a ) )for i=3:k %是通过劳斯的判断标准判断系统的稳定性的IFA (I-1，1 )=0%判断第一列的元素是否小于或等于0q=1；break结束结束PS=1disp (系统不稳定)else输出disp (系统稳定) %系统稳定性判定结果结束8。