matlab求解微分方程问题-ode45

guanchangbin

建立f1.m
function f=f1(s,p1)
f=175.534*(7251.3067*sin(s)+1.13832*sqrt(6*10^7-2*p1)*s^2-5*10^7*p1);

在命令窗口输入： [s,p1]=ode45('f1',[0 0.1],150000)

为什么求不出来啊？？？

[ 本帖最后由 ljelly 于 2009-2-11 16:31 编辑 ]

ljelly

出现什么提示？
把问题说的详细一些

guanchangbin

就是一直busy，好像是一直在计算，但是没有结果

luyang_11

function f=f1(s,p1),里面的第一个量一般都用t(时间)来表示,一定是标量形式的自变量.
我们一般求解的微分方程不止一个,所以第二个量一般用y来表示,必须是列向量.
f=175.534*(7251.3067*sin(s)+1.13832*sqrt(6*10^7-2*p1)*s^2-5*10^7*p1);而这个方程我看不懂,
我不知道为什么还有s和p1,如果有两个未知数的话必须有两个方程才能求解.
显然,我觉得你这个方程列的似乎有些问题.

shlshl

试试ode23s可以解不

babysoul

原微分方程模型是什么样子的啊？lz贴出来

guanchangbin

谢谢大家的帮助，我的问题解决了哦

ljelly

谢谢大家的帮助，我的问题解决了哦
guanchangbin 发表于 2009-3-1 18:42

解决了，能不能分享一下，什么原因，如何解决的？
谢谢！