请解释一下这个拉格朗日插值的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
随机阅读
- 中山横栏去中山影视城远吗 请问,丛中山影视城站到横栏镇六沙村怎么去?
- 文登有多少个镇? 威海市文登区张家产镇地址
- 气压撑杆安装方法 气压杆尺寸
- 形容学生沉着冷静的四字词语 形容沉着冷静的词语
- 郁金香.玫瑰,水仙花,菊花,月季花分别在一天中的哪个时间开放? 水仙郁金香要剥皮吗
- 求异形庇护所 谁有异形庇护所第四季第四集第五集的中文字幕~~~~跪求啊~~等得我想死!~~~
- 永川金科公园王府怎么样?好不好?值不值得买? 金科王府 业主论坛
- 二氧化碳的性质和用途 自然氧化塘工艺
- PID控制中PI,PID分别代表什么意思 pid调节装置中的调节过程
- 三木四火一土零金零水 五行中4金4木1水4火5土是什么意思啊、求大师解释!
- 六安儿童水乐园 六安室内儿童乐园哪家比较好?
- 杉杉来了赵丽颖那么多吻戏陈晓不吃醋吗 陈晓赵丽颖吻戏
- 为什么人血白蛋白很难买到? 贵州民博会志愿者
- 头孢替唑钠32个月小孩一天用两只多不多 头孢替唑钠口服效果
- 我刚买了兰蔻的水分缘眼霜,大牌爱用者都进来分享一下经验吧 水份缘舒缓眼霜好不好
- 耐旱的蔬菜有哪些? 什么耐寒又耐旱蔬菜
- 横道河子网络 哈尔滨到横道河子的火车多少钱?
- 岩板的优点和缺点分别是什么?专业人士给讲讲? 石材全自动水冲击
- 不孕不育专科医院是哪家? 哪家是不孕不育专科医院
- 扬正气比奉献名言警句 谁有有关“扬正气 树新风”的名言或诗句!!!