求助，做FFT只有第一次计算结果正确，后面的结果越来越大

cao651474007 · 发表于 2016-10-20 13:06:37

我的程序在C编译器上都很正常，但是移植到STM32上会出现问题。具体情况如下：
C编译器执行FFT程序很正常，结果也都正确；但是移植到STM32循环执行FFT，就只有复位后第一次的计算结果是正确的，后面的数据都是越来越大，我排查程序发现其中两条语句屏蔽掉后就不会出现计算结果越来越大的情况，但是这两条语句中也没出现数据没清零的情况。

cao651474007 · 发表于 2016-10-20 13:06:38

这个问题解决了，还是程序的问题：计算过程中，我将得到的相位放在了需补数组中，造成后面循环中的虚部部位零

cao651474007 · 发表于 2016-10-20 13:18:06

[mw_shl_code=c,true]void data_count( float re[64][12],float DATA_OUT[5][6])
{


   #define pi 3.1415926
   #define n 64
   #define m 6
   #define Ntab 4
   unsigned int i,j,k,p,q,s,x,a;
// u8 a=0;
   float tr,ti,vr,vi,rms,wr,wi;

for (a=0;a<1;a++)
//----------------------------------------------------------------------------------------------------
{

   i=n;
   for (j=0;j<n;j++)
   {
   p=j;
   q=0;
   for (i=0;i<m;i++)
   {
   k=(p>>1);
   q=(q<<1)+(p-(k<<1));
   p=k;
   }
   if (q>j)
   {
   vr=re[j][a];
   re[j][a]=re[q][a];
   re[q][a]=vr;
   vi=im[j][a];
   im[j][a]=im[q][a];
   im[q][a]=vi;
   }
   }

//----------------------------------?????????----------------------------------------------
   for (i=0;i<m;i++)
   {
   p=(1<<(i+1));
   q=(p>>1);
   for (j=0;j<q;j++)
   {
   x=(1<<(m-i-1))*j*Ntab;
   vr=costab[x];
   vi=sintab[x];
   for (k=j;k<n;k+=p)
   {
   s=k+q;
   tr=re[a]*vr-im[a]*vi;
   ti=re[a]*vi+im[a]*vr;

   wr=re[k][a]-tr;
   re[a]=wr;

   wi=im[k][a]-ti;
                                                im[a]=wi;

//-------------------------------------------------------------------------------------
      wr=re[k][a]+tr;
      re[k][a]=wr;
               //这四行代码屏蔽掉就不出现结果越来越大的情况
      wi=im[k][a]+ti;
      im[k][a]=wi;

//--------------------------------------------------------------------------------------
   }
   }
   }
}[/mw_shl_code]