同时从多个 MATLAB 文件运行 ODE

Question

我在两个 .m 文件中的两个函数中创建了 ode。我正在尝试结合这两个函数的 ode 并同时运行它。

以下是函数，

主要函数：runtoy - 使用 ODE45 两次分别解决名为 toy 和 yot 的两个函数中的 ode。

function runtoy()
%% toy
Nnodes = 5;
vc = 50;
vz0 = [280 100 150 180 110];
vz = vz0(2:4)';
tspan = [0 2];
[tz,dvzdt] = ode45(@(t,y) toy(t,y,vc), tspan,vz)
plot(tz,dvzdt)

hold all
plot(tz,vc,'*')
%% yot

[tc,dvcdt] = ode45(@(t,z) yot(t,z,vz(2)), tspan,vc)
hold all
plot(tc,dvcdt)
hold all
plot(tc,vz(2),'o')
legend({'toy-vz1','toy-vz2','toy-vz3','toy-vc','yot-vc','yot-vz2'})

%%EDIT

init = [vz' vc]; 
Combined = ode45(@(t,v) combined_ODEs(t,v), tspan ,init')
end

函数toy

中的第一组ODE

function dvz = toy(t,y,vc)
vz = y;
Matrix = [-2 0 0;1 -3 1;0 -4 1];
Connect = [0;conn(vz(1),vc);0]; 
dvz = Matrix*vz -Connect;
end

ODE 存在于函数 yot

function dvc=yot(t,z,vz)
dvc = conn(vz,z);
end

函数 yot 和 toy

中调用的函数 conn

function c= conn(vz,vc)
c = vz - vc/3;
end

有人可以就如何同时解决两个函数中的 ode 问题提出建议吗？我可以在一个文件中拥有这些功能。但是，我将其用作测试用例来实现我的真实系统。

编辑：评论中建议的方法试用

function dV = combined_ODEs(t,V)
vz = V(1:3);
vc = V(4);
dV(1:3,1) = toy(t,vz,vc);
dV(4,1) = yot(t,vc,vz(2)); % not sure what z is here??
end

Answer 1

由于这两个 ODE 是耦合的，我的建议是编写（但）另一个函数来引用定义 ODE 的两个（或更多）函数并解决这个问题，即类似的东西：

function dV = combined_ODEs(t,V)
vz = V(1);
vc = V(2);
dV(1) = toy(t,vz,vc);
dV(2) = yot(t,z,vz); % not sure what z is here??
end