Discussion :

[Lhaneslo]

Bonjour à tous,

Ma question risque de vous paraître un peu triviale et je m'excuse d'avance si elle a déjà été posée.

Je voudrais simplement savoir comment faire le produit de convolution de deux fonctions dans le programme R. J'ai essayé plein de choses, mais je n'ai rien de concluant, et je galère pas mal avec la fonction 'convolve'.

Mon idée (qui ne marche pas en l'état) était la suivante: si [a,b] est la réunion (ensemble fini) des domaines de mes deux fonctions f et g.
x = seq(from=a,to=b,length=2*n+1) #2*n est le nombre d'intervalles de longueur (b-a)/(2*n) entre a et b

Code :

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

1
2
3
4
#je vais approcher mon produit de convolution par une somme de Riemann
y1=f(x[-1])
y2=g(x[-1])
P = convolve(y1,y2,type='open')[1: (2*n)]*(b-a)/(2*n)

Si vous avez une réponse à mon problème (ou à mon idée), merci d'avance!

Cordialement
Lhaneslo

[Meseira]

Bonjour,

une solution simple consisterait à revenir à la définition du produit de convolution et d'utiliser la fonction integrate pour calculer l'intégrale. Par exemple, voici comment calculer (f*g)(x) sur le domaine [a,b],

Code :

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

integrate(function(t) { f(x-t)*g(t) }, a, b)