음성,혹은 영상 신호 처리에사용하는 fourier transform 으로 fft 보다 느리지만 구현이 쉽다는 이유로 만들었다. ㅡㅡ;
/************************************
제목 : 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;
}
}
