ljiewding.m

详细介绍了Mathlab的使用方法与技巧

function K=ljiewding(Gsys,Hsys)
%系统临界稳定K值的判定，由函数ljiewding和chdihshu组成。Gsys：前向通道传递函数；
%Hsys：反馈通道传递函数；Gsys、Hsys：TF模型
k=0.0001:1:1000;nn=length(k);        %K取值以1为步长，初步判定临界不稳定K的值
for i=1:1:800
   p=chdihshu(Gsys,Hsys,k(i));    %当系统增益为k(i)时，求系统特征方程的根向量
   ii=find(real(p)>=0);n=length(ii);
   if(n>0),break,end                                %若找到K值，退出循环
end
if n==0                                        %若没有找到位于S右半平面的根，则n=0
   disp('系统稳定K的取值范围为K>0');
else if k(i)==0.0001                 %若K=0.0001时，系统不稳定，判定系统不稳定
      disp('系统在K>0取值范围不稳定');
else if k(i)==1.0001                 %若K=0.0001时，系统不稳定，在[0.001，1.001]
   k=0.0001:0.001:k(i);nn=length(k);            %区间搜索使系统临界稳定的K值
   for i=1:nn
      p=chdihshu(Gsys,Hsys,k(i));  
      ii=find(real(p)>=0);n=length(ii);
      if(n>0),disp('系统稳定临界值');k(i);break
   end;end
else                                      %若K>0.001时，系统才不稳定，在[K-1,K+1]
   k=k(i)-1:0.001:k(i)+1;nn=length(k);          %区间搜索使系统临界稳定的K值 
   for i=1:nn
      p=chdihshu(Gsys,Hsys,k(i));    
      ii=find(real(p)>=0);n=length(ii);
      if(n>0),disp('系统稳定临界值');k(i);break
   end,end
end;end;end
K=k(i);                                              %返回系统稳定临界值K