FFT汇编程序代码.docVIP

FFT汇编程序代码 FFT算法汇编原程序与大家共享 ;----------------------------------------------------------- ; 快速付里叶变换子程序 ; ; 入口 : 一维四字节浮点数数组的首地址 ; ; ; 以 2 为基β的值,信号抽样的点数 = 2^β; ; ; 出口 : 在原数组的位置保存结果的值; ; ; 补充说明 : ; ; 1. 该子程序所需 RAM 的容量为 2^β*12 字节,另外需要 ; ; 少量的堆栈作为临时变量的存储空间; ; ; 2. 所需 RAM 空间以输入的首地址为基址,向增加的方向 ; ; 扩展; ; ; 应用举例 : ; ; PUSH #0A000H ; 数组首地址压栈 ; PUSH #6 ; β值压栈 ; CALL FFT ; ;----------------------------------------------------------- PROC FFT FFT: LD FFT_CX,2[SP] LD FFT_AX,#0001H ; SHL FFT_AX,FFT_CL ; 计算采样点数 PUSH FFT_AX ; 将采样点数压栈 SHL FFT_AX,#1 ; LD FFT_BX,6[SP] ; 根据采样点数, ADD FFT_BX,FFT_AX ; 计算出其余 PUSH FFT_BX ; 三个被占用 ADD FFT_BX,FFT_AX ; 空间的首地址; PUSH FFT_BX ; 并依次将首地址 SHL FFT_AX,#1 ; ADD FFT_BX,FFT_AX ; 压栈; PUSH FFT_BX ; LD FFT_CX,6[SP] LD FFT_FX,#2 LOOP1: CLR FFT_AX SUB FFT_EX,FFT_CX,#1 LD FFT_BX,FFT_EX LD FFT_DX,10[SP] LOOP2: SHL FFT_AX,#1 ; SHR FFT_EX,#1 ; 产生倒序数 ADDC FFT_AX,0 ; DJNZW FFT_DX,LOOP2 ; SHL FFT_AX,#2 SHL FFT_BX,#2 ; FFT_AX : 产生源地址 ADD FFT_AX,12[SP] ; FFT_BX : 产生目的地址 ADD FFT_BX,2[SP] ; 输入数据按倒序排列 BMOV FFT_AX,FFT_FX ; DJNZW FFT_CX,LOOP1 ; LD FFT_AX,12[SP] LD FFT_BX,4[SP] LD FFT_DX,#0080H ST FFT_DX,[FFT_AX]+ ST 0,[FFT_AX] ST 0,[FFT_BX]+ ST 0,[FFT_BX] LD FFT_AX,#0DB82H ; 计算分度值 SUB FFT_AX,10[SP] ; 即,(2*π)/(2^β) PUSH #490FH ; 的值 PUSH FFT_AX ;