数值分析实验报告Newton插值法
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
山东师范大学数学科学学院实验报告
2) 选取插值节点,计算插值节点画出cosx 在[0,1.2]上图像,如图3-1,并标出插值算法所计算的点。
>> X=[0 0.2 0.4 0.6 0.8 1.0 1.2;1 0.980067 0.921061 0.825336 0.696707 0.540302 0.362358]; >> x1=0:0.1:1.2; >> format long
>> [A,Pn]=Newton(X,x1) >> x2=0:0.01:1.2; >> plot(x2,cos(x2)) >> hold on
>> plot(x1,Pn,'*')
>> legend('cosx','Newton 插值');
用Newton 插值法计算出x1=0:0.1:1.2处的值为表3-2:
表3-2 x
0 0.1 0.2
0.3 0.4 0.5 0.6 Pn 1 0.995005 0.980067 0.955337 0.921061 0.877583 0.825336
x 0.7 0.8 0.9 1 1.1 1.2
Pn 0.764843 0.696707 0.62161
0.540302 0.453596 0.362358
计算插值节点画出cosx 在[0,1.2]上图像,并标出插值算法所计算的点。
图3-1
总结
(1)使用Newton 插值法编程,在计算差商表时应注意把),,,(10k x x x f 到
)
,,,(1n k n k n x x x f +--),,2,1(n k =都计算出来; (2)Newton 插值法的误差大约为-7101⨯。