我需要编写一个程序来生成Cauchy分布的随机实现
{1美元^
具有空位置和单位比例。在
我还需要做一个介于-5和5个格之间的直方图,对于1000个点的随机实现,和理论曲线确保它们具有相同的单位。在
我计算了Cauchy分布的累积分布函数
:
我写了下面的python代码:
在
在from __future__ import division
import scipy
import random
import matplotlib.pyplot as plt
import numpy as np
import math as m
valuesX = []
for q in range(1000):
R = random.random()
x = m.tan(m.pi*(R-0.5)) #Cumulative Function
valuesX.append(x)
z = np.linspace(-10,10,1000)
y = 1/(m.pi*(1+z**2)) #Theoretical Cauchy
plt.plot(y,z)
plt.hist(valuesX, bins = 50, range = [-5,5], normed=True)
我的结果是:
我不知道这是否可以接受,因为我正在根据概率密度函数绘制离散值(随机实现)。我怎么能比较他们两个呢?因为我需要找到上面曲线的分数差,并确定两条曲线之间的整体均方根偏差,作为随机实现大小的函数。在