梯形求积公式 和 复合梯形求积公式 Matlab 实现

梯形求积公式 和 复合梯形求积公式 Matlab 实现

梯形求积公式

仅使用区间两点 x1,f(x1),x2,f(x2)
组成的梯形面积S代替 x2x1f(x)dx 的近似方法
求积

x2x1f(x)dxS=x2x12×(f(x2)f(x1))

复合梯形求积公式

将求积区间[a,b]分为n个区间,每个区间步长为h( h=ban )在每个区间求梯形积分
复合梯形求积
Si为第i个梯形的面积

x2x1f(x)dxi=0nSi=h2(f(a)+k=1n1f(xk)+f(b))

说明

上述公式是我用mathjax写的,如有错误请联系我修正
敬请指正
概述省略了部分推导过程,请查阅详细推导资料

Matlab 实现代码

梯形求积公式

将该函数存为m文件

function res = Trapezium(f,a,b)
    format long;
    if b < a
        c = b;
        b = a;
        a = c;
    end
    res = (b-a)*(f(a)+f(b))/2;

调用下面语句测试函数

f = inline('sin(x)','x')
Trapezium(f,0,pi/2)

复合梯形求积公式

将该函数存为m文件

function res = ComTrapezium(f,n,a,b)
    format long;
    if b < a
        c = b;
        b = a;
        a = c;
    end
    h = (b-a)/n;
    d = f(a);
    for i = a+h:h:b-h 
        d = d + (2 * f(i));
    end
    d = d + f(b);
    res = (d * h / 2);

调用下面语句测试函数

f = inline('sin(x)','x')
ComTrapezium(f,4,0,pi/2)

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
©️2022 CSDN 皮肤主题:大白 设计师:CSDN官方博客 返回首页
评论 6
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值