欢迎来到淘文阁 - 分享文档赚钱的网站! | 帮助中心 好文档才是您的得力助手!
淘文阁 - 分享文档赚钱的网站
全部分类
  • 研究报告>
  • 管理文献>
  • 标准材料>
  • 技术资料>
  • 教育专区>
  • 应用文书>
  • 生活休闲>
  • 考试试题>
  • pptx模板>
  • 工商注册>
  • 期刊短文>
  • 图片设计>
  • ImageVerifierCode 换一换

    数值分析实验报告2——Runge现象(共9页).docx

    • 资源ID:15114697       资源大小:134.57KB        全文页数:9页
    • 资源格式: DOCX        下载积分:20金币
    快捷下载 游客一键下载
    会员登录下载
    微信登录下载
    三方登录下载: 微信开放平台登录   QQ登录  
    二维码
    微信扫一扫登录
    下载资源需要20金币
    邮箱/手机:
    温馨提示:
    快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。
    如填写123,账号就是123,密码也是123。
    支付方式: 支付宝    微信支付   
    验证码:   换一换

     
    账号:
    密码:
    验证码:   换一换
      忘记密码?
        
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,就可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰。
    5、试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。

    数值分析实验报告2——Runge现象(共9页).docx

    精选优质文档-倾情为你奉上数值分析课程实验报告插值逼近题目一.Runge函数的插值1. Runge函数Runge函数的表达式为:其在-1,1区间上的函数图像如图1.1。在课程学习中我们知道,对Runge函数进行高次插值时有可能在两端出现不收敛的情况,即Runge现象。下面将分别用四种不同的插值方法在-1,1区间上对Runge函数进行插值,并分析是否产生Runge现象,比较插值效果。图1.1.Runge函数在-1,1区间的函数图像2.Newton插值首先根据课本上的Newton插值算法进行编程(代码略)。核心思想就是用符号变量进行中间运算,以便将最终的插值函数用符号表达式表示出来,并进一步生成图像。此处插值节点选择为等距插值节点,即:其中h=0.1。插值曲线与原曲线的对比如图1.2(蓝色为原曲线,红色为插值曲线)。从图中看出,在区间中部,二者吻合较好;但在区间两端二者则产生了明显偏差,甚至可以达到一个非常大的数值(e20量级)。因此,在等距节点的20次Newton插值下,产生了明显的Runge现象。图1.2.Newton插值曲线与原曲线对比3. Lagrange插值此处同样是根据Lagrange插值的具体算法进行编程。但插值节点不再是等距分布,而是如下形式:插值曲线与原曲线的对比如图1.3(蓝色为原曲线,红色为插值曲线)。从图中看出,插值曲线与原曲线吻合的很好,没有产生明显的Runge现象。对比产生了明显Runge现象的20次Newton插值,Lagrange插值的最高次数虽然也是20,但由于此处的插值节点不是等距分布的(事实上,此处采用的插值节点正是Chebyshev多项式的零点),而是中间疏两边密,因此两侧较密的节点很好地抑制了Runge现象。图1.3. Lagrange插值曲线与原曲线对比4. 分段线性插值分段线性插值是这几种插值方法中最容易处理的一个,只需要将每个节点对应的函数值求出再将相邻的数据点两两用直线相连即可。此处采用了等距节点,所得插值曲线与原曲线对比如图1.4(蓝色为原曲线,红色为插值曲线)。从图中图1.4. 分段线性插值曲线与原曲线对比看出,此处分段线性插值的效果也还是不错的,二者只在区间中部略微存在一些偏差,而在其他区域整体上吻合的很好,并且不存在Runge现象。这是由于分段线性插值通过对插值区间分段的方法将插值函数的次数有效降低,因而即使是等距节点分布,也很好地避免了出现Runge现象的倾向。5. 三次样条插值三次样条插值是这四种插值方法中编程最麻烦的,但并不是说存在多大的技术难度,只是因为插值过程中的步骤比较繁琐,因而代码也显得较为冗长。此处依然采用等距节点,所得插值曲线与原曲线对比如图1.5(蓝色为原曲线,红色为插值曲线)。从图中看出,三次样条插值的效果比分段线性插值更胜一筹,三次样条插值曲线和原曲线在整个插值区间都基本处于重合状态,几乎没有肉眼可见的偏差。同样,由于三次样条插值的插值函数最高次数只有3,在等距节点下也没有产生Runge现象。图1.5.三次样条插值曲线与原曲线对比题目二.分段函数的插值1. 分段函数定义在-1,1区间的分段函数的函数表达式为:其函数图像如图2.1。分段函数最大的特点就是在个别点上函数值或导数值存在突变,因此可以预计,除了可能出现的Runge现象外,在那些突变点附近的插值结果也可能会出现较大的偏差。下面将分别采用之前的四种插值方法在该函数的-1,1定义域内对其进行插值。图2.1.分段函数图像2.Newton插值首先根据课本上的Newton插值算法进行编程。此处插值节点选择为等距插值节点,即:其中h=0.1。插值曲线与原曲线的对比如图2.2(蓝色为原曲线,红色为插值曲线)。从图中看出,与Newton法对Runge函数的插值结果相比,Newton法对于该分段函数的插值效果显得更加糟糕:不仅在区间两端产生了极强烈的震荡(即Runge现象),就连区间中部也存在较小的上下震荡。因此,从整体来看,几乎所有距插值节点稍远的点都存在较大的偏差,这表明该分段函数在等距节点下的20次Newton插值效果非常不理想。图2.2.Newton插值曲线与原曲线对比3. Lagrange插值此处同样是根据Lagrange插值的具体算法进行编程。但插值节点不再是等距分布,而是如下形式:插值曲线与原曲线的对比如图2.3(蓝色为原曲线,红色为插值曲线)。从图中看出,与同样次数的Newton法相比,Lagrange法所得的插值曲线虽然在区间中部的分布与其相似,但在区间两端较好地收敛到了原曲线上,即较好地消除了Runge现象。这同样是因为此处的插值节点不是等距分布的(事实上,此处采用的插值节点正是Chebyshev多项式的零点),而是中间疏两边密,因此两侧较密的节点很好地抑制了Runge现象。图2.3. Lagrange插值曲线与原曲线对比4. 分段线性插值分段线性插值是这几种插值方法中最容易处理的一个,只需要将每个节点对应的函数值求出再将相邻的数据点两两用直线相连即可。此处采用了等距节点, 图2.4. 分段线性插值曲线与原曲线对比所得插值曲线与原曲线对比如图2.4(蓝色为原曲线,红色为插值曲线)。从图中看出,此处分段线性插值的效果较好,二者只在区间中部函数值的突变点附近存在一些偏差,而在其他区域整体上吻合的很好,不存在Runge现象。这是由于分段线性插值通过对插值区间分段的方法将插值函数的次数有效降低,因而即使是等距节点分布,也很好地避免了出现Runge现象的倾向。5. 三次样条插值三次样条插值是这四种插值方法中编程最麻烦的,但并不是说存在多大的技术难度,只是因为插值过程中的步骤比较繁琐,因而代码也显得较为冗长。此处依然采用等距节点,所得插值曲线与原曲线对比如图2.5(蓝色为原曲线,红色为插值曲线)。从图中看出,三次样条插值的效果与分段线性插值相近,也是在区间中部的函数值突变处有一定的偏差,而其他区域都吻合较好,也没有产生Runge现象。同样,这也是由于三次样条插值的插值函数最高次数只有3,因此在等距节点下进行插值也没有产生Runge现象。图2.5.三次样条插值曲线与原曲线对比三.总结本文通过matlab编程分别采用Newton插值法、Lagrange插值法、分段线性插值法以及三次样条插值法对Runge函数和一个分段函数进行了插值逼近,插值区间-1,1,插值节点21个,并通过分析计算结果主要得到了以下结论:1.插值多项式次数过高时会产生严重的Runge现象。本实验中,无论是Runge函数还是分段函数的20次Newton插值多项式都产生了严重的Runge现象,区间两端处的插值出现剧烈震荡,严重失真。2.同样是高次插值多项式,若适当选取插值节点能够在一定程度上抑制Runge现象。本实验中,Runge函数和分段函数的20次Lagrange插值多项式由于采用了中间疏两边密的非等距结点,而不是Newton插值多项式所用的等距节点,有效地抑制或消除了本应出现的Runge现象(这里说“本应出现”是经过计算验证的,若取等距节点,则20次的Lagrange插值多项式也会出现严重Runge现象)。3.降低插值多项式的次数能有效避免Runge现象。本实验中,分段线性插值法(各区间上均为1次)和三次样条插值法(最高次数为3)都取得了较为理想的差值逼近效果,没有出现Runge现象,且在整个插值区间都与原函数的图像吻合的很好。4.与连续函数相比,存在不连续点的分段函数的插值逼近误差更大,且更加不稳定。本实验中,对连续的Runge函数进行插值逼近时,除了等距节点的高次Newton多项式出现严重Runge现象,其余三种方法基本都收敛到了原曲线上,取得了不错的插值逼近效果;而对分段函数进行插值逼近时,除了等距节点的高次Newton多项式的逼近效果非常糟糕外(巨大偏差,严重震荡),其余三种方法虽然没有出现Runge现象,但在不连续点(x=0)的附近区域都存在一定的误差,整体逼近效果逊色于对连续Runge函数的插值逼近。综上,在实际运用中,为了取得较好的插值逼近效果,应尽量保证以下几点:不采用次数过高的插值多项式;适当选取插值节点;避免函数值突变,若不得已对存在不连续点的函数进行插值逼近,可以尝试分段插值,并将不连续点都处理到子区间的端点上,从而原函数在各子区间内分段连续,以便提高插值逼近的效果。专心-专注-专业

    注意事项

    本文(数值分析实验报告2——Runge现象(共9页).docx)为本站会员(飞****2)主动上传,淘文阁 - 分享文档赚钱的网站仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知淘文阁 - 分享文档赚钱的网站(点击联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    关于淘文阁 - 版权申诉 - 用户使用规则 - 积分规则 - 联系我们

    本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

    工信部备案号:黑ICP备15003705号 © 2020-2023 www.taowenge.com 淘文阁 

    收起
    展开