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]
)