echo off;
clear
clc
close all
disp('[step 1]から[step 5]において読み取った振幅比と位相差の値をボード線図上にプロットします．')
disp('[step 1]から[step 5]で読み取った振幅比[dB]と位相差[deg]の値を入力してください');

amp1=input('[step 1]の振幅比=');
phase1=input('[step 1]の位相差=');
amp2=input('[step 2]の振幅比=');
phase2=input('[step 2]の位相差=');
amp3=input('[step 3]の振幅比=');
phase3=input('[step 3]の位相差=');
amp4=input('[step 4]の振幅比=');
phase4=input('[step 4]の位相差=');
amp5=input('[step 5]の振幅比=');
phase5=input('[step 5]の位相差=');

echo on;
% --------------------------------------------------------
% 読み取った振幅比と位相差をボード線図上にプロットします．
% --------------------------------------------------------

pause	% press any key to continue.

figure(1)
semilogx(0.01,amp1,'o')
set(gca,'Ylim',[-40 20]);
xlabel('\omega [rad/s]','fontsize',12)
ylabel('Gain [dB]','fontsize',12)
title('ゲイン線図','fontsize',12)
grid on
hold on
semilogx(0.1,amp2,'o')
semilogx(1.0,amp3,'o')
semilogx(10,amp4,'o')
semilogx(100,amp5,'o')

pause	% press any key to continue.

figure(2)
semilogx(0.01,phase1,'o')
set(gca,'Ylim',[-180 90]);
xlabel('\omega [rad/s]','fontsize',12)
ylabel('Phase [deg]','fontsize',12)
title('位相線図','fontsize',12)
grid on
hold on
semilogx(0.1,phase2,'o')
semilogx(1,phase3,'o')
semilogx(10,phase4,'o')
semilogx(100,phase5,'o')

pause	% press any key to continue.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 読み取った振幅比と位相差の値がボード線図上にプロットされました．
% 読み取った値(点)が，振幅比と位相差ともに右下がりになっていること
% を確認してください．
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

pause	% press any key to continue.

% --------------------------------------------------------
% ボード線図のシミュレーション応答を重ねてプロットし，
% 読み取った値と比較します．
% --------------------------------------------------------

pause	% press any key to continue.

figure(1)
w=logspace(-2,2,100);
[g,p]=bode([1],[1 1],w);
semilogx(w,20*log10(g));
hold on

pause	% press any key to continue.

figure(2)
semilogx(w,p);

pause	% press any key to continue.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 読み取った値(点)が，シミュレーション応答による曲線に沿って
% いるかどうかを確認してください．
%
% 以上の結果から，1次系では，周波数が低い範囲では，ゲインはほぼ 0 
% [dB] となり，周波数が高い範囲では -20[dB/dec] で減少すること
% を確認してください．また，1次系では，位相が最大で -90[deg] 
% 遅れることも確認してください． 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

pause	% press any key to continue.

