基于MATLAB的快速傅里叶变换

一个基于MATLAB的快速傅里叶变换，给新手练练手，真是很不错

load woke.txt;
l_s=length（woke）;
n=ceil（l_s/2）;
for i=1:l_s
    t1（i）=woke（i1）;
    f1（i）=woke（i2）;
end
%自己设置采样频率（这里认为500）
Fs=500;



subplot（211）;
plot（t1f1）;
ylabel（‘原始图‘）;




ps=fft（f1l_s）;
%!取幅值大小
%Re = real（ps）;
%Im = imag（ps）;

%for （ i= 1:1:2048 ）
%	Re2（i） = Re（i） * Re（i）;
%	Im2（i） = Im（i） * Im（i）;
%	A（i） = 8*sqrt（ Re2（i） + Im2（i） ） /2048;
%end
A=2*abs（ps）/l_s;
A（1）=A（1）/2;

%!画出频域图
tt=Fs*（0:n-1）/l_s;
subplot（212）;plot（ttA（1:n））;

ylabel（‘频谱图‘）;
%可以将tt和A的数据保存，他们就是频谱分析的结果，分别是频率和幅值

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        602  2005-12-04 16:51  fft.m

----------- ---------  ---------- -----  ----

                  820                    2