function?y=FFT（xorder）;
%disp（‘2的階數：‘）;
%order=input（‘‘）;??????????????????????????%獲取蝶形運算的級數??
num=power（2order）;
%t=0:0.001:（（num-1）*0.001）;
%x=sin（2*pi*50*t）;????????????????????????%正弦函數離散化

n=0:num-1;?????????????????????????????%碼位倒置算法
for?h=1:length（n）
????temp2=dec2bin（n（h）order）;
????c=uint8（temp2）-48;?????????????????%char型轉二進制
????temp1=c;
???for?j=order:-1:1
???????temp1（j）?=?c（order-j+1）;??????
???end
???c=temp1+48;
???b（h）=bin2dec（char（c））+1;?????????????%二進制轉十進制?????????????
end
a=b;
for?m=1:num???????????????????????????????????????????????????????????
????y（m）=x（a（m））;
end
x=y;

w0=exp（-i*2*pi/num）;??????????%FFT算法???????????????????????????????????
pow=1;
count=num;
for?m=1:order;?????????????????????????????????????????
????pow=pow*2;
????count=count/2;
????e=count;
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????65024??2010-01-27?22:08??FIR.doc

?????文件????????707??2010-01-27?22:07??FIR.m

?????文件???????1335??2009-12-31?16:00??FFT.m

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

????????????????67066????????????????????3