filename = '普通歩き止まりなし.csv';%ファイル入れる

start = 47;%変換する範囲の最初
endpoint = 1466; %終わるポイント
axis = 2;%0がx1がｙ2がz軸の値を求める
M = csvread(filename,start,axis,[start,axis,endpoint,axis]);%[a,b,c,d]のa,cで時間指定、b,dで軸指定
Fs = 5;            % Sampling frequency                    
T = 1/Fs;             % Sampling period       
L = length(M);             % Length of signal
t = (0:L-1)*T;        % Time vector

Y = fft(M,L);
Y(1) = 0;

F = (0 : 1/L : 1/2-1/L)*Fs;


plot(F,abs(Y(1:L/2)))

title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('Frequency')
ylabel('magnitude')
