[求助]用Routh判据求系统稳定性,用MATLAB怎么编程

zhaoqinquan

用Routh判据求系统稳定性,涉及到不等式的求解

如 A<B-C*D,  已知B,D, 且是常数, C是变量, C=1:20:400,    画出A和C 的关系图

用MATLAB怎么编

[ 本帖最后由 xmwhit 于 2007-10-21 21:18 编辑 ]

vib

给出个实例吧

bainhome

你给的是系统特征方程式吗？
请教这跟routh判据有什么关系？

zhengbc

好像有一本书编了这个程序.......,具体忘记了,唉........是关于Matlab在控制中的用法的书.......

realone

function [rtab,info]=routh0(den)
info=[];
vec1=den(1:2:length(den)); nrT=length(vec1);
vec2=den(2:2:length(den)-1);
rtab=[vec1; vec2, zeros(1,nrT-length(vec2))];
for k=1:length(den)-2,
   alpha(k)=vec1(1)/vec2(1);
   for i=1:length(vec2),
       a3(i)=rtab(k,i+1)-alpha(k)*rtab(k+1,i+1);
   end
   if sum(abs(a3))==0
      a3=polyder(vec2);
      info=[info,'All elements in row ',...
            int2str(k+2) ' are zeros;'];
   elseif abs(a3(1))<eps
      a3(1)=1e-6;
      info=[info,'Replaced first element;'];
   end
   rtab=[rtab; a3, zeros(1,nrT-length(a3))];
   vec1=vec2; vec2=a3;
end

邓庆文

回复 5 # realone 的帖子

劳驾帮忙逐句解析下，谢谢