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
%自己設(shè)置采樣頻率（這里認(rèn)為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;

%!畫(huà)出頻域圖
tt=Fs*（0:n-1）/l_s;
subplot（212）;plot（ttA（1:n））;

ylabel（‘頻譜圖‘）;
%可以將tt和A的數(shù)據(jù)保存，他們就是頻譜分析的結(jié)果，分別是頻率和幅值

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

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

-----------?---------??----------?-----??----

??????????????????820????????????????????2