基于dsp数字信号处理器的fft实现

基于dsp数字信号处理器的fft实现内容摘要：

F MVMM DATA_PROC_BUF,REORDERED STM K_FFT_SIZE1,BRC RPTBD bit_rev_end1 STM K_FFT_SIZE,AR0 MVDD *ORIGINAL_INPUT+,*REORDERED+ MVDD *ORIGINAL_INPUT,*REORDERED+ MAR *ORIGINAL_INPUT+0B bit_rev_end: .asg AR1,GROUP_COUNTER .asg AR2,PX .asg AR3,QX .asg AR4,WR .asg AR5,WI .asg AR6,BUTTERFLY_COUNTER .asg AR7,STAGE_COUNTER STM 0,BK LD 1,ASM STM fft_data,PX STM fft_data+K_DATA_IDX_1,QX STM K_FFT_SIZE/21,BRC LD *PX,16,A RPTBD stage1end1 STM K_DATA_IDX_1+1,AR0 SUB *QX,16,A,B ADD *QX,16,A STH A,ASM,*PX+ ST B,*QX+ ||LD *PX,A SUB *QX,16,A,B ADD *QX,16,A STH A,ASM,*PX+0% ST B,*QX+0% ||LD *PX,A stage1end: STM fft_data,PX STM fft_data+K_DATA_IDX_2,QX STM K_FFT_SIZE/41,BRC LD *PX,16,A RPTBD stage2end1 STM K_DATA_IDX_2+1,AR0 SUB *QX,16,A,B ADD *QX,16,A STH A,ASM,*PX+ ST B,*QX+ ||LD *PX,A SUB *QX,16,A,B ADD *QX,16,A STH A,ASM,*PX+ STH B,ASM,*QX+ MAR *QX+ ADD *PX,*QX,A SUB *PX,*QX,B STH A,ASM,*PX+ SUB *PX,*QX,A ST B,*QX ||LD *QX+,B ST A,*PX ||ADD *PX+0%,A ST A,*QX+0% ||LD *PX,A stage2end: STM K_TWID_TBL_SIZE,BK ST K_TWID_IDX_3,d_twid_idx STM K_TWID_IDX_3,AR0 STM cosine,WR STM sine,WI STM K_LOGN21,STAGE_COUNTER ST K_FFT_SIZE/81,d_grps_t STM K_FLY_COUNT_31,BUTTERFLY_COUNTER ST K_DATA_IDX_3,d_data_idx stage: STM fft_data,PX LD d_data_idx,A ADD *(PX),A STLM A,QX MVDK d_grps_t,GROUP_COUNTER group: MVMD BUTTERFLY_COUNTER,BRC RPTBD butterflyend1 LD *WR,T MPY *QX+,A MAC *WI+0%,*QX,A ADD PX,16,A,B。 B:=(QR*WR+QI*WI)+PR ST B,*PX。 PR39。 :=((QR*WR+QI*WI)+PR)/2 ||SUB *PX+,B ST B,*QX ||MPY *QX+,A MAS *QX,*WR+0%,A ADD *PX,16,A,B ST B,*QX+ ||SUB *PX,B LD *WR,T ST B,*PX+ ||MPY *QX+,A butterflyend: PSHM AR0 MVDK d_data_idx,AR0 MAR *PX+0 MAR *QX+0 BANZD group,*GROUP_COUNTER POPM AR0 MAR *QX LD d_data_idx,A SUB 1,A,B STLM B,BUTTERFLY_COUNTER STL A,1,d_data_idx LD d_grps_t,A。