clear
clc
close all
%輸入典型二階系統(tǒng)參數(shù)，確定系統(tǒng)傳遞函數(shù)模型
alph=input('please input alph=');
wn=input('please input wn=');
num=wn^2;
den=[1 2*alph*wn wn^2];
%判斷系統(tǒng)是否穩(wěn)定
[z,p,k]=tf2zp(num,den);
ii=find(real(z)>0);
n1=length(ii);
jj=find(real(p)>0);
n2=length(jj);
if(n2>0)
   disp('the system is unstable')
   disp('it is no use for getting 動態(tài)參數(shù)')
   else
      %調(diào)用求取二階系統(tǒng)階躍響應動態(tài)參數(shù)的函數(shù)文件
      [y,x,t]=step(num,den);
      plot(t,y)
      [pos,tp,tr,ts2]=stepchar(y,t)
      
end