Aide en Matlab

[invitec32a9991]

Bonjour ,
je veux tracer cette courbe Cs=((Ka*[A]*[B0])/(Ka*[A]+Kd))*(1-exp(-(Ka*[A]+Kd)*t)) (une courbe exponentielle)
a différentes valeurs de Ka (1000:1000:10000) de telle sorte j’obtiens une figure à plusieurs courbes, j'ai essayée avec une boucle for mais ca marche pas ( voila le code .


clear all
close all
clc
% A + B <===> AB
Cs=0;
%concentration initiale du complexe [AB]
[A]= 0.001;
%concentration initiale d'analyte[A] en [mol/m^3]
Kd= 0.05 ; %constante de dissociaton en s^-1
[B0]=4.7e-9;
%concentration initiale des recepteurs [B] en [mol/m^2]
Ka=0;
%constante d'association en [m^3/(mol*s)]

t=0:2:20 ;
% temps en [s]


for Ka =1000:1000:10000;

Cs=((Ka*[A]*[B0])/(Ka*[A]+Kd))*(1-exp(-(Ka*[A]+Kd)*t));


end



plot(t,Cs)

xlabel('temps,s');

ylabel('variable dependant Cs');
hold on
-----

[fiatlux]

Salut

(Oublie pas d'utiliser la banière CODE stp)
c'est normal car dans ta boucle tu "écrases" à chaque fois ta précédente courbe en la sauvant dans Cs !
Donc soit tu sauves toutes tes courbes dans une matrice (et non un seul vecteur Cs comme tu le fais maintenant), soit tu les affiches directement depuis la boucle.
Par exemple tu remplaces ta fin de code par:

figure, hold on; %hold on: à mettre AVANT de tracer les courbes, ou après avoir tracé la 1ère, afin que quand tu traces la 2e, la 1e ne s'efface pas du graphe
colors = jet(length(1000:1000:10000)); %juste pour avoir des courbes de différentes couleurs
for Ka =1000:1000:10000;
    Cs=((Ka*[A]*[B0])/(Ka*[A]+Kd))*(1-exp(-(Ka*[A]+Kd)*t));
    plot(t,Cs,'color',colors(Ka/1000,:))
end
xlabel('temps,s');
ylabel('variable dependant Cs');

[invitec32a9991]

merci beaucoup pour votre réponse et aussi votre conseil

[fiatlux]

de rien

[A voir en vidéo sur Futura]