CEP
CEP
System (a)
The general transfer function is:
C(s) = b / s(s 2 + as + b)
C(s) = b / s(s 2 + as + b)
= A / s +( Bs + C / s 2 + as + b )
b = A(s 2 + as + b) + (Bs + C)s
Equate coefficients:
s 2 : A + B = 0 =⇒ B = −A
S : Aa + C = 0 =⇒ C = −Aa
Constant: Ab = b =⇒ A = 1
A = 1, B = −1, C = −a
C(s) = 1 / s − (s + a / s 2 + as + b )
Complete the square in the denominator:
s 2 + as + b = ( s + a / 2 )2 + ( b − a 2 / 4 )
Assuming (b − a 2 / 4 > 0) (underdamped case), we rewrite:
C(s) = 1 / s − s + a / (( s + a / 2 )2 + ( b − a 2 / 4 ))
Split the second term:
s + a / (( s + a / 2 )2 + ω2 )
= ( s + a / 2 )/ (( s + a / 2 )2 + ω2 ) + (a − a / 2) / (( s + a / 2 )2 + ω2 )
= ( s + a / 2 )/ (( s + a / 2 )2 + ω2 ) + ( a / 2) / (( s + a / 2 )2 + ω2 )
c(t) = 1 − e − (a /2) t cos(ωt) − a/ 2ω e − (a/2) t sin(ωt), ω = √ b − a 2 / 4
System (b):
C(s) = 9 / s(s 2 + 9 s+9)
Rewrite the denominator:
s 2 + 9s + 9 = (s + √ 9 − √ 6)(s + √ 9 + √ 6)
= (s + 1.468)(s + 7.532)
C(s) = 9 / s(s + 1.468)(s + 7.532)
= A / s + (B / s + 1.468) + (C / s + 7.532)
9 = A(s + 1.468)(s + 7.532) + Bs(s + 7.532) + Cs(s + 1.468)
Evaluate at roots:
At s = 0: 9 = A(1.468)(7.532) =
A = 9 / 1468 × 7.532 = 0.814
At s = −1.468: 9 = B(−1.468 + 7.532) =
B = 9/ 6.064 = 1.484
At s = −7.532: 9 = C(−7.532 + 1.468) =
C = 9 / -6.604 = -1.484
C(s) = 0.814 / s +( 1.484 / (s + 1.468) )− (1.484 / (s + 7.532 ))
c(t) = 0.814 + 1.484 e −1.468t − 1.484 e −7.532t
System (c):
Solve:
9 = A(s 2 + 9) + (Bs + C)s
s2:A=1
s: B = 0
Constant: 9A = 9
A=1
Thus:
C(s) = 1 / s − (s / s 2 + 9)
c(t) = 1 − cos(3t)
System (e):
C(s) = 9 / s(s 2+6s+9)
s 2 + 6s + 9 = (s + 3)2
C(s) = 9 / s(s + 3)2
= A / s + B / s + 3 + C / (s + 3)2
Solve:
9 = A(s + 3)2 + Bs(s + 3) + Cs
s2:A+B=0=
B = −A
s: 6A + 3B + C = 0
Constant: 9A = 9 =
A=1
Thus: A = 1, B = −1, C = −3
C(s) = 1 / s − (1 /s + 3) − (3 / (s + 3)2 )
c(t) = 1 − e −3t − 3te−3t
System B
Matlab code for step response..
% Define time vector from 0 to 5 seconds
t = 0:0.01:5;
% Output in Laplace
C_s = G * R;
% Time vector
t_vals = 0:0.01:5;
c_vals = c_func(t_vals);
% Plotting
figure;
plot(t_vals, c_vals, 'LineWidth', 2);
grid on;
xlabel('Time (s)');
ylabel('c(t)');
title('Step Response of System (c) with Time-Domain Expression');
% Display equation as text on plot
eqn_str = ['c(t) = ', char(vpa(c_t, 3))];
text(0.5, 0.8, eqn_str, 'Units', 'normalized', 'FontSize', 10, 'BackgroundColor', 'white');
System D
Poles and zeros :
% Define transfer function variable
s = tf('s');
Response :
syms s t
% Plotting
figure;
plot(t_vals, c_vals, 'LineWidth', 2);
grid on;
xlabel('Time (s)');
ylabel('c(t)');
title('Step Response of System (d) with Time-Domain Expression');
Response:
syms s t
Q2:
Part a:
% Given second-order transfer function parameters
wn = sqrt(361); % Natural frequency
zeta = 16 / (2 * wn); % Damping ratio
% Calculations
Ts = 4 / (zeta * wn); % Settling time (2% criterion)
Tp = pi / (wn * sqrt(1 - zeta^2)); % Peak time
PO = exp(-zeta * pi / sqrt(1 - zeta^2)) * 100; % Percent Overshoot
Tr = (1.76 * zeta^3 - 0.417 * zeta^2 + 1.039 * zeta + 1) / wn; % Rise time approx
% Display results
fprintf('Natural Frequency (wn): %.2f rad/s\n', wn);
fprintf('Damping Ratio (zeta): %.4f\n', zeta);
fprintf('Settling Time (2%%): %.4f seconds\n', Ts);
fprintf('Peak Time: %.4f seconds\n', Tp);
fprintf('Percent Overshoot: %.2f%%\n', PO);
fprintf('Rise Time: %.4f seconds\n', Tr);
Part b :
% Symbolic solution for step response
syms s t