离散傅里叶变换（DFT） C代码

是《数字信号处理C语言程序集》的DFT的C语言源代码，已测试！

#include“math.h“
void dft（xyabnsign）
	/*x-双精度实型一维数组，长度为n。存放要变换数据的实部。
	   y-双精度实型一维数组，长度为n。存放要变换数据的虚部。
	   a-双精度实型一维数组，长度为n。存放变换结果的实部。
	   b-双精度实型一维数组，长度为n。存放变换结果的虚部。
	   n-整型变量。数据长度。
	   sign-整型变量。当sign=1时，子函数dft（）计算离散傅里叶正变换；当sign=-1时，dft（）计
	   算离散傅里叶反变换。
	   	*/
    int  nsign;
    double x[]y[]a[]b[];
{
	
    int ik;
    double cdqws;
    q = 6.28318530715/n;
    for（k=0;k    {
        
        w = k*q;
        a[k] = b[k] = 0.0;
        for（i=0;i        {
            d = i*w;
            c = cos（d）;
            s = sin（d）*sign;
            a[k] += c*x[i] + s*y[i];
            b[k] += c*y[i] - s*x[i];
        }
        
    }
    if（sign == -1）
    {
        c=1.0/n;
        for（k=0;k        {
            a[k] = 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1002  2009-07-21 17:13  dft.c

     文件       2002  2009-07-21 17:13  dft.m

     文件       3252  2009-07-21 17:12  DFT.OBJ

     文件      27590  2009-07-21 17:12  DFT.EXE

----------- ---------  ---------- -----  ----

                33846                    4