我正在尝试查找[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]中得到该函数的零,而不必通过这种漫长而幼稚的方法?