我正在尝试运行以下代码,以在上述模型的DSP处理器中实现FFT。该代码是用C语言编写的

//All codes submitted by  Srivatsa Sinha
/*Include Library Files*/
#include <stdio.h>
#include<math.h>

/*Templates*/
//#define ll long long
#define cmp struct complex
#define print_cmp(z) printf("%lf + %lfj\n",z.re,z.im)
#define M_PI 3.14159265358979

struct complex
{
    float re;
    float im;
}res,w,w_var,t,u;

cmp signal[4];
cmp result[4];
int j,s,k,n,z,x,no_of_bit,rev,pow2[10],m,a;
void calc()
{
    pow2[0]=1;
    for(k=1; k<10; k++) { pow2[k] = pow2[k-1]*2;}
}

cmp add(cmp a, cmp b)
{
    res.re = a.re + b.re;
    res.im = a.im + b.im;
    return res;
}

cmp sub(cmp a, cmp b)
{
    res.re = a.re - b.re;
    res.im = a.im - b.im;
    return res;
}

cmp mul(cmp a, cmp b)
{
    res.re = a.re*b.re - a.im*b.im;
    res.im = a.re*b.im + a.im*b.re;
    return res;
}

cmp cmp_pow(cmp base, int exp)
{
    res.re = 1.0; res.im = 0.0;
    for(a=1; a<=exp; a++)
        res = mul(res,base);
    return res;
}

void bit_rev()
{
    no_of_bit = log2(n);
    for(z=0; z<n; z++)
    {
      rev = 0;
      for(x=0; x<no_of_bit; x++)
      {
        if((z>>x)&1)
        {
          rev = rev|(1<<(no_of_bit-1-x));
        }
      }
      result[rev] = signal[z];
    }
}

void fft()
{
  bit_rev();
  for(s=1; s<=(int)log2(n); s++)
  {
     m = pow2[s];
     w.re = (float)cos(((float)2.0*(float)M_PI)/(float)m); w.im = (float)-1.0*sin(((float)2.0*(float)M_PI)/(float)m);
     for(k=0; k<n; k+=m)
     {
       w_var.re = 1; w_var.im = 0;
       for(j=0; j<=(m/2 - 1); j++)
       {
         t = mul(w_var, result[k+j+(int)(m/2)]);
         u = result[k+j];
         result[k+j] = add(u,t);
         result[k+j+(int)(m/2)] = sub(u,t);
         w_var = mul(w_var, w);
       }
     }

  }
}

void dft()
{
    w.re = (float)cos(((float)2.0*(float)M_PI)/(float)n); w.im = (float)-1.0*sin(((float)2.0*(float)M_PI)/(float)n);
    for(k=0; k<n; k++)
    {
        result[k].re = 0.0; result[k].im = 0.0;
        for(j=0; j<n; j++)
        {
            result[k] = add(result[k],mul(signal[j],cmp_pow(w,k*j)));
        }
    }
}

int main()
{
   calc();
   printf("Enter length of the signal: ");
   scanf("%d",&n);
   printf("Enter the signal sequence: \n");
   for(k=0;k<n; k++)
   {
       scanf("%f",&signal[k].re);
       scanf("%f",&signal[k].im);
   }
   fft();
   printf("The Result of FFT is: \n");
   for(k=0; k<n; k++)
   {
     print_cmp(result[k]);
   }
   dft();
   printf("The Result of DFT is: \n");
   for(k=0; k<n; k++)
   {
     print_cmp(result[k]);
   }
   return 0;
}


在编译之前添加了以下文件:

hello.cmd The Linker file
rts6700.lib The Library file


虽然代码可以成功编译,但是在编译时会出现以下错误:

can't allocate .stack size 00000400 (page 0) in IRAM (avail: 000003ec)
can't allocate .sysmem size 00000400 (page 0) in IRAM (avail: 000003ec)
can't allocate .cio size 00000120 (page 0) in IRAM (avail: 0000010a)
errors in input - ./Debug/FFTa.out not built


请帮忙??

最佳答案

谢谢您的回忆。自从我使用TI工具以来已经有两年了。我不再有权访问其文档,但似乎您需要设置链接器控制文件。如果使用CodeComposer,则需要创建一个目标(或指定一个仿真模型),该目标将定义RAM的位置。如果需要从处理器中获得高性能,则需要挖掘各种内存类型,以最大程度地减少总线延迟。

10-07 18:58
查看更多