Fortran学习21:数值计算方法4

本文介绍了数值计算中的两种基本插值方法——拉格朗日插值和牛顿插值,并通过Fortran代码实现了这两种方法,同时利用Python进行了结果可视化。

今天要学习的就是插值方法了,这一节量很大的,我感觉今天一天搞不完的,能学到那就算哪吧。

数值计算方法我建议可以取读一本书,高尚、秦斌合写的计算方法,算是对母校的书籍进行一个推广吧。

1、拉格朗日插值方法:

经过n个数据点,生成最高项为n-1的多项式,

假设:又(X1,Y1)、(X2,Y2)、(X3,Y3)......(XN,YN)个数据连,经过拉格朗日插值法获得的多项式就是:

g(x)=\frac{(x-x_{2})(x-x_{3})(x-x_{4})...(x-x_{n})}{(x_{1}-x_{2})(x_{1}-x_{3})(x_{1}-x_{4})...(x_{1}-x_{n})}y_{1}+ \frac{(x-x_{1})(x-x_{3})(x-x_{4})...(x-x_{n})}{(x_{2}-x_{1})(x_{2}-x_{3})(x_{2}-x_{4})...(x_{2}-x_{n})}y_{2}+....+ \frac{(x-x_{1})(x-x_{2})(x-x_{3})...(x-x_{n})}{(x_{n}-x_{1})(x_{n}-x_{2})(x_{n}-x_{3})...(x_{n}-x_{n-1})}y_{n}

我无语的地方是这个用的是作者自己搞的sgl库来 绘制图形,emm...咋说呢,我没有光驱,那就算了,直接写代码吧,后面用matlab来作图。

module interpolate_utility
    implicit none
    type point
        
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值