




已阅读5页,还剩1页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
数值分析实验报告实验名称: 观察龙格(Runge)现象实验 班级:12级信息与计算科学(1)班姓名: 学号:33 16 17 59实验日期: 2014.10.11 周次: 6 实验地点: A14-504 实验目 的. 观察拉格朗日插值的龙格(Runge)现象.实验内容对于函数进行拉格朗日插值,取不同的节点数n,在区间-5,5上取等距间隔的节点为插值点,把f(x)和插值多项式的曲线画在同一张图上进行比较。(a可以取任意值)具体步骤:1、 a=1时,1)取n=4,作出f(x)和插值多项式的曲线图; 2)取n=10,作出f(x)和插值多项式的曲线图;2、a=0.25时,1)取n=4,作出f(x)和插值多项式的曲线图; 2)取n=10,作出f(x)和插值多项式的曲线图;3、分析上述曲线图,你可以得出什么结论?实验步骤拉格朗日插值的调用函数function y=lagrange(x0,y0,x)n=length(x0);m=length(x);for i=1:m z=x(i); L=0.0; for j=1:n T=1.0; for k=1:n if k=j T=T*(z-x0(k)/(x0(j)-x0(k); end end L=T*y0(j)+L; end y(i)=L;enda=1时,1) 取a=1;n=4x0=-5:2:5;y0=5./(1+x0.2);x=-5:0.1:5;y=lagrange(x0,y0,x);y1=5./(1+x.2);plot(x,y,-r)hold onplot(x,y1,-b)hold off 2)取n=10x0=-5:1:5;y0=5./(1+x0.2);x=-5:0.1:5;y=lagrange(x0,y0,x);y1=5./(1+x.2);plot(x,y,-r)hold onplot(x,y1,-b)hold off2、a=0.25时,1)取n=4x0=-5:2:5;y0=5./(0.25*0.25+x0.2);x=-5:0.1:5;y=lagrange(x0,y0,x);y1=5./(1+x.2);plot(x,y,-r)hold onplot(x,y1,-b)hold off2) 取3)4) a=0.25n=10x0=-5:1:5;y0=5./(0.25*0.25+x0.2);x=-5:0.1:5;y=lagrange(x0,y0,x);y1=5./(1+x.2);plot(x,y,-r)hold onplot(x,y1,-b)hold off实验结果分析a=1时,1)取n=4a=1时,2)取n=10a=0.25时,1)取n=4a=0.25时,2)取n=10上述现象告诉我们,并不是插值多项式的次数越高(即插值节点越多)精度越高,从数值计算上可解释为高次插值多项式的计算会带来舍入误差的增大,从而引起计算失真。因此,实际应用做插值时一般只用一次、二次最多用三次插值多项式。那么如何提高插值精度呢?采用分段插值是一种办法。教师评语多种插值对比function y=fun(x); y=5./(1+x.2);endhours=-5:5;h=-5:0.1:5;temps=fun(hours);t1=interp1(hours,temps,h,spline); %()t2=interp1(hours,temps,h, nearest);t3=interp1(hours,temps,h, linear);t4=interp1(hours,temps,h, cubic);subplot(2,2,1); plot(hours,temps, bo,h,t1, r); title(spline ); %subplot(2,2,2); plot(hours,temps, bo,h,t2, r); title( nearest); subplot(2,2,3); plot(hours,temps, bo,h,t3, r); title
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- JJF 2313-2025多参数在线激光粒度仪校准规范
- 2025年福州东盟海产品交易所有限公司总经理职位职业经理人市场化选聘考前自测高频考点模拟试题及参考答案详解一套
- 2025福建三明大田县公开招聘紧缺急需专业教师7人考前自测高频考点模拟试题及1套完整答案详解
- 小学安全培训收费标准表课件
- 2025年临沂兰陵县教育系统部分事业单位公开招聘教师(5人)考前自测高频考点模拟试题有答案详解
- 2025江苏连云港市海州湾发展集团有限公司及子公司招聘20人考前自测高频考点模拟试题及完整答案详解
- 2025广东广州市中山大学孙逸仙纪念医院超声科医教研岗位招聘模拟试卷及答案详解(历年真题)
- 安全培训教学课件制作
- 2025江西吉安市直三家公立医院编外招聘33人考前自测高频考点模拟试题附答案详解(黄金题型)
- 2025福建省高速公路集团有限公司招聘43人考前自测高频考点模拟试题及参考答案详解
- 临床医师定期考核必刷题库及答案(一)
- 职业本科《大学英语》课程标准
- 2024年承包建设工程合同
- 英语语法课程教学大纲
- 《陆上风电场工程概算定额》NBT 31010-2019
- 水平四初中羽毛球大单元教学教案(18课时)
- 2024年河北石家庄市高速公路集团限公司面向社会公开招聘收费人员150名公开引进高层次人才和急需紧缺人才笔试参考题库(共500题)答案详解版
- 酒店住宿抵款协议书
- 【基于WBS分解图的工程项目施工进度管理与优化案例探析22000字(论文)】
- 配电箱安全专项教育培训课件
- 智慧医保监管一体化平台建设方案
评论
0/150
提交评论