matrix - 离散勒让德多项式的 Octave 矩阵

Question

我需要得到 N x 列 (L) 的勒让德多项式矩阵，在 L 上对任意 N 进行评估。

除了显式计算每一行的多项式向量之外，是否有更好的计算矩阵的方法？这种方法 (N = 4) 的代码片段在这里:

L = linspace(-1,1,800);

# How to do this in a better way?
G = [legendre_Pl(0,L); legendre_Pl(1,L); legendre_Pl(2,L); legendre_Pl(3,L)];

谢谢，沃伊塔

Answer 1

创建一个匿名函数。文档位于 http://www.gnu.org/software/octave/doc/interpreter/Anonymous-Functions.html

f = @(x) legendre_Pl(x,L);

然后使用 arrayfun 将函数 f 应用于数组 [1:N] 文档位于 http://www.gnu.org/software/octave/doc/interpreter/Function-Application.html

CellArray = arrayfun(f, [1:N], "UniformOutput", false);

这为您提供了一个元胞数组。如果您想要矩阵中的答案，请使用 cell2mat

G = cell2mat(CellArray);