由于MATLAB中没有现成的拉格朗日插值命令,因此我们可以根据Lagrange插值法的定义编写一个Lagrange插值命令。

定义:

给定n个插值节点x1,x2,...,xn和各个节点对应的函数值y1,y2,...yn,利用拉格朗日插值多项式公式(注:k从0到n累加),其中,  ,则可以得到在插值区间内的任意x的函数值y为。从公式中可以看出,生成的多项式与用来插值的数据密切相关,数据变化则函数就要重新计算。

以下为M语言编写的函数文件lagrange.m

%拉格朗日插值命令
function yy=lagrange(x,y,xx)
%Lagrange 插值,求数据(x,y)所表达式的函数在插值点xx处的插值
m=length(x);
n=length(y);
if m~=nerror('向量x,y长度必须一致');
end
s=0;
for i=1:nt=ones(1,length(xx));for j=1:nif j~=it=t.*(xx-x(j))/(x(i)-x(j));endends=s+t*y(i);
end
yy=s;

读者可以自己找一个题目自行实验,通过画出的图像可以发现,拉格朗日插值的一个明显的特点为——拟合出的多项式通过原x,y的每一个对应的点。

MATLAB——数值插值之拉格朗日(Lagrange)插值法相关推荐

  1. matlab 数值插值,MATLAB数值实验一(数据的插值运算及其应用完整版)

    佛山科学技术学院 实 验 报 告 课程名称 数值分析 实验项目 插值法与数据拟合 专业班级 机械工程 姓 名 余红杰 学 号 10 指导教师 陈剑 成 绩 日 期 月 日 一.实验目的 1.学会Lag ...

  2. matlab中怎么写拉格朗日函数,MATLAB编辑n次拉格朗日函数插值法的程序

    n次拉格朗日函数插值 程序: function LagrangesNs() %用于求过n点的拉格朗日n-1次插值多项式options={'Name of data file'}; title='Lag ...

  3. 实验Matlab数值运算,MATLAB数值实验一(数据的插值运算及其应用完整版

    <MATLAB数值实验一(数据的插值运算及其应用完整版>由会员分享,可在线阅读,更多相关<MATLAB数值实验一(数据的插值运算及其应用完整版(6页珍藏版)>请在人人文库网上搜 ...

  4. matlab全域基函数,多项式函数插值:全域多项式插值(一)单项式基插值、拉格朗日插值、牛顿插值 [MATLAB]...

    全域多项式插值指的是在整个插值区域内形成一个多项式函数作为插值函数.关于多项式插值的基本知识,见"计算基本理论". 在单项式基插值和牛顿插值形成的表达式中,求该表达式在某一点处的值 ...

  5. matlab实现查值,[数学建模(六)]使用MATLAB实现插值

    常用的插值:拉格朗日多项式插值.牛顿插值.分段线性插值.Hermite 插值和三次样条插值. 1.拉格朗日插值法 function y=lagrange(x0,y0,x); n=length(x0); ...

  6. 数学建模准备 插值(拉格朗日多项式插值,牛顿多项式插值,分段线性插值,分段三次样条插值,分段三次Hermite插值)

    文章目录 摘要(必看) 0 基础概念 什么是插值 插值用途 什么是拟合 插值和拟合的相同点 插值和拟合的不同点 1 常用的基本插值方法 1.1 多项式插值法 1.1.1 拉格朗日多项式插值法 多项式插 ...

  7. MATLAB数值分析学习笔记:多项式拟合

    在很多情况下,我们需要估计精确值之间的可能取值,常用的方法是多项式插值法.多项式插值法的原理是n个数据点可以确定n-1次多项式:例如三点确定一条抛物线. 工程上常用的多项式插值的方法有两种:牛顿插值法 ...

  8. 拉格朗日插值_拉格朗日插值定理的理论基础

    缺失,几乎是不可避免的.只要做数据处理,不可避免的工作就是插值.而插值里面比较常用的方法之一就是拉格朗日插值法,这篇文章就跟大家讲讲拉格朗日插值的理论基础. 为什么需要进行插值 我们进行数据处理的理想 ...

  9. 拉格朗日多项式插值法 思路详解+详细代码

    在数学建模过程中大家经常会使用插值法对数据进行处理,而其中拉格朗日多项式插值法是较为常用到的. 以下是我在学习拉格朗日插值法时通过阅读许多大佬博主的文章时发现,要么只有代码,要么只有理论讲解或者例题, ...

最新文章

  1. 高德地图SDK大致使用
  2. codeforces 483B Friends and Presents
  3. 让你IDA Pro 快速上手使用
  4. Linux的实际操作:文件目录类的实用指令(ln history)
  5. artcam 9.0英文版本下载_tomcat 下载安装
  6. IoC(控制反转)的主要组件和注入的两种方式
  7. python爬虫怎么挣钱_python爬虫12 | 爸爸,他使坏,用动态的 Json 数据,我要怎么搞?...
  8. LeetCode4. 寻找两个有序数组的中位数
  9. [线性相关] 皮尔森相关系数的计算及假设检验
  10. P6615 Kruskal + 构造
  11. 什么人适合学平面设计?
  12. 不喜欢现在的领导,怎么办?不懂跟领导相处,你到哪都混不好
  13. 【Java】多人协同在线编辑Excel源代码及部署包
  14. 部署本地thinkphp6(iis+php7)
  15. Niushop 分销推广员管理
  16. ubuntu 生成桌面快捷方式(有图)
  17. 获取frame 标签里#documents里面的内容
  18. 【算法】求解钱币兑换问题
  19. html5+文本间距,html怎么设置文字的间距
  20. 【手机投影】安卓手机投影到WIN10

热门文章

  1. 中化CIO彭劲松:IT治理让我明明白白做事
  2. 吃鸡游戏(18.11.24)
  3. 齐齐哈尔市太美家政公司管理系统
  4. 点到反比例函数最短距离怎么求_数学原来靠“背”的,这几首顺口溜瞬间帮你记住数学重点公式和法则!...
  5. springboot项目打成公共jar包被其他项目单独引用(包含mybatis-plus)
  6. 【C语言】exit()和_exit()的区别
  7. removebg图片换背景_Python版
  8. 苹果电脑提示打不开 因为它来自身份不明的开发者 不能安装下载的苹果软件解决方法
  9. Mysql:max函数参数为字符类型需要通过CAST函数转换为数值型
  10. 第56章、复制网站、发布网站区别?