【正文】 我的關心和幫助。附錄I附錄II/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////主函數：從主函數中可以看出FFT的基本過程和步驟，以及碼位倒置、產生三角因子、蝶形運算、計算功率譜。 .title ″zxy_fft5k_program″ .mmregs .global_c_int00 .include ″″ .include ″″ .include ″″STACK .usect ″stack″，stack_sizeSYSTEM_STACK .set STACK+stack_size .sect ″vectors″reset: BD _c_int00 STM SYSTEM_STACK，SP .space 31*4*16 .text_c_int00: NOP CALL c54x_init 。Initialize CPU and clear A and B CALL bit_rev 。Here AR7→fft_data，others is free NOP CALL fft_task NOP CALL power NOPmain_start: B main_start .end///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////初始化CPUstack_size .set 200hK_FFT_SIZE .set 256 .textc54x_init: 。CPU initialize BC TstStatu，ANEQ STM 9FF7h，CLKMD 。CPU Clock Feg.:100MHz STM 00A0h，PMST STM 8FFFh，SWWSR 。IO wait for 7 clk STM 8802h，BSCR 。Clockout=50MHz LD 0h，DP 。Data Page Points at LD 0h，A LD 0h，B RET////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////，其中，，，，分別為 =cos()jsin()=cos()j0=1j0 =cos()jsin()=cos()jsin()=0j1 =cos()jsin()=cos()jsin()= =cos()jsin()=cos()jsin()= =cos()jsin()=cos()jsin()= =cos()jsin()=cos()jsin()= =cos()jsin()=cos()jsin()= =cos()jsin()=cos()jsin()=把上述系數的實部和虛部分別乘以32768（即）得到如下數值：32768，0，0，32768，23170，23171，30273，12540，32128，6393，32610，3212，32728，1608，32758，805。 .data d_input_circ 。輸入數據.word 0，10125，19260，26509，31164，32767，31164，26509，19260，10125 .word 0，10126，19261，26510，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 0，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 0，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 0，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 0，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 1，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 0，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 0，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 1，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 1，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 1，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 0，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 1，10126，19261，25610，31165，32768，31165，25610，19261，10126.word 1，10125，19260，26509，31164，32767，31164，26509，19260，10125.word 0，10126，19261，25610，31165，32768temp1 .space 7*256*16 。為了碼位倒置結果正確而設的一段空間d_input_addr 。存放插零值后的數據，256個數隔一個數差一個零值 space 512*16fft_data 。存放碼位倒置后的數據結果 .space 512*16weight_fact 。系數，按實部、虛部、實部、虛部排列 .word 32768，0，0，32768，23170，23171 。，的值 .word 30273，12540，32138，6393，32610，3212 。，的值 .word 32728，1608，32758，805 。，的值L_t .word 0hLE_v .word 2，4，8，16，32，64，128，256LE1_v .word 1，2，4，8，16，32，64，128U_v .word 0h，0h 。存放的值T_v .word 0h，0hJ_v .word 0hI_v .word 0hfft_data_tst 。存放FFT后的值 .space 256*16Temp .word 0h//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// .textbit_rev: 。碼位倒置 　NOP 　 RSBX CMPT 　 SSBX FRCT 。Fractional mode is on 　 SSBX SXM 　 LD 1，ASM 　 STM d_input_addr，AR3 　STM d_input_circ，AR2 。modify new data 　STM K_FFT_SIZE1，BRC 　 RPTBD data_move1 　 STM 256，BK 　 MVDD *AR2，*AR3+ 　 MAR *AR2+% 　 ST 0h，*AR3+ 。Real input data 256→plex data 512data_move: 　 STM d_input_addr，AR3 。AR3→1 st original input 　 STM fft_data，AR7 。AR7→data processing buffer 　 MVMM AR7，AR2 。AR2→1 st bitreversed data 　 STM K_FFT_SIZE1，BRC 　 RPTBD bit_rev_end1 　 STM K_FFT_SIZE，AR0 。AR0=1/2 size of circ buffer 　MVDD *AR3+，*AR2+ 　 MVDD *AR3，*AR2+ 　 MAR *AR3+0Bbit_rev_end: 　NOP 　 RET 。AR7 is reserved for fft_datapower:　　　　　　　　　　　　　　　　　　 。計算功率譜 STM fft_data，AR2 STM fft_data，AR3 MAR