Matrix 离散勒让德多项式的倍频程矩阵

我需要得到任意N的勒让德多项式在L上的nx列（L）矩阵

有没有更好的方法来计算矩阵，而不仅仅是显式地计算每行的多项式向量？此方法（N=4）的代码段如下：

谢谢，

---

Vojta创建一个匿名函数。文件在

然后使用arrayfun将函数，

f

应用于位于

这将为您提供一个单元格数组。如果您想在矩阵中找到答案，请使用

cell2mat

G = cell2mat(CellArray);

---

首先，您如何定义“最佳”以及第二，是否可以接受任何不准确之处？例如，你只需在每一秒点上计算L并进行插值，我只需要在1到n（n是结果矩阵的行数）之间找到L上legendre_Pl（n，L）的行，并将它们作为Nxcols（L）的矩阵。使用n=1:n作为legendre_Pl的参数在这里不起作用。我所说的“最佳”是指生成N个变量，因此我不必手动键入所有legendre_Pl。此外，最好不要使用for循环。不准确不是这个问题的主题

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

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

G = cell2mat(CellArray);