各位老师，有一个gumbel分布，概率分布函数为： y = exp( -exp( -a * ( x - b ) ) )

已经估计出了参数a和b。

请问，它的概率密度函数应该怎么画？

谢谢

做一个数值微分就可以了。

  import numpy as np
import pylab as pl

def diff(f, x, dx=1e-10):
    return (f(x) - f(x-dx))/dx
    
def gumbel(x):
    a = 1
    b = 1
    return np.exp(-np.exp(-a*(x-b)))
    
if __name__ == "__main__":    
    x = np.linspace(-2, 2, 1000)    
    y1 = gumbel(x)
    y2 = diff(gumbel, x)
    
    pl.plot(x,y1)
    pl.plot(x,y2)
    
    pl.show()
 

有道理。

每次提问，RY老师都能及时回复，非常感谢