问题描述
我正在创建尝试从我创建的Pareto分布创建新的密度函数:
import numpy as np
import matplotlib.pyplot as plt
x_m,alpha = 2+(1/n_i),3.
#drawing samples from distribution
samples_2 = (np.random.pareto(alpha,10000) + 1) * x_m
count,bins,_ = plt.hist(samples_2,10000,normed=True)
fit = alpha*x_m**alpha / bins**(alpha+1)
plt.plot(bins,max(count)*fit/max(fit),linewidth=2,color='r')
plt.xlabel('bins',fontsize=15)
plt.ylabel('probability density',fontsize=15)
plt.title('Probability Density Function',fontsize=15)
plt.grid(b=True,color='grey',alpha=0.3,linestyle='-.',linewidth=2)
plt.rcParams['figure.figsize'] = [8,8]
plt.show()
samples_2
This is the Pareto Distribution i simulate
其中n_1:n_i = random.randint(0,len(t))
,t
是离散随机变量的样本数组。
您注意到sample_2
也是此帕累托样本的集合,但仍然不需要。
我需要创建一个新数组,其中该数组的成分必须是len(samples_2)
的随机n samples_2处的和,从数学上讲:
enter image description here
很明显,我知道这需要一个循环来迭代总和以创建新数组,但是我仍然不能。
我该如何循环以获得所需的?
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)