拉格朗日插值逐次线性插值法.pptVIP

下载本文档

9
0
约1.97千字
约 16页
2017-07-05 发布于四川
举报
版权申诉

拉格朗日插值逐次线性插值法.ppt

1、本文档共16页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

第二章插值与拟合例2.1 已知 , , 解: 这里x0=100，y0=10，x1=121，y1=11, 利用线性插值利用线性插值求例2.1* 已知 , , 解: 这里x0=100，y0=10，x1=121，y1=11, x2=144，y2=12，利用抛物线插值公式利用抛物线插值求例1：已知x0=100, x1=121, x2=144, 求在x=115时的近似值。则（与利用抛物线插值得到的结果一样） 2.1.3 逐次线性插值法为克服这一缺点，通常可用逐次线性插值方法求得高次插值。例如在例2.1-2.1*中：则当k=0时为线性插值。 k=1时插值节点为三点，公式为计算时可由k=0到k=n-1逐次求得所需的插值多项式。计算过程如下公式(**)也可以改成下面的计算公式称之为NEVILLE（列维尔）算法，计算过程如下例：已知分别利用 sin x 的1次、2次 Lagrange 插值和牛顿插值计算 sin 50?。 Function y=chazhi1(x0,y0,x1,y1,x)%拉格朗日插值法 y=(x-x1)/(x0-x1)*y0+( x-x0)/(x1-x0)*y1 Function y=chazhi1(x0,y0,x1,y1,x)%牛顿插值法 y=y0+(y1-y0)/(x1-x0)*( x-x0) function chazhi1 %一次插值 x0=pi/6;y0=1/2; x1=pi/4;y1=1/sqrt(2); x=pi/180*50; y=chazhi12(x0,y0,x1,y1,x) error=y-sin(x) function y=chazhi11(x0,y0,x1,y1,x)%拉格朗日插值法 y=(x-x1)/(x0-x1)*y0+( x-x0)/(x1-x0)*y1; function y=chazhi12(x0,y0,x1,y1,x)%牛顿法 y=y0+(y1-y0)/(x1-x0)*( x-x0); function chazhi2 %二次插值 x0=pi/6;y0=1/2; x1=pi/4;y1=1/sqrt(2); x2=pi/3;y2=sqrt(3)/2; x=pi/180*50; y=chazhi22(x0,y0,x1,y1,x2,y2,x)%选择插值函数 error=y-sin(x) function y=chazhi21(x0,y0,x1,y1,x2,y2,x)%拉格朗日插值法 y=((x-x1)*(x-x2))/((x0-x1)*(x0-x2))*y0+((x-x0)*(x-x2))/... ((x1-x0)*(x1-x2))*y1+((x-x0)*(x-x1))/((x2-x0)*(x2-x1))*y2; function y=chazhi22(x0,y0,x1,y1,x2,y2,x)%牛顿插值法 y=y0+(y1-y0)/(x1-x0)*( x-x0)+(x-x0)*(x-x1)*... ((y2-y1)/(x2-x1)-(y1-y0)/(x1-x0))/(x2-x0); Function y=chazhi1(x0,y0,x1,y1,x2,y2,x)%拉格朗日插值法 y=((x-x1)*(x-x2))/((x0-x1)*(x0-x2))*y0+ ((x-x0)*(x-x2))/((x1-x0)*(x1-x2))*y1+… ((x-x0)*(x-x1))/((x2-x0)*(x2-x1))*y2 Function y=chazhi1(x0,y0,x1,y1, x2,y2,x)%牛顿插值法 y=y0+(y1-y0)/(x1-x0)*( x-x0)+… (x-x0)*(x-x1)*((y2-y1)/(x2-x1)-(y1-y0)/(x1-x0))/(x2-x0) 例1:给定f(x)=lnx的数据表 xi 2.20 2.40 2.60 2.80 3.00 f(xi) 0.78846 0.87547 0.95551 1.02962 1.09861 已经测得在某热敏电阻在不同温度下的电阻如下：电阻（欧姆） 600 300 100 50 10 5 温度（摄氏度）10 30 50 70 90 100 根据这些数据，希望合理地估计出其它温度（如25摄氏度，40摄氏度…）时的电阻第二章插值与拟合