C++

Une fonction qui ne fonctionne qu'une fois !!

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : Une fonction qui ne fonctionne qu'une fois !!

pierrot7777

Bonjour,
Je suis en train d'ecrire un programme qui effectue entre autre du traitement de donnees.
Je dois calculer la magnitude de deux vecteurs, j'ai pour cela ecrit une fonction Magnitude, que j'utilise de la maniere suivante :

Code :

MagnitudeBgT= Magnitude(BgTLp,BgT_size,RedFactor,Nfft);
MagnitudeTemplate = Magnitude(TemplateLp,Template_size,RedFactor,Nfft);

Je ne pense pas que cela soit tres important mais ci-dessous vous trouverez l'implementation de la fonction.

Peu importe les parametres, toujours est-il que la fonction renvoit des resultats corrects la premiere fois (ici donc le tableau MagnitudeBgT sera rempli correctement), tandis que le second tableau ne le sera pas (premiere case du tableau non nulle puis que des zeros).

Si j'ordonne mes lignes dans ce sens:

Code :

MagnitudeTemplate = Magnitude(TemplateLp,Template_size,RedFactor,Nfft);
MagnitudeBgT= Magnitude(BgTLp,BgT_size,RedFactor,Nfft);

C'est exactement l'inverse ! MagnitudeTemplate est correctement rempli mais pas MagnitudeBgT...
Est ce que l'appel d'une fonction peut la modifier ???

Voici le code de la fonction

Code :

double* Magnitude(double* input,int len_input,float fact_red,int nfft)
{
complex_t * datafft = new complex_t[nfft];
complex_t * resfft = new complex_t[nfft];
double* magnitude = new double[nfft/2]; //we keep only half of the spectrum since it is symetric
double E=1; //total energy
datafft[0].im = 0;
datafft[0].re = input[0]/nfft;
if(nfft > len_input)
{
for(int n = 1;n<len_input;n++) //filter to increase high frequencies. fact_red must be between 0 and 1. 0.7 gives good results.
{
datafft[n].im = 0;
datafft[n].re = (input[n] - fact_red*input[n-1])/nfft;
}
for(int n = len_input;n<nfft;n++) //filter to increase high frequencies. fact_red must be between 0 and 1. 0.7 gives good results.
{
datafft[n].im = 0;
datafft[n].re = 0;
}
}
else
{
for(int n = 1;n<nfft;n++) //filter to increase high frequencies. fact_red must be between 0 and 1. 0.7 gives good results.
{
datafft[n].im = 0;
datafft[n].re = (input[n] - fact_red*input[n-1])/nfft;
}
}
resfft = FFT_simple(datafft,nfft);
for(int n = 0;n<nfft/2;n++)
{
magnitude[n] = module(resfft[n]);
}
E = sum(magnitude,0,nfft/2);
for(int n = 1;n<nfft/2;n++)
{
magnitude[n] = magnitude[n]/E;
}
return magnitude;
}

Merci pour votre aide !

Publicité

Taz

bisounours-codeur

Vu qu'en C++ les indices de tableaux vont de 0 à n-1, t'as du chemin à faire ...

pierrot7777

Bingo,
en effet dans la derniere ligne je vais de 1 a nfft/2, au lieu de 0 a nfft/2-1 ... l'habitude de matlab !
Merci Taz, ca n'explique qu'a moitie le probleme mais au moins maintenant ca marche !

Joel F

Real men use unique_ptr

ca l'expique totalement. En sortant du tableau, tu niaque les données de celui d'aprés qui a de forte change d'etre allouées de manières contigues.

Apres t'as aussi le droit d'utiliser des libs de FFT plutot que de reinventer mal la roue : goto FFTW

Taz

bisounours-codeur

C'est pas qu'à la dernière ligne, c'est partout ...

pierrot7777

Non non c'est juste la derniere ligne parce que j'ai implemente une fonction sum qui fait que je me suis plante dans les indices. Pour les boucles, je fais i<nfft/2 (et pas i<=nfft/2 )

Joel, j'aimerais bien utiliser FFTW mais lorsque je suis alle sur leur site, impossible de telecharger une version pour windows (je travaille avec msvc), le lien est mort, et j'ai jamais reussi a trouver un autre miroir ... si l'un de vous a une version qui fonctionne avec windows et s'il pouvait m'envoyer ca par mail je suis plus que preneur !

Taz

bisounours-codeur

OK, autant pour moi, en C++ les indices de tableaux ne vont de 0 à n-1 que sur la dernière ligne.

FORUM HardWare.fr

Programmation

C++

Une fonction qui ne fonctionne qu'une fois !!

Sujets relatifs
AJAX fonction globale?	inline et pointeur sur fonction
fonction exit et CSS	Scheme : fonction intersection
[Besoin d'aide]: php dans requête SQL qui ne fonctionne pas	erreur dans ma fonction
[GTK+] Ajouter un boutton a une fenetre grace à une fonction	Fonction substring dans un swf
pointeur vs reference en retour de fonction	Apelle de la fonction Click de composant crée ou cour de l'execution
Plus de sujets relatifs à : Une fonction qui ne fonctionne qu'une fois !!

Page générée en 0.037 secondes