自动控制原理实验用Matlab软件编制劳斯判据程序并解题(《学习辅导》例4.3.5)

上 海 电 力 学 院实 验 报 告自 动 控 制 原 理 实 验 课 程题目：用Matlab软件编制劳斯判据程序并解题(学习辅导例4.3.5)班级： 姓名： 学号： 时间： 2012年11月4日 自动化工程学院自动化（电站自动化）专业实 验 报 告 目 录一、 问题描述3二、 理论方法分析3三、 实验设计与实现3四、 实验结果与分析5五、 结论与讨论6六、 实验心得体会6七、 参考文献7八 附录7一、问题描述用MATLAB编制劳斯判据列出其劳斯矩阵并判断相对应系统的稳定性二、理论方法分析采用M文件实现Matlab编程。1)M文件的建立与调用 从Matlab操作桌面的“File”菜单中选择“New”菜单项，再选择 “M-file”命令，屏幕将出现Matlab文本编辑器的窗口。 在Matlab命令窗口的“File”菜单中选择“Open”命令，则屏幕出现“Open”对话框，在文件名对话框中选中所需打开的M文件名。2)M文件的调试在文件编辑器窗口菜单栏和工具栏的下面有三个区域，右侧的大区域是程序窗口，用于编写程序；最左面区域显示的是行号，每行都有数字，包括空行，行号是自动出现的，随着命令行的增加而增加；在行号和程序窗口之间的区域上有一些小横线，这些横线只有在可执行行上才有，而空行、注释行、函数定义行等非执行行的前面都没有。在进行程序调试时，可以直接在这些程序上点击鼠标以设置或去掉断点。三、实验设计与实现（1）程序%RouthMatrix*劳斯矩阵(带参数的特征多项式)并判断对应系统稳定性*clear;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); if 2*(i+1)n a(2,i+1)=0; break end a(2,i+1)=p(2*(i+1);endfor k=3:n %计算从第三行开始劳斯矩阵内容 for j=1:ceil(n-k+1)/2) if a(k-1,1)=0 %判断是否有共轭虚根 disp(系统有共轭虚根) breaksign=1; break end a(k,j)=(a(k-1,1)*a(k-2,j+1)-a(k-1,j+1)*a(k-2,1)/a(k-1,1); endenddisp(劳斯矩阵) %输出对应的劳斯矩阵disp(double(a)for i=3:k %用劳斯判据判断系统的稳定性 if a(i-1,1)n a(2,i+1)=0; break end a(2,i+1)=p(2*(i+1);endfor k=3:n %计算从第三行开始劳斯矩阵内容 for j=1:ceil(n-k+1)/2) if a(k-1,1)=0 %判断是否有共轭虚根 disp(系统有共轭虚根) breaksign=1; break end a(k,j)=(a(k-1,1)*a(k-2,j+1)-a(k-1,j+1)*a(k-2,1)/a(k-1,1); endenddisp(劳斯矩阵) %输出对应的劳斯矩阵disp(double(a)for i=3:k %用劳斯判据判断系统的稳定性 if a(i-1,1)=0 %判断第一列元素是