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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

最近在模拟一个实验，最近归结为要解一个方程组，我就打算用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 的帖子

谢谢，我已经用你这个方法解决了

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 关于ode45中方程无法输入的困难

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。