我试图在 [0,1] 范围内找到函数的零点。我的函数通常看起来像这样:
F(t)=2*sin(8*Pi*t)+2*sin(2*Pi*t)-1
到目前为止,我尝试过的天真的方法是将函数的非常粗略的图绘制到 Pari,plot(t=0,1,F(t))
中,然后从那里“猜测”一个小范围 [a1,a2] 第一个零的位置。然后我一直在使用solve(t=a1,a2,F(t))
找到那个零。然后找到第二个小范围 [a2,a3],并重复,直到我找到所有零。
显然不是一个好方法,但它完成了工作,到目前为止我看过的许多例子只有几个零。但是,我现在正在查看的示例(例如我在上面定义的 F(t))有更多的零,所以我正在寻找一种更快的方法。
我知道solvestep
,但我尝试过使用它,但它永远不会返回任何解决方案。我猜这种类型的函数并不适合它使用的拆分方法(或者可能我用错了..)。
有没有什么好方法可以让我在 [0,1] 中返回这个函数的零,而不必经历这个漫长而幼稚的方法?