Discussion :

[Casimirs]

Bonjour à tous,

J'ai un petit problème pour coder une somme qui dépend d'une variable ('w'). Voici ce que je souhaite calculer

\sum_n An/(w-sn) pour 1<n<10 et plotter ensuite cette somme en fonction de w. Voici les codes que j'ai écris pour y arriver :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
w = linspace(1,10,100); % valeurs de w
 
A = [10,22,35,46,55,86,78,88,98,10]; % ma suite de An
s = [15,22,3,24,5,62,75,8,29,10];      % ma suite de sn
 
% la somme pour tous les w
somme = 0;
for i=1:10
somme = somme+A(i)./(w-s(i))
end

ou

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
w = linspace(1,10,100); % valeurs de w
 
A = [10,22,35,46,55,86,78,88,98,10]; % ma suite de An
s = [15,22,3,24,5,62,75,8,29,10];      % ma suite de sn
 
% la somme pour tous les w
somme = 0;
for j=1:100
  for i=1:10
somme(j) = somme+A(i)./(w(j)-s(i))
  end
end

Je m'attends à ce que le résultat final ('somme') soit un vecteur de 100 éléments dont chacun correspond à la valeur de la somme précédente calculée pour chacune des valeurs de w.
Je souhaite enfin plotter cette somme en fonction de w.

Les deux codes précédents me donnent n'importe quoi !

Quelqu'un a t-il une idée pour résoudre ce problème ?

Merci de votre aide !

[Invité]

Bonsoir,

ta deuxième proposition est presque bonne

Code :

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

1
2
3
4
5
6
somme = zeros(1,length(w)); % Voir FAQ
for j=1:100
    for i=1:10
        somme(j) = somme(j)+A(i)./(w(j)-s(i)) % le point '.' n'est pas nécessaire : tu travailles ici avec des scalaires
    end
end

Ou un peu plus optimisé :

Code :

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

1
2
3
for j=1:100
    somme(j) = sum( A./(w(j)-s) ); % ici par contre oui
end