2、拉格朗日插值的matlab实现2、拉格朗日插值实例：根据下面的数据点求出其拉格朗日插值格式,并计算当x=1。6时y的值。x00。51.01.52。02.53。0y00。47940。81450。99750.90930.59850。1411function[f,f0］=Language(x，y,x0)%求已知数据点的拉格朗日插值多项式%已知数据点的x坐标向量:x％已知数据点的y坐标向量：y％插值点的x坐标：x0%求得的拉格朗日插值多项式:f%x0处的插值：f0symst;if（length(x)==length（...
2、拉格朗日插值实例：根据下面的数据点求出其拉格朗日插值格式,并计算当x=1。6时y的值。x00。51.01.52。02.53。0y00。47940。81450。99750.90930.59850。1411function[f,f0］=Language(x，y,x0)%求已知数据点的拉格朗日插值多项式%已知数据点的x坐标向量:x％已知数据点的y坐标向量：y％插值点的x坐标：x0%求得的拉格朗日插值多项式:f%x0处的插值：f0symst;if（length(x)==length（y））n=length（x）;elsedisp(’x和y的维数不相等!’）；return；end %检错f=0。0；for（i=1:n）l=y(i）;for(j=1:i—1）l=l*(t—x（j）)/(x(i)—x(j）)；end;for(j=i1：n)l=l*（t-x(j))/(x(i)-x(j））； %计算拉格朗日基函数end；f=fl； ％计算拉格朗日插值函数simplify（f）； %化简endf0=subs（f,'t’,x0）； %计算插值点的函数值运行程序；x=0:0.5:3；y=［00.47940。84150。99750。90930.59850。1411];[f,f0]=Language（x，y，1。6) ％计算输出的拉格朗日插值多项式％计算x=1.6时的插值输出值f0运行结果f=-799/3125*t*(t—1)*(t—3/2)*（t-2)*(t—5/2）*(t-3）561/500*t＊（t-1/2）*(t—3/2）＊(t—2）＊（t—5/2）＊（t-3)—133/75*t*（t—1/2)＊（t-1)＊（t-2）*（t—5/2)＊（t-3)3031/2500＊t*(t—1/2)＊（t—1)*（t-3/2）＊（t—5/2)*（t-3)-399/1250*t＊(t—1/2)＊(t-1)＊(t-3/2）＊(t—2）＊（t—3）1411/112500*t*(t-1/2）*(t-1）*（t—3/2)＊（t-2）＊（t-5/2)f0=0.9996
本文档为【2、拉格朗日插值的matlab实现】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483