IDFT的公式及C语言实现.doc
百度空间?|?百度首页?
???????????????
sszhou's blog
track, true view, thinking
?
主页博客相册|个人档案 |好友
?查看文章?IDFT的公式及C语言实现
2007-08-16 16:42
IDFT的公式及C语言实现
该公式与DFT的计算公式极其类似,只是符号发生了变化,而且多了一个1/N。可见IDFT的代码核DFT的代码很相似,在此只给出核心的代码:
void my_idft_01(
data_t *xr, // input real part A(k)
data_t *xi, // input image part: B(k)
data_t *yr, // output real part: a(n)
data_t *yi, // output image part: b(n)
int N // N
)
{
int n, k;
data_t Q;
Q = 2*PI/N;
for (k=0; k
{
yr[k] = 0.0;
yi[k] = 0.0;
for (n=0; n
{
// Different with DFT
yr[k] += xr[n]*cos(Q*n*k) - xi[n]*sin(Q*n*k);
yi[k] += xi[n]*cos(Q*n*k) + xr[n]*sin(Q*n*k);
}
}
// Different with DFT * (1/N)
for (k=0; k
{
yr[k] /= N;
yi[k] /= N;
}
}