用 mathemtical 完成下例子.pptVIP

用mathemtical完成下例子;例5. 赶火车过程仿真（随机问题的计算机仿真）; 这个人到达B站时刻的频率分布为; 我们用计算机仿真的方法来解决这个问题，就是通过多次计算机实验来分析此人能赶上火车的概率有多大，这需要我们建立一个仿真模型利用计算机仿真这个实际过程。 引入以下变量： T1 火车从A站开出的时刻； T2 火车从A站运行到B站所需要的时间； T3 此人到达B站的时刻； 此人能赶上火车的条件是 T3 T1+T2 T1, T2, T3 是随机变量，其概率分布为;T2服从N（30，4）; T1 =0, if 0 x =0.7, T1 =5, if 0.7x =0.9, T1 =10, if 0.9x =1, ;; 正方形的四个顶点各有一人.在某一时刻,四人同时出发以匀速v按顺时针方向追赶下一个人.如果他们始终保持对准目标,试求出每个人的行进路线.;解: 建立平面直角坐标系xoy ; ;For[k = 1, k = 12, k++, (t = k; dt = 0.02; v = 1; n = t/dt; robit = {{{0, 10}}, {{10, 10}}, {{10, 0}}, {{0, 0}}}; For[j = 1, j = n, j++, For[i = 1, i = 4, i++, xx1 = robit[[i, j, 1]];yy1 = robit[[i, j, 2]]; If[i != 4, xx2 = robit[[i + 1, j, 1]]; yy2 = robit[[i + 1, j, 2]], xx2 = robit[[1, j, 1]];yy2 = robit[[1, j, 2]]]; dd = Sqrt[(xx2 - xx1)^2 + (yy2 - yy1)^2] // N; xx1 = xx1 + v*dt*(xx2 - xx1)/dd; yy1 = yy1 + v*dt*(yy2 - yy1)/dd; robit[[i]] = Append[robit[[i]], {xx1, yy1}]]]; g1 = ListPlot[robit[[1]], PlotJoined - True, DisplayFunction - Identity]; g2 = ListPlot[robit[[2]], PlotJoined - True, DisplayFunction - Identity]; g3 = ListPlot[robit[[3]], PlotJoined - True, DisplayFunction - Identity]; g4 = ListPlot[robit[[4]], PlotJoined - True, DisplayFunction - Identity]; g5 = ListPlot[{{0, 10}, {10, 10}, {10, 0}, {0, 0}}, PlotJoined - True, DisplayFunction - Identity]; Show[{g1, g2, g3, g4, g5}, DisplayFunction - $DisplayFunction, AspectRatio - Automatic])]