请解释一下这个拉格朗日插值的matlab函数实现,希望详细些。。各个参数的意义也解释一下啊。。一头雾 教材上学过的,建议跟着教材看就能明白了。做了一个测试,希望有所帮助。代码:%用matlab编写拉格朗日插值算法的程序,并以下面给出的函数表为数据基础,在整个插值区间上采用拉格朗日插值法计算f(0.6),写出程序源代码,输出计算结果x-2.15-1.00 0.01 1.02 2.03 3.25y 17.03 7.24 1.05 2.03 17.06 23.05function main()clc;x=[-2.15-1.00 0.01 1.02 2.03 3.25];y=[17.03 7.24 1.05 2.03 17.06 23.05];x0=0.6;f=Language(x,y,x0)function f=Language(x,y,x0)求已知数据点的拉格朗日插值多项式已知数据点的x坐标向量:x已知数据点的y坐标向量:y插值点的x坐标:x0求得的拉格朗日插值多项式或在x0处的插值:fsyms t l;if(length(x)=length(y))n=length(x);elsedisp('x和y的维数不相等!');return;检错endh=sym(0);for(i=1:n)l=sym(y(i));for(j=1:i-1)l=l*(t-x(j))/(x(i)-x(j));end;for(j=i+1:n)l=l*(t-x(j))/(x(i)-x(j));end;h=h+l;endsimplify(h);if(nargin=3)f=subs(h,'t',x0);计算插值点的函数值elsef=collect(h);f=vpa(f,6);将插值多项式的系数化成6位精度的小数end结果:f=0.0201>;>;拉格朗日插值法matlab程序 syms t;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=i+1:n)l=l*(t-x(j))/(x(i)-x(j));计算拉格朗日基函数end;f=f+l;计算拉格朗日插值函数simplify(f);化简if(i=n)if(nargin=3)f=subs(f,'t',x0);计算插值点的函数值elsef=collect(f);将插值多项式展开f=vpa(f,6);将插值多项式的系数化成6位精度的小数endend
随机阅读
- 请问新化县城有图书馆吗 新化南门药店大桥店
- 从上沙到彩田南路坐什么车 下沙到彩田路海鹰大厦
- 永不退缩下载链接 永不退缩2双语版百度云
- 横店? 横店四海中学还招生吗
- 个人所得税的专项扣除项目及标准是什么 综合所得税扣除标准
- 胜利雅苑隔壁的汇丰 孟晚舟27号会被释放吗?
- 52路(北京西站-平乐园),线路图 北京西站到平乐园是做52路吗
- 蝙蝠侠前传黑暗骑士崛起猫女 蝙蝠侠前传3黑暗骑士崛起结局什么意思
- 早睡早起对青少年的发育有哪些好处 小孩早睡早起的好处
- 邹城最繁华的街道是哪里?商业圈在哪里?女人逛街的地方? 邹城购物街
- 音响喇叭的小球瘪了怎么办? 调音台防尘罩
- 新华联国际大厦在哪里 西安市南门外领先国际大厦
- 杞菊地黄丸浓缩丸是什么 杞菊地黄丸大蜜丸好还是浓缩丸好
- 朝阳旺角板楼 住房公积金可以代替现金买房子吗?
- 谁能告诉我东莞东城哪家海鲜酒楼既好吃价格又合理呢? 东莞食为先海鲜大酒楼
- 有哪些 100 元以下,实用且贴心,用了之后相见恨晚甚至感动出泪的东西? 花粥格子衫牛仔裤
- 哪些演员演谁都像自己? 我要当八路军演员表小
- 唢呐吹不响是为什么? 赣榆话木头
- 老年服务中心的存在 如何办理老年人养老服务中心的手续
- 宜宾天池一带要修中山小学吗? 天池路小学