main_fft_2_4.c

来自「用汇编实现的FFT 程序。基于ADI的DSP 实时处理程序」· C语言 代码 · 共 56 行

#include<filter.h>
#include<fract.h>
#include<complex.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;

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++)
    {
        for(i=0;i<L;i++)
        {
			x[i]=In1[k*L+i];
        }
 	   rfft_fr16(x,temp,x_fft,twiddle_table,twiddle_stride,N,block_exponent,scale_method);
 	   for(i=0;i<N;i++)
 	   {
	       yr_fft[i] = cmlt_fr16(x_fft[i],h_fft[i]);
 	       yr_fft[i].re <<=8;
 	       yr_fft[i].im <<=8;
 	   } 	   
 	   ifft_fr16(yr_fft,temp,yr,twiddle_table,twiddle_stride,N,block_exponent,scale_method);
 	   for(i=0;i<N;i++)
 	   {
 	       y[k*L+i] += yr[i].re<<8;
 	   }
    }
    while(1)    
    {
    }
}