قلم

برای به دست آوردن مشخصات حالت گذرای یک سیستم مرتبه دوم دلخواه با استفاده از برنامه متلب می توانید از کد زیر استفاده نمایید :

zita=input('please enter zita :');
wn=input('please enter wn :');
T=1./(zita.*wn);
disp('T : ')
disp(T)
Ts=4.*T;
disp('Ts : ')
disp(Ts)
t=0:0.001:10;
h=1-(1./sqrt(1-zita.^2)).*exp(-zita.*wn.*t).*sin(wn.*t.*sqrt(1-zita.^2)+atan((sqrt(1-zita.^2))./zita));
Tr=t(min(find(h==1)));
disp('Tr : ')
disp(Tr)
Mpt=1+exp((-pi.*zita)/(sqrt(1-zita.^2)));
disp('Mpt : ')
disp(Mpt)
Tp=pi./(wn.*sqrt(1-zita.^2));
disp('Tp : ')
disp(Tp)
PO1=100.*(Mpt-1);
disp('PO1 : ')
disp(PO1)
PO2=100.*exp((-pi.*zita)/(sqrt(1-zita.^2)));
disp('PO2 : ')
disp(PO2)
plot(t,h)