Tracer une fonction sous Matlab

Résolu/Fermé
atiffa Messages postés 54 Date d'inscription lundi 4 février 2008 Statut Membre Dernière intervention 10 juin 2009 - 27 mars 2008 à 15:22
Sacabouffe Messages postés 9427 Date d'inscription dimanche 19 août 2007 Statut Membre Dernière intervention 29 mai 2009 - 27 mars 2009 à 16:26
Bonjour,

je n'arrive pas a tracer la fonction (NaKpump) suivante:

function [JNa, JK, JNH4, Ja] = NaKpump(Ja_max, Cnac, Cks, CNH4s, KNa, Kk, KNH4)


Ja = Ja_max*1./(1+(KNa./Cnac)^3);

JNa = -3.*Ja;

JNH4 = 2*Ja./(((Cks/Kk)./(CNH4s/KNH4))+1);

Jk = 2*Ja*(1-1./(((Cks/Kk)./(CNH4s/KNH4))+1));

end

Je veux tracer JNa en fonction de différentes valeurs de Cnac

pourriez vous m'aider?
merci
A voir également:

3 réponses

Sacabouffe Messages postés 9427 Date d'inscription dimanche 19 août 2007 Statut Membre Dernière intervention 29 mai 2009 1 832
27 mars 2009 à 15:31
Salut

Je sais que le sujet date d'un mois mais au cas où t'aurais pas trouvé et/ou quelqu'un aurait un souci de ce genre.
Une des solutions, plutôt que de faire une boucle pour définir un vecteur des valeurs de JNa en fonction de valeurs de Cnac en appelant la fonction plusieurs fois, c'est de définir ta fonction comme agissant sur un vecteur pour Cnac. Donc comme d'hab, faut rajouter des . pour les opérations élémentaires quand elles agissent sur un vecteur. Je suppose que toutes les autres variables sont scalaires.
function [JNa, JK, JNH4, Ja] = NaKpump(Ja_max, Cnac, Cks, CNH4s, KNa, Kk, KNH4)
Ja = Ja_max./(1+(KNa./Cnac).^3);
JNa = -3*Ja;
JNH4 = 2*Ja/(((Cks/Kk)/(CNH4s/KNH4))+1);
Jk = 2*Ja*(1-1/(((Cks/Kk)/(CNH4s/KNH4))+1));
end
Les variables de sortie JNa, JK, JNH4, Ja sont alors vectorielles.
plot(Cnac,JNa);
Et les carottes sont cuites ! ;-)

Ciao
1
atiffa Messages postés 54 Date d'inscription lundi 4 février 2008 Statut Membre Dernière intervention 10 juin 2009 12
27 mars 2009 à 16:22
:-) Merci sacabouffe de ta réponse, mais je me suis débrouiller pour le résoudre.
Par contre est ce que tu connais la méthode de newton Raphson, j'ai besoin de savoir comment on utilise cette méthode pour trouver les racines d'une fonctions . (j'ai pas compris grand chose sur internet)

Merci infiniment sacabouffe

((Et les carottes sont cuites ! ;-) hum ça l'air so delicious ;-) hiiiiiiiiiiiiii


Merci encore
0
Sacabouffe Messages postés 9427 Date d'inscription dimanche 19 août 2007 Statut Membre Dernière intervention 29 mai 2009 1 832
27 mars 2009 à 16:26
:-D
Je m'en doutais un peu mais bon, ça pourra peut-être resservir...
'Fin bon... y a pas de quoi du coup ! ;-)

J'ai déjà répondu dans ton autre sujet pour la méthode de Newton-Raphson ;-)
http://www.commentcamarche.net/forum/affich 11720941 matlab methode de newton raphson
0