function?snr=MY_FFT（w）
%PART_1
%-------------測試函數激勵-------------------------
%用正弦函數作為信噪比的測試函數，頻率為120HZ，采到整-
%頻點數，然后計數FFT結果中峰值與噪聲最大值的比值-----
%--------------------------------------------------
x=Sin_Gen（w）;
%PART_2
%-------------------------------------------------
%-------------蝶形運算模塊-------------------------
%-------------------------------------------------

%-------------------------------------------------
%第一級運算，沒有乘法器----------------------------
%[ABCD]=Step_One（abcd）----------------------
%輸入：第一級蝶形運算的數據，數據位寬為14-----------
%輸出：經過第一級蝶形運算產生的數據，數據位寬為16---
%-------------------------------------------------
for?m=1:256
????[ABCD]=Step_One（x（m）x（m+256）x（m+512）x（m+768））;
????X_1（m）=A;
????X_1（m+256）=B;
????X_1（m+512）=C;
????X_1（m+768）=D;
end
%-------------------------------------------------
%第二級運算
%[ABCD]=Step_Two（abcdm）--------------------
%輸入：第二級蝶形運算的數據，數據位寬為16-----------
%??????m為在本級中的循環級數-----------------------
%輸出：經過第二級蝶形運算產生的數據，數據位寬為18---
%注：1?該級為防溢出結果縮減為原來的1/2--------------
%????2?存儲蝶形因子采用了17位，實際擴大2^15---------
%-------------------------------------------------
for?m=1:4
????for?n=1:64
????????[ABCD]=Step_Two（X_1（256*（m-1）+n）X_1（256*（m-1）+n+64）X_1（256*（m-1）+n+128）X_1（256*（m-1）+n+192）m）;
????????X_2（256*（m-1）+n）=A;
????????X_2（256*（m-1）+n+64）=B;
????????X_2（256*（m-1）+n+128）=C;
????????X_2（256*（m-1）+n+192）=D;
????end
end
%-------------------------------------------------
%第三級運算
%[ABCD]=Step_Three（abcdm）------------------
%輸入：第三級蝶形運算的數據，數據位寬為18-----------
%??????m為在本級中的循環級數-----------------------
%輸出：經過第三級蝶形運算產生的數據，數據位寬為18---
%注：1?該級為防溢出結果縮減為原來的1/4--------------
%????2?存儲蝶形因子采用了18位，實際擴大2^16---------
%-------------------------------------------------
for?m=1:16
????for?n=1:16
????????[ABCD]=Step_Three（X_2（64*（m-1）+n）X_2（64*（m-1）+n+16）X_2（64*（m-1）+n+32）X_2（64*（m-1）+n+48）m-1）;
????????X_3（64*（m-1）+n）=A;
????????X_3（64*（m-1）+n+16）=B;
????????X_3（64*（m-1）+n+32）=C;
????????X_3（64*（m-1）+n+48）=D;
????end
end
%--------------------------------------------------
%第四級運算
%[ABCD]=Step_Four（abcdm）--------------------
%輸入：第四級蝶形運算的數據，數據位寬為18-----------
%??????m為在本級中的循環級數-----------------------
%輸出：經過第四級蝶形運算產生的數據，數據位寬為18---
%注：1?該級為防溢出結果縮減為原來的1/8--------------
%????2?存儲蝶形因子采用了18位，實際擴大2^16---------
%-------------------------------------------------
for?m=1:64
????for?n=1:4
????????[ABCD]=Step_Four（X_3（16*（m-1）+n）X_3（16*（m-1）+n+4）X_3（16*（m-1）+n+8）X_3（16*（m-1）+n+12）m-1）;
????????X_4（16*（m-1）+n）=A;
????????X_4（16*（m-1）+n+4）=B;
????????X_4（16*（m-1）+n+8）=C;
????????X_4（16*（m-1）+n+12）=D;
????end
end
%---------------------------------------------------------------
%第五級運算
%[ABCD]=Step_Five（abcdm）--------------------
%輸入：第五級蝶形運算的數據，數據位寬為18-----------
%??????m為在本級中的循環級數-----------------------
%輸出：經過第五級蝶形運算產生的數據，數據位寬為16---
%注：1?本級不需要防溢出--------------
%????2?存儲蝶形因子采用了18位，實際擴大2^16---------
%????3?數據輸出除以四是為了截斷到輸出需要的16位
%---------

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

?????文件???????4182??2010-08-20?13:36??FFT14_16\MY_FFT.m

?????文件??????????0??2010-08-30?17:01??FFT14_16\MY_FFT輸入的數據是頻率，自動生成該頻率的正弦函數.txt

?????文件????????161??2010-08-20?13:33??FFT14_16\Sin_Gen.m

?????文件??????29091??2010-08-30?17:06??FFT14_16\snr.jpg

?????文件??????22528??2010-08-30?19:28??FFT14_16\snr.xls

?????文件????????130??2010-08-20?13:47??FFT14_16\snr_cal.m

?????文件???????4395??2010-08-15?15:14??FFT14_16\Step_Five.m

?????文件???????5604??2010-08-20?13:05??FFT14_16\Step_Four.m

?????文件????????483??2010-08-15?10:23??FFT14_16\Step_One.m

?????文件???????5292??2010-08-15?15:12??FFT14_16\Step_Three.m

?????文件???????4085??2010-08-15?15:11??FFT14_16\Step_Two.m

?????目錄??????????0??2010-08-30?19:28??FFT14_16

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

????????????????75951????????????????????12