[問題]fitzhugh-nagumo 模型 以I為參數的繪圖
大家好 最近再重做前人的模擬 遇到苦思不解的麻煩
fitzhugh-nagumo 模型 (神經元活動的數學模型)
想要畫出 I(x軸)-V(y軸) 大部分的文章都是 t-V圖
最後找到幾乎接近的 code 還是改不出來 XD
希望板上的大大們能指引1下 小弟程式不太行 :(
我想要模擬的 fitzhugh-nagumo 模型 I(電流)=0.05+(0.0005*t)^p p=1 or 2 or
0.5
(steven baer 2008的文章 Slow acceleration and deacceleration through a Hopf
bifurcation: Power ramps,
target nucleation, and elliptic bursting. 中的 fig 1& 2) 畫出其 I-V圖.
幾乎可行的code (參數已改為要模擬的)
global I; % Need this so that function "f" knows about variable I
for I = [0.05:0.001:1]; % Loop over different values of applied current I
f = @(t,y) [-y(1)*(y(1)-0.2)*(y(1)-1)- y(2) + I; 0.05*(y(1) -0.4*y(2)) ];
g = @(y) f(0,y);
fp = fsolve(g,[0 0]); % Find the fixed point
Vss = fp(1); Wss = fp(2); % Get the steady-state V and W values from "fp"
hold on;
[T Y] = ode45(f, [0 100], Y(length(T),:) ); % Run some more in order
[T Y] = ode45(f, [0 100], Y(length(T),:) ); % to settle at equilibrium
plot(I, min(Y(:,1)), 'm.'); % minimum of V(t): a magenta point
plot(I, max(Y(:,1)), 'm.'); % maximum of V(t): a magenta point
plot(I,Vss,'.','color',color); % Fixed point: a red or blue point
title('Bifurcation diagram (V_{min}, V_{max}, V*)', 'fontsize', 12);
xlabel('I'); ylabel('V equilibrium');
axis tight;
drawnow;
end;
因為試過迴圈 改I=y(3) 副程式等方式還是試不出來 希望有大大能幫忙 感謝 :)
能完整執行的code 在連結中 我的列子 要先執行過原檔 (我只抓右下角的IT圖)
http://web.njit.edu/~matveev/UBM/matlab.html 的 FN_Hopf.m
(Matlab code for the Fitzhugh-Nagumo model, showing the solution with initial
condition close to the fixed point, as the applied current I is varied.
Subthreshold Hopf bifurcation arises as the fixed point loses stability and
gives way to periodic spiking solution.)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 120.107.155.218
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
