用
MATLAB
模拟掷硬币过程
我们掷一枚硬币,
它出现的情况只可能是要么出现正面,
要么出现反面,
这一随机事件的样
本点有限且等可能。
所以掷硬币这一随机事件为古典概型,
它出现的样本点是有限的且等可能。
为了模拟掷
硬币出现正面或者反面,规定随机数小于
0.5
时为反面,否则为正面。
在
MATLAB
中提供了一个在
[0
,
1]
区间上均匀分布的随机函数
rand()
。可用
round()
函
数将其变成
0
—
1
阵,然后将整个矩阵的各元素值加起来再除以总的元素个数即为出现正面
的概率。
一、连续掷
100
次硬币,运行
1000
次,程序如下:
for i=1:100
a(i)=sum(sum(round(rand(1000))))/100
end
以下值是截取部分运行结果,发现正面出现的概率总是在
0.5
左右浮动。
a =
Columns 1 through 6
0.5080 0.4935 0.5059 0.5035 0.5040 0.4969
Columns 7 through 12
0.4945 0.5010 0.4875 0.4948 0.5042 0.5022
Columns 13 through 18
0.4894 0.4965 0.4977 0.4969 0.5008 0.5154
Columns 19 through 24
0.5013 0.5022 0.4969 0.5006 0.4980 0.4949
Columns 25 through 30
0.4878 0.4994 0.4964 0.4944 0.4995 0.4950
Columns 31 through 36
0.5066 0.4989 0.4940 0.4988 0.4939 0.4909
Columns 37 through 42
0.4977 0.5025 0.4927 0.4977 0.4904 0.5040
Columns 43 through 48
0.5080 0.4935 0.5059 0.5035 0.5040 0.4969
Columns 49 through 54
0.5029 0.4977 0.5024 0.4956 0.4857 0.5035
Columns 55 through 60
0.5006 0.5116 0.5035 0.4953 0.4974 0.5012
Columns 61 through 66
0.4997 0.5039 0.5009 0.5012 0.5037 0.5021
Columns 67 through 72
1