I have the following function:
function ypdiff = ypdiff(t,y)
a = 0.01;
b = 0.1;
ypdiff(1) = -a*y(1)*y(2);
ypdiff(2) = b*y(1)*y(2)-b*y(2);
ypdiff(3) = b*y(2);
ypdiff = [ypdiff(1) ypdiff(2) ypdiff(3)]';
如果我想解决这个问题,我会调用ode45函数,如下所示:
[t y] = ode45(@ypdiff, [to tf], yo);
但是如果我想向这个函数传递参数,我该如何使用ode45函数?具体来说,我尝试以下公式:
function ypdiff = ypdiff(t,y,u)
a = 0.01;
b = 0.1;
ypdiff(1) = -a*u*y(1)*y(2);
ypdiff(2) = b*u*y(1)*y(2)-b*y(2);
ypdiff(3) = b*u*y(2);
ypdiff = [ypdiff(1) ypdiff(2) ypdiff(3)]';
ode45(@(t, y) ypdiff(t, y, u)
,而不是只使用ode45(@ypdiff(u))
。为什么需要t
和y
? - Ka Wa Yip