clc;
clear all;
num1 = 1;
den1 = [9e-6 9.5e-3 1]
sys1 = tf(num1,den1)
step(sys1,1.2);
hold on;
num2 = 1;
den2 = [4.2e-5 0.037 1]
sys2 = tf(num2,den2)
step(sys2,1.2);
hold on;
num3 = 1;
den3 = [3.92e-4 0.074 1]
sys3 = tf(num3,den3)
step(sys3,1.2);
hold on;
info = stepinfo(sys1);
disp(info);
set_Time = info.SettlingTime;
disp(['Settling Time TF1: ' num2str(set_Time) 'seconds']);
deviation = 1-0.98;
err = deviation*100;
disp(['Steady-State error TF1: ' num2str(err) '%']);
max_OS=((info.Peak-0.98)/0.98)*100;
disp(['Maximum Overshoot TF1: ' num2str(max_OS) '%']);
info2 = stepinfo(sys2);
disp(info2);
set_Time2 = info2.SettlingTime;
disp(['Settling Time TF2: ' num2str(set_Time2) 'seconds']);
deviation2 = 1-0.977;
err2 = deviation2*100;
disp(['Steady-State error TF2: ' num2str(err2) '%']);
max_OS2=((info2.Peak-0.977)/0.977)*100;
disp(['Maximum Overshoot TF2: ' num2str(max_OS2) '%']);
info3 = stepinfo(sys3);
disp(info3);
set_Time3 = info3.SettlingTime;
disp(['Settling Time TF3: ' num2str(set_Time3) 'seconds']);
deviation3 = 1-0.982;
err3 = deviation3*100;
disp(['Steady-State error TF3: ' num2str(err3) '%']);
max_OS3=((info3.Peak-0.982)/0.982)*100;
disp(['Maximum Overshoot TF3: ' num2str(max_OS3) '%']);