赶上火车的概率.ppt

赶上火车的概率 【问题 】如图，一列火车从A站途经B站开往C站，某人每天从B站上这趟火车赶往C站. 赶上火车的问题 问题分析 问题分析 基本假设 建立模型 建立模拟模型 模拟算法设计 模拟结果 课内练习 请求解模型P{T3＜T1＋T2 }的解析解。 小结 数学建模 * 数学建模 Yong Zhang A B 火车运行方向 某人 他已了解到： 1）火车从A站到B站的运行时间是均值为30分钟，标准差为2分钟的随机变量； 思考：请研究他能否赶上这趟火车。 C 0.7 0.2 0.1 频 率 午后1:00 午后1:05 午后1:10 出发时刻 他到达B 站的时刻的频率分布为 0.3 0.4 0.2 0.1 频率 午后1:28 午后1:30 午后1:32 午后1:34 时刻 他能否及时赶上火车？ 2）火车在下午大约1点离开A站，离开时刻 的频率分布如下： 明确问题：他能及时赶上火车的概率是多少？ 分析法：用概率统计知识建立分析模型求解析解。（思考） 模拟法：用概率统计知识建立模型，通过模拟求解近似解。 即先建立模拟模型，然后通过计算机模拟得到问题的近似解。在同样条件下多次试验，计算他能及时赶上火车的频率。 建模方向（思路）： 能及时赶上火车的充要条件是： 其中 T1—火车从A站出发的时刻； T2—火车的运行时间； T3—他到达B站的时刻。 是什么变量？如何模拟？ 2.将午后1时记为t=0，设火车运行时间T2 服从 正态分布： T2 ~N（30，22）。 1.假设T1，T2，T3都是相互独立的随机变量; 火车出发时刻T1和人到达B站时刻T3的分布律分别为： 0.7 0.2 0.1 0 5 10 P(t) T1（分） 0.3 0.4 0.2 0.1 28 30 32 34 P(t) T3（分） 为了简化计算，将下午1点记为初始时刻。得到随机变量T1和T3的分布律如下： 能及时赶上火车的概率p=P{T3＜T1＋T2 } 车A站出发 人到达B站 如果r为在(0,1)均匀分布的随机数，为了模拟随机变量T1和T3，可以通过如下方法： 则t1和t3可分别用来模拟随机变量T1和T3。 主要变量说明： n 模拟次数 k 临时变量，存储当前累计模拟次数 count 存储赶上火车的次数 两种不同风格的算法描述 输入：赶火车次数（天数） 输出：赶上火车的频率 第 1 步 输入模拟次数n 第 2 步 k=1，count=0 第 3 步 当k=n，执行Step 4－9，否则执行第11步 第 4 步 生成均匀分布随机数赋给r 第 5 步 由r及公式确定T1模拟火车出发时刻 第 6 步 生成均匀分布随机数赋给r； 第 7 步 由r及公式确定T3模拟人达到时刻 第 8 步 生成正态分布随机数T2模拟火车运行时间 第 9 步 IF T1+T2 T3, count=count+1，END 第 10 步 执行第3步 第 11 步 输出赶上火车频率p=count/n 算法1（分步骤描述） 1.初始化： 输入模拟次数n; count=0; 2.模拟n次 for i=1 to n, 模拟随机变量T1,T2,T3， 分别赋给t1,t2,t3; if t1+t2 t3, count=count+1 end if end for app_prb=count / n; 算法2（伪代码描述） n=input(输入模拟次数：); count=0; for i=1:n, rt1=rand; %模拟随机变量t1（火车从A站出发的时刻） if rt10.7 T1=0; elseif rt1=0.7 rt10.9 T1=5; else T1=10; end T2=30+randn*2; %模拟随机变量t2（火车的运行时间） %模拟随机变量t3（他到达B站的时刻） rt3=rand; if rt30.3 T3=28; elseif rt3=0.3 rt30.7 T3