Problème avec transformée de Fourier

**scarabeous** · 06/06/2009, 15h28

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**Goten** · 06/06/2009, 15h56

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 double<vector>

c'est pas plutôt :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 std::vector<double>

**scarabeous** · 06/06/2009, 16h26

Envoyé par Goten

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 double<vector>

c'est pas plutôt :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**Dreambeliever** · 06/06/2009, 18h47

Salut,

Avec un débugger,

on voit que l'erreur est à la ligne:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

Res= (voir->reel());

Puis dans la fonction reel():

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

           x2[i] = 0;

Parce que tu n'as pas initialisé x2!

utilise le débugger

**scarabeous** · 07/06/2009, 14h37

Envoyé par Dreambeliever

Salut,

Avec un débugger,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

           x2[i] = 0;

Parce que tu n'as pas initialisé x2!

utilise le débugger

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
vector<double> x2(n);

En tout cas là ça marche bien, merci pour votre aide!!

Problème avec transformée de Fourier

C++

Discussions similaires

Partager

Partager