matlab实现lagrange插值和分段线性插值

数值分析作业

姓名:虞驰程

题目:

函数:f x =1+x2在[-5,5]上,取n=10,对其进行分段线性插值和拉格朗日插值,在Matlab 中实现且绘图。

Matlab 实现:

首先定义函数f ,在Matlab 中用function.m 文件编写, 具体代码如图1所示: 1

图1 f(x)函数

定义分段线性插值的基本函数,用function.m 文件编写,具体代码如图2所示:

图2 分段线性插值基本函数

定义拉格朗日插值的基本函数,用function.m 文件编写,具体代码如图3所示:

图3 拉格朗日插值的基本函数

进行分段线性插值并绘图和原函数进行对比的Matlab 实现代码如图4所示:

图4 分段线性插值函数绘制

其结果如图5所示,其中红色代表分段线性插值结果,蓝色代表原函数:

图5 分段线性插值和原函数对比

同理可以进行拉格朗日插值并绘图,其Matlab 实现代码如图6所示,其结果如图7所示:

图6 拉格朗日插值函数绘制

图7 拉格朗日插值和原函数对比

最后我们可以将分段线性插值、拉格朗日插值和原函数进行对比,其实现代码如图8所示,最终结果如图9所示(黑色代表原函数,蓝色是分段线性插值,红色是拉格朗日插值):

图8 两种插值方法和原函数对比实现

图9 两种插值方法和原函数对比

数值分析作业

姓名:虞驰程

题目:

函数:f x =1+x2在[-5,5]上,取n=10,对其进行分段线性插值和拉格朗日插值,在Matlab 中实现且绘图。

Matlab 实现:

首先定义函数f ,在Matlab 中用function.m 文件编写, 具体代码如图1所示: 1

图1 f(x)函数

定义分段线性插值的基本函数,用function.m 文件编写,具体代码如图2所示:

图2 分段线性插值基本函数

定义拉格朗日插值的基本函数,用function.m 文件编写,具体代码如图3所示:

图3 拉格朗日插值的基本函数

进行分段线性插值并绘图和原函数进行对比的Matlab 实现代码如图4所示:

图4 分段线性插值函数绘制

其结果如图5所示,其中红色代表分段线性插值结果,蓝色代表原函数:

图5 分段线性插值和原函数对比

同理可以进行拉格朗日插值并绘图,其Matlab 实现代码如图6所示,其结果如图7所示:

图6 拉格朗日插值函数绘制

图7 拉格朗日插值和原函数对比

最后我们可以将分段线性插值、拉格朗日插值和原函数进行对比,其实现代码如图8所示,最终结果如图9所示(黑色代表原函数,蓝色是分段线性插值,红色是拉格朗日插值):

图8 两种插值方法和原函数对比实现

图9 两种插值方法和原函数对比


相关文章

  • 线性拟合法
  • 摘 要 摘 要 插值法和曲线拟合是两种来源于实际,同时又广泛应用于实际的重要的数值 计算方法.随着计算机技术的不断发展以及人类计算机水平的逐步提高,他们在国民经济和科学研究中占据了越来越重要的地位.插值法与曲线拟合结合计算机技术例如MATL ...查看


  • 数值分析课程设计(最终版)
  • 本文主要通过Matlab 软件,对数值分析中的LU 分解法.最小二乘法.复化Simpon 积分.Runge-Kutta 方法进行编程,并利用这些方法在MATLAB 中对一些问题进行求解,并得出结论. 实验一线性方程组数值解法中,本文选取LU ...查看


  • 什么叫拟合?什么叫插值?二者的区别是什么?
  • 什么叫拟合?什么叫插值?二者的区别是什么? 插值和拟合都是函数逼近或者数值逼近的重要组成部分 他们的共同点都是通过已知一些离散点集M上的约束,求取一个定义 在连续集合S(M包含于S)的未知连续函数,从而达到获取整体规律的 目的,即通过&qu ...查看


  • 数值分析学习心得体会
  • 数值分析学习感想 一个学期的数值分析,在老师的带领下,让我对这门课程有了深刻的理解和感悟.这门 课程是一个十分重视算法和原理的学科,同时它能够将人的思维引入数学思考的模式,在处 理问题的时候,可以合理适当的提出方案和假设.他的内容贴近实际, ...查看


  • 计算方法拉格朗日插值
  • 第二章 插值法 知识点:拉格朗日插值法,牛顿插值法,误差,龙格现象,分段插值. 1.背景 实践活动中,表现事务变化的信息往往只是一些离散点值,例如 每个6小时记录一次温度,以此反映一天的气温变化状况,如下表图 温度(.C) 34 32 30 ...查看


  • MATLAB实现拉格朗日插值
  • 数值分析上机报告 题目:插值法学号:姓名:靳会有 201014924 一.调用MATLAB内带函数插值 1.MATLAB内带插值函数列举如下: 2.取其中的一维数据内插函数()为例,程序如下: 其调用格式为: yi=interp1(x, y ...查看


  • 插值算法与matlab代码
  • Matlab 中插值函数汇总和使用说明 MATLAB 中的插值函数为interp1,其调用格式为: yi= interp1(x,y,xi,'method') 其中x ,y 为插值点,yi 为在被插值点xi 处的插值结果:x,y 为向量, ' ...查看


  • 分段线性插值 2
  • 1.4分段插值 一. 分段线性插值 在每个区间 i i + 1 上,用1阶多项式 (直线) 逼近 f (x): for x ∈[x i , x i +1] 即用折线代替曲线.设f (x )连续 一f (x ) P 1h (x ) →记 h ...查看


  • 基于MATLAB的发动机万有特性曲线绘制方法
  • 2009年第2期(总第110期) 内燃机与动力装置I.c.E&Powerplant 2009年4月 [设计研究] 基于MATLAB的发动机万有特性曲线绘制方法 周广猛1,郝志刚2,刘瑞林1,陈东3,管金发1,张春海4 (1.军事交通 ...查看


热门内容