complexe dans matlab

Version imprimable

03/05/2011, 10h58
autocadNUL

complexe dans matlab

Bonjour!
voilà mon code me génère des complexe alors que je n'en veux pas. Ceux ci on une partie imaginaire ridicule comparé à la partie réelle (50 ordre de grandeur de différence!) j'ai donc pensé à une erreur numérique. Est ce possible avec matlab?
03/05/2011, 11h24
Jean Dumoncel

Bonjour,

difficile de te répondre sans connaître le calcul effectué et donc le code...

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
Ttampon=0;
Tdelta=Tg;
while abs(Tdelta-Ttampon)>1e-6
    %calcul du libre parcours moyen
    L=k*Tdelta/(sqrt(2)*Pref*pi*(3.15e-10)^2);
 
    %détermination de l'épaisseur de la couche limite
    gamma1=1+L/r;
    gamma2=1+(L/r)^2;
    delta=r/(((r/L)^2)*(0.2*gamma1^5-(gamma1^3)*gamma2/3+2/15*gamma2^2.5)-1);
    somme2=(calc_cv(T)-calc_cv(Tdelta));
    Ttampon=Tdelta;
    tampon2=alpha*(R_GP/(2*M_v)+somme2);
    Tdelta=T-(r+delta)*k0*((Tdelta-Tg)^(w+1))/((w+1)*(Tg^w)*mu*tampon2);
end
%calcul des fonction des équation différencielles
tampon=1/(3*r*rho_s*Cp_p);
dT=tampon*((eps*q/4+sigma*(Te^4-T^4))-mu*(delta_H/M_s+tampon2*(T-Tdelta)))

voila une partie du code. La puissance 2.5 utilisée peut expliquer une partie de l'erreur, mais je n'y rentre pas une valeur négative...

03/05/2011, 11h34
Jean Dumoncel

Et quel est la ligne de code qui renvoie un nombre complexe? As-tu essayé de mettre des points d'arrêt pour analyser ton code?
03/05/2011, 11h37
autocadNUL

A vrai dire je ne sais pas trop quelle ligne est concernée... Le truc c'est que ce bout de code est dans une fonction. Je ne peut donc pas mettre de point d'arret!
03/05/2011, 11h40
Jean Dumoncel

Bien sur que tu peux mettre des points d'arrêt : clic droit sur la ligne sur lesquelles tu soupçonnes un problème, puis "Set/Clear breakpoint". Un premier clic ajoute un point d'arrêt, un deuxième le retire. Tu peux en mettre plusieurs pour vérifier ainsi à chaque étape l'état de tes variables.
03/05/2011, 11h41
autocadNUL

La ligne 5 renvoye un nombre complexe de temps en temps! pourtant elle n'a que des réels en argument non?
03/05/2011, 12h10
autocadNUL

rebonjour
je crois que j'ai compris! à la ligne 14, je fait une différence à la puissance w+1. Est il possible que matlab me produise un complexe si cette différence est négative. Car dans la résolution de cette équation même si cette différence est normalement positive, la récursivité fait que occasionnellement elle peut être négative.

merci
03/05/2011, 12h15
autocadNUL

c'est bien la raison, vérifié et approuvé...