数值分析实验3报告_第1页
数值分析实验3报告_第2页
数值分析实验3报告_第3页
数值分析实验3报告_第4页
数值分析实验3报告_第5页
全文预览已结束

下载本文档

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

文档简介

中南民族大学实验报告实验课名称： 数值分析 导教师： 学生姓名： 年级专业： 2011信息 学号： 实验名称： Runge现象 实验日期： 2014 年 3 月 31 日 实验成绩： 实验目的1、观察多项式插值问题的Runge现象，了解数值不稳定现象；2、借助Lagrange插值，从图形中观察准确曲线与插值函数对应的曲线；3、熟练掌握Matlab中的画图工具。实验内容对于函数进行Lagrange插值，取不同的节点数（如取：），在区间上取等距间隔的节点为插值点。把准确曲线和插值多项式的曲线画在同一张图上进行比较。实验内容实验步骤 1、编写算法实现题目要求；2、观察图形，并调整节点个数，看精度有哪些变化；如，或 ，分别观察准确曲线与插值多项式的曲线的误差，回答：是否插值节点越多，即插值多项式次数越高，精度越高？3、验证分析实验结果，并整理报告；实验注意事项本程序的关键，需要调用实验1的Lagrange插值程序。如果区间划分为10等份：x0=-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 off实验内容在Matlab中建立M文件lagrange.m：function f=lagrange(x,y,t)f=zeros(size(t);n=length(y);for i=1:n z=ones(size(t); for j=1:n if i=j z=z.*(t-x(j)/(x(i)-x(j); end end f=f+z*y(i);endreturnn=10的时候，在命令窗口输入：x0=-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 off实验内容N=5的时候，输入：x0=-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(
人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

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

最新文档

评论

0/150

提交评论