【数值分析】区间折半法，matlab实现

最新推荐文章于 2026-08-21 12:36:00 发布

原创最新推荐文章于 2026-08-21 12:36:00 发布 · 818 阅读

7 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

收录于

数值分析

当前文章被收录于：

数值分析

21 篇文章 19 人学习

订阅专栏查看详情

当前文章被以下社区和专栏收录：

文章介绍了如何使用区间折半法求解定积分问题，通过逐步减小步长并计算梯形和辛普森法则的近似值，直至达到预设的精度。Matlab代码示例展示了该方法的具体应用。

区间折半法

从梯形公式出发，上一步步长为 ${h}$ ，则有步长折半后的积分
$T_{2n}= \frac{1}{2}T_n+ \frac{h}{2} \sum_{i=0}^{ n-1}f(x_{i+0.5})$
matlab实现

%% 区间折半法例子
format long
[I i] = halfStep(@f,0,1,1e-7,10)

function r = f(x)
    if x == 0
        r = 1;
    else
        r = sin(x)./x;
    end
end

%% 区间折半法
% 输入函数，范围，精度，最大迭代次数
% 输出积分值，迭代次数
function [I,i] = halfStep(f,a,b,eps,max_iter)
    h = b-a;
    T0 = h/2*(f(a)+f(b))
    for i = 1:max_iter
        x_n = a+h/2:h:b-h/2;
        T1 = 0.5*T0+0.5*h*sum(f(x_n))
        if abs(T1-T0)<eps
            I = T1;
            break
        end
        T0 = T1;
        h = h/2;
    end
end