MATLAB实现变步长梯形和变步辛普森求积法

在MATLAB中实现变步长梯形法和变步长辛普森求积法可以使用以下步骤。我们将使用函数f(x) = x^2，x范围从0到1，作为示例。
首先，定义函数f(x)和x的范围：

f = @(x) x.^2;  % 定义函数f(x) = x^2
x = 0:0.01:1;  % 定义x的范围从0到1，步长为0.01

接下来，使用变步长梯形法进行数值积分：

T_approx = zeros(size(x));  % 初始化近似值数组
% 循环计算每个区间的梯形面积
for i = 1:(length(x)-1)
h = x(i+1) - x(i);  % 计算当前区间的步长
a = f(x(i));  % 计算当前区间左端点的函数值
b = f(x(i+1));  % 计算当前区间右端点的函数值
T_approx(i) = h * (a + b) / 2;  % 计算梯形面积并存储到近似值数组中
end
% 将最后一个区间处理为特殊情况
T_approx(end) = T_approx(end-1);  % 处理最后一个区间
T_approx = trapz(x, T_approx);  % 使用trapz函数进行数值积分，得到最终的近似值

现在我们可以使用同样的方法来实现变步长辛普森求积法：

S_approx = zeros(size(x));  % 初始化近似值数组
% 循环计算每个区间的辛普森面积
for i = 1:(length(x)-1)
h = x(i+1) - x(i);  % 计算当前区间的步长
a = f(x(i));  % 计算当前区间左端点的函数值
b = f(x(i+1));  % 计算当前区间右端点的函数值
c = f((x(i) + x(i+1)) / 2);  % 计算区间中点的函数值
S_approx(i) = h * (a + b + c) / 3;  % 计算辛普森面积并存储到近似值数组中
end
% 将最后一个区间处理为特殊情况
S_approx(end) = S_approx(end-1);  % 处理最后一个区间
S_approx = simpson(x, S_approx);  % 使用simpson函数进行数值积分，得到最终的近似值

以上代码演示了如何使用MATLAB实现变步长梯形法和变步长辛普森求积法。请注意，这里使用了trapz和simpson函数来进行数值积分，它们分别对应于变步长梯形法和变步长辛普森求积法的数值积分公式。通过循环计算每个区间的梯形或辛普森面积，并将它们累积起来，我们可以得到最终的近似值。