% Define the transfer function coefficients
numerator = 8.008;
denominator = [1 0.54 7.95]; % Coefficients of s^2, s, and the constant term
% Create the transfer function
sys = tf(numerator, denominator);
% Create the Bode plot
figure;
bode(sys);
grid on;
title('Bode plot of the system');
% Evaluate the frequency response at specific frequencies
f1 = 0.1; % frequency in Hz
f2 = 7.5; % frequency in Hz
w1 = 2 * pi * f1; % Convert to rad/s
w2 = 2 * pi * f2; % Convert to rad/s
response1 = evalfr(sys, 1j*w1);
response2 = evalfr(sys, 1j*w2);
% Display the magnitude and phase of the frequency responses
disp(['Frequency response at ' num2str(f1) ' Hz:']);
disp(['Magnitude: ' num2str(abs(response1)) ', Phase: ' num2str(angle(response1) * 180/pi) ' degrees']);
disp(['Frequency response at ' num2str(f2) ' Hz:']);
disp(['Magnitude: ' num2str(abs(response2)) ', Phase: ' num2str(angle(response2) * 180/pi) ' degrees']);
