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

2019-10-16 02:26发布

我的程序在C编译器上都很正常,但是移植到STM32上会出现问题。具体情况如下:
C编译器执行FFT程序很正常,结果也都正确;但是移植到STM32循环执行FFT,就只有复位后第一次的计算结果是正确的,后面的数据都是越来越大,我排查程序发现其中两条语句屏蔽掉后就不会出现计算结果越来越大的情况,但是这两条语句中也没出现数据没清零的情况。
友情提示: 此问题已得到解决,问题已经关闭,关闭后问题禁止继续编辑,回答。
该问题目前已经被作者或者管理员关闭, 无法添加新回复
15条回答
cao651474007
1楼-- · 2019-10-16 05:15
 精彩回答 2  元偷偷看……
cao651474007
2楼-- · 2019-10-16 05:22
[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]
cao651474007
3楼-- · 2019-10-16 10:26
二楼是代码,麻烦哪位可以帮忙看一下,找找原因
shikihane
4楼-- · 2019-10-16 10:43
太难了
cao651474007
5楼-- · 2019-10-16 11:07
刚刚一直运行发现,到后面会变成一个固定数值左右变化的情况,但是C编译器运行输出一直都是固定的值
不吃鱼的老鼠
6楼-- · 2019-10-16 11:51
 精彩回答 2  元偷偷看……

一周热门 更多>