数值积分——高斯求积

高斯求积推导(以三次精度为例)

将a~b之间的积分表示成如下形式:

I = c_0f(x_0) + c_1f(x_1)

其中c0,c1,x0,x1均为未知量,假设上面的式子能够精确计算常数、线性函数、抛物线和三次函数的积分,这样我们能够得到4个方程(a=1,b=1):

c_0+c_1 = \int_{-1}^11dx = 2

c_0x_0+c_1x_1 = \int_{-1}^1xdx = 0

c_0x_0^2+c1x_1^2 = \int_{-1}^1x^2dx = \frac{2}{3}

c_0x_0^3+c_1x_1^3 = \int_{-1}^1x^3dx=0

解得:

c_0 = c_1 = 1;x_0 = -\frac{1}{\sqrt{3}},x_1 = \frac{1}{\sqrt{3}}

通用性解析

        假设一个新的变量xd,其与x呈线性关系x = a_1+a_2x_d,下限x=a对应xd = -1,上限x = b对应xd = 1,得到:

x = \frac{(b+a)+(b-a)x_d}{2}

dx = \frac{b-a}{2}dx_d

代入原积分公式,可转换为通用形式。

  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

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

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值