python – 在数值上解决集成限制?
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了python – 在数值上解决集成限制?,小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含1066字,纯文字阅读大概需要2分钟。
内容图文
![python – 在数值上解决集成限制?](/upload/InfoBanner/zyjiaocheng/777/035eebe4d4174d6198243749e33c52af.jpg)
我试图在两个函数的曲线下找到一个区域,使得该区域的面积等于给定的数量m.例如,通过集成两个功能:
f(x) = 3x - 3x^2
g(x) = x
象征性地从a到b,并将每个等于m(例如m = 1/3).可以求解系统以找到a和b的值:
1.5 b^2 - b^3 - 1.5 a^2 + a^3 = 1/3
0.5 b^2 - 0.5 a^2 = 1/3
我得到的正解是:a = 0.364823,b = 0.894294
现在,我的问题是:如果f和g作为向量(或python中的数组)而不是符号函数给我,我怎么能找到a和b?我知道如何使用trapz()来查找向量的区域,但我不知道是否有办法用它来在数字上找到上述问题中的a和b.
谢谢!
解决方法:
你可以尝试sympy
>>> from sympy import *
>>> x,a,b = symbols('x a b')
>>> fx = "3*x - 3*x**3"
>>> gx = "x"
>>> m = 1/3
>>> int_fx = integrate(fx, (x,a,b))
3*a**4/4 - 3*a**2/2 - 3*b**4/4 + 3*b**2/2
>>> int_gx = integrate(gx, (x,a,b))
-a**2/2 + b**2/2
>>> solve([Eq(int_fx, m), Eq(int_gx, m)],(a, b))
[(-0.577350269189626, -1.00000000000000), (-0.577350269189626, 1.00000000000000), (0.577350269189626, -1.00000000000000), (0.577350269189626, 1.00000000000000)]
顺便说一下,我很难理解你是如何设法获得一个明确的整体解决方案的
1.5 b^2 - b^3 - 1.5 a^2 + a^3 = 1/3
对于
f(x) = 3x - 3x^3
内
[a,b]
内容总结
以上是互联网集市为您收集整理的python – 在数值上解决集成限制?全部内容,希望文章能够帮你解决python – 在数值上解决集成限制?所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。