快速FFT变换C语言程序

快速FFT变换源程序

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. typedef struct                                                                                        //定义复数结构体变量
   
5. {
   
6.         double real;
   
7.         double imag;
   
8. }complex;
   
9. 
   
10. #define N 5
    
11. 
    
12. #define PI 3.1514926535897932384626433832795028841971        //定义圆周率
    
13. 
    
14. void c_plus(complex, complex, complex *);                                //复数加运算
    
15. void c_sub(complex, complex, complex *);                                //复数减运算o
    
16. void c_mul(complex, complex, complex *);                                //复数乘运算
    
17. void Wn_i(int, int, complex *);                                                    //旋转因子
    
18. 
    
19. int main()
    
20. {
    
21.         complex f[N] = { 1, 2, 2, 3 };
    
22.         complex wn;
    
23.         Wn_i(8, 3, &wn);
    
24.         c_mul(f[0], wn, &f[2]);
    
25.         c_sub(f[1], f[2], &(f[3]));
    
26.         c_plus(f[1], f[2], &(f[4]));
    
27.         printf("f[%d]为：%lf + %lf        \n", 2, f[2].real, f[2].imag);
    
28.         printf("f[%d]为：%lf + %lf        \n", 3, f[3].real, f[3].imag);
    
29.         printf("f[%d]为：%lf + %lf        \n", 4, f[4].real, f[4].imag);
    
30.         return 0;
    
31. }
    
32. 
    
33. void c_plus(complex a, complex b, complex *c)                        //复数加法
    
34. {
    
35.         c->real = a.real + b.real;
    
36.         c->imag = a.imag + b.imag;
    
37. }
    
38. void c_sub(complex a, complex b, complex *c)                        //复数减法
    
39. {
    
40.         c->real = a.real - b.real;
    
41.         c->imag = a.imag - b.imag;
    
42. }
    
43. void c_mul(complex a, complex b, complex *c)                        //复数乘法
    
44. {
    
45.         c->real = a.real*b.real - a.imag*b.imag;
    
46.         c->imag = a.real*b.imag + a.imag*b.real;
    
47. }
    
48. 
    
49. void Wn_i(int n1, int i, complex *Wn)                                        //定义旋转因子
    
50. {
    
51.         Wn->real = cos(2 * PI*i / n1);
    
52.         Wn->imag = -sin(2 * PI*i / n1);
    
53. }
    

复制代码

全部资料51hei下载地址：
fft.zip (699 Bytes, 下载次数: 42)

2019-5-31 15:02 上传

点击文件名下载附件
源程序
下载积分: 黑币 -5