V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
movq
V2EX  ›  程序员

Python sympy 里面解方程怎么会出现这种情况?要窒息了

  •  
  •   movq · 2022-04-10 09:19:51 +08:00 · 1759 次点击
    这是一个创建于 951 天前的主题,其中的信息可能已经有所发展或是发生改变。

    T9LXPyhj23W4e1E

    这个一元三次方程,sympy 给出的结果竟然带虚数,我还以为这个是这个方程本身有问题

    前面的代码检查了好几遍,都没问题

    后来我把这个方程放到网上解,直接解出来了。。。

    obNWZLu7dfQnMSX

    太窒息了

    这种情况如何解决呢?这么复杂的虚数结果,我也没法看出来是多少啊

    7 条回复    2022-04-10 11:23:35 +08:00
    movq
        2
    movq  
    OP
       2022-04-10 09:32:50 +08:00
    @ericls 这样的结果是:[10*CRootOf(x**3 - 14*x**2 + 6*x + 192, 0), 10*CRootOf(x**3 - 14*x**2 + 6*x + 192, 1), 10*CRootOf(x**3 - 14*x**2 + 6*x + 192, 2)]
    ericls
        3
    ericls  
       2022-04-10 09:38:33 +08:00 via iPhone
    你 600 是不是少乘以了一个 x
    movq
        4
    movq  
    OP
       2022-04-10 09:48:24 +08:00
    @ericls

    我一开始的第二个截图上面是少乘了一个 x ,后来加上去了,然后更新了截图

    print(real_roots(x ** 3 - 140 * x ** 2 + 600*x + 19200))

    [2*CRootOf(x**3 - 70*x**2 + 150*x + 2400, 0), 2*CRootOf(x**3 - 70*x**2 + 150*x + 2400, 1), 2*CRootOf(x**3 - 70*x**2 + 150*x + 2400, 2)]

    这个有办法化简吗?

    我现在找到的解决办法是,nsimplify(x1,tolerance=0.0000001)这样得到实根
    ericls
        5
    ericls  
       2022-04-10 10:08:26 +08:00
    polys.polytools.nroots
    dontmind
        6
    dontmind  
       2022-04-10 10:55:12 +08:00
    Eureka0
        7
    Eureka0  
       2022-04-10 11:23:35 +08:00   ❤️ 1
    https://s2.loli.net/2022/04/10/EdsFUrufGaTt1kK.png

    如图,这个数量级的虚部大概是浮点数运算精度引起的吧,完全可以再写个规则过滤一下
    关于   ·   帮助文档   ·   博客   ·   API   ·   FAQ   ·   实用小工具   ·   959 人在线   最高记录 6679   ·     Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 · 29ms · UTC 21:38 · PVG 05:38 · LAX 13:38 · JFK 16:38
    Developed with CodeLauncher
    ♥ Do have faith in what you're doing.