% 此範例展示一個簡單正弦波的傅立葉轉換，以雙邊頻譜來顯示
% 此正弦波的頻率恰巧是 freqStep 的整數倍，所以雙邊頻譜應該只有兩個非零點

N = 256;                        % 點數
fs = 8000;                      % 取樣頻率
freqStep = fs/N;              % 頻域的頻率的解析度
f = 10*freqStep;                % 正弦波的頻率，恰是 freqStep 的整數倍
time = (0:N-1)/fs;                % 時域的時間刻度
y = cos(2*pi*f*time);           % Signal to analyze
Y = fft(y);                     % Spectrum
Y = fftshift(Y);                % 將頻率軸的零點置中

% Plot time data
subplot(3,1,1);
plot(time, y, '.-');
title('Sinusoidal signals');
xlabel('Time (seconds)'); ylabel('Amplitude');
axis tight

% Plot spectral magnitude
freq = freqStep*(-N/2:N/2-1);   % 頻域的頻率刻度
subplot(3,1,2);
plot(freq, abs(Y), '.-b'); grid on
xlabel('Frequency)'); 
ylabel('Magnitude (Linear)');

% Plot phase
subplot(3,1,3);
plot(freq, angle(Y), '.-b'); grid on
xlabel('Frequency)'); 
ylabel('Phase (Radian)');