文档详情

用matlab编写拉格朗日插值算法地程序.doc

发布:2018-07-04约7.7千字共9页下载文档
文本预览下载声明
用matlab编写拉格朗日插值算法的程序10[?标签: HYPERLINK /z/Search.e?sp=Smatlabch=w.tag.yjj \t _blank matlab, HYPERLINK /z/Search.e?sp=S%E6%8F%92%E5%80%BC%E7%AE%97%E6%B3%95ch=w.tag.yjj \t _blank 插值算法, HYPERLINK /z/Search.e?sp=S%E7%A8%8B%E5%BA%8Fch=w.tag.yjj \t _blank 程序?] 用matlab编写拉格朗日插值算法的程序,并以下面给出的函数表为数据基础,在整个插值区间上采用拉格朗日插值法计算f(0.6),写出程序源代码,输出计算结果 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 匿名 回答:1 人气:6 解决时间:2011-05-24 19:58 满意答案好评率:83% 做了一个测试,希望有所帮助。代码:% 用matlab编写 HYPERLINK /z/Search.e?sp=S%E6%8B%89%E6%A0%BC%E6%9C%97%E6%97%A5ch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 拉格朗日 HYPERLINK /z/Search.e?sp=S%E6%8F%92%E5%80%BC%E7%AE%97%E6%B3%95ch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 插值算法的程序,并以下面给出的 HYPERLINK /z/Search.e?sp=S%E5%87%BD%E6%95%B0%E8%A1%A8ch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 函数表为数据基础, % 在整个 HYPERLINK /z/Search.e?sp=S%E6%8F%92%E5%80%BCch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 插值区间上采用拉格朗日 HYPERLINK /z/Search.e?sp=S%E6%8F%92%E5%80%BC%E6%B3%95ch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 插值法计算f(0.6),写出程序 HYPERLINK /z/Search.e?sp=S%E6%BA%90%E4%BB%A3%E7%A0%81ch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 源代码,输出计算结果 % 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 function 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) %求已知数据点的拉格朗日插值 HYPERLINK /z/Search.e?sp=S%E5%A4%9A%E9%A1%B9%E5%BC%8Fch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 多项式 %已知数据点的x坐标 HYPERLINK /z/Search.e?sp=S%E5%90%91%E9%87%8Fch=w.search.yjjlinkcid=w.search.yjjlink \t _blank 向量: x %已知数据点的y坐标向量: y %插值点的x坐标: x0 %求得的拉格朗日插值多项式或在x0处的插值: fsyms t l; if(length(x) == length(y)) ??? n = length(x); else ??? disp(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)); ?
显示全部
相似文档