1 pièce(s) jointe(s)
Produit de convolution par une gaussienne
Bonjour à tous,
Dans le cadre de mes études je suis amené (à un moment) à devoir filtrer un signal de sortie par une courbe gaussienne.
Après pas mal de recherches, il s'avère que c'est un produit de convolution que je dois faire.
En gros j'ai un spectre discrétisé dans un tableau de double "x" [0 à 1024] avec les valeurs pour chacun des x[i] dans un tableau de double "y".
Je dois filtrer ce signal par une gaussienne d'écart type 0.07 et mettre les nouvelles ordonnées dans un tableau de double z.
Voici la représentation de mon spectre grâce à JFreeChart que j'ai obtenu à l'aide de mes tableau x et y :
Pièce jointe 169386
Cela fait deux jours que je cherche, et je commence vraiment à désespérer :(
Merci d'avance pour votre aide,
Monax
Edit: j'ai aussi ce code là:
Code:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
| public static double[] convolute( double[] v , double[] m )
{
final int mid = m.length / 2;
double[] r = new double[v.length];
double sum;
for (int i=0, ie=v.length; i<ie; i++)
{
sum = 0;
for (int jv=i-mid, j=0, je=m.length; j<je; j++, jv++)
{
if (jv >= 0 && jv < v.length)
{
sum += m[j];
r[i] += v[jv] * m[j];
}
}
r[i] /= sum;
}
return r;
} |
mais la il me faudrait une discrétisation d'une gaussienne d’écart type 0.07 sur 1024 canaux dans un tableau de double, quelqu'un sait comment faire svp?
