彥霖 實驗筆記

由此程式碼驗證得知 "快速傅立葉轉換" 再進行 "快速傅立葉反轉換" 後的數據 (音波) 將保持不變。

y = wavread('E:\2.音樂\onj001.wav');

s = fft(y); % 快速傅立葉轉換
f = abs(s); % 取絕對值

y2 = ifft(s); % 快速傅立葉反轉換
f2 = abs(fft(y2)); % 快速傅立葉反轉換 後的結果 再進行一次快速傅立葉轉換

subplot(2,2,1), plot(f); % 原檔案的頻譜
subplot(2,2,2), plot(y); % 原檔案的音波
subplot(2,2,3), plot(f2); % 快速傅立葉反轉換後的頻譜
subplot(2,2,4), plot(y2); % 快速傅立葉反轉換後的音波

sound(y2, 44100); % 撥放快速傅立葉反轉換的結果