Python sympy 里面解方程怎么会出现这种情况?要窒息了 - V2EX
V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
请不要在回答技术问题时复制粘贴 AI 生成的内容
movq
V2EX    程序员

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

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

    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     Solana     2650 人在线   最高记录 6679       Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 32ms UTC 02:26 PVG 10:26 LAX 19:26 JFK 22:26
    Do have faith in what you're doing.
    ubao msn snddm index pchome yahoo rakuten mypaper meadowduck bidyahoo youbao zxmzxm asda bnvcg cvbfg dfscv mmhjk xxddc yybgb zznbn ccubao uaitu acv GXCV ET GDG YH FG BCVB FJFH CBRE CBC GDG ET54 WRWR RWER WREW WRWER RWER SDG EW SF DSFSF fbbs ubao fhd dfg ewr dg df ewwr ewwr et ruyut utut dfg fgd gdfgt etg dfgt dfgd ert4 gd fgg wr 235 wer3 we vsdf sdf gdf ert xcv sdf rwer hfd dfg cvb rwf afb dfh jgh bmn lgh rty gfds cxv xcv xcs vdas fdf fgd cv sdf tert sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf shasha9178 shasha9178 shasha9178 shasha9178 shasha9178 liflif2 liflif2 liflif2 liflif2 liflif2 liblib3 liblib3 liblib3 liblib3 liblib3 zhazha444 zhazha444 zhazha444 zhazha444 zhazha444 dende5 dende denden denden2 denden21 fenfen9 fenf619 fen619 fenfe9 fe619 sdf sdf sdf sdf sdf zhazh90 zhazh0 zhaa50 zha90 zh590 zho zhoz zhozh zhozho zhozho2 lislis lls95 lili95 lils5 liss9 sdf0ty987 sdft876 sdft9876 sdf09876 sd0t9876 sdf0ty98 sdf0976 sdf0ty986 sdf0ty96 sdf0t76 sdf0876 df0ty98 sf0t876 sd0ty76 sdy76 sdf76 sdf0t76 sdf0ty9 sdf0ty98 sdf0ty987 sdf0ty98 sdf6676 sdf876 sd876 sd876 sdf6 sdf6 sdf9876 sdf0t sdf06 sdf0ty9776 sdf0ty9776 sdf0ty76 sdf8876 sdf0t sd6 sdf06 s688876 sd688 sdf86