Discussion :

[mewiine]

Bonsoir à tous !

dans le cadre de mon cours d'info, je dois réaliser un programme de traitement audio sous Lazarus.

Pour ce faire, je récupère tout d'abord les amplitudes dans un fichier texte grâce à Matlab. Ensuite, afin de tracer le spectre en fréquences, je dois calculer les transformées de Fourier. Ne souhaitant pas utiliser FFT, je dois calculer chaque fréquence manuellement. J’ai trouvé ceci après quelques recherches : http://www.google.fr/url?sa=t&rct=j&...7RvOSW3qO_AZMg
J'utilise alors les deux formules de la page 5.

Cependant, j’ai quelques problèmes :
- Pour les destructeurs, j’obtiens un « illegal qualifer » pour les free, je ne comprends pas pourquoi..
- Après un certain temps de calcul, Lazarus me renvoie un message d'erreur "external SIGSEGV". Mais même en utilisant le debugger je ne vois pas où ça coince..

Qu'en pensez-vous ?
Je vous remercie d'avance pour votre aide, et bonnes soirée !

[Paul TOTH]

Projet
Options du projet...
Options du compilateur
Code
Intervalle (-Cr)

tu coches la case Intervalle ci-dessus et tu relances le projet

[mewiine]

Merci de ta rapide réponse !

Dur dur la doc Lazarus.. A quoi correspond cette option ?

J'ai cocher cette case, j'ai maintenant une erreur "RunError(58336)" avant même de rentrer dans la boucle.. =/

[Paul TOTH]

Merci de ta rapide réponse !

Dur dur la doc Lazarus.. A quoi correspond cette option ?

J'ai cocher cette case, j'ai maintenant une erreur "RunError(58336)" avant même de rentrer dans la boucle.. =/

cette option ajoute du code à la compilation qui s'assure que tu n'utilises pas un index hors limite, exemple :

Code :

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

1
2
3
4
5
6
7
8
 
type
 TMonTableau = array[0..2] of Byte;
var
  T: TMonTableau;
begin
  T[12] := 1; // Big ! erreur
end;

Dans ton code tu fixes des dimensions arbitraires à tes tableaux, et les valeurs choisie sont TRES inférieures aux dimensions requises.

Tu devrais utiliser des tableaux dynamiques

Code :

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

1
2
3
4
5
6
7
 
var
  T: array of Byte;
begin
  SetLength(T, 13); // donc de 0 à 12
  T[12] := 1; // pas de problème
end;

et pour avoir une aide digne de ce nom il faut taper dans Delphi 6 ou 7

[mewiine]

Dac!

J'ai corrigé les tableaux, ainsi qu'une petite erreur dans la boucle de somme. Si je n'ai dès lors plus de message d'erreur, j'ai maintenant un calcul interminable, l'appui sur le bouton tableau ne donnant rien.. D'où le problème peut-il venir ?

Edit : Ah si, autant pour moi, après un certain temps j'obtiens un "range check error COL[1] ROW[100]".. pourquoi ici ?
De plus, et oh miracle!, des valeurs sont présentes dans le tableau (en espérant qu'elles correspondent bien à celles voulues). D'ailleurs, pour ce qui est du côté mathématiques de la FFT, mon "adaptation" semble-t-elle correcte ?

Merci !

[Paul TOTH]

pour le range check c'est que ton tableau est mal dimensionné.

pour le temps de cacul, tu devrais sortir la création des Ccomplexe de la boucle car c'est très consommateur de temps, je n'ai plus ton code ne tête mais en gros

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
var
  C, Cn: Ccomplexe;
begin
  C := Ccomplexe.create;
  Cn := Ccomplexe.Create;
  for ...
  begin
    c.init(...) // au lieu d'un constructor, une methode qui  initialise le composant
   ...
  end;
  C.Free;
  Cn.Free;
end;

quand à la validité mathématique ....

[mewiine]

J'ai ajouté une méthode Init :

Code :

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

1
2
3
4
5
6
procedure Ccomplexe.init;
begin
  fRho:=0;
  fTheta:=0;
  SetCartesien(fRho,fTheta);
end;

et changé mon code comme suit :

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
constructor CTabSpectre.create(T:CTableau);
var Theta:real; n,k:integer; C,Cn:Ccomplexe;
begin
  SetLength(fFreq,T.fNe);
  SetLength(fAmp,T.fNe);
  C:=CComplexe.create_polaire(0,0);
  Cn:=CComplexe.create_polaire(0,0);
  for k:=0 to T.fNe-1 do
  begin
  fFreq[k]:=(k*(1/T.fTe))/T.fNe;
     for n:=0 to T.fNe-1 do
     begin
       Theta:=(-2*k*pi)/T.fNe;
       Cn.init;
       Cn.SetCartesien(T.fX[n],Theta*n);
       C.somme(Cn);
     end;
  fAmp[k]:=1/T.fNe*C.GetRho;
  end;
  Cn.destroy;
  C.destroy;
end;

Cependant, le temps de réponse est toujours aussi long.. Était-ce bien ce qu'il fallait faire ?

[Paul TOTH]

pas exactement...

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
 
constructor CTabSpectre.create(T:CTableau);
var Theta:single; n,k:integer; C,Cn:Ccomplexe;
begin
  SetLength(fFreq,T.fNe);
  SetLength(fAmp,T.fNe);
  C:=CComplexe.create;
  Cn:=CComplexe.create;
  for k:=1 to T.fNe do
  begin
  fFreq[k]:=(k*(1/T.fTe))/T.fNe;
   C.init_polaire(0,0);
     for n:=0 to T.fNe-1 do
     begin
       Theta:=(-2*k*pi)/T.fNe;
       Cn.init_polaire(T.fX[n],Theta*n);
       C.somme(Cn);
     end;
  fAmp[k]:=1/T.fNe*C.Rho;
  end;
  Cn.destroy;
  C.destroy;
end;

où la méthode init_polaire reprend le code du constructor.

Code :

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

1
2
3
4
5
6
7
 
procedure cComplexe.init_polaire (aRho,aTheta:real);
begin
  fRho:=aRho;
  fTheta:=aTheta;
  SetCartesien(aRho,aTheta);
end;

ensuite il faudrait utiliser des Single ou Double au lieu de Real

ceci dit avec 5345.656.996 itérations ça prend du temps

[mewiine]

Okay! En fait je remplace mes constructeurs de Ccomplexe par les Init avec paramètres, c'est ça ? Du coup, si je ne redéfinit pas les Create ce seront ceux par défaut de TObject qui seront utilisés ?

Tu remplaces les Real par des Single, c'est pour optimiser la mémoire ? Nous n'avons jamais utilisé que des real en cours, du coup c'est un peu flou..

[mewiine]

Une autre question :

Afin de tracer mon histogramme (Amplitudes en fonction des Fréquences), je pense devoir utiliser un TChart, mais je ne trouve pas comment.. J'ai bien trouvé la méthode AddSerie mais je ne vois pas quel TBasicChartSerie mettre en paramètre. Dois-je créer autant de série que de barres voulues ? (celles-ci devraient se trouver tous les k*fFe)

[Paul TOTH]

Okay! En fait je remplace mes constructeurs de Ccomplexe par les Init avec paramètres, c'est ça ? Du coup, si je ne redéfinit pas les Create ce seront ceux par défaut de TObject qui seront utilisés ?

oui

Tu remplaces les Real par des Single, c'est pour optimiser la mémoire ? Nous n'avons jamais utilisé que des real en cours, du coup c'est un peu flou..

non c'est que Single et Double sont natif au processeur, pas Real

Une autre question :

Afin de tracer mon histogramme (Amplitudes en fonction des Fréquences), je pense devoir utiliser un TChart, mais je ne trouve pas comment.. J'ai bien trouvé la méthode AddSerie mais je ne vois pas quel TBasicChartSerie mettre en paramètre. Dois-je créer autant de série que de barres voulues ? (celles-ci devraient se trouver tous les k*fFe)

pourquoi ne pas simplement dessiner un graphe 2D dans un Bitmap ? une simple règle de trois et tu pourras déterminer les positions x,y des points

[mewiine]

Je peux tracer mon tableau, merci !

Je m'occupe donc maintenant de mon histo. Le composant TChart n'est pas adapté à ton avis ? Je trouvais pourtant ça cool, à condition de savoir s'en servir ^^
Concernant le Bitmap je ne vois pas trop comment procéder, avec cette histoire de règle de trois..

[mewiine]

En cherchant un peu d'aide sur le composant TChart j'ai codé ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
procedure TForm1.Button1Click(Sender: TObject);
var T:CTableau; S:CtabSpectre; i:integer;
begin
  T:=CTableau.create;
  S:=CTabSpectre.create(T);
  for i:=0 to T.fNe-1 do
      SpectreBarSeries1.AddXY(round(i*T.fNe),round(S.fAmp[i]*1000000000));
  T.destroy;
  S.destroy;
end;

en ayant crée auparavant un TBarSerie : SpectreBarSeries1

Cependant, à l’exécution lorsque je clique sur le bouton rien ne se passe.. Que manque-t-il ?

[Paul TOTH]

la dernière fois que j'ai utilisé ce composant c'était en 2002 pour l'écriture de "Delphi 7 Studio" aux éditions Eyrolles ^^

je dirais (d'après la page 115 dudit ouvrage)...

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
procedure TForm1.Button1Click(Sender: TObject);
var T:CTableau; S:CtabSpectre; i:integer;
begin
T:=CTableau.create;
S:=CTabSpectre.create(T);
for i:=0 to T.fNe-1 do
SpectreBarSeries1.Add({XY(round(i*T.fNe),round(}S.fAmp[i]{*1000000000)});
T.destroy;
S.destroy;
end;