定点浮点fft.docVIP

浮点FFT#include math.h#define PI 3.14159void my_float_fft(float xr[], float xi[], short n){ int i,t; int N; // 总点数 int cnt; // 计数器 int rcnt; // 位反计数器 int step; // 位反计数器步长 float temp; // 交换用临时变量 int stage; // 变换级数 int group_start; // 组起始地址 int p_pos,q_pos ; // 组内蝶形P,Q的位置 float wr,wi; // 旋转因子实部、虚部 float qwr,qwi; // QW的实部、虚部 float pr,pi; // P的实部、虚部 float ph; // 旋转因子每次的旋转量 N=1n; // 得到总点数 rcnt = 0; // 位反计数器=0 // 倒序 for (cnt=1; cntN; cnt++) // 从1计到N { step = N; do // 搜索位反计数器的step { step = 1; } while (rcnt+step = N); rcnt = (rcnt (step-1)) + step; // 位反计数器+1 if (rcnt = cnt) continue; // 位反计数器小于等于计数器，则不交换 temp=xr[cnt]; // 交换 xr[cnt]=xr[rcnt]; xr[rcnt]=temp; temp=xi[cnt]; xi[cnt]=xi[rcnt]; xi[rcnt]=temp; } for(stage=1;stage=n;stage++) { t=1(stage-1); // 设置组内循环次数，第一级为1，每级增加一倍 ph=2*PI/N*(1(n-stage)); // 设置相位增量，第一级为PI，每级减小一半 for(group_start=0;group_startN;group_start+=(1stage)) { for(i=0;it;i++) { wr=cos(ph*i); // W的实部 wi=sin(ph*i); // W的虚部 p_pos=group_start+i; // 得到P的位置 q_pos=group_start+(1(stage-1))+i; // 得到Q的位置 qwr = wr*xr[q_pos] + wi*xi[q_pos]; // Q*W的实部 qwi = wr*xi[q_pos] - wi*xr[q_pos]; // Q*W的虚部 pr = xr[p_pos]; // P的实部 pi = xi[p_pos]; // P的虚部 xr[q_pos] = pr - qwr; // Q=P-Q*W的实部 xi[q_pos] = pi - qwi; // Q=P-Q*W的虚部 xr[p_pos] = pr + qwr; // P=P+Q*W的实部 xi[p_pos] = pi + qwi; // P=P+Q*W的虚部 } } } return;}#include tables.h /* inline short FIX_MPY(short a, short b){ long c = ((long)a * (long)b) 14; // (int16 * int16)-int32, int3214 a = (c 1) + (c 0x01); //判断最低位, 四舍五入 return a;}*/int_smpyr(int a, in