trying to plot a phase portrait ODE

8 次查看（过去 30 天）
Jguyson
Jguyson 2016-12-18
回答： Sally Al Khamees 2016-12-22
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?

请先登录，再回答此问题。

回答（1 个）

Sally Al Khamees
Sally Al Khamees 2016-12-22
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

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by