Discussion :

[arhemza]

bonjour

j'ai un problème avec la fonction fsolve, malgré mes recherches dans les problèmes qui sont déjà cité dans le forum ainsi que sur l'internet, je n'ai pas arrivé à trouver la solution, donc je besoin de vous.


j'ai un système d'équations non linaire 5*5 (5 équations, 5 inconnues).

quand j'ai fait l'appelle de la fonction fsolver il me donne toujours cet erreur.
merci.

[Beltharion]

Sans le code qui va avec, nous allons avoir du mal a t'aider.

[arhemza]

voici le système de 5 équations.

Code :

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

1
2
3
4
5
F =[Np*Iirrf - Np*I0rf*(exp(q*Vocrf/Ns*nrf*K*Trf)-1) - (Vocrf/(Ns/Np)*RPrf);
    Np*Iirrf - Iscrf - Np*I0rf*(exp(q*Iscrf*RSrf/Np*nrf*K*Trf) - 1) - (Iscrf*(Ns/Np)*RSrf)/((Ns/Np)*RPrf);
    Np*Iirrf - Imrf - Np*I0rf *(exp(q*(Vmrf + Imrf*(Ns/Np)*RSrf)/Ns*nrf*K*Trf) - 1) - (Vmrf + Imrf*(Ns/Np)*RSrf)/((Ns/Np)*RPrf);
    ((q*Np*I0rf/Ns*nrf*K*Trf)*exp(q*(Vmrf + Imrf*(Ns/Np)*RSrf)/Ns*nrf*K*Trf) +( 1/(Ns/Np)*RPrf))/(1 + (q*I0rf*RSrf/nrf*K*Trf)*exp(q*(Vmrf + Imrf*(Ns/Np)*RSrf)/Ns*nrf*K*Trf) + RSrf/RPrf) - (Imrf/Vmrf);
    Np*Iirrf*(G/Grf)*(1+(alphaT/Iirrf)*(T-Trf)) - Np*I0rf*((T/Trf)^3)*exp((Egref/Keg*Trf)-(Eg/Keg*T))*exp(q*((Vocrf + betaT*(T-Trf))/Ns*nrf*K*T)-1) - ((Vocrf + betaT*(T-Trf))/(Ns/Np)*((G/Grf)*RPrf))];

voici les variables et les constantes.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
%unknown parameters
Iirrf = Y(1);
I0rf = Y(2);
RSrf = Y(3);
RPrf = Y(4);
nrf = Y(5);
 
%constants
 
Vocrf = 21.6; Iscrf = 3.1; Np = 1; Ns = 36; Imrf =  2.9; Vmrf = 17.2;
q =1.6*(10^(-19)); Keg = 8.61*(10^(-5)); Trf = 298; Egref = 1.16; Grf = 1000;
 alphaT = 0.00095; betaT = -0.079; T = 298; G = 1000; K = 1.38*(10^(-23));

code de commande.

Code :

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

1
2
3
4
5
6
7
8
9
% script in the command window
Iirrf0 = 2.9999; I0rf0 = 1*10^(-9); RSrf0 = 0.002; RPrf0 = 95;  nrf0= 1.2; 
 
 
Y0 = [Iirrf0,I0rf0,RSrf0,RPrf0,nrf0]; % initial guess 
% call fsolve
 
options=optimset('Display','iter');   
[Y,fval] = fsolve(@myfun,Y0,options);