Discussion :

[Por Qué]

Bonjour,

je viens de calculer un coefficient avec le programme suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
 
x=0.1*10^-6;
w=0.08*10^-6;
Lambda=400e-9:10e-9:800e-9; % en nm
c=3e8; 
h= 4.135e-15;
E=h.*c./Lambda;
C = [5030.02, 483.916;1634.30,79.4079]; % en cm-1
C3= 1046.08;  % en cm-1
Ek = 1.09969; % en eV
E3 = 1.40985; % en eV
dE = [0.0583148,0.0220161]; % en eV
N = 0.394122; 
dN = 1.23084; % en eV-1.
a=x+w;
 
for w = 1:2
for j=1:2
                 Alpha= C(j,w)* power(E-Ek+power(-1,j)*dE(w)+abs(E-Ek+power(-1,j)*dE(w))/2,2)./E+C3*power((E-E3)+abs((E-E3)),N+dN*E)./E
 
end
end
 
 exp (-Alpha.*(x+w));

le probleme qui ce pose que la fonction exp donne toujours la valeurs 0 alors que si on prend une valeur du vecteur alpha nous tombe sur une valeur et pas 0.
par exemple:
si on prend la première valeur de Alpha = 5.8517*10^4 alors
exp(5.8517*10^4*(w+x)) = 0.9
alors que le programme me donne un 0
j'ai pas compris ou est le probleme!

Merci d'avance

[shaiHulud]

Bonjour,

Tout d'abord, comme Alpha est un vecteur que tu affiche en un coup, l'affichage matlab va arrondir et donc les valeurs relativement trop faibles seront affichées comme des 0 même si elles sont pas nulles.

Ensuite, ta première valeur de Alpha ne vaut pas 5.8517, mais 5.8517*10^4, et Alpha(1).*(x+w) est de l'ordre de 10^5. Du coup exp(-10^5) vaut bien 0 au niveau de la précision machine

[Por Qué]

d'abord j'ai oublier de mettre le 10^4 mais j'ai ajouté dans l' exp
deuxiement 5.8517*10^4*0.18*10^-6 donne bien un chiffre de l'ordre de 10^-2 et non pas 10^5

[Cheetor]

Bonjour,

En fait le problème vient du fait que tu ré-affectes une nouvelle valeur à w pendant ton programme. En effet , ligne 2:

Code :

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

w=0.08*10^-6;

et ligne 16 :

Code :

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

for w = 1:2

Ton calcul en dernière ligne fait plutôt :
5.8517*10^4*(2+0.18*10^-6), où tu repasses à une puissance de 10^5. Et donc à 0.

Le mieux serait donc de renommer ta valeur en entrée de boucle.