当前位置：文江博客话题详情

四次方程求解器的实现不起作用

发布于 2024-10-12 08:35:57 字数 434 浏览 3 评论 0原文

我的四次方程求解器的实现出了什么问题？
这是我的代码 GitHub 上。我按照以下步骤操作：http://www.1728.com/quartic2.htm
事实上，真正的实现从第 271 行开始，我在其中创建了单子多项式。

如果我尝试使用具有 4 个实根的多项式，则效果很好（例如 3x^4 + 6x^3 - 123x^2 - 126x + 1,080），否则给出错误的根。

谢谢，
鲁比克

我将函数命名为 __quartic 因为它仍在开发中

原文

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

瞎闹 2024-10-19 08:35:57

如果这不是你的三次解算器的问题（正如phkahler建议的那样），那么我打赌这是这些行中int除法的情况：

f = float(c - (3*b**2 / 8))
g = float(d + (b**3/ 8) - (b*c / 2))
h = e - (3*b**4 / 256) + (b**2 * c / 16) - (b*d / 4)

如果你使用的Python版本低于3.0并且没有导入与 __future__ 相除，那么这里可能会遇到重大问题。例如，如果 b = 3，则 (3*b**2/8) 通过 int 除法等于 3，而不是正确的浮点除法值 3.375。纠正这个问题的最佳方法可能是在开始计算之前将 a、b、c、d 和 e 转换为浮点数。在进行此更正后，在对 f 和 g 进行所有计算后，您实际上没有必要将其转换为浮点型，而且它们现在对您并没有多大帮助。

有关此划分问题及其历史的更多信息，请参阅：

http://www.python。 org/dev/peps/pep-0238/

http://python-history.blogspot.com/2009/03/problem-with-integer-division.html

If it isn't a problem with your cubic solver (as phkahler suggested), then I'm betting that it is a case of int division in these lines:

f = float(c - (3*b**2 / 8))
g = float(d + (b**3/ 8) - (b*c / 2))
h = e - (3*b**4 / 256) + (b**2 * c / 16) - (b*d / 4)

If you are using a version of Python that is less than 3.0 and not importing division from __future__, then you can have major problems here. For instance, if b = 3, then (3*b**2/8) is equal to 3 by int division rather than the correct float division value of 3.375. Probably the best way to correct this is to convert a,b,c,d and e to floats before you start doing calculations. Your conversion to float after all of the calculations for f and g are really unnecessary after you've made this correction and they aren't really helping you much now.

For more information about this division issue and its history see:

http://www.python.org/dev/peps/pep-0238/

http://python-history.blogspot.com/2009/03/problem-with-integer-division.html

回复收藏 0 原文