数学建模之计算机仿真.pptVIP

图3 :接受服务事件子程序 服务台空 置服务台”忙” 产生服务结束时刻 登记到事件表 排队人数+1 否 是 图4 :寻找服务事件子程序 排队人数≥1 置服务台”忙” 产生服务结束时刻 排队人数-1 排队人数+1 否 是 五：蒙特卡洛方法 在用传统方法难以解决的问题中，有很大一部分可以用概率模型进行描述．由于这类模型含有不确定的随机因素，分析起来通常比确定性的模型困难．有的模型难以作定量分析，得不到解析的结果，或者是虽有解析结果，但计算代价太大以至不能使用．在这种情况下，可以考虑采用Monte Carlo方法。Monte Carlo方法是计算机模拟的基础，它的名字来源于世界著名的赌城——摩纳哥的蒙特卡洛，其历史起源于1777年法国科学家蒲丰提出的一种计算圆周π的方法——随机投针法，即著名的蒲丰投针问题。 18世纪，法国数学家蒲丰和勒可莱尔提出的“投针问题”，记载于蒲丰1777年出版的著作中：“在平面上画有一组间距为a的平行线，将一根长度为L（L=a/2）的针任意掷在这个平面上，求此针与平行线中任一条相交的概率。”蒲丰本人证明了，这个概率是： p=2L/（πa）π为圆周率 利用这个公式可以用概率的方法得到圆周率的近似值。下面是一些资料 试验者 时间 投掷次数 相交次数 圆周率估计值 Wolf 1850年 5000 2532 3.1596 Smith 1855年 3204 1218.5 3.1554 C.De Morgan 1860年 600 382.5 3.137 Fox 1884年 1030 489 3.1595 Lazzerini 1901年 3408 1808 3.1415929 Reina 1925年 2520 859 3.1795 五：蒙特卡洛方法 Monte Carlo方法的基本思想是首先建立一个概率模型，使所求问题的解正好是该模型的参数或其他有关的特征量．然后通过模拟一统计试验，即多次随机抽样试验（确定m和n），统计出某事件发生的百分比．只要试验次数很大，该百分比便近似于事件发生的概率．这实际上就是概率的统计定义．利用建立的概率模型，求出要估计的参数．蒙特卡洛方法属于试验数学的一个分支． 五：蒙特卡洛方法 蒙特卡洛方法适用范围很广泛，它既能求解确定性的问题，也能求解随机性的问题以及科学研究中的理论问题．例如利用蒙特卡洛方法可以近似地计算定积分，即产生数值积分问题． 蒙特卡洛法求圆周率 clear n=50000 X=rand(n,1); Y=rand(n,1); k=0; for i=1:n; if X(i)^2+Y(i)^2=1 k=k+1; end end 4*k/n %1/4*PI*r^2==k/n 蒲丰投针求圆周率的M程序 function pi_value=pinpi(k,a,l) %k投针次数 %a线距 %l针长（必须小于等于a） %pi_value返回pi值 m=0; if la; fprintf(error：针长必须小于等于%d\n,a); fprintf(请重新调用函数pinpi(k,d,l) \n); pi_value=0; else for i=1:k a*rand(1)=l*sin(pi*rand(1)) %判断投针下沿与最近平行线的距离与投针垂直高度的关系，距离小于等于高度就相交 if m=m+1; end end p=m/k; pi_value=2*l/(a*p); fprintf(投针法求得pi=%d\n,pi_value); end format long g pinpi(100000,4,3) 投针法求得pi=3.143797e+000 ans = 3.14379728795087 任意曲边梯形面积的近似计算 一个古老的问题：用一堆石头测量一个水塘的面积．应该怎样做呢？测量方法如下：假定水塘位于一块面积已知的矩形农田之中．如图8．2所示．随机地向这块农田扔石头使得它们都落在农田内．被扔到农田中的石头可能溅上了水，也可能没有溅上水，估计被“溅上水的”石头量占总的石头量的百分比．试想如何利用这估计的百分比去近似计算该水塘面积？ 任意曲边梯形面积的近似计算 任意曲边梯形面积的近似计算 结合图8．2中的图形(1)分析，只要已知各种参数及函数（a，b，H，f(x)），有以下两种方法可近似计算水塘面积． 1．随机投点法 1）赋初值：试验次数n=0，成功次数m=0；规定投点试验的总次数N； 2）随机选择m个数对xi,yi,1im，其中axib, 0yiH，置 n＝n＋l； 3）判断