Discussion :

[Newenda]

Bonjour,

Je débute sur python et je cherche à passer un code de Matlab en ce dernier.
s est un signal de longueur 228352
N = 11025
et je veux créer sif tel que:

voici en matlab:

Code :

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

1
2
3
for ii = 1:round(size(s,1)/N)-1
    sif(ii,:) = s(1+(ii-1)*N:ii*N);
end

que je traduis en python

Code :

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

1
2
3
4
5
6
7
 
iend=np.arange(round(len(s)/N))
iend=len(iend)-1
sif=np.ones((iend,len(s)))
 
for ii in range(iend):
   sif[ii,:]=s[1+ii*N:(ii+1)*N]

avec comme erreur:

ValueError: could not broadcast input array from shape (11024) into shape (228352)

je comprends pas pourquoi il parle des shapes de N et de s là..

Merci pour toute aide.

[wiztricks]

Salut,

Quand vous écrivez:

Code :

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

       sif[ii,:]=s[1+ii*N:(ii+1)*N]

Vous voulez copier les N valeurs successives (1+ii*Nii+1)*N) de s.
Le message d'erreur vous dit seulement que ce N là ne correspond pas au nombre de cases de sif[ii,:]

La vraie question serait plutôt "pourquoi écrire sif=np.ones((iend,len(s))) plutôt que sif=np.ones((iend,N)) ?"

- W

[Newenda]

Merci pour votre réponse,

effectivement, la taille de sif est bien iend, N
Mais j'image que si on se trompe sur la taille de sif avant la boucle ce n'est pas très grave non? juste comme Matlab, ça ne sert qu'à gagner du temps de calcul?

Mais en changeant cela:

Code :

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

sif=np.ones((iend,N))

j'ai toujours le même msg d'erreur

ValueError: could not broadcast input array from shape (0) into shape (206)

.. je ne comprends pas d'où vient le "from shape (0)"

[wiztricks]

.. je ne comprends pas d'où vient le "from shape (0)"

Ca vient toujours de la dimension de la séquence qu'on veut recopier ("from shape").
Par curiosité, vous pourriez lancer la console Python pour voir les différents messages d'erreurs:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>>> import numpy as np
>>> ar = np.ones((3, 5))
>>> ar[0,:] = np.ones(0)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ValueError: could not broadcast input array from shape (0) into shape (5)
>>> ar[0,:] = np.ones(6)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ValueError: could not broadcast input array from shape (6) into shape (5)
>>>

- W

[Newenda]

ha ma première erreur a été d'oublier qu'on incrémente à partir de 1 dans matlab et 0 dans python (pourquoi d'ailleurs, étonnant de commencer à 0 mais bon..)

donc j'ai changé un peu ma boucle

Code :

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

1
2
3
 
  for ii in range(iend):
         sif[ii,:]=s[1+(ii)*N:(ii+1)*N]

mais toujours:

ValueError: could not broadcast input array from shape (11024) into shape (11025)

Désolé, dans ce cas je ne vois pas d'où vient le 11024 et donc je vois pas trop bien le corriger.. le sif.shape me retourne bien (206,11025)

j'essaye naivement:

Code :

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

1
2
3
In [28]: for ii in range(iend):
    ...:     sif[ii,:]=s[1+(ii)*(N+1):(ii+1)*(N+1)]
    ...:

mais me retourne:

/home/anaconda2/bin/ipython:2: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
if __name__ == '__main__':

C'est grave?

[wiztricks]

Désolé, dans ce cas je ne vois pas d'où vient le 11024 et donc je vois pas trop bien le corriger.. le sif.shape me retourne bien (206,11025)

Ce 11024 vient de 1+(ii)*N:(ii+1)*N.
Et comme vous l'avez remarqué Python commence à 0.

- W