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

区间折半法

从梯形公式出发,上一步步长为 h {h} h ,则有步长折半后的积分
T 2 n = 1 2 T n + h 2 ∑ i = 0 n − 1 f ( x i + 0.5 ) T_{2n}= \frac{1}{2}T_n+ \frac{h}{2} \sum_{i=0}^{ n-1}f(x_{i+0.5}) T2n=21Tn+2hi=0∑n−1f(xi+0.5)

matlab实现

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
相关推荐
yongui478341 天前
基于MATLAB的8QAM调制解调仿真与BER性能分析
开发语言·matlab
电力程序小学童1 天前
【预测】基于CNN-BiLSTM-Attention的光伏预测模型
matlab·预测·光伏预测
机器学习之心1 天前
MATLAB基于加速遗传算法投影寻踪模型的企业可持续发展能力评价研究
matlab·加速遗传算法·投影寻踪模型·企业可持续发展能力评价
Matlab仿真实验室2 天前
基于Matlab实现图像栅格化处理
图像处理·计算机视觉·matlab·图像栅格化处理
jllllyuz2 天前
matlab裂纹检测与延展分析系统
人工智能·计算机视觉·matlab
jghhh012 天前
针对大尺度L1范数优化问题的MATLAB工具箱推荐与实现
开发语言·算法·matlab
不枯石2 天前
Matlab通过GUI实现点云的统计滤波(附最简版)
开发语言·图像处理·算法·计算机视觉·matlab
不枯石2 天前
Matlab通过GUI实现点云的导向(引导)滤波(附最简版)
开发语言·图像处理·算法·计算机视觉·matlab
熊猫_豆豆2 天前
MATLAB 九大行星太阳系运行程序
开发语言·nginx·matlab
xiaoningaijishu2 天前
MATLAB中的Excel文件操作:从入门到精通
其他·算法·matlab·excel