Question

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

以下是功能，

主要功能： runtoy -分别使用ODE45两次解决名为 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

函数还

中存在ODE

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

在函数 yot 和玩具

中调用的函数 conn

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

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

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

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

同时从多个MATLAB文件运行ODE

1 个答案: