用mathematics实现拉格朗日插值及其图形化.docx

下载文档

27
0
约2.4千字
约 5页
2017-12-17 发布于河南
举报
版权申诉
保障服务

用mathematics实现拉格朗日插值及其图形化.docx

1、本文档共5页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

两点的线性插值直接运用公式给出x0=0.5;y0=-0.693147;x1=0.6;y1=-0.510826;f[x_]:=(x-x1)/(x0-x1)*y0+(x-x0)/(x1-x0)*y1;f[0.54]-0.620219p1=Plot[f[x],{x,-1,1}];p2=ListPlot[{{0.54,-0.620219}},PlotStyle?{RGBColor[1,0,1],PointSize[0.02]}];Show[p1,p2]Clear[x0,y0];运用for循环进行两点几多点的拉格朗日插值并将其图形化m={{0.5,-0.693147},{0.6,-0.510826}};y0=0;For[i=1,i?Length[m],i++,y=1;For[j=1,j?Length[m],j++, If[i?j,y=y*(x0-m[[j,1]])/(m[[i,1]]-m[[j,1]])]];y0=y0+m[[i,2]]*y];Print[y0];y0/.x0?0.546.93147 (-0.6+x0)-5.10826 (-0.5+x0)-0.620219p1=ListPlot[m,PlotStyle?{RGBColor[1,0,0],PointSize[0.02]}];p2=Plot[y0,{x0,-1,1},PlotStyle?{RGBColor[1,1,0],Thickness[0.01]}];Show[p1,p2]m={{1,3},{2,6},{3,9}};x0=1.5;y0=0;For[i=1,i?Length[m],i++,y=1;For[j=1,j?Length[m],j++, If[j?i,y=y*(x0-m[[j,1]])/(m[[i,1]]-m[[j,1]])]];y0=y0+m[[i,2]]*y];y04.5 m={{1,1.2},{3,2.6},{5,3.8},{6,5.5},{7,6.7},{8,9.8},{10,12}};x0=9;y0=0;For[i=1,i?Length[m],i++,y=1;For[j=1,j?Length[m],j++, If[j?i,y=y*(x0-m[[j,1]])/(m[[i,1]]-m[[j,1]])]];y0=y0+m[[i,2]]*y];y015.2343Clear[x0,y0]m={{1,3},{2,6},{3,9}};y0=0;For[i=1,i?Length[m],i++,y=1;For[j=1,j?Length[m],j++, If[j?i,y=y*(x0-m[[j,1]])/(m[[i,1]]-m[[j,1]])]];y0=y0+m[[i,2]]*y];Print[y0];y0/.x0?1.53/2 (2-x0) (3-x0)+6 (3-x0) (-1+x0)+9/2 (-2+x0) (-1+x0)4.5p3=ListPlot[m,PlotStyle?{RGBColor[1,0,0],PointSize[0.02]}];p4=Plot[y0,{x0,-10,10},PlotStyle?{RGBColor[1,1,0],Thickness[0.01]}];Show[p3,p4]Clear[x0,y0]m={{1,4},{3,10},{5,116},{7,150},{9,201}};y0=0;For[i=1,i?Length[m],i++,y=1;For[j=1,j?Length[m],j++, If[j?i,y=y*(x0-m[[j,1]])/(m[[i,1]]-m[[j,1]])]];y0=y0+m[[i,2]]*y];Print[y0];y0/.x0?101/96 (3-x0) (5-x0) (7-x0) (9-x0)+5/48 (5-x0) (7-x0) (9-x0) (-1+x0)+29/16 (7-x0) (9-x0) (-3+x0) (-1+x0)+25/16 (9-x0) (-5+x0) (-3+x0) (-1+x0)+67/128 (-7+x0) (-5+x0) (-3+x0) (-1+x0)42503/128p5=ListPlot[m,PlotStyle?{RGBColor[1,0,0],PointSize[0.02]}];p6=Plot[y0,{x0,-10,10},PlotStyle?{RGBColor[1,1,0],Thickness[0.01]}];p9=ListPlot[{{x0,y0}},PlotStyle?{RGBColor[1,0,0],PointSize[0.02]}];Show[p5,p6]这是分段线性插值Clear[x0,y0]m={{1,4},{3,10},{5,116},{7,150},{9,201}};y