【发布时间】:2016-10-27 11:42:47
【问题描述】:
我正在尝试运行以下代码以在上述型号的 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
请帮忙??
【问题讨论】:
-
在包含这个 rts6700.lib 之前,是否已经为您的目标 DSP 构建了一个简单的“hello.cmd”?
-
您的 IRAM 快用完了。要么更改内存配置以增加 IRAM,要么将这些部分移动到其他部分,如 SDRAM 或其他部分。
-
您能否解释一下 IRAM 和 SDRAM 的术语,以及我们如何将部分传输到 SDRAM?我是一名本科生,对 DSP 处理器很陌生,因此不知道。它是我每天在大学实验室试验技术的一部分而编写的代码
标签: c signal-processing fft processor