1 pièce(s) jointe(s)
Méthode iterative par matlab
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:
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:
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
