Blame view

fft/fftcow.cpp 1.11 KB
6c227d28   rsimonin   premier push
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
  #include <iostream>
  #include "wavdata.h"
  #include "fft.h"
  #include <math.h>
  
  
  #define FREQ 22400
  #define AMPLITUDE 10
  #define FREQDO 261
  #define FREQLA	440
  #define SIZE FREQ*4
  
  
  int main(int argc, char **argv)
  {
6a961b68   rsimonin   modifications
16
17
  //Definition variables
      int i;
6c227d28   rsimonin   premier push
18
19
  	float pi=3.141592;
  	WavData w;
6a961b68   rsimonin   modifications
20
21
22
23
24
25
26
27
  	w.load("COW.WAV");
  	char *data = w.data();
  	printf("%d\n",w.datasize());
  	char *data2 = new char[w.datasize()];
  
      double test[w.datasize()*2][2];
  	double test1[w.datasize()*2][2];
  	double test2[w.datasize()*2][2];
6c227d28   rsimonin   premier push
28
  	
6a961b68   rsimonin   modifications
29
30
31
  //Creation de la donnée
  	
  	for(i=0;i<w.datasize();i++){
6c227d28   rsimonin   premier push
32
  		test[i][0]=(double)data[i];
6a961b68   rsimonin   modifications
33
  		test[i][1]=0;
6c227d28   rsimonin   premier push
34
      }
6c227d28   rsimonin   premier push
35
  	
6a961b68   rsimonin   modifications
36
37
38
39
  //FFT
  	fft(32768*2,test,test2);
  
  //FFT inverse
6c227d28   rsimonin   premier push
40
41
  	ifft(32768*2,test1,test2);
  
6a961b68   rsimonin   modifications
42
43
  //Récupération des données
  	for(i=0;i<w.datasize();i++){
3e4da888   rsimonin   samedi 16h
44
45
46
  		//printf("%lf\n",test1[i][0]);
  		//if(test1[i][0]>50)test1[i][0]=50;
  		//if(test1[i][0]<-50)test1[i][0]=-50;
6c227d28   rsimonin   premier push
47
48
49
  		data2[i]=test1[i][0];
  	}
  
6a961b68   rsimonin   modifications
50
  //Preapation format final
6c227d28   rsimonin   premier push
51
52
53
54
55
56
57
58
  	w.setAudioFormat(1);
  	w.setNbrChanel(1);
  	w.setFrequency(FREQ);
  	w.setBytePerBloc(4);
  	w.setBytePerSec(FREQ);
  	w.setBitsPerSample(8);
  	w.clearData();
  	
6a961b68   rsimonin   modifications
59
  	w.setDatasize(w.datasize());
6c227d28   rsimonin   premier push
60
61
62
63
  	w.setData(data2);
  	w.save("cowfft.WAV");
  
  }