V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
推荐学习书目
Learn Python the Hard Way
Python Sites
PyPI - Python Package Index
http://diveintopython.org/toc/index.html
Pocoo
值得关注的项目
PyPy
Celery
Jinja2
Read the Docs
gevent
pyenv
virtualenv
Stackless Python
Beautiful Soup
结巴中文分词
Green Unicorn
Sentry
Shovel
Pyflakes
pytest
Python 编程
pep8 Checker
Styles
PEP 8
Google Python Style Guide
Code Style from The Hitchhiker's Guide
flyzy2005
V2EX  ›  Python

Python 怎么解方程啊?

  •  
  •   flyzy2005 · 2018-04-23 10:03:34 +08:00 · 5275 次点击
    这是一个创建于 2412 天前的主题,其中的信息可能已经有所发展或是发生改变。

    用的 sympy 中的 solve,这个库是不是只能解有解的方程?

    p1 = 0.8
    p2 = 1
    r1 = 500
    r2 = 3500
    x = symbols('x')
    print('in')
    s=solve(p1*exp(-x*r2)-p2*exp(-x*r1)+p2-p1, x)
    print(s)
    

    一直输出不了 s。。。 像这种方程用 python 怎么解?

    17 条回复    2018-04-24 06:00:27 +08:00
    est
        1
    est  
       2018-04-23 10:12:13 +08:00
    又来一个廖学峰课后作业系列?
    flyzy2005
        2
    flyzy2005  
    OP
       2018-04-23 10:18:05 +08:00 via iPhone
    @est。。不是 = = 要算个方程
    l1093178
        3
    l1093178  
       2018-04-23 10:31:50 +08:00
    WolframAlpha 了解一下?
    flyzy2005
        4
    flyzy2005  
    OP
       2018-04-23 10:53:12 +08:00 via iPhone
    @l1093178 谢谢 算出来了
    princelai
        5
    princelai  
       2018-04-23 12:14:08 +08:00 via Android
    前两天不是刚问完,sympy
    flyzy2005
        6
    flyzy2005  
    OP
       2018-04-23 12:18:11 +08:00 via iPhone
    @princelai 我用了那个 没解算不出来好像 换一个 x-2 很快算出来
    asdf123101
        7
    asdf123101  
       2018-04-23 12:23:43 +08:00 via Android
    自己写个二分法?知道导数可以用牛顿。这种方程自己写一个很快的。
    sciooga
        8
    sciooga  
       2018-04-23 12:47:07 +08:00
    这样行不行?

    import random

    while True:
    x = random.random()
    然后将 x 带入方程两边如果相等就 break
    reself
        9
    reself  
       2018-04-23 12:54:51 +08:00 via Android
    @sciooga 哈哈猴子算法
    flyzy2005
        10
    flyzy2005  
    OP
       2018-04-23 14:00:04 +08:00 via iPhone   ❤️ 1
    @sciooga 😂😂😂 用了那个 wolframalpha
    pandachow
        11
    pandachow  
       2018-04-23 14:05:28 +08:00
    sympy 啊。。。
    flyzy2005
        12
    flyzy2005  
    OP
       2018-04-23 14:08:04 +08:00 via iPhone
    @pandachow 我就用的这个库啊……你看我代码 卡在 solve 算不出来
    idlerlestat
        13
    idlerlestat  
       2018-04-23 14:22:01 +08:00 via iPhone
    sympy 是给解析解,你这个方程没有解析解啊。
    flyzy2005
        14
    flyzy2005  
    OP
       2018-04-23 15:39:09 +08:00 via iPhone
    @idlerlestat 嗯。。所以我来问应该用什么
    idlerlestat
        15
    idlerlestat  
       2018-04-23 16:58:59 +08:00
    @flyzy2005 numpy 和 scipy.optimize.fsolve 就行
    inflationaaron
        16
    inflationaaron  
       2018-04-23 20:11:15 +08:00 via iPhone
    这种需求用 Mathematica 啊
    20150517
        17
    20150517  
       2018-04-24 06:00:27 +08:00 via Android
    不是我说,这种最好用 tensorflow,做个 variable 然后用模型能帮你算出系数
    关于   ·   帮助文档   ·   博客   ·   API   ·   FAQ   ·   实用小工具   ·   1343 人在线   最高记录 6679   ·     Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 · 25ms · UTC 17:42 · PVG 01:42 · LAX 09:42 · JFK 12:42
    Developed with CodeLauncher
    ♥ Do have faith in what you're doing.