Problème fsolve MATLAB résolution système non linéaire

[invite9b92ae3d]

Bonjour à tous !

J'essaye de résoudre un système de 4 équations à inconnues non linéaire à l'aide de fsolve de Matlab, mais j'ai toujours une erreur.
J'ai pourtant suivi l'aide en ligne de Matlab.

C'est assez urgent donc si vous avez la solution merci beaucoup !!

J'ai tout d'abord créé un fichier myfun.m comprenant mon système de 4 équations (cpc1, Tc1,etc. sont des constantes définies dans mon prog principal) :

function F = myfun(x)
F = [cpc1*(Tc1-x(1)) - cpf1*(x(2) - x(3));
    C*(x(2) - x(3)) - E*(Tc1 - x(2) - x(1) + x(3))/(log((Tc1 - x(2))/(x(1) - x(3))));
    A*(x(1) - x(3))-B*(x(4) - Tje);
    B*(x(4) - Tje) - D*(x(1) - x(4) - x(3) + Tje)/(log((x(1) - x(4))/(x(3) - Tje)))];

Et ensuite j'essaye de résoudre ça à la fin de mon programme principal :

x0 = [100; 100; 100; 100];           % Make a starting guess at the solution
options=optimset('Display','iter');   % Option to display output
[x,fval] = fsolve(@myfun,x0,options)  % Call solver

Mais j'ai le message d'erreur suivant :

??? Input argument "cpc1" is undefined.

Error in ==> myfun at 2
F = [cpc1*(Tc1-x(1)) - cpf1*(x(2) - x(3));

Error in ==> fsolve at 254
            fuser = feval(funfcn{3},x,varargin{:});

Error in ==> monprog at 44
[x,fval] = fsolve(@myfun,x0,options)  % Call solver
Caused by:
    Failure in initial user-supplied objective function evaluation. FSOLVE cannot
    continue.

Il ne semble pas reconnaître les constantes que j'ai pourtant calculées dans mon programme principal....

Merci d'avance !
-----

[invitea29b3af3]

Salut, bienvenue sur le forum.

Tu as donné toi-même la réponse à ta question:

cpc1, Tc1,etc. sont des constantes définies dans mon prog principal

La portée des variables définies dans un fichier matlab s'arrête à ce fichier. Si tu veux les utiliser dans ta fonction il faut les passer en arguments de la fonction. Oui je sais c'est embêtant mais c'est comme ça.

[invite9b92ae3d]

Merci pour ta réponse !
Que veux-tu dire par "les passer en arguments" ? Est-ce que tu pourrais me donner un exemple ? (désolé je suis novice !)

[invite79d10163]

Tu peux aussi ajouter le terme "global" devant les variables que tu souhaitent partager dans différentes fonctions.

[A voir en vidéo sur Futura]

[invitea29b3af3]

Tu peux aussi ajouter le terme "global" devant les variables que tu souhaitent partager dans différentes fonctions

Eh ben j'ai appris quelque chose moi aussi Je me disais bien que quelque chose comme ça devait exister, mais je ne connaissais pas.

Ce que j'entendais par "les passer en arguments" c'est simplement au lieu d'appeler ta fonction avec fsolve(@myfun,x0,options), tu l'appeles avec fsolve(@myfun,x0,options,cpf1, Tc1,....etc.).

[invite9b92ae3d]

OK, merci beaucoup à tous les deux !
J'ai essayé de mettre toutes mes constantes en argument et ça marche !!