第三天
使用Pygal模拟掷骰子
0.安装Pygal
在Linux和OS X系统中,应执行的命令类似于下面这样:
pip install --user pygal
在Windows系统中,命令类似于下面这样:
python -m pip install --user pygal
同样也可以:
pip install pygal --user
--user可以不加,但是如果系统提示没有权限,则一定需要加上
Pygal画廊
要了解使用Pygal可创建什么样的图表,请查看图表类型画廊:访问http://www.pygal.org/ ,单击Documentation,再单击Chart types。每个示例都包含源代码,让你知道这些图表是如何 生成的。
创建Die类
下面的类是在模拟掷骰(tou)子:
die.py
from random import randint
class Die():
"""表示一个骰子的类"""
def __init__(self,num_sides=6):
"""骰子默认为六面"""
self.num_sides = num_sides
def roll(self):
"""返回一个位于1和骰子面数之间的随机值"""
return randint(1,self.num_sides)
方法__init__()接受一个可选参数,创建这个类的实例时,如果没有指定任何实参时,面数则会默认为6,如果指定了实参,这个值用于设置骰子的面数。骰子是根据面数命名的,6面的骰子名为D6,8面的骰子名为D8,以此类推。
方法roll()使用函数randint()来返回一个1和面数之间的随机数。这个函数可能返回起始值1、终止值num_sides或这两个数之间的任何整数。
掷骰子
使用这个类来创建图表前,先来掷骰子,将结果打印出来,并检查结果是否合理:
die_visual.py
from die import Die
#创建一个D6
die = Die()
#掷几次骰子,并将结果存储在一个列表中
results = []
for roll_num in range(100): #这里的roll_num可以换成任意一个为使用过的变量,就是是 _ 也可以。(_表示占位符)
result = die.roll()
results.append(result)
print(results)
我们创建了实例die(而且没有给予参数,所以num_sides默认为6),我们用使用for roll_num in range(100)
,使骰子掷100次。并且每一次的结果都存储在了列表results中。
结果大致如下(运算的结果相同的概率是1/6**100那么多)
[4, 6, 5, 6, 1, 5, 6, 3, 5, 3, 5, 3, 2, 2, 1,
3, 1, 5, 3, 6, 3, 6, 5, 4, 1, 1, 4, 2, 3, 6,
4, 2, 6, 4, 1, 3, 2, 5, 6, 3, 6, 2, 1, 1, 3,
4, 1, 4, 3, 5, 1, 4, 5, 5, 2, 3, 3, 1, 2, 3,
5, 6, 2, 5, 6, 1, 3, 2, 1, 1, 1, 6, 5, 5, 2,
2, 6, 4, 1, 4, 5, 1, 1, 1, 4, 5, 3, 3, 1, 3,
5, 4, 5, 6, 5, 4, 1, 5, 1, 2]
从结果我们可以知道,每个出现的数字都在1~6之间,而且每个数字都出现过,这表面每个数字的出现都是随机的。
分析几个
为了分析掷一个D6骰子的结果,我们计算每个点出现的次数:
die_visual.py
from die import Die
#创建一个D6
die = Die()
#掷几次骰子,并将结果存储在一个列表中
results = []
for roll_num in range(1000): #这里的roll_num可以换成任意一个为使用过的变量,就是是 _ 也可以。(_表示占位符)
result = die.roll()
results.append(result)
#分析结果
frequencies = []
for value in range(1,die.num_sides+1):
frequency = results.count(value)
frequencies.append(frequency)
print(frequencies)
这样的话,我们没有显示每一次掷出的结果,但是我们根据统计,把每一个点数出现的次数进行了统计。为了使统计更精确我们这次运行1000次,而且为了存储每个点数出现的次数,我们创建了另一个空列表frequencies,用于存储它们。 通过下面这一部分:
for value in range(1,die.num_sides+1):
frequency = results.count(value)
frequencies.append(frequency)
我们遍历可能的点数(这里为1~6),计算每个点数在results中出现了多少次,将其传递给变量frequency,然后运用for循环和append()函数,将每个值都加到列表frequencies的末尾,最后把这个列表打印出来。
绘制直方图
有了频率后,我们就可以绘制结果的直方图(条形统计图)了。
die_visual.py
import pygal
from die import Die
#创建一个D6
die = Die()
#掷几次骰子,并将结果存储在一个列表中
results = []
for roll_num in range(1000): #这里的roll_num可以换成任意一个为使用过的变量,就是是 _ 也可以。(_表示占位符)
result = die.roll()
results.append(result)
#分析结果
frequencies = []
for value in range(1,die.num_sides+1):
frequency = results.count(value)
frequencies.append(frequency)
#对结果进行可视化
hist = pygal.Bar()
hist.title = "Results of rolling one D6 1000 times."
hist.x_labels = ['1', '2', '3', '4', '5', '6']
hist.x_title = "Result"
hist.y_title = "Frequency of Result"
hist.add("D6",frequencies)
hist.render_to_file("die_visual.svg")
为了创建条形图,我们创建了一个pygal.Bar()实例,并将其存储在变量hist中,接下来,我们设置hist的属性:
-
hist.title(用于表示直方图名称的的字符串)
-
hist.x_labels(用于表示x轴的标签)在这里我们通过
hist.x_labels = ['1', '2', '3', '4', '5', '6']
,用掷D6骰子的可能的结果做为x轴的标签。 -
hist.x_title 与 hist.y_title(分别表示x轴与y轴的标题)。
然后,我们使用hist.add()
将一系列值添加 到图表中(向它们传递要给添加的值指定的标签,还有一个列表,其中包含将出现在图表中的值)。
最后,我们使用hist.render_to_file()
将这个图表渲染成一个名为 “die_visual.svg” 的SYG文件,而且这种危机的扩展名必须为.svg的格式。die_visual.svg默认存储在die_visual.py所在的文件夹里。
要想查看生成的直方图,最简单的方式就是用Wed浏览器。
同时掷两个面数相同的骰子
同时掷两个骰子,得到的点数更多,结果分布情况也不同。所以我们需要创建两个D6骰子,且将两个骰子的点数相加。并将结果存储在列表results中:
dice_visual.py
import pygal
from die import Die
#创建一个D6
die_1 = Die()
die_2 = Die()
#掷几次骰子,并将结果存储在一个列表中
results = []
for roll_num in range(1000): #这里的roll_num可以换成任意一个为使用过的变量,就是是 _ 也可以。(_表示占位符)
result = die_1.roll() + die_2.roll()
results.append(result)
#分析结果
frequencies = []
max_result = die_1.num_sides + die_2.num_sides
for value in range(1,max_result+1):
frequency = results.count(value)
frequencies.append(frequency)
#对结果进行可视化
hist = pygal.Bar()
hist.title = "Results of rolling one D6 1000 times."
hist.x_labels = ['2', '3', '4', '5', '6', '7', '8', '9', '10', '11', '12']
hist.x_title = "Result"
hist.y_title = "Frequency of Result"
hist.add("D6 + D6",frequencies)
hist.render_to_file("die_visual.svg")
创建两个Die实例后,我们掷骰子多次,并计算每次的总点数,可能出现的最大点数位12为两个骰子的最大点数之和,我们定义了一个变量max_result,并将最大点数之和存储在其中。尽管我们计算2到max_result的各种点数出现的次数(仅仅使用range(2,13)就可以了),但是这个只适用于数据较小时,但在模拟现实世界情况,最后编写轻松地模拟各种情形的代码。因此,不管骰子有多少面,我们都可以计算出。
同时掷两个面数不同的骰子
比如我们创建一个6面和10面的骰子。看看同时运行这两个骰子50000次的结果。
different_dice.py
import pygal
from die import Die
#创建一个D6
die_1 = Die()
die_2 = Die(10)
#掷几次骰子,并将结果存储在一个列表中
results = []
for roll_num in range(50000): #这里的roll_num可以换成任意一个为使用过的变量,就是是 _ 也可以。(_表示占位符)
result = die_1.roll() + die_2.roll()
results.append(result)
#分析结果
frequencies = []
max_result = die_1.num_sides + die_2.num_sides
for value in range(1,max_result+1):
frequency = results.count(value)
frequencies.append(frequency)
#对结果进行可视化
hist = pygal.Bar()
hist.title = "Results of rolling a D6 and a D10 50000 times."
hist.x_labels = ['2', '3', '4', '5', '6', '7', '8', '9', '10', '11', '12', '13', '14', '15', '16']
hist.x_title = "Result"
hist.y_title = "Frequency of Result"
hist.add("D6 + D10",frequencies)
hist.render_to_file("die_visual.svg")