trying to plot a phase portrait ODE
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hi
I'm attempting to plot a phase portrait of a neuron model. The equation is;
C v' = k(v−vr)(v−vt)−u + I
u' = a{b(v−vr)−u}
if v ≥ vpeak, then
v ← c, u ← u + d
How would you go about doing this?
0 commentaires
Réponses (1)
Sally Al Khamees
le 22 Déc 2016
Phase Portraits are created by generating Phase Plots for different starting points. Phase plot can be created using ODEPHAS2 or ODEPHAS3 in matlab as a value for 'OUTPUTFCN' property in ODESET while creating options object. Following is an example.
function pdecaller()
options = odeset('OutputFcn',@odephas2);
for i=-4:4
for j= -4:4
[t,y] = ode45(@vdp1,[0 20],[i;j],options); hold on;
end
end
end
function dydt = vdp1(t,y)
dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];
end
0 commentaires
Voir également
Catégories
En savoir plus sur Ordinary Differential Equations dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!