数值分析(第2版)课件6.4 高斯型求积公式_第1页
数值分析(第2版)课件6.4 高斯型求积公式_第2页
数值分析(第2版)课件6.4 高斯型求积公式_第3页
数值分析(第2版)课件6.4 高斯型求积公式_第4页
数值分析(第2版)课件6.4 高斯型求积公式_第5页
已阅读5页,还剩54页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

6.4高斯型求积公式(1)曲绍波12正交多项式是满足正交性的多项式序列,具备许多良好的性质,被广泛应用于数学物理、工程技术、科学计算等诸多领域,在数学理论研究与工程实践中都起着非常重要的作用。6.4.1

正交多项式

3

4

5内积的性质:

6例1

三角函数族

验证:

7

8

9正交化方法:

10正交多项式的性质:

11

(3)相邻的正交多项式之间存在递推关系

12

证明:

13

证明:

于是

14(一)勒让德多项式

Rodrigue(1814)给出了勒让德多项式的简单表达式

15

展开

16勒让德多项式的性质:(1)正交性

(2)奇偶性

17(4)递推公式

利用递推公式可得

1819(二)切比雪夫多项式

(三)第二类切比雪夫多项式

20(四)拉盖尔多项式

(五)埃尔米特多项式

6.4高斯型求积公式(2)曲绍波2122考虑如下形式积分的计算6.4.2

Gauss型求积公式的一般理论

构造如下形式的求积公式

23

插值型求积公式选取求积节点;构造插值多项式;对插值多项式积分,即得插值型求积公式(求积系数);24

插值型求积公式若节点取为积分区间的等分点,则得到牛顿-柯特斯公式。问题:

25例1

构造形如

(梯形公式)26解:令求积公式

由第二式和第四式得

再结合第一式和第三式得

取代入方程组,解得

27于是得到求积公式

可以验证,该求积公式具有3次代数精度。

问题:28

证明:

29

证明:

30定理1插值型求积公式

31证明必要性假设求积公式

32

证明充分性33

所以有

34

定理1插值型求积公式

35

构造高斯型求积公式:

36例2

构造如下形式的高斯型求积公式

37于是,得

解得

于是,得高斯型求积公式38求积余项:高斯型求积公式的求积余项为

证明考察Lagrange插值基函数

39

从而有

推论2高斯型求积公式是数值稳定的。40

6.4高斯型求积公式(3)曲绍波4142高斯型求积公式6.4.3

常用的Gauss型求积公式

43

(一)高斯-勒让德求积公式

44此时,节点为2次勒让德多项式(1)两点高斯-勒让德求积公式

的两个零点,即

45

于是,有

46在两点高斯-勒让德求积公式中,

高斯型求积公式的求积余项为

于是得两点高斯-勒让德求积公式的余项

两点高斯-勒让德求积公式具有三次代数精度。47此时,节点为3次勒让德多项式(2)三点高斯-勒让德求积公式

的三个零点,即

48

于是,有

从而得到三点高斯-勒让德求积公式

解得

三点高斯-勒让德求积公式具有5次代数精度。49表1高斯-勒让德求积公式的节点和系数00.000000002.000000004±0.90617985±0.538469310.000000000.236926890.478628670.568888891±0.577350271.000000002±0.774596670.000000000.555555560.888888895±0.93246951±0.66120939±0.238619190.171324490.360761570.467913933±0.86113631±0.339981040.347854850.6521451550作变换

对于等号右端的积分可以使用高斯-勒让德求积公式进行计算。51解:作变量代换例1

用高斯-勒让德求积公式计算积分

分别使用两点、三点和五点高斯-勒让德求积公式。52经计算得

同理

53高斯-勒让德求积公式牛顿-柯特斯公式两点梯形公式三点辛普森公式五点柯特斯公式表2计算结果比较

54(二)高斯-切比雪夫求积公式

55

经过计算,得求积系数

于是得高斯-切比雪夫求积公式

求积余项为

56

于是得两点高斯-切比雪夫求积公式

于是得三点高斯-切比雪夫求积公式

57例2

用高斯-切比雪夫求积公式计算积分

解作变量替换

则有

权函数被

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论