Discussion :

[vivi78]

Bonjour !

Je cherche simplement à entrer la somme suivante sur Matlab :

pour i>=5

N(K(2))*(K(2)-K(1)) + N(K(3))*(K(3)-K(2)) + N(K(3))*(K(3)-K(2)) + ... + N(K(i-1))*(K(i-1)-K(i-2))

où N et K sont des vecteurs donc K(j) renvoie l'élément j de K Et ainsi N(K(j)) renvoie le K(j) ème élément de N.

J'ai essayé la forme suivante :

Code :

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

1
2
3
for i = [3:1:numel(K)]
    sum(Nr(K(2))*(K(2)-K(1)) : Nr(K(i-1))*(K(i-1)-K(i-2)))
end

Donc très logiquement ça marche pour i=3 mais à partir de 4 ça ne renvoie pas les résultats attendus. Le problème vient donc la manipulation sum/ produit de matrice dans sum...

Je vous remercie par avance de vos réponses !

[Myrne]

Je ferais quelque chose comme ça:

Code :

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

1
2
3
4
5
6
 
L=length(k);
 
N2=N(K(2:L).*(K(2:L)-K(1:L-1)));
 
somme=sum(N2);

Ce que tu fais dans ton code ne peut pas marcher. Regarde la doc de sum, cette fonction ne fait que sommer les différents termes de ton vecteur, tu ne devrais pas avoir besoin de faire une boucle for si tu as défini ton vecteur correctement.

Actuellement, ce que tu demandes de faire est un truc du genre :

Code :

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

1
2
3
 
pour i=3:L
<div style="margin-left:40px">somme( N(f(3):f(i)) )</div>fin pour

où f(i) = K(i-1))*(K(i-1)-K(i-2))

Donc ça marche en effet pour la première boucle puisque tu ne sommes qu'un seul terme (!), mais ça n'a aucune raison de marcher pour i=4 si f(4) se retrouve inférieur à f(3). Et si par malheur ta fonction f avait été croissante, tu aurais eu quelque chose qui marche, mais qui n'aurai rien eu à voir avec ce que tu voulais obtenir.

[vivi78]

Merci pour ces explications.
La boucle for ne me sert effectivement à rien pour ce que je demande, si je l'ai mise c'est en fait parce que je veux évaluer la valeur de cette somme de i =3 à i = numel(K).

Cependant quelque chose me surprend d'emblée dans ton code :

Code :

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

1
2
3
4
5
L=length(k);
 
N2=N(K(2:L).*(K(2:L)-K(1:L-1)));
 
somme=sum(N2);

ce qui m'intéresse est : N(K(2:L)).*(K(2:L)-K(1:L-1))
et non pas N(K(2:L).*(K(2:L)-K(1:L-1)))

[Myrne]

Ah ben oui j'ai lu vite fait et j'm'est planté dans le parenthèsage chef.

Cependant, la méthode reste la même ce me semble =)

EDIT: et dans ce cas, pour alléger l'écriture, tu peux même faire ça:

Code :

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

N2=N(K(2:L)).*diff(K);

[vivi78]

Merci beaucoup ça marche super !!!