s10fig9c.m

Industrial Mathematics

% fig 10.9a Bode Plots
%
fmin = 0;           % min freq (in Hz) to view
fmax = .5;           % max freq (in Hz) to view
df=0.01;             % step size
f=fmin:df:fmax;     % frequency vector
w = 2*pi*f;         % radian freq vector
s = w*i;            %  
N= s.*s  + 1;       % numerator of transfer fn
D = s.*s  + s + 1;  % denominator
H= N./D;            % transfer fn
r = abs(H);         % gain at freq w
phi = angle(H);     % phase lead
% plot(f,r,'k-')    % gain vs freq in Hz
%plot(f,phi*180/pi,'k-')    %for phase plot in degrees
 plot(f,20*log10(r),'k-')     %  for decibels 
% replace with   plot(w,r) for       gain vs radian freq