请解释一下这个拉格朗日插值的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
随机阅读
- 无氧运动肌肉量减少 想减少身体的脂肪含量 同时增强肌肉力量 应该以有氧运动为主还是无氧运动? 运动强度 时间 多少合适?
- 侠客风云传武林任清 侠客风云传情圣路线能攻略夜叉吗
- 怎样看建筑总平面图 乡村道路工程总平面布置图
- 青海省民和县巴州镇的经度,纬度是多少 青海省民和县有几个镇
- 怎样去除渗水井里得油 渗水井不渗水了,怎样清洗民用渗水井里的油污啊
- 禁毒工作主要成效 为实现禁毒工作目标,我国应加强哪些方面的工作?
- 赵欣怡市西 如何评价 2016 清华招生宣传片 iTsinghua?
- 六氟化硫微泡输卵管造影的危害? 六氟化硫造影
- 常州遥观镇碧水云庭房价 常州遥观观澜云庭附近有没有好的英语辅导班
- 北半球中纬度西风分布纬度范围
- 大海战2 法国防空炮 是炮手么???? 那我的 防空兵 没用了???? 大海战 法国 防空炮手 转职
- 在广州找了租屋 广州市天河区东圃联合社区9号公寓
- 不同螺旋性的椭圆偏振光 线偏振光垂直入射到四分之一波片,偏振方向于四分之一波片所成角度的不同,如何影响形成圆偏振光还是椭圆偏振光的呢?是方向还是振幅在影响呢?
- 路从今夜白的大结局是什么 魏子健和路悠言是什么电视剧
- 哈尔滨民南街 哈尔滨桥南街二十号在哪
- 衡量电力系统电能质量的三个重要指标是什么? 供电质量的主要指标
- 取保候审意味着什么? 取保候审是因为什么
- 走河堤路能到兴平吗 从周至到兴平走河堤路从兴户大桥能走吗?
- 墅外商铺周边环境怎么样?生活便利吗? 康桥永乐街商铺如何
- 山药和豌豆可以同吃吗 豌豆花生牛肉可以一起煮粥吗