Discussion :

[cacatozor]

Bonsoir,

Je souhaite faire le produit de convolution de deux fonctions, pour cela, j'ai utilisé le module numpy et son outil convolve.
J'ai testé sur la fonction porte en faisant la convolution de cette fonction par elle-même, mais je n'obtient pas le résultat escompté.

Voici le code:

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
 
import numpy as np
import matplotlib.pyplot as plt
 
def porte(t):
      a=[]
      for i in t:
           if -0.5<=i<=0.5:
               a.append(1.)
           else:
               a.append(0)
       return a
 
t=np.linspace(-1,1,1000)
c=np.convolve(porte(t),porte(t))
plt.plot(np.linspace(-1,1,len(c)),c)

Le graphe ne me donne pas le même graphe que sur wikipédia, je trouve la même forme mais avec une largeur plus petite (entre -0.5 et 0.5 au lieu de -1 et 1).
Savez-vous d'où peut venir ce problème ?

[lg_53]

Primo ton code est super lent alors que tu ne fais quasiment rien. Pardis, utilise la puissance de numpy !

Code :

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

1
2
3
4
5
 
def porte(t):
      res=np.zeros(t.shape)
      res[ (-0.5<=t)*(t<=0.5) ] = 1
      return res

(le and devant être remplacé par un symbole multiplier, le or serait remplacé par un +)

Ensuite une lecture de la documentation de cette fonction de convolution
https://docs.scipy.org/doc/numpy-1.1....convolve.html
t'aurais permit peut-être d'essayer qqch comme ceci :

Code :

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

c=np.convolve(porte(t),porte(t),'same')

[cacatozor]

Merci Ig_53,


En effet, c'est beaucoup plus rapide ainsi et le "same" me donne les bonnes valeurs pour les abscisses. Cependant, la valeur de la fonction en 0 est de 25 alors qu'elle devrait être de 1. En divisant par max(c) on trouve la bonne courbe.