计算相关勒让德多项式的精确导数可能很棘手。 即使在高级文本中,它们通常也被写成递归关系和/或包含阶乘的(归一化)因子。 因此,一个简单的软件实现将很快遇到 IEEE754 双精度的限制,导致 NaN/inf 或显着的精度损失,已经在非常低的 N 度。 LEGENDRE_DERIVATIVE 是导数计算的完全矢量化、数值稳定且经过稳健验证的实现。 它允许快速准确地计算任何 N 次的导数。
它的工作原理与 MATLAB 自己的 LEGENDRE 相同,不同之处在于它不计算多项式值,而是计算导数的值。
2021-09-08 15:43:14
6KB
matlab
1