mathematica 解微分方程组 30
求高手指点:解微分方程组equ1和equ2,公式没写错,个人认为程序应该也没错,但每次计算都出错,并且提示:NDSolve::ntdvdae:Cannotsolvetof...
求高手指点:解微分方程组equ1和equ2,公式没写错,个人认为程序应该也没错,但每次计算都出错,并且提示:NDSolve::ntdvdae: Cannot solve to find an explicit formula for the derivatives. NDSolve will try solving the system as differential-algebraic equations. >>和NDSolve::ndsv: Cannot find starting value for the variable H. >>实在不晓得该怎么办了,求高手指教。(下面是我写的程序)
H0 = 0.000003;
\[Phi]0 = 10;
\[Phi]f = 0.000741;
equ1 = 2/3*\[Phi]'[t]^3 -
Sqrt[4*\[Phi]'[t]^6/9 + \[Phi]'[t]^4/3 -
Exp[2*\[Phi][t]]*\[Phi]'[t]^2/3] == H[t];
equ2 = 4*\[Phi]'[t]^2*\[Phi]''[t] - 3*Exp[2*\[Phi][t]]*\[Phi]'[t]^2 +
2*\[Phi]'[t]^4 + 4*\[Phi]'[t]^6 -
6*\[Phi]'[t]^3*(2*\[Phi]'[t]^3/3 -
H[t]) - (4*\[Phi]'[t]^5*\[Phi]''[t] + 3*H[t]*H'[t] -
6*\[Phi]'[t]^2*\[Phi]''[t]*H[t] - 2*\[Phi]'[t]^3*H'[t])/
Sqrt[2*\[Phi]'[t]^3/3 - H[t]] == 0;
sol = NDSolve[{equ1, equ2, H[t0] == H0, \[Phi][t0] == \[Phi]0, \[Phi]'[t0] == [Phi]f}, {H, \[Phi]}, {t, -10, 10} ]
大仙们有没运行一下?其实我觉得语法没错啊,为毛就不能算? 展开
H0 = 0.000003;
\[Phi]0 = 10;
\[Phi]f = 0.000741;
equ1 = 2/3*\[Phi]'[t]^3 -
Sqrt[4*\[Phi]'[t]^6/9 + \[Phi]'[t]^4/3 -
Exp[2*\[Phi][t]]*\[Phi]'[t]^2/3] == H[t];
equ2 = 4*\[Phi]'[t]^2*\[Phi]''[t] - 3*Exp[2*\[Phi][t]]*\[Phi]'[t]^2 +
2*\[Phi]'[t]^4 + 4*\[Phi]'[t]^6 -
6*\[Phi]'[t]^3*(2*\[Phi]'[t]^3/3 -
H[t]) - (4*\[Phi]'[t]^5*\[Phi]''[t] + 3*H[t]*H'[t] -
6*\[Phi]'[t]^2*\[Phi]''[t]*H[t] - 2*\[Phi]'[t]^3*H'[t])/
Sqrt[2*\[Phi]'[t]^3/3 - H[t]] == 0;
sol = NDSolve[{equ1, equ2, H[t0] == H0, \[Phi][t0] == \[Phi]0, \[Phi]'[t0] == [Phi]f}, {H, \[Phi]}, {t, -10, 10} ]
大仙们有没运行一下?其实我觉得语法没错啊,为毛就不能算? 展开
2个回答
展开全部
别的不说,你t0 是 t=0 还是 t=-10.。。
而且 ,不是所有的微分方程都有解。看看初值条件是不是给少了。
而且 ,不是所有的微分方程都有解。看看初值条件是不是给少了。
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
Sievers分析仪
2024-10-13 广告
2024-10-13 广告
是的。传统上,对于符合要求的内毒素检测,最终用户必须从标准内毒素库存瓶中构建至少一式两份三点标准曲线;必须有重复的阴性控制;每个样品和PPC必须一式两份。有了Sievers Eclipse内毒素检测仪,这些步骤可以通过使用预嵌入的内毒素标准...
点击进入详情页
本回答由Sievers分析仪提供
展开全部
好像没有。你可以转写成y',y''型的再用NSolve。当然,也可以用函数D表示导数(包括偏导数),我想dy,dx形式的转写成用函数D表示要比转写成y',y''型方便。函数D可以在帮助查到。
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询