快速傅里叶算法C语言完整实现.docVIP

快速傅里叶算法C语言完整实现 #includestdio.h #includemath.h #define uint unsigned int #define uchar unsigned char typedef struct Complex//复数结构体 { double Re; double Im; }Complex; Complex Mul(Complex a,Complex b)//复数相乘 { Complex c; c.Re=a.Re*b.Re-a.Im*b.Im; c.Im=a.Re*b.Im+a.Im*b.Re; return c; } Complex Plus(Complex a,Complex b)//复数相加 { Complex c; c.Re=a.Re+b.Re; c.Im=a.Im+b.Im; return c; } Complex Sub(Complex a,Complex b)//复数相减 { Complex c; c.Re=a.Re-b.Re; c.Im=a.Im-b.Im; return c; } /****************************************************************************************** 函数说明：倒序运算 参数说明：N点数，等于2的M次方，*x序列的首地址 返回值：无 ******************************************************************************************/ void Invert_order(uint N,Complex *x,uchar M) { uint power,i,j,m,tempp;//m是i的倒位序数 Complex temp;//是临时变量 for(i=0;iN;i++) { power=1; m=0; for(j=0;jM;j++) { tempp=i; tempp=(M-1-j); if(tempp0x01) { m=m+power; } power=power*2; } if(im) { temp.Im=(x+i)-Im; temp.Re=(x+i)-Re; (x+i)-Re=(x+m)-Re; (x+i)-Im=(x+m)-Im; (x+m)-Im=temp.Im; (x+m)-Re=temp.Re; } } } /****************************************************************************************** Invert_order函数结束： ******************************************************************************************/ /****************************************************************************************** 函数说明：指数运算 参数说明：BaseNumber底数,IndexNumber指数 返回值：无符号整数 ******************************************************************************************/ uint Pow(uchar BaseNumber,uint IndexNumber) { uint m=1;//指数函数值 uchar i; for(i=0;iIndexNumber;i++) { m=BaseNumber*m; } return m; } /****************************************************************************************** Pow函数结束： ******************************************************************************************/ /****************************************************************************************** 函数说明：对数运算 参数说明：BaseNumber底数,An