拯救女神-北太天元学习23

一个女神挽救了很多人的生命，因为她是O型血，万能输血者。
天有不测风云，人有旦夕祸福。
女神生病了，要做手术，需要3管O型血，现在知道人群中的O型血
的比例是 41%,  要抽够3管血时停止，设抽到第X个人恰好抽了3管
O型血， 那么X的期望是多少？X>14的概率是多大？

% 北太天元模拟献血问题
% 中国O型血的人比例是41%， 来自网上搜索
% 我们希望得到3管O型血，假设每个献血者仅仅捐献一管血
% 有一天可能我们采集3个人就得到了3个O型血，有一天
% 可能采集了20个人也没有得到3个O型血
% 恰好得到3个O型血时采集的捐献者的人的个数是一个随机数.
% 我们用数值模拟尝试求这个随机数的期望值(也就是平均值)
clc
clear all
clf
close all

模拟总次数 = 2000;
num_O型 = zeros(模拟总次数, 1); %采集O型血的管数
num_献血者 = zeros(模拟总次数,1); %献血者的人数
for k = 1:模拟总次数 %for 每次模拟
     while num_O型(k) < 3 % 持续循环直到O型血个数等于3个
            if rand < 0.41 %if 一个献血者是O型血
                num_O型(k) = num_O型(k) + 1; %num_O型增加1
            end
            num_献血者(k) = num_献血者(k) + 1; % 献血者增加1
     end
end
 histogram(num_献血者) % 画出达到3个O型血采集的献血者的个数
 xlabel('献血者的个数')
 ylabel('频率')

 mean_献血者 = mean(num_献血者) %平均需要的献血者
 标准差 = std(num_献血者) %计算标准差
 mean_超过10献血者的比例 = mean(num_献血者 > 10) % 超过10献血者的比例
mean_超过11献血者的比例 = mean(num_献血者 > 11) % 超过11献血者的比例
mean_超过12献血者的比例 = mean(num_献血者 > 12) % 超过12献血者的比例
mean_超过13献血者的比例 = mean(num_献血者 > 13) % 超过13献血者的比例
mean_超过14献血者的比例 = mean(num_献血者 > 14) % 超过14献血者的比例