% Bode線図の[Hz]での表示
% 2006/06/02    MATLAB用
% 2006/06/09	f,wをf_bode,w_bodeに変更，bode入力部修正

f_bode = [1:0.1:100];
w_bode = 2*pi*f_bode;
NUM = [ 0 1];
DEN = [ 1/(2*pi*1) 1];

[MAG,PHASE]=Bode(NUM,DEN,w_bode);

figure(1)
subplot(2,1,1)
    semilogx(f_bode,20*log10(MAG))
    axis([min(f_bode) max(f_bode) min(20*log10(MAG)) max(20*log10(MAG))])
    title('')
    xlabel('Frequency[Hz]')
    ylabel('Gain [dB]')
    grid on
subplot(2,1,2)
    semilogx(f_bode,PHASE)
    axis([min(f_bode) max(f_bode) min(PHASE) max(PHASE)])
    xlabel('Frequency[Hz]')
    ylabel('Phase[deg]')
    grid on
    