wolfram-mathematica - 在 Mathematica 中找到微分方程解的零点

Question

给定以下代码：

s := NDSolve[{x''[t] == -x[t], x[0] == 1, x'[0] == 1}, x, {t, 0, 5 }]
Plot[Evaluate[{x[t]} /. s], {t, 0, 3}]

这绘制了微分方程的解。我将如何数值求解 x[t] 的零，其中 t 介于 0 和 3 之间？

Answer 1

@rcollyer 回答了最初的问题。我正在回答您在对 rcollyer 的回答的第一条评论中发布的问题：

但是如果我们的 s 是 "s := NDSolve[{x'[t]^2 == -x[t]^3 - x[t] + 1, x[0] == 0.5}, x, { t, 0, 5}]" 然后 FindRoot 函数只返回一个错误，而绘图显示在 0.6 左右有一个零。

所以：

s = NDSolve[{x'[t]^2 == -x[t]^3 - x[t] + 1, x[0] == 0.5}, 
             x, {t, 0, 1}, Method -> "StiffnessSwitching"];
Plot[Evaluate[{x[t]} /. s], {t, 0, 1}]
FindRoot[x[t] /. s[[1]], {t, 0, 1}]

{t -> 0.60527}

编辑

回答 rcollyer 的评论，“第二行”来自平方导数，如：

s = NDSolve[{x'[t]^2 == Sin[t], x[0] == 0.5}, x[t], {t, 0, Pi}];
Plot[Evaluate[{x[t]} /. s], {t, 0, Pi}]

来自（哪里：

DSolve[{x'[t]^2 == Sin[t]}, x[t], t]
(*
{{x[t] -> C[1] - 2 EllipticE[1/2 (Pi/2 - t), 2]}, 
 {x[t] -> C[1] + 2 EllipticE[1/2 (Pi/2 - t), 2]}}
*)

Answer 2

FindRoot作品

In[1]:=  FindRoot[x[t] /. s, {t, 0, 3}]
Out[1]:= {t -> 2.35619}