Matlab编程实现拉格朗日插值法计算插值
在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。
操作方法
(01)理解多项式插值的内涵,明确拉格朗日插值多项式的性质和计算方法。
(02)打开Matlab软件。
(03)新建M文件。1:在Matlab中点击左上角的New加号按钮。2:在弹出的下拉菜单中选择Script。(Ps:可直接使用快捷键【Ctrl+N)】软件就会弹出M函数编辑界面。)
(04)编写拉格朗日插值函数程序: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:nxp = 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
(05)保存文件,需要特别注意的是文件名应与函数名相同。
(06)在Matlab中输入给定数据x和y,以及需要计算插值点的值xh。如图输入后按【Enter】键,载入数据。
(07)调用编写的程序,并运行。在Matlab的命令窗口输入【lagrange (x,y,xh)】按【Enter】键即可得到拉格朗日插值函数计算的插值。
特别提示
本文笔者所用的Matkab版本为2014a,其他版本均适用,只是2012a及以前的版本界面差别较大。
-
电脑上怎样申请QQ号
很多人在打开电脑的时候,都会习惯登录QQ跟别人聊天,如果没有QQ账号的话,可以在电脑上免费申请QQ号,或者之前的QQ号是别人给的,自己不喜欢,也可以重新申请一个新的QQ号,在这里,跟大家分享一下在电脑上如何免费申请QQ号的方法。电脑上免费申请QQ号的步骤(01)在电脑上下载和...
-
QQ浏览器如何安装应用插件
QQ浏览器如何安装应用插件,浏览器的应用插件的安装对我们来说是非常方便的。下面以安装护眼导航应用为例,介绍如何安装浏览器应用插件。操作方法(01)打开QQ浏览器,点击右上角的“”+“”如图所示(02)进入了安装应用的首页,如图所示(03)比如安装一个护眼的应用,鼠标停留在...
-
EPSON R230打印机清零软件无法使用的解决方法
当EPSONR230打印机两个红灯闪烁,并提示“打印机的墨管已到它们的使用寿命”时就需要使用清零软件进行处理,可是在我们使用清零软件的过程中,总是出现这样或那样的问题。那么我们在遇到这些问题时该如何解决呢?操作方法(01)首先当打开清零软件可执行文件时,系统会弹出...
-
在Word文档中怎么样进行简繁体转换?
字体有繁体和简体之分,而Word文档主要的功能是记录保存文本文字,可能由于区域不一样,有的喜欢用简体有的喜欢用繁体,为了方便查看,就需要用到简繁体的相互转换了。操作方法(01)方法一:比较常用的,就是用Word文档自带的简繁体转换功能,找到审阅——简转繁——便可将选择的...