Music is my life

  제목 : DFT , InverseDFT
  용도 : 음성 신호 처리용
  학번 : 1999108145
  이름 : 이영곤

************************************/

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define PI 3.141592 // PI

typedef struct
{
double real,img;
}COMPLEX; // Complex 구조체 선언

/*****
Module : N-point InverseDFT
input : COMPLEX ARRAY
output : COMPLEX ARRAY
N : Size of DFT
*****/
void idft(COMPLEX *input , COMPLEX *output , int N)
{
int i,k;
for (k = 0 ; k < N ; k++)
{
output[k].real = 0.0; // 초기화
output[k].img = 0.0; // 초기화

for ( i = 0 ; i < N ; i++)
{
   output[k].real += input[i].real*cos(2*PI*i*k/(double)N);
   output[k].img += input[i].real*sin(2*PI*i*k/(double)N);
}

output[k].real = 1/(double)N*output[k].real;// 1/N 으로 Scaling
output[k].img = 1/(double)N*output[k].img; // 1/N 으로 Scaling
}

}

/*****
Module : N-point InverseDFT
input : COMPLEX ARRAY
output : COMPLEX ARRAY
N : Size of DFT
*****/
void dft(COMPLEX *input , COMPLEX *output , int N)
{
int i,k;

for (k = 0 ; k < N ; k++)
{
output[k].real = 0.0; // 초기화
output[k].img = 0.0; // 초기화

for ( i = 0 ; i < N ; i++)
{
   output[k].real += input[i].real*cos(2*PI*i*k/(double)N);
   output[k].img += input[i].real*sin(2*PI*i*k/(double)N);
}
output[k].real = output[k].real;
output[k].img = -output[k].img;
}
}