matlab 拉格朗日插值 Matlab编程实现拉格朗日插值法计算插值

Matlab编程实现拉格朗日插值法计算插值――简介

在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫・拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。

Matlab编程实现拉格朗日插值法计算插值――工具/原料
Matlab2014aMatlab编程实现拉格朗日插值法计算插值――方法/步骤

Matlab编程实现拉格朗日插值法计算插值 1、

理解多项式插值的内涵,明确拉格朗日插值多项式的性质和计算方法。


Matlab编程实现拉格朗日插值法计算插值 2、

打开Matlab软件。


Matlab编程实现拉格朗日插值法计算插值 3、

新建M文件。

1:在Matlab中点击左上角的New加号按钮。

2:在弹出的下拉菜单中选择Script。

(Ps:可直接使用快捷键【Ctrl+N)】

matlab 拉格朗日插值 Matlab编程实现拉格朗日插值法计算插值

软件就会弹出M函数编辑界面。)

Matlab编程实现拉格朗日插值法计算插值_约瑟夫・拉格朗日


Matlab编程实现拉格朗日插值法计算插值 4、

编写拉格朗日插值函数程序:

function yh=lagrange (x,y,xh)

n = length(x);

m = length(xh);

yh = zeros(1,m);

c1 = ones(n-1,1);

c2 = ones(1,m);

for i=1:n

xp = x([1:i-1 i+1:n]);

yh = yh + y(i)*prod((c1*xh-xp'*c2)./(x(i)-xp'*c2));

end

注:该程序只可一次计算实现一个插值计算。可实现多个插值计算的程序如下:

function yh=lagrange(x,y,xh)

n = length(x);

m = length(xh);

x = x(:);

y = y(:);

xh = xh(:);

yh = zeros(m,1);

c1 = ones(1,n-1);

c2 = ones(m,1);

for i=1:n,

xp = x([1:i-1 i+1:n]);

yh = yh + y(i) * prod((xh*c1-c2*xp')./(c2*(x(i)*c1-xp')),2);

end


Matlab编程实现拉格朗日插值法计算插值 5、

保存文件,需要特别注意的是文件名应与函数名相同。


Matlab编程实现拉格朗日插值法计算插值 6、

在Matlab中输入给定数据x和y,以及需要计算插值点的值xh。如图输入后按【Enter】键,载入数据。

Matlab编程实现拉格朗日插值法计算插值_约瑟夫・拉格朗日


Matlab编程实现拉格朗日插值法计算插值 7、

调用编写的程序,并运行。在Matlab的命令窗口输入【lagrange (x,y,xh)】按【Enter】键即可得到拉格朗日插值函数计算的插值。

Matlab编程实现拉格朗日插值法计算插值――注意事项
本文笔者所用的Matkab版本为2014a,其他版本均适用,只是2012a及以前的版本界面差别较大。

  

爱华网本文地址 » http://www.aihuau.com/a/8104200103/207928.html

更多阅读

插值法内插法、试误法 简要说明 技术权益的简要说明

学财务管理的时候,很多地方都会用到插值法,基本每本财务管理的书都会有,比较常用的有求适当的折现率使某项目的现金净流量为零、求债券的到期收益率、内含报酬率的计算等。其基本原理为在抛物线f(x)在很短的一段里为近似直线y=a+bx,再

ArcGIS插值图分类面积统计_阳光温热 arcgis反距离权重插值

在用ArcGIS进行地统计分析时,遇到一种情况,在插值完成后,想根据插值等级进行面积的统计分析,通过查资料解决了这个问题,现在分享一下。1.根据已有数据进行克里金插值。具体的探索性数据分析、数据变换、交叉验证之类的步骤,可以查找网上其

2.2 完整系的拉格朗日方程 欧拉 拉格朗日方程

用广义坐标表出的动力学方程称为拉格朗日方程,可以直接由牛顿第二定律导出。(1)达朗贝尔方程设受约束的质点系中质点i所受的主动力和约束力分别为和 ,位矢为 ,由牛顿第二定律有给质点i以虚位移 ,得对整个质点系在理想约束条件下,有(2.6)

声明:《matlab 拉格朗日插值 Matlab编程实现拉格朗日插值法计算插值》为网友超能分享!如侵犯到您的合法权益请联系我们删除