Effacer les filtres
Effacer les filtres

trying to plot a phase portrait ODE

7 vues (au cours des 30 derniers jours)
Jguyson
Jguyson le 18 Déc 2016
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?

Réponses (1)

Sally Al Khamees
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

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by