CHART4 带权的插值型求积公式_第1页
CHART4 带权的插值型求积公式_第2页
CHART4 带权的插值型求积公式_第3页
CHART4 带权的插值型求积公式_第4页
CHART4 带权的插值型求积公式_第5页
已阅读5页,还剩27页未读 继续免费阅读

下载本文档

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

文档简介

1、带权的插值型求积公式其中 为a,b上的权函数。8/7/202218/7/20222带权的插值型求积公式(2)代数精确度至少为n次8/7/20223求积公式 含有2n+2个待定参数 能否通过节点以及求积系数的选择将代数精确度提的更高,即超过n次,最高能达到多少次?5 Gauss型求积公式8/7/20224也就是说无论节点及其系数如何选择,求积公式的精度不可能达到2n+2次。 讨论续无论怎么选择节点,总存在多项式8/7/20225若 是 上的一组互异节点,且求积公式达到2n+1次代数精度,则称该求积公式为Guass型求积公式,其求积节点 (k=0,1,n)称为高斯点,系数 称为高斯系数。 Gaus

2、s型求积公式定义8/7/20226高斯型求积公式一定是插值型求积公式,其系数由高斯点唯一确定。 直接利用代数精确度的定义求得高斯点以及高斯系数需要解非线性方程组,很困难。 当高斯点确定后,可以用基函数插值的方法或者解线性方程组的方法求得。高斯型求积公式是精度最高的求积公式。 结论8/7/20227高斯点确定以后,高斯系数也可以由如下插值型求积公式中的系数公式确定: 确定.即可由线性方程组8/7/20228求积公式 至少具有n次代数精确度的充要条件是它是插值型的 充分性 :如果求积公式为插值型,利用截断误差知对于任意次数n的多项式 f(x), 有Rf=0,故求积公式至少具有n次精度。引理8/7/

3、20229必要性 :设求积公式 具有n次代数精度。 用n次插值函数, 仍有,根据插值型求积公式定义知,其求积公式为插值型求积公式。 必要性证明8/7/202210定理:插值型求积公式中的节点 是高斯点的充要条件是,在a,b上,以这些点为零点的n+1次多项式 与任意次数不超过n的多项式P(x)带权 正交,即高斯点的选取定理8/7/202211必要性 :证明: 设 是高斯点,于是对任意次数不超过n的多项式P(x) ,的次数不超过2n+1 充分性 :对任意次数不超过2n+1的多项式 f(x)用 除的商为p(x),余项为q(x)。对任意次数不超过n的多项式P(x) 有8/7/202212充分性所给的求

4、积公式是插值型的,其代数精度至少为n。 即求积公式具有2n+1次代数精度,从而 是一组高斯点。8/7/202213当 为正交多项式系中的n+1次多项式取 ,则 有n+1个互异的零点,且对任意次数不超过n的多项式有 a,b上带权 正交的n+1次多项式的零点就是高斯型求积公式的一组高斯点。由正交多项式的性质知它在开区间上存在n+1个互不相同的零点。 Remark8/7/202214.高斯型求积公式是收敛的。.高斯型求积公式是稳定的。(j=0,1,n) 故高斯求积系数Aj一定为正。高斯公式是稳定的。高斯型求积公式的收敛性和稳定性注:收敛性论证需用Weierstrass定理。8/7/202215高斯型

5、求积公式的截断误差定理: 设 在 内只有2n+2阶导数,则高斯型求积公式的余项为:证明: 设 为满足的Hermite插值多项式,则 次数 。8/7/202216由于高斯型求积公式的代数精度为2n+1,故 高斯型求积公式具有代数精度高、且总是收敛、稳定的优点。也可构造复化高斯求积公式。8/7/2022171.高斯勒让德求积公式 几种特殊的高斯型求积公式8/7/202218当积分区间为 时,可通过变换 将 变换为 高斯点 为n+1次切比雪夫多项式的零点: 高斯切比雪夫求积公式 8/7/202219高斯拉盖尔求积公式8/7/202220高斯-埃尔米特求积公式8/7/202221例:求高斯型求积公式的系数 及节点解:对函数类f(x)=1, 积分公式精确成立的。高斯型求积公式的构造举例 8/7/202222求解方法设高斯点是二次函数的零点.8/7/202223 6 数值微分6.1 插值法建立求导公式:插值型求导公式以离散数据 近似表达 插值多项式 的导数作为未知函数的导数近似8/7/202224 两点公式8/7/202225 三点公式(等距节点)8/7/2022268/7/2022278/7/202228如二阶三点公式高阶导数数值微分公式8/7/202

温馨提示

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

评论

0/150

提交评论