我正在尝试计算数组的二阶导数,尽管在使用样条拟合和微分时观察噪声,或者在使用 np.gradient 直接微分时观察噪声。
通过阅读有关此主题的信息,我相信龙格现象可能会导致这种噪音。一种建议的解决方案是不等的网格间距——例如使用 Chebyshev 节点。
我在 Python 中看不到任何这样的函数,所以想知道是否有人在 Python 中有过这种数学实现的经验?
此外,直接计算我的值数组以在这个新的 Chebyshev 网格上进行微分更好,还是我可以在微分之前简单地插值到新网格上以提高准确性?
任何提示或参考都非常受欢迎,谢谢。