%　[Hz]単位のBode線図，実験結果と理論式の比較．
%　実験結果はsin_tr.mを用いてテキスト形式データとして保管されているものを用いる．
%　すなわち，周波数，Gain（dB)，位相(Phase）の順のデータで保管されているものを用いる．
%　2005.08.10 小川 Octave用．

clear;
subplot(2,1,1);
clearplot();
subplot(2,1,2);
clearplot();

%理論式の周波数特性求める．

num = [ 0 1];
den = [ 1 1 1];

SYS = tf2sys(num,den,0);

[Gain,Phase,w] = bode(SYS);

%実験結果の周波数特性求める．

loadfile='D:\\05kadai_03\\98Test\\05_09\\05_09_09ASS050_020RV_f_response\\Sin_low\\vpos2xfreq.csv';

exdata = dlmread(loadfile,',');

f_ex = exdata(:,1);
Gain_ex = exdata(:,2);
Phase_ex = exdata(:,3);

subplot(2,1,1);
%figure(1);
grid('on')
%xlabel('[Hz]')
ylabel('Gain[dB]')

%範囲指定
axis([ 0.01 10 -20 0],'auto');
%axis([ 0.01 10 -20 0]);

semilogx(w/(2*pi()),20*log10(Gain),'-;theory;',f_ex,Gain_ex,'-;experiment;');


subplot(2,1,2);
%figure(2);
xlabel('[Hz]')
ylabel('Phase[deg]')
grid('on')

%範囲指定
axis([ 0.01 10 -180 180],'auto');
%axis([ 0.01 10 -180 180]);

semilogx(w/(2*pi()),Phase,'-;theory;',f_ex,Phase_ex,'-;experiment;');
