trying to plot a phase portrait ODE

Question

Jguyson 2016-12-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/317344-trying-to-plot-a-phase-portrait-ode

回答： 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?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Sally Al Khamees 2016-12-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/317344-trying-to-plot-a-phase-portrait-ode#answer_248211

在 MATLAB Online 中打开

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