《数值分析》-- 高斯求积公式

本文深入探讨了高斯型求积公式,包括高斯-勒让德和高斯-切比雪夫公式,阐述了它们的理论基础和特性。在节点个数固定的情况下,讨论了能否通过选择节点位置提高求积公式的精度。同时,介绍了利用正交多项式构造求积公式的方法,并给出了不同节点数量下的示例和求积公式,如两个点的高斯-勒让德公式。文章还回顾了复合梯形公式和Simpson公式等经典积分方法。
摘要由CSDN通过智能技术生成


概述

  • 问题
    那么,在节点个数一定的情况下,是否可以在[a,b]上自由选择节点的位置,使求积公式的精度提得更高?
    在这里插入图片描述
  • 例题
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

一、高斯型求积公式的一般理论

  • 一般理论
    在这里插入图片描述

1.1 高斯型求积公式和高斯点

在这里插入图片描述
习题

1.2 高斯点的特征

在这里插入图片描述

  • 利用正交多项式构造高斯求积公式
    在这里插入图片描述
    习题

二、常用的高斯求积公式⭐

2.1 高斯-勒让德求积公式( Gauss-Legendre )

  1. Legendre 多项式族:
    在这里插入图片描述
    低阶Legendre多项式
    在这里插入图片描述
  2. 高斯-勒让德求积公式(G-L求积公式)⭐
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
  1. 一般有界区间[a, b]上的高斯-勒让德求积公式
    (G-L求积公式)
    目的转换区间到 [-1,1] :
    在这里插入图片描述
    习题

2.2 高斯-切比雪夫求积公式( Gauss-Chebyshe)

切比雪夫多项式
在这里插入图片描述
在这里插入图片描述

  • 截断误差
    在这里插入图片描述

总结

  • 例题
  1. 在这里插入图片描述
    待定系数法构造高斯求积公式:
    在这里插入图片描述
  • 问题
    为什么2点Gauss公式有应该有三次代数精度

一般n+1个节点的求积公式的代数精度最高为2n+1次

在这里插入图片描述
在这里插入图片描述

  1. 在这里插入图片描述
    在这里插入图片描述
  • 详细过程
    在这里插入图片描述
  1. 在这里插入图片描述
    在这里插入图片描述
  • 详细过程
    在这里插入图片描述
    在这里插入图片描述
  • 例题
  1. 在这里插入图片描述
    利用正交多项式构造高斯求积公式:
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
  • 例题
  1. 在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
  • 公式回忆
  1. 复合梯形公式
    在这里插入图片描述
    三个点: T 2 T2 T2
    五个点: T 4 T4 T4
  2. Simpson公式
    在这里插入图片描述
  3. 三个点的高斯-勒让德求积公式
    在这里插入图片描述
  • 例题
  1. 在这里插入图片描述
    在这里插入图片描述
  • 公式
    在这里插入图片描述
  1. 在这里插入图片描述
    在这里插入图片描述
  • 两点高斯-勒让德求积公式
    在这里插入图片描述

  • 例题
  1. 在这里插入图片描述
  • 42
    点赞
  • 214
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
高斯比雪夫求积是一种数值积分方法,用于求解积分。该方法利用比夫多项式的零点和系数逼近被积函数,从而得到分的近似值。 在MATLAB中,可以使用"Gauss-Chebyshev quadrature"函数来实现高斯比雪夫求积。该函数使用比雪夫多项式的零点和系数,以及被积函数的值来计算积分的近似值。 以下为使用MATLAB的高斯比雪夫求积的示例代码: ```matlab % 定义被积函数 f = @(x) x.^2 + 2.*x + 1; % 求解高斯比雪夫求积 n = 5; % 求积点数 [a, b = [-1, 1]; % 积分区间 x = cos(pi*(0:n-1)/(n-1)); % 比雪夫零点 w = pi/(n-1)*ones(size(x)); % 比雪夫系数 approximation = sum(w.*f((b-a)/2*x+(b+a)/2)); % 近似积分值 % 显示结果 disp(approximation); ``` 在上述示例代码中,我们定义了一个被积函数f(x) = x^2 + 2x + 1。然后,我们选择了5个积分点,并计算了比雪夫零点和系数。最后,我们使用近似求积公式计算了积分的近似值。 请注意,以上代码仅为示例,如果你有具体的被积函数,需要根据函数的特点和积分区间来调整代码。另外,MATLAB中还提供了其他的数值积分函数,你可以根据实际需求选择合适的方法进行求解。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* [Matlab数值积分法汇总_牛顿辛普森_Matlab数值积分法汇总_高斯求积_高斯埃尔米特_matlab数值积分](https://download.csdn.net/download/weixin_42696271/25541913)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 50%"] - *2* *3* [比雪夫求积分法(附全过程代码)](https://blog.csdn.net/m0_37772174/article/details/88542702)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

胜天半月子

打不打商的无所谓,能帮到你就好

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值