main_fft_2_4_asm.c

用汇编实现的FFT 程序。基于ADI的DSP 实时处理程序

#include<filter.h>
#include<fract.h>
#include<complex.h>
#include<math.h>

#define Len_In 1182
#define P 119
#define N 256
#define L 138
#define t 9

fract16 h[N] = {
    #include"FIR_LP_coeff.dat"
}; 
fract16 In1[t*L] = {
    #include"Indata.dat"
};
fract16 y[t*L+P-1];  //output

//调用库函数rfft_fr16和ifft_fr16所需参数的声明
complex_fract16 temp[N]; 
complex_fract16 twiddle_table[N];
int twiddle_stride=1; 
int block_exponent=0; 
int scale_method =0;

//circular 
extern void fftasm1(int ,fract16 In[],fract16 x[]);
extern void fftasm2(int ,complex_fract16 x_fft[],complex_fract16 h_fft[],complex_fract16 yr_fft[]);
extern void fftasm3(int ,complex_fract16 yr[], fract16 y[]);

void main(void)
{
	complex_fract16 h_fft[N],x_fft[N],yr_fft[N],yr[N];
	fract16 x[N];
	
    twidfftrad2_fr16(twiddle_table,N);
    rfft_fr16(h,temp,h_fft,twiddle_table,twiddle_stride,N,block_exponent,scale_method);
    int k=0,i=0;
    for(k=0;k<t;k++)
    {
       fftasm1(L,&In1[k*L],x);
 	   rfft_fr16(x,temp,x_fft,twiddle_table,twiddle_stride,N,block_exponent,scale_method);
 	   fftasm2(N,x_fft,h_fft,yr_fft); 
 	   ifft_fr16(yr_fft,temp,yr,twiddle_table,twiddle_stride,N,block_exponent,scale_method);
 	   fftasm3(N,yr,&y[k*L]); 	   
    }
    while(1)    
    {
    }
}