计算方法实验,matla实现拉格朗日插值法。

我觉得自己这一篇小报告写的还可以,记录一下。

一.实验内容:
用matlab对应的命令以及提供的数组实现lagrange插值法。

二.实验步骤:1.编写函数,函数名为ldj(我的名字简写)。2.输入已知点并在图中用红色+将它们标注出来,为了后面将计算出的点也标注出来,与它们做比对,键入hold on。3.很明显,已知的数组为3,对应于lagrange插值法中的n=2,故使用抛物插值公式。

三.实验结果以及实验改进:从图1可以看出由数组xi[95 1 8] 通过sqrt函数并不能准确预测我的学号47对应的sqrt近似值,结果误差过于大(9.5480的平方整数部分为91,与47相去甚远)。我于是主要误差来源,通过翻查课本26页以及思考发现:结论1.当已知数组的sqrt函数不能较为准确算出的话,误差会越来越大(可能往过大也可能往过小发展)。于是将已知数组改为xi[100 121 144],可以发现预测值得到了很大的改善,见图2(7.1071平方的整数部分为50),得到该结果之后总觉得这个结果差强人意。于是经过一番思考以及实践之后得到了结论2.当满足结论1之后,若已知数组与待求数距离越近,lagrange插值估计数也就越接近真实值。见图3,若6.8548平方的整数部分为46,且保留四位小数的结果竟然达到了46.988!这与47之间的差距已经如此之小。

四.实验收获:通过以上实验,我至少可以得到以下两个结论:对于lagrange抛物插值来说,结论1.当已知数组的函数不能较为准确算出的话,预测值误差会越来越大(可能往过大也可能往过小发展)。结论2.当满足结论1之后,若已知数组与待求数距离越近,lagrange插值估计数也就越接近真实值。

不知道是不是代码写错了出现这些问题。
图片和m文件后面我上传,也可以问我要,虽然脚步写的挺拉胯。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值