Discussion :

[cocanouar]

Bonjour,
J'ai une équation itérative et je suis débutant dans la programmation et je cherche un aide pour résoudre cette équation ( dans les pièces jointes):
avec

Code :

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

1
2
3
4
k0=1;
a=1;
k=0.5;
p=0.1

avec Re c'est la partie réal et Im c'est la partie imaginaire

et Voilà le petit programme que j'ai essayé :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
k0=1;
a=1;
k=0.5;
p=0.1;
for j=1:50,
    Re(j)=(1+p/(1+(pi^2/16)*(k0*a)^4*(j/k)^4))/((1+p/(1+(pi^2/16)*(k0*a)^4*(j/k)^4))^2+p^2*((j/k)+1/((pi/4)*(k0*a)^2*(j/k)^2))^2);
    Im(j)=-p*((j/k)+1/((pi/4)*(k0*a)^2*(j/k)^2))/((1+1/(1+(pi^2/16)*(k0*a)^4*(j/k)^4))^2+p^2*((j/k)+1/((pi/4)*(k0*a)^2*(j/k)^2))^2); 
end
plot(Re,'r');
 
pause;
print -djpeg figure 1
end

Merci d'avance

Bonne année

[Caro-Line]

Et quel est le problème exactement ?

[cocanouar]

Merci de m'avoir répondu,

Le problème ce que dans l'algorithme que j'ai pas respecter lé réglé du méthode altération et précisément la condition de convergence |k(j+1)-K(j)|<epsilon

Bonne année

[Caro-Line]

J'avoue que je ne comprends pas :
Quel critère de convergence ?
Les variables k(j) n'apparaissent nulle part ni dans ton code, ni dans l'équation.

[cocanouar]

Le k = Re(k) + Im(K)

le k apparu dans la faormule qui est dans la piéce jointe et même dans le code j'ai écris la partie réal et la partie imaginaire de k .
avec Re(j) signifié Re(K)(j)
Merci d'avance

[Caro-Line]

Peut-être mais n'empêche qu'il nous manque des éléments pour pouvoir t'aider.
Quel est exactement l'énoncé du problème ?

Là tu nous a d'abord donné une équation puis ensuite tu parles d'un critère de convergence.

Je ne comprends pas en plus quel genre de problème tu rencontres : est-ce réellement un problème de programmation MATLAB ou plus un souci de mise en place de l'algorithme ?