離散傅裡葉變換matlab程式

離散傅裡葉變化可以将一個序列展開成不同頻率的正弦信号的和，matlab中采用fft函數即可實作。x信号為2000Hz和3000Hz的正弦信号之和。

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

%%離散時間傅裡葉變換-正弦信号

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

clc;clear all;close all;

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

f1 = 2000;

f2 = 3000;

Fs = 10000;%采樣頻率為10KHz

dt = 1/Fs;%采樣時間

L = 1000;%采樣點數

t = (0:L-1)*dt;%時間向量

x = cos(2*pi*f1*t)+cos(2*pi*f2*t);

w = (0:L-1)*2*pi/L-pi;

F = fft(x)/L;

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

stem(W,abs(F(1:L/2)));%stem用于繪制莖狀圖