Bonjour!
Voici mon probleme:
J'ai une matrice M(n,m), ainsi que trois vecteurs A(m), B(m) et C(m). Je souhaite construire la matrice N(n,m), où, pour tout couple (n,m), j'ai:
N_{m,n} = \sum_{i=1}^m \frac{A_i}{B_i - M_{m,n}} C_i
Evidemment, avec des boucles imbriquées, c'est facile, mais c'est horriblement lent. Actuellement, mon code passe 50% de son temps à calculer ça, alors que ce n'est pas la partie "essentielle".
Je pense que la solution passe par la vectorisation, mais je ne sais pas comment la mettre en place, notamment en raison du B_i - M_{m,n}.
Actuellement, le code a cette tête:
Merci d'avance pour votre aide!
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7 for k = 1:n for j=1:m Zkj = Z(k,j); N(k,j) = sum( (A./(B - Zkj )).*C); end end
Partager