% 保存してあった周波数応答結果を重ねがきでBode線図（[Hz]）
% 2006/06/04    MATLAB用

f = [1 200];
MAG = [-20 20];
PHASE = [-180 180];

w = 2*pi*f;

f_res1 = csvread('D:\06kadai_01\00_simulation\アクティブダンサ\060603_01frq_FrefToFact_sim.csv');

f_res2 = csvread('D:\06kadai_01\98Test\06_06\06_06_03アクティブダンサ試験装置300Nシリンダ単体F特\060603_02f_res_ActFrcRefToActFrcAct.csv');

f1 = f_res1(:,1);       %[Hz]
MAG1 = f_res1(:,2);     %[dB]
PHASE1 = f_res1(:,3);   %[deg]

f2 = f_res2(:,1);       %[Hz]
MAG2 = f_res2(:,2);     %[dB]
PHASE2 = f_res2(:,3);   %[deg]

figure(1)
subplot(2,1,1)

    semilogx(f1,MAG1)
    hold on
    semilogx(f2,MAG2,'ro')
    hold off
%    axis([min(f) max(f) min([MAG1;MAG2]) max([MAG1;MAG2])])
    axis([min(f) max(f) min(MAG) max(MAG)])
    title('')
    xlabel('Frequency[Hz]')
    ylabel('Gain [dB]')
    grid on
    
    
subplot(2,1,2)
    semilogx(f1,PHASE1)
    hold on
    semilogx(f2,PHASE2,'ro')
    hold off
%    axis([min(f) max(f) min([PHASE1;PHASE2]) max([PHASE1;PHASE2])])
    axis([min(f) max(f) min(PHASE) max(PHASE)])
    xlabel('Frequency[Hz]')
    ylabel('Phase[Deg]')
    grid on
    