针对特定正弦函数的最快求根算法咨询

先直接回应你的问题：首先得明确，你的函数 f(x) = sin²(πx) + sin²(π·323/x) 是两个平方项的和，所以它的值永远非负——这就是为什么你在隔离的区间里看不到函数值变号，二分法这类依赖符号变化的算法自然用不了。

接下来拆解根的本质：f(x) = 0 的充要条件是两个平方项各自为0，也就是：

• sin(πx) = 0，这意味着 x 必须是整数（x = k, k ∈ ℤ）
• sin(π·323/x) = 0，这意味着 323/x 必须是整数（也就是 x 是323的整数因数）

所以 x 得同时满足这两个条件：既是整数，又是323的因数。先对323做质因数分解：323 = 17 × 19，它的所有整数因数是：±1, ±17, ±19, ±323。你说排除1和323，那剩下的非平凡根就是 ±17, ±19。

如果你的目标就是找这些精确根，那完全不需要数值迭代算法——直接分解因数是最快的方法，比任何数值算法都高效。

不过我猜你可能是不是混淆了“根”和“极小值点”？比如你可能想找函数值接近0的点？如果是这样的话，因为函数在这些非根的区间里是正的，你可以用以下数值算法：

• 牛顿法：如果能算出函数的导数，牛顿法收敛速度很快，不过需要初始值足够接近极小值点。你的函数导数化简后是 f’(x) = π sin(2πx) - (646π / x²) sin(2π·323/x)，只要你有一个合适的初始值（比如你隔离的区间里的某个点），牛顿法能快速收敛到极小值点。
• 黄金分割法：如果你不确定导数，或者不想算导数，黄金分割法是个不错的选择——它是单峰区间上找极小值的经典方法，不需要导数，只需要不断缩小区间，收敛速度虽然比牛顿法慢，但胜在稳定。
• 布伦特法：这是黄金分割法和抛物线插值的结合，收敛速度比黄金分割法快，同样不需要导数，也是很常用的无导数优化算法。

最后再敲个黑板：如果真的是找f(x)=0的根，因数分解秒杀一切数值方法；如果是找极小值，就根据自己的情况选上面的迭代算法就行。

备注：内容来源于stack exchange，提问作者Canbach