% Program for Bode Plot
clc;
clear;
num=[50];
den=[1 9 30 40];
sys=tf(num,den);
w=logspace(-1,5,1000); % Specify Frequency range
bode(sys,w);
margin(sys);
grid on;
xlabel('Frequency')
ylabel('Phase angle in degrees')
title('Bode Plot of system (50/s^3+9s^2+30s+40)')
[Gm, Pm, Wcp, Wcg]=margin(sys) % Finding Gain and Phase Margins and
%Gain & Phase Crossover Frequencies
OUTPUT MATLAB:
OUTPUT OCTAVE: