离散傅里叶变换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用于绘制茎状图