0

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

以下是功能，

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

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

功能玩具中出现的第一组 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

存在于函数yot中的 ODE

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
4

1 回答 1

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
于 2018-10-30T12:21:43.220 回答