Legendre、Laguerre 和 Hermite 的高斯正交节点和权重的 m 文件,推荐你看下这个打包好的资源。三个 .m 文件用的是对称伴随矩阵来搞定特征值计算,节点权重直接就出来了。嗯,方法比较稳,不用担心高阶时候根变复数那种玄学问题,挺省心。

对称矩阵带来的好处,除了稳定性,还有一点挺关键:所有特征值都是实数,你就放心大胆用。之前你用过些开源代码,跑高阶就不稳定?这个资源刚好可以替代,用起来还简单。

你要是经常在 数值积分信号 或是一些 有限元计算 场景下跑高斯正交,真建议下一个备用。支持任何 n 阶,扩展性也不错。

顺手推荐几篇相关内容,看完你会更有感觉:Hermite 权重的自定义Gramschmidt 正交矩阵构造,都挺实用。

如果你正好在做 正交多项式 方向的项目,可以把这资源收藏了,用起来会比较顺手。