如何在python中使用for循环创建NXN矩阵

Question

我是python的新手，并试图为CLT（中心限制理论）创建一个骰子示例如下

import numpy as np
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt

SampleSize=1000000
NumberofDice=2

Dice=np.zeros([SampleSize,NumberofDice])

for i in range(NumberofDice+1):
 Dice[:,i]=np.random.randint(1,7,SampleSize)

但我说错误

  

追踪（最近一次调用最后一次）：文件“”，第11行，中    IndexError：索引2超出轴1的大小为2

的范围

我不知道如何处理它，基本上我想创建一个2 X 1000000矩阵，每行都是一个独立的滚动骰子过程。

有人可以帮忙吗？

Answer 1

range(NumberofDice+1)只删除+1，numpy索引为0到n-1，而不是1到n，range也遵循该惯例。

顺便说一句，你也可以直接这样做：

Dice = np.random.randint(1,7,size=(SampleSize,NumberofDice))