数值分析实验插值与拟合.doc_第1页
数值分析实验插值与拟合.doc_第2页
数值分析实验插值与拟合.doc_第3页
数值分析实验插值与拟合.doc_第4页
数值分析实验插值与拟合.doc_第5页
已阅读5页,还剩1页未读 继续免费阅读

下载本文档

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

文档简介

数值分析课程实验一:插值与拟合一、实验目的1. 理解插值的基本原理,掌握多项式插值的概念、存在唯一性;2. 编写MATLAB程序实现Lagrange插值和Newton插值,验证Runge现象;3. 通过比较不同次数的多项式拟合效果,理解多项式拟合的基本原理;4. 编写MATLAB程序实现最小二乘多项式曲线拟合。二、实验内容1. 用Lagrange插值和Newton插值找经过点(-3, -1), (0, 2), (3, -2), (6, 10)的三次插值公式,并编写MATLAB程序绘制出三次插值公式的图形。2. 设如果用等距节点xi = -5 + 10i/n (i = 0, 1, 2, , n)上的Lagrange插值多项式Ln(x)去逼近它。不妨取n = 5和n = 10,编写MATLAB程序绘制出L5(x)和L10(x)的图像。3. 在某冶炼过程中,根据统计数据的含碳量与时间关系如下表,试求含碳量与时间t的拟合曲线。t (min)0510152025303540455055y (10-5)01.272.162.863.443.874.154.374.514.584.024.64(1) 用最小二乘法进行曲线拟合;(2) 编写MATLAB程序绘制出曲线拟合图。三、实验步骤1. (1) Lagrange插值法:在线性空间Pn中找到满足条件:的一组基函数,li(x)的表达式为有了基函数,n次插值多项式就可表示为(2) Newton插值法:设x0, x1, , xn是一组互异的节点,yi = f(xi) (i = 0, 1, 2, , n),f(x)在处的n阶差商定义为则n次多项式差商表的构造过程:xif(xi)一阶差商二阶差商三阶差商四阶差商x0f(x0)x1f(x1)fx0, x1x2f(x2)fx1, x2fx0, x1,x2x3f(x3)fx2, x3fx1, x2,x3fx0, x1,x2,x3x4f(x4)fx3, x4fx2, x3,x4fx1, x2,x3,x4fx0, x1,x2,x3,x4MATLAB程序实现:试验结果:2. MATLAB程序实现:试验结果:3. 多项式拟合的一般方法可归纳为以下几步:(1)由已知数据画出函数粗略的图形散点图,确定拟合多项式的次数n;

温馨提示

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

评论

0/150

提交评论