Discussion :

[nakor]

Bonjour!
Apres la FFT, je m'attaque a la convolution, toujours dans NRiC.
(dans ce programme on cherche a realiser la convolution d'une serie de N donnees "Data[N]", par une fonction reponse "response[N]"

J'ai un peu du mal a comprendre comment est rempli le tableau response :

The response function must be stored in wraparound order: The first half of the array respns contains the impulse response function at positive times, while the second half of the array contains the impulse response function at negative times, counting down from the highest element respns[m-1]

Pourquoi doit-on positionne les elements negatifs a la fin du tableau ?
Sachant qu'en suite on va realiser une operation du type

Code :

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

1
2
no2=n>>1;
resul=data*respons/no2;

avec dataout et responseout les transformees de fouriers de nos deux series de donnes, data et response.

En faisant confiance a NRiC sur la repartition des points, meme si je la comprend pas pour l'instant, j'ai rempli mon tableau response comme suit (ma fonction response est la fonction de Lorentzian) :

Code :

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

1
2
 
data=0,5*x/(PI*(-jo)*(-jo)+(0,5*x)*(0,5*x));

Merci de votre aide !

[Jerome Briot]

C'est (pour le moment) plus un problème d'algorithme que de C

Je déplace donc dans un forum plus adéquat

[nakor]

Le code venant d'un livre de C (numerical recipie) je pensai que quelqu'un pourrait l'avoir deja utilise. Mais c'est vrai que c'est plutôt mathématiques comme problème.

D'ailleurs, quelqu'un aurait il une idée pour vérifiée si mon algorithme fonctionne en utilisant la convolution de deux fonctions qui donnerait un résultat connu, ou facile a vérifier ?