songqiqqq 发表于 2012-8-22 22:18

关于ode45中方程无法输入的困难

最近在模拟一个实验，最近归结为要解一个方程组，我就打算用ode45函数来解，可以在方程输入的时候遇到了困难。具体讲，方程的形式是这样的dx/dt=int f(x,theta,u)du,因为目前关于theta我只能得出一个数值解（也就是对于u是离散的取值），所以我打算把定积分int近似为求和来解。所以问题变为dx/dt=xigema: f(x,theta,u)du,可是这样一来，左边就有几千项公式要我输入。
    我本采用匿名函数，ode45('odefun',[],[]),odefun=@(t,y)........    可是这样无法输入啊，急啊，怎么把方程右边一大串弄进去

louislwm 发表于 2012-8-24 00:14

右边单独编一个m文件，
function dy=odefun(t,y)
    dy=......
然后再用ode45('odefun',[],[])不行么？

songqiqqq 发表于 2012-8-27 10:20

回复 2 # louislwm 的帖子

谢谢，我已经用你这个方法解决了
页: [1]
查看完整版本: 关于ode45中方程无法输入的困难