移植codec2的fft出现与电脑运行结果有区别

科技骏马

我把codec2移植到STM32F409+CubeMX+freeRTOS+fatfs，开发板不是原子的，是杂牌的，从SD卡上读取raw语音，用450bps编码后再写到SD卡，在电脑上解码后听不清，用电脑编码的同一段声音再用电脑解码还能听清。我调试中发现最初与电脑运行过程有分歧的是make_analysis_window函数中codec2_fft的结果W数组，可能是arm_cfft_f32与kiss_fft有差异造成的，但是执行analyse_one_frame中dft_speech的codec2_fft_inplace里面也有arm_cfft_f32，结果都与电脑执行结果相同。我贴一下核心代码，在sine.c中

1. void make_analysis_window(C2CONST *c2const, codec2_fft_cfg fft_fwd_cfg, float w[], COMP W[])
   
2. {
   
3.   float m;
   
4.   COMP  wshift[FFT_ENC];
   
5.   COMP  temp;
   
6.   int   i,j;
   
7.   int   m_pitch = c2const->m_pitch;
   
8.   int   nw      = c2const->nw;
   
9.   //这里删除了一些注释
   
10. 
    
11.   m = 0.0;
    
12.   for(i=0; i<m_pitch/2-nw/2; i++)
    
13.   w[i] = 0.0;
    
14.   for(i=m_pitch/2-nw/2,j=0; i<m_pitch/2+nw/2; i++,j++) {
    
15.   w[i] = 0.5 - 0.5*cosf(TWO_PI*j/(nw-1));
    
16.   m += w[i]*w[i];
    
17.   }
    
18.   for(i=m_pitch/2+nw/2; i<m_pitch; i++)
    
19.   w[i] = 0.0;
    
20. 
    
21.   /* Normalise - makes freq domain amplitude estimation straight
    
22.      forward */
    
23. 
    
24.   m = 1.0/sqrtf(m*FFT_ENC);
    
25.   for(i=0; i<m_pitch; i++)
    
26.   w[i] *= m;
    
27.   //这里删除了一些注释
    
28.   for(i=0; i<FFT_ENC; i++) {
    
29.   wshift[i].real = 0.0;
    
30.   wshift[i].imag = 0.0;
    
31.   }
    
32.   for(i=0; i<nw/2; i++)
    
33.   wshift[i].real = w[i+m_pitch/2];
    
34.   for(i=FFT_ENC-nw/2,j=m_pitch/2-nw/2; i<FFT_ENC; i++,j++)
    
35.   wshift[i].real = w[j];
    
36.   codec2_fft(fft_fwd_cfg, wshift, W);//这里开始进入codec2_fft.h的内联函数，这里出错了
    
37.   //这里删除了一些注释
    
38.   for(i=0; i<FFT_ENC/2; i++) {
    
39.   temp.real = W[i].real;
    
40.   temp.imag = W[i].imag;
    
41.   W[i].real = W[i+FFT_ENC/2].real;
    
42.   W[i].imag = W[i+FFT_ENC/2].imag;
    
43.   W[i+FFT_ENC/2].real = temp.real;
    
44.   W[i+FFT_ENC/2].imag = temp.imag;
    
45. }
    
46. }

复制代码

在codec2_fft.h中，USE_KISS_FFT未定义，执行宏定义else中的语句

1. static inline void codec2_fft(codec2_fft_cfg cfg, codec2_fft_cpx* in, codec2_fft_cpx* out)
   
2. {
   
3. #ifdef USE_KISS_FFT
   
4.   kiss_fft(cfg, (kiss_fft_cpx*)in, (kiss_fft_cpx*)out);
   
5. #else
   
6.   memcpy(out,in,cfg->instance->fftLen*2*sizeof(float));
   
7.   arm_cfft_f32(cfg->instance,(float*)out,cfg->inverse, 1);
   
8.   //这里删除了一些注释
   
9.   if (cfg->inverse)
   
10.   {
    
11.     arm_scale_f32((float*)out,cfg->instance->fftLen,(float*)out,cfg->instance->fftLen*2);
    
12.   }
    
13. #endif
    
14. }

复制代码

需要完整代码的，请查看附件

gpf88

怎么样,问题解决了吗