Control Stystem PDF
Control Stystem PDF
Control Stystem PDF
num1=[1 2];
den1=[3 1 0];
G1=tf(num1,den1);
G2=tf([2],[1 7]);
G3=tf([1 5],[1 6 3]);
G4=tf([1],[1 0]);
T1=parallel(G1,G2);
T2=series(T1,G3);
T=feedback(T2,G4,-1)
Example - 2 - Laplace Transforms
clc
clear all
syms t n a b;
%% example-1
disp('(a) laplace of t is:');
laplace(t)
%% example-2
disp('(b) laplace of t^2 is:');
laplace(t^2)
%% example-3
disp('(c) laplace of t^3 is:');
laplace(t^3)
%% example-4
disp('(d) laplace of t^n is:');
laplace(t^n)
%% example-5
disp('(e) laplace of exp(a*t) is:');
laplace(exp(a*t))
%% example-6
disp('(f) laplace of t*exp(a*t) is:');
laplace(t*exp(a*t))
%% example-7
disp('(g) laplace of t^2*exp(a*t) is:');
laplace(t^2*exp(a*t))
%% example-8
disp('(h) laplace of t^n*exp(a*t) is:');
laplace(t^n*exp(a*t))
%% example-9
disp('(i) laplace of cos(a*t) is:');
laplace(cos(a*t))
%% example-10
disp('(j) laplace of sin(a*t) is:');
laplace(sin(a*t))
%% example-11
disp('(k) laplace of t*cos(a*t) is:');
laplace(t*cos(a*t))
%% example-12
disp('(l) laplace of t*sin(a*t) is:');
laplace(t*sin(a*t))
%% example-13
disp('(m) laplace of t^2*cos(a*t) is:');
laplace(t^2*cos(a*t))
%% example-14
disp('(n) laplace of t^2*sin(a*t) is:');
laplace(t^2*sin(a*t))
%% example-15
disp('(o) laplace of exp(a*t)*cos(b*t) is:');
laplace(exp(a*t)*cos(b*t))
%% example-16
disp('(p) laplace of exp(a*t)*sin(b*t) is:');
laplace(exp(a*t)*cos(b*t))
%% example-17
disp('(q) laplace of cosh(a*t) is:');
laplace(cosh(a*t))
%% example-18
disp('(r) laplace of sinh(a*t) is:');
laplace(sinh(a*t))
%% example-19
disp('(s) laplace of t*cosh(a*t) is:');
laplace(t*cosh(a*t))
%% example-20
disp('(t) laplace of t*sinh(a*t) is:');
laplace(t*sinh(a*t))
NOTE:
Use ‘ilaplace’ keyword to take inverse laplace transform
NOTE:
n!=Γ(n+1)
Example – 3
clc
clear all
p =
-6.0000
-6.0000
-4.0000
-2.0000
k =
[]
From the MATLAB result, the output expression can be written as,
Example – 4 – Time Domain Specifications
1) Peak time
2) Peak overshoot
3) Rise time
4) Settling time
>> num=130;
>> den=[1 15 130];
>> T=tf(num,den)
T =
130
----------------
s^2 + 15 s + 130
omega_n =
11.4018
>> zeta=den(2)/(2*omega_n)
zeta =
0.6578
>> ts=4/(zeta*omega_n)
ts =
0.5333
>> tp=pi/(omega_n*sqrt(1-zeta^2))
tp =
0.3658
>> tr=(1.76*zeta^3-0.417*zeta^2+1.039*zeta+1)
/omega_n
tr =
0.1758
>> mp=100*exp(-zeta*pi/sqrt(1-zeta^2))
mp =
6.4335
>> step(T)
Method-2
>> sys = tf([130],[1 15 130]);
>> % default 2% settling time
>> % default 10% to 90% rise time
>> s = stepinfo(sys)
s =
RiseTime: 0.1758
SettlingTime: 0.5272
SettlingMin: 0.9037
SettlingMax: 1.0643
Overshoot: 6.4307
Undershoot: 0
Peak: 1.0643
PeakTime: 0.3684
Example – 5 – Bode Plot
Create bode plot of the dynamic system.
>> bode(h)
Example – 6 – Root locus plot
Create root locus plot of the dynamic system.
>> rlocus(h)
Example – 7 – Nyquist plot
Create nyquist plot of the dynamic system.
>> nyquist(h)
Example – 8 – Nicholas plot
Create nicholas plot of the dynamic system.
>> nichols(h)
Example – 9 – Types of signal responses
i) Step response
>> step(h)
i) Impulse response
>> impulse(h)