ode45 时定义function时,用A B作变量时为什么需要先设置A和B初始状态
% Drug absorption model - basic
k0 = 0.01;
k1 = 0.035;
[Tout, Cout] = ode45(@drugAbsorption, [0, 200], [200, 0, 0],[],k0, k1);
function Cout = drugAbsorption(t, ICs, k0, k1)
A = ICs(1);
B = ICs(2);
L = ICs(3);
dAdt = -k0*A;
dBdt = k0*A - k1*B;
dLdt = k1*B;
Cout = [dAdt; dBdt; dLdt];
end
如题,ode45 时定义function时,用A B作变量时为什么需要先设置A和B初始状态,就是code里面的A = ICs(1);B = ICs(2)这两部分,我看官方ode45里都没有这一步,官方的例子在下面。
function dydt = vdp1(t,y)dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];[t,y] = ode45(@vdp1,[0 20],[2; 0]);第一次发帖,谢谢各位解答啦!