matlab编程题-三点法和五点法Gauss-Legendre

Gauss-Legendre积分是一种数值积分方法，用于近似一个函数在一定区间内的积分。这种方法特别适用于当函数的解析积分很难或无法直接求得时。Gauss-Legendre积分使用特定的点（称为Gauss点）和相应的权重来近似积分，这些点和权重是通过Legendre多项式的根和相应的积分公式计算得出的。

三点Gauss-Legendre积分

在三点Gauss-Legendre积分中，使用三个点和三个相应的权重来近似区间 [-1, 1] 上的积分。这些点是Legendre多项式（在这种情况下是三阶多项式）的根，而权重则是用于确保最高可能的精确度。对于三点Gauss-Legendre积分，点和权重分别为：

• 点:
• 权重: w1​=w3​=5/9​,w2​=8/9​

五点Gauss-Legendre积分

五点Gauss-Legendre积分使用五个点和相应的权重来进行类似的操作，但这次是基于五阶Legendre多项式。这提供了更高的精度。

应用到积分 

要应用这些公式计算给定区间 [0, 1] 的积分，需要进行一个线性变换，将积分区间从 [-1, 1] 映射到 [0, 1]。这可以通过一个简单的变换实现，例如对于任何在 [-1, 1] 区间的点 x，变换为 ​ 使之适应 [0, 1] 区间。 

 MATLAB代码

% 定义被积函数
f = @(x) (x .* exp(x) ./ (1 + x).^2);

% 三点Gauss-Legendre积分
x3 = [-sqrt(3/5), 0, sqrt(3/5)];  % 三点的高斯点
w3 = [5/9, 8/9, 5/9];            % 三点的权重
I3 = 0.5 * sum(w3 .* f(0.5 * x3 + 0.5)); % 计算积分

% 五点Gauss-Legendre积分
x5 = [-sqrt(5 + 2*sqrt(10/7)) / 3, -sqrt(5 - 2*sqrt(10/7)) / 3, 0, sqrt(5 - 2*sqrt(10/7)) / 3, sqrt(5 + 2*sqrt(10/7)) / 3]; % 五点的高斯点
w5 = [(322 - 13*sqrt(70)) / 900, (322 + 13*sqrt(70)) / 900, 128/225, (322 + 13*sqrt(70)) / 900, (322 - 13*sqrt(70)) / 900]; % 五点的权重
I5 = 0.5 * sum(w5 .* f(0.5 * x5 + 0.5)); % 计算积分

% 显示结果
fprintf('三点Gauss-Legendre积分结果: %f\n', I3);
fprintf('五点Gauss-Legendre积分结果: %f\n', I5);

运行结果