[求助]用Routh判据求系统稳定性,用MATLAB怎么编程
用Routh判据求系统稳定性,涉及到不等式的求解如 A<B-C*D,已知B,D, 且是常数, C是变量, C=1:20:400, 画出A和C 的关系图
用MATLAB怎么编
[ 本帖最后由 xmwhit 于 2007-10-21 21:18 编辑 ] 给出个实例吧 你给的是系统特征方程式吗?
请教这跟routh判据有什么关系? 好像有一本书编了这个程序.......,具体忘记了,唉........是关于Matlab在控制中的用法的书....... function =routh0(den)
info=[];
vec1=den(1:2:length(den)); nrT=length(vec1);
vec2=den(2:2:length(den)-1);
rtab=;
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=;
end
rtab=;
vec1=vec2; vec2=a3;
end 回复 5 # realone 的帖子
劳驾帮忙逐句解析下,谢谢
页:
[1]