標題: 快速FFT變換C語言程序 [打印本頁]
作者: IQ_kuo    時間: 2019-5-31 15:03
標題: 快速FFT變換C語言程序
快速FFT變換源程序
  1. #include <stdio.h>
  2. #include <math.h>

  4. typedef struct                                                                                        //定義復數(shù)結(jié)構(gòu)體變量
  5. {
  6.         double real;
  7.         double imag;
  8. }complex;

  10. #define N 5

  12. #define PI 3.1514926535897932384626433832795028841971        //定義圓周率

  14. void c_plus(complex, complex, complex *);                                //復數(shù)加運算
  15. void c_sub(complex, complex, complex *);                                //復數(shù)減運算o
  16. void c_mul(complex, complex, complex *);                                //復數(shù)乘運算
  17. void Wn_i(int, int, complex *);                                                    //旋轉(zhuǎn)因子

  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. }

  33. void c_plus(complex a, complex b, complex *c)                        //復數(shù)加法
  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)                        //復數(shù)減法
  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)                        //復數(shù)乘法
  44. {
  45.         c->real = a.real*b.real - a.imag*b.imag;
  46.         c->imag = a.real*b.imag + a.imag*b.real;
  47. }

  49. void Wn_i(int n1, int i, complex *Wn)                                        //定義旋轉(zhuǎn)因子
  50. {
  51.         Wn->real = cos(2 * PI*i / n1);
  52.         Wn->imag = -sin(2 * PI*i / n1);
  53. }
復制代碼

全部資料51hei下載地址：
fft.zip (699 Bytes, 下載次數(shù): 43)




歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1