Problème avec transformée de Fourier

Version imprimable

Bonjour,

Je fais mes débuts en C++ et j'essaie de programmer la transformée de Fourier Discrète. Je veux en gros que ma classe TFD dispose de deux méthodes, reel et imaginaire, qui m'affiche la partie réelle et imaginaire de la transformée de Fourier de mon échantillon de départ. Quand je compile mon programme, il m'affiche 'double' is not a template et invalid types 'double[int]' for array subscript. Que dois-je faire svp pour régler ce problème, ça me prend la tête depuis quelques heures?

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
 
 
#include "TFD.hh"
#include <cmath>
#include <vector>
 
using namespace std;
 
 
 TFD::TFD(const int& _dir, const int& _taille, const double<vector>& _x1, const double<vector>& _y1)
                                          :dir(_dir),taille(_taille),x1(_x1),y1(_y1){};
 
 double<vector> TFD::reel() const{
 
       double<vector> x2;
       double arg;
       double cosarg,sinarg;
       int i,k;
 
       for (int i=0;i<taille;i++) {
           x2[i] = 0;
           arg = - dir * 2.0 * 3.141592654 * (double)i / (double)taille;
           for (int k=0;k<taille;k++) {
               cosarg = cos(k * arg);
               sinarg = sin(k * arg);
               x2[i] += (x1[k] * cosarg - y1[k] * sinarg);
               }
       }
 
       if (dir == -1) {
          for (int i=0;i<taille;i++) {
              x2[i] = x2[i] / (double)taille;
              }
          }
 
       return x2;
 }
 
 double<vector> TFD::imaginaire() const{
 
       double<vector> y2;
       double arg;
       double cosarg,sinarg;
       int i,k;
 
 
       for (int i=0;i<taille;i++) {
           y2[i] = 0;
           arg = - dir * 2.0 * 3.141592654 * (double)i / (double)taille;
           for (int k=0;k<taille;k++) {
               cosarg = cos(k * arg);
               sinarg = sin(k * arg);
               y2[i] += (x1[k] * sinarg + y1[k] * cosarg);
               }
       }
 
       if (dir == -1) {
          for (int i=0;i<taille;i++) {
              y2[i] = y2[i] / (double)taille;
              }
          }
 
 
       return y2;
 }
 
 TFD::~TFD(){}

Merci d'avance pour votre aide

06/06/2009, 15h56
Goten

Code:

double<vector>

c'est pas plutôt :

Code:

std::vector<double>

Citation:

Envoyé par Goten

Code:

double<vector>

c'est pas plutôt :

Code:

std::vector<double>

Oui effectivement, merci beaucoup, je n'avais pas vu cette erreur bête. Mais maintenant pas d'erreurs quand je compile mais je n'arrive pas à faire des tests avec mon main.

Voici le fichier .hh

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
 
#ifndef _aze_H_
#define _aze_H_
#include <vector>
using std::vector;
 
class TFD {
 
private:
        int dir;      // prend comme valeur 1 ou -1 qui nous permet de spécifier si on veut une TF directe ou inverse
        int taille;   // Taille de l'échantillon
        vector<double> x1;  // Vecteur représentant la partie réelle de la série de points (xi) dont on cherche la TF
        vector<double> y1;  // Vecteur représentant la partie imaginaire de la série de points (xi) dont on cherche la TF
 
 
 
public:
       TFD(const int& _dir, const int& _taille, const vector<double>& _x1, const vector<double>& _y1); //constructeur par défaut
       vector<double> reel() const; //Méthode renvoyant la partie réelle du résultat de la TF
       vector<double> imaginaire() const; //Méthode renvoyant la partie imaginaire du résultat de la TF
       virtual ~TFD();
};
#endif

Le fichier .cpp

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
 
#include "TFD.hh"
#include <cmath>
#include <vector>
 
using namespace std;
using std::vector;
 
 
 TFD::TFD(const int& _dir, const int& _taille, const vector<double>& _x1, const vector<double>& _y1):dir(_dir),taille(_taille),x1(_x1),y1(_y1){};
 
 vector<double> TFD::reel() const{
 
       vector<double> x2;
       double arg;
       double cosarg,sinarg;
       int i,k;
 
       for (int i=0;i<taille;i++) {
           x2[i] = 0;
           arg = - dir * 2.0 * 3.141592654 * (double)i / (double)taille;
           for (int k=0;k<taille;k++) {
               cosarg = cos(k * arg);
               sinarg = sin(k * arg);
               x2[i] += (x1[k] * cosarg - y1[k] * sinarg);
               }
       }
 
       if (dir == -1) {
          for (int i=0;i<taille;i++) {
              x2[i] = x2[i] / (double)taille;
              }
          }
 
       return x2;
 }
 
 vector<double> TFD::imaginaire() const{
 
       vector<double> y2;
       double arg;
       double cosarg,sinarg;
       int i,k;
 
 
       for (int i=0;i<taille;i++) {
           y2[i] = 0;
           arg = - dir * 2.0 * 3.141592654 * (double)i / (double)taille;
           for (int k=0;k<taille;k++) {
               cosarg = cos(k * arg);
               sinarg = sin(k * arg);
               y2[i] += (x1[k] * sinarg + y1[k] * cosarg);
               }
       }
 
       if (dir == -1) {
          for (int i=0;i<taille;i++) {
              y2[i] = y2[i] / (double)taille;
              }
          }
 
 
       return y2;
 }
 
 TFD::~TFD(){}

Et mon main:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
 
#include <iostream>
#include <cmath>
#include <vector>
#include "TFD.hh"
 
using namespace std;
 
int main(){
    vector<double> v1(4);
    vector<double> v2(4);
    vector<double> Res(4);
    double test;
 
    v1[0]=1.0;
    v1[1]=7.0;
    v1[2]=5.0;
    v1[3]=8.0;
 
    v2[0]=3.0;
    v2[1]=3.0;
    v2[2]=4.0;
    v2[3]=9.0;
 
TFD* voir=new TFD(1, 4, v1, v2);
Res= (voir->reel());
test=Res[0]+Res[1]+Res[2]+Res[3];
 
 
        cout << "La partie réelle de la TFD est R=" << test << endl;
 
 
        system("PAUSE");
delete voir;
voir=0;
}

06/06/2009, 18h47
Dreambeliever

Salut,

Avec un débugger,

on voit que l'erreur est à la ligne:

Code:

Res= (voir->reel());

Puis dans la fonction reel():

Code:

x2[i] = 0;

Parce que tu n'as pas initialisé x2!

utilise le débugger :calim2:
07/06/2009, 14h37
scarabeous
Citation:

Envoyé par Dreambeliever

Salut,

Avec un débugger,

Code:

x2[i] = 0;

Parce que tu n'as pas initialisé x2!

utilise le débugger :calim2:

J'ai utilisé le debugger et je savais que c'était de la méthode réel que venait le problème et effectivement il fallait initialiser le vecteur. Il y aussi un truc qui marche, c'est de mettre la taille du vecteur quand on le déclare, ça résout le problème:
Code:

1 2 vector<double> x2(n);
En tout cas là ça marche bien, merci pour votre aide!!