设为首页收藏本站如何选择编译器？在线F77转F90 繁体浏览

切换到窄版开启辅助访问

Fortran Coder»论坛 › 研讨区 › Fortran语法讨论 › 拉格朗日插值法，跑Aitken算法时出现NaN与Infinity ...

返回列表发新帖

查看: 7299|回复: 3

上一主题

下一主题

[数值问题] 拉格朗日插值法，跑Aitken算法时出现NaN与Infinity

已解决

3 帖子	1 主题	0 精华

F 币: 25 元
贡献: 14 点

发消息

电梯直达

跳转到指定楼层

楼主

发表于 2020-2-22 13:50:52 | 只看该作者 |只看大图回帖奖励

回帖奖励

|倒序浏览 |阅读模式

aitken.txt (388 Bytes, 下载次数: 2)

  如附件和图1，我编了aitken算法后，人工验证后觉得原理没有问题，运行也没有bug。
  但如图2，跑完后发现，无论n，l，x，y怎么取值， i=2 k=1 时算出的 y(k+i) 总会出现nan或infinity 。
  图2是拿了个简单的线性插值验证下，发现没有除零问题。到底是为什么呢？实在想不通，请问有大佬解惑吗？

分享到: 微信

收藏0 点赞0 点踩0

回复

使用道具举报

实名认证

2033 帖子	12 主题	5 精华

臭石头雪球

F 币: 1643 元
贡献: 709 点

发消息

沙发

发表于 2020-2-22 14:27:14 | 只看该作者

y(k+i)=(l-x(k+i))*y(i)/(x(i)-x(k+i))+(l-x(i))*y(k+i)/(x(k+1)-x(i))
^
确定这里是1？？
如果是1的话。当 k=1,i=2 时，分母 (x(k+1)-x(i)) 为 0

回复支持反对

使用道具举报

3 帖子	1 主题	0 精华

F 币: 25 元
贡献: 14 点

发消息

板凳

楼主| 发表于 2020-2-22 16:03:55 | 只看该作者

fcode 发表于 2020-2-22 14:27
y(k+i)=(l-x(k+i))*y(i)/(x(i)-x(k+i))+(l-x(i))*y(k+i)/(x(k+1)-x(i))
...

非常感谢！是我看久眼花没看出来了的问题了，诶人工验证时都是算成k+i，没想到我这个怎么手残写成了k+1.

回复支持反对

使用道具举报

3 帖子	1 主题	0 精华

F 币: 25 元
贡献: 14 点

发消息

地板

楼主| 发表于 2020-2-22 16:14:22 | 只看该作者

幼红月发表于 2020-2-22 16:03
非常感谢！是我看久眼花没看出来了的问题了，诶人工验证时都是算成k+i，没想到我这个怎么手残写成了k+1. ...

问题已解决

回复支持反对

使用道具举报

返回列表发新帖

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-27 15:43

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表