%本程序是用小干擾法判斷系統(tǒng)的靜態(tài)穩(wěn)定
S0=input('請(qǐng)輸入初始功率:S0=');
V0=input('請(qǐng)輸入無(wú)限大系統(tǒng)母線電壓:V0=');
Xd=input('請(qǐng)輸入系統(tǒng)直軸等值電抗:Xd=');
w0=input('請(qǐng)輸入同步電角速度：w0=');
Tj=input('請(qǐng)輸入慣性時(shí)間常數(shù):Tj=');
zn=input('請(qǐng)輸入綜合阻尼系數(shù):D=');
Eq=sqrt((V0+imag(S0)*Xd./V0)^2+(real(S0)*Xd./V0)^2);
dtj0=atan(real(S0)*Xd./(V0*(V0+imag(S0)*Xd./V0)));
Psl=Eq*V0./Xd;
C=w0;
D=-1./Tj*Eq*cos(dtj0)*V0./Xd;
E=w0*zn./(2*Tj);
Kp=(Psl-real(S0))./real(S0);
root1=-E+sqrt(E^2+C*D);
root2=-E-sqrt(E^2+C*D);
disp('該系統(tǒng)的靜態(tài)穩(wěn)定極限:Psl=');
disp(Psl);
disp('該系統(tǒng)的靜態(tài)儲(chǔ)備系數(shù)為Kp=');
disp(Kp);
disp('該系統(tǒng)線性微分方程的特征根為:');
disp(root1);disp(root2);
if zn>=0
    if (root1==conj(root2))|((real(root1)<0)&(real(root2)<0));
        disp('該系統(tǒng)是靜態(tài)穩(wěn)定的。');
    else
         disp('該系統(tǒng)是不穩(wěn)定的。');
     end
 elseif zn<0
      disp('該系統(tǒng)是不穩(wěn)定的。');
  end