如何解决python scipy.integrate反向ppf
你好,我想给我一个特定的数值整数的x值是什么。
到目前为止,我正在使用scipy.integrate并且工作正常,但是给定一个整数值,有一种方法可以知道给出该结果的x是什么?
让我说我有一个函数f(x)= | 2-x | 1
我想知道给我0.25(第一个四分位数)的x值是什么。
scipy.stats中的正态分布是norm.ppf(),在这种特殊情况下,我将其用于PDF(概率密度函数),但它可以是任何整数。
感谢和问候
解决方法
我使用binary search在logarithmic time中找到答案(非常快)。您也可以运行code below online here。
import math,scipy.integrate,scipy.optimize
def f(x):
return math.sin(x)
a,b = 0,10
integral_value = 0.5
res_x = scipy.optimize.bisect(
lambda x: scipy.integrate.quad(f,a,x)[0] - integral_value,b
)
print(
'found point x',res_x,',integral value at this point',scipy.integrate.quad(f,res_x)[0]
)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。