Matlab:
N = 200; % 采样点数量
A = 3; % 振幅
fv = 10; % 波的重复次数
f = fv / 200; % 频率
t = [1:200]; % 时间
y = A * sin(t*2*pi*f); % 生成波形采样数据
plot(y); % 显示时域图
figure;
Y = fft(y, N); % FFT变换
plot(abs(Y)); % 显示频域图
--------------------------------------------
傅立叶变换入门:http://news.cnblogs.com/n/209014/
FFT是基于复数运算的,而实际采样点的数据一般是实数,因此要转换成复数形式。标准的做法是复数的实数部分为采样点的实数数据,虚数部分全部填0。(但这样做会有一定运算效率的问题,因为填0的部分也参与运算,对此也有一些优化算法。)