请教个用循环调用ode函数解带变量的微分方程的问题

morrissey22

问题是这样，用ode45解下面图片中这样一个微分方程组
方程组中的k从1变化到1000，就是要循环调用1000次ode函数算一千遍这样的方程组，每次方程组中的k的值都不一样
原来调用ode都是将原方程组先编写为m文件再用ode调用即可
现在这种情况我不知道程序该怎么写，知道的朋友帮个忙想想办法吧，万分感谢！！！

sigma665

最外面加个循环不可以吗

morrissey22

不知道怎么加循环量k啊
原来调用都是直接[t,x]=ode45('shg',[0,1],x0); 其中shg就是方程组的m文件名
现在方程组中有个k需要不断变化，调用ode的时候该怎么写呢？

梅梅

我记得好像有人回答过这个问题：[t,x]=ode45('shg',[0,1],x0，k); 并且shg定义时写成dydt=shg(t,y,flag,k).我写没试过，你试试看。希望对你有帮助

sogooda

原帖由 梅梅 于 2008-3-19 12:27 发表

                               
登录/注册后可看大图

我记得好像有人回答过这个问题：[t,x]=ode45('shg',[0,1],x0，k); 并且shg定义时写成dydt=shg(t,y,flag,k).我写没试过，你试试看。希望对你有帮助

差不多是这样，不过我记得应该是
[t,x]=ode45('shg',[0,1],x0,[],k);

定义shg同意4楼的。

morrissey22

楼上朋友是对的，问题已经解决，非常感谢！