我有一个2d pdf函数，它是一个双变量正态分布，

enter image description here

我绘制它 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 与原始数据集进行比较
enter image description here

问题是， density plot 是不够的，我也想做一个 q-q plot . 但是，由于pdf是集成的结果，是否可以使用Scipy来实现？

我不知道 . 也许一个可能的解决方案是诉诸 Sympy . 由于 bivariate pdf 的参数都是已知的，我们可以生成角度方向的pdf .

python scipy