请解释一下这个拉格朗日插值的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
随机阅读
- 求教Wittig-Horner反应中氢化钠使用注意事项 用氢化钠的注意事项
- 娄桥街道下斜城中村F06 温州有哪些镇
- 行人止步急切盼望(打一个字) 行人止步(打一字)字谜
- 当兵政审 我要去当兵 武装部要学校盖章 盖什么章呢 找那个部门盖呢 当兵需要章学校没人怎么办
- 请问口吃真的可以看好吗,我儿子现在有口吃现象,我好急 口吃有看好的吗
- 葡萄皮肉分离是不是超新鲜了 巨峰葡萄的皮肉不易分离是怎么回事?
- 上海什么景点有玻璃栈道 眼镜男走玻璃栈道
- 前天去医院,医生给开了个叫空气波治疗的,说是压力疗法,有没有谁知道什么是空气波啊?我看那个机 医用空气压力波治疗
- 支付宝交电费显示此账单已经缴纳 你好用支付宝交电费为什么显示此帐单已缴纳怎么回事啊
- 全国十大小商品批发市场都有哪呀? 中国轻纺城副食品市场
- 不足 对生产经营工作涉猎太少 如何提升自身职业素养
- 孕妇梦见各式各样的she 孕妇梦见各式各样的围巾
- 我国政府为了保障食品安全采取了哪些措施 运动会食品安全保障工作情况
- 商城县武桥村派出所 我家住在信阳商城武桥!我想问广东佛山这边还有谁也住在那里?如题 谢谢了
- 北京永乐进京检查站电话 去北京办进京证的地方在哪
- 英雄联盟日服的周免什么时候更新 请问,英雄联盟日服什么时候开
- 有人喜欢笙这个乐器的么? 丰巢快递柜 水语人家
- 气化炉灰渣的含碳率 麻烦哪位老司机!谁能帮忙!!煤炭灰熔融性是什么意思?
- 为什么说资本家剥削了工人? 企业家雇佣资本理论
- 对照检查材料 三讲一比对照材料