% Octave用   2005/06/20
%
% txtファイル(TAB区切り）のデータ（1列目：時間,2列目：データ）
% を読取り,周波数分析する.

clear
format long e

testdata =  load 030206c1_oct.txt;  %周波数分析するデータファイル

% pr0827

start_number = 1;
end_number = size(testdata,1);

time_testdata = testdata(1:end_number,1);

data1 = testdata(1:(end_number-start_number+1),3);

F = 1/100e-6;     %サンプリング周波数[Hz]
MAG_max=.0001;    %最大値(出力用)
MAG_min=0;        %最小値(出力用)  

frq_max=100;%F/2; %最大周波数(出力用)
frq_min=0;        %最小周波数(出力用)

y=fft(data1);

[Row,Col]=size(data1);

MAG=abs(y)/(Row/2);
% PHASE=angle(y);
f = (0:length(y)-1)'*F/length(y);

    axis([frq_min frq_max MAG_min MAG_max]);
    xlabel('Frequency[Hz]')
    ylabel('Rod side Pressure (abs) [MPa]')
    grid "on"
    plot(f,MAG)

