% xfft.m

%

% ＜　Fourier Tansform ＞

%

Fs = 200; % サンプリング周波数 (Hz)

T = 1; % 観測時間（1sec）

dt = 1/Fs; % 時間分解能（サンプリング間隔）

N = T/dt; % サンプル数

t = (0:dt:T-dt); % time vector

%

% ----------------------signal---------------------------

x = cos(2*pi*10*t) + sin(2*pi*50*t); % 信号（時系列データ）

figure(1)

subplot(411); plot(x); legend('signal'); % 観測信号のプロット

%

% ----------------fast fourier tansform------------------

w = fft(x)/N;

subplot(412); plot(real(w)); % フーリエ係数（実部）

legend('fourier coefficient (real)');

subplot(413); plot(imag(w)); % フーリエ係数（虚部）

legend('fourier coefficient (imagenary)');

subplot(414); plot(abs(w)); % フーリエ係数（絶対値）

legend('fourier coefficient (absolute value)');

%

% ------------------Perseval theorem---------------------

pwrt = sum(x.^2)/N

pwrf = sum(abs(w).^2)

% -------------------------------------------------------

% 2002.4.29 by K.Tsukada