%　[Hz]単位のBode線図，実験結果の伝達関数の計算結果を図示．
%　実験結果はCSV形式データとして
%　周波数，Gain（dB)，位相(Phase）の順のデータで保管されているものを用いる．
%　2005.09.29 小川 Octave用．

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

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

loadfile='D:\\05kadai_03\\98Test\\as320_002freq_sn002_250HZ.csv';
exdata = dlmread(loadfile,',');

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

%位相を0～-360[deg]にあわせる．
%for ii = 1:size(Phase_ex,2)
%	if Phase_ex(1,ii) >= 0
%		Phase_ex(1,ii) = Phase_ex(1,ii)-360;
%	end
%end

gset('lmargin 13')

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

%範囲指定
%axis([ 1 100 -70 20],'auto');
axis([10 1000 -20 5]);

semilogx(f_ex,Gain_ex,'1o;;');

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

%範囲指定
%axis([ 1 100 -300 100],'auto');
axis([ 10 1000 -360 0]);

semilogx(f_ex,Phase_ex,'1o;;');
