Scipy：有没有可能用自定义概率密度函数制作qq图？

Question

我有一个二维pdf函数，它是一个二元正态分布

我绘制它angle-wisely，与我得到的原始数据进行比较。也就是说，我将(x,y)转换为(r, theta)，然后进行一些集成。使用Scipy.integrate可以很容易地做到这一点

# Let f be the pdf
x_vals = np.arange(0, max_speed, 0.5)
y_vals_ =[integrate.nquad(f, 
                          [[x_val-0.5, x_val+0.5],[angle_radian-incre_radian/2, angle_radian+incre_radian/2]]) 
          for x_val in x_vals]

我可以将pdf与原始数据集进行比较

问题是density plot不够，我还想做一个q-q plot。但是，既然pdf是集成的结果，那么可以用Scipy来实现吗？在

我对此一无所知。也许一个可能的解决办法是求助于Sympy。由于bivariate pdf的参数都是已知的，我们可以生成相应的角度pdf。在