用MATLAB求解一个ODEs系统？

Question

这是我第一次发帖，所以如果你需要更多的信息或者我没有做正确的事情，请让我知道！

F=@(t,x) [
.the RHS of my first order ODES.
];


[t x]=ode45(F,[Range], [Initial conditions]);

其中f(a)可以是另一个基于常数的函数a。

我已经按照以下顺序设置了我的代码；

常数，然后是f(a)形式的方程，然后是我的微分方程，后面是ode45求解器和plot命令。然而，我得到了几个错误“未定义的函数或变量”，因为早期的方程依赖于变量/方程，直到后来才被定义。

非常感谢您的帮助：）

正如建议的，这里是我的代码类型的示例；

`%Constants
a=34
b=31
c=20
%Equations
A=b*cos(2*pi)
B=a*EQ1
C=c*x(2)
%DifferentialEquations
EQ1=x(1)*A
EQ2=(EQ3-EQ1)*(B-C)
EQ3=x(2)*x(3)
F=@(t,x) [EQ1;EQ2;EQ3;];[t x]=ode45(F,[0 10], [0 0 0 ]);

Answer 1

根据注释，您试图在定义变量之前使用它。在MATLAB中，你需要定义你的变量和函数，然后再使用它们。因为方程是彼此的函数，所以您可以使用匿名函数来定义f。

示例：

a=34;
b=31;
c=20;
A = b*cos(2*pi);
EQ1 = @(x) x(1)*A;
B = @(x) a*EQ1(x);
C = @(x) c*x(2);
EQ3 = @(x) x(2)*x(3);
EQ2 = @(x) (EQ3(x) - EQ1(x))*(B(x)-C(x));
F = @(t,x) [EQ1(x);EQ2(x);EQ3(x)];
[t,x] = ode45(F,[0 10], [0 0 0]);