Discussion :

[aziz1015]

Bonjour, j'ai besoin de votre aide afin de modéliser corrrectement les nullclines du modèle de Morris-Lecar :
Vous remarquerez les différences entre les deux images, j'obtiens l'image avec les graphiques en couleur sur octave qui est très loin du résultat comme vous pouvez le voir.

Le modèle de Morris-Lecar est le suivant :



J'obtiens les fonctions des nullclines en annulant les dérivées, ce qui me donne :



d'où le code en octave :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
g_L=2;
g_K=8;
g_Ca=4;
V_L=-50;
V_K=-70;
V_Ca=100;
V_1=10;
V_2=-15;
V_3=-1;
V_4=14.5;
 
u=linspace(-0.2,1,10);
 
h=linspace(-50,50,10);
N_inf=0.5*(1+tanh((h-V_3)/V_4));
M_inf=0.5*(1+tanh((h-V_1)/V_2));
 
 
I = [100 200 300 400 500];
for i=1:5
   for j=1:10
   	 V(j,i) = (I(i)+g_Ca*V_Ca*M_inf(j)+g_K*V_K*u(j))./(g_L+g_Ca*M_inf(j)+g_K*u(j));
   end
end
 
figure;
hold on;
plot(V(:,1),u,"b",'LineWidth',1.5);
plot(V(:,2),u,'r','LineWidth',1.5);
plot(V(:,3),u,'g','LineWidth',1.5);
plot(V(:,4),u,'c','LineWidth',1.5);
plot(V(:,5),u,'m','LineWidth',1.5);
plot(h,N_inf,"k",'LineWidth',2);
hold off;
axis([-50 50 0 1]);
 
 
legend("I=100","I=200","I=300","I=400","I=500")
xlabel("V(MV)");
ylabel("N");