编写拉盖尔多项式相关MATLAB代码遇错，求系数代码修正

1. 先修正拉盖尔多项式函数的返回值问题 #

你的Laguerre_Polynomial函数定义了输出参数M，但最后没有将计算结果赋值给它，导致调用时无法正确获取多项式数组。同时优化代码结构，预分配符号数组提升效率：

function M = Laguerre_Polynomial(N)
    syms t k
    L = sym(zeros(N+1, 1)); % 预分配符号数组
    for i = 0:N
        x = factorial(i);
        % 直接计算第i阶拉盖尔多项式
        L(i+1) = symsum( ((-1)^k * x)/(factorial(k)^2 * factorial(i-k)) * t^k, k, 0, i );
    end
    M = transpose(L); % 将结果赋值给输出参数M
end

2. 修正系数计算函数的语法错误与逻辑 #

原Coeff函数的核心错误是积分函数的调用格式错误，同时需要明确拉盖尔多项式的正交内积计算规则（权重为exp(-t)）。修改后的代码：

function D = Coeff(N)
    syms t
    f = t^3 + t^2 + 1;
    Lag = Laguerre_Polynomial(N); % 正确调用拉盖尔多项式函数
    C = sym(zeros(N+1, 1)); % 预分配系数数组
    for j = 0:N
        % 计算正交系数：内积<f, L_j> = ∫₀^∞ f(t)L_j(t)exp(-t)dt
        integrand = f * Lag(j+1) * exp(-t);
        C(j+1) = int(integrand, t, 0, inf);
    end
    D = C; % 赋值给输出参数D
end

错误点说明 #

• 原Coeff函数第7行的语法错误：你错误拆分了被积函数，int(f*exp(-t))*Lag(j+1)不符合MATLAB的int函数调用规则，正确的做法是将整个内积项作为int的第一个参数，同时指定积分变量t和上下限。
• 原Laguerre_Polynomial函数未返回输出参数：函数定义了M作为输出，但最后没有将计算结果赋值给它，导致调用时无法获取有效数据。

内容的提问来源于stack exchange，提问作者Natasha Azmat