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

图1到来事件子程序系统人数+1产生下一个顾客到来时刻调用接受服务事件的程序第37页,共61页，星期六，2024年，5月图2离去事件子程序系统人数-1置服务台”闲”已服务人数+1调用寻找服务事件子程序第38页,共61页，星期六，2024年，5月图3:接受服务事件子程序服务台空置服务台”忙”产生服务结束时刻登记到事件表排队人数+1否是第39页,共61页，星期六，2024年，5月图4:寻找服务事件子程序排队人数≥1置服务台”忙”产生服务结束时刻排队人数-1排队人数+1否是第40页,共61页，星期六，2024年，5月五：蒙特卡洛方法在用传统方法难以解决的问题中，有很大一部分可以用概率模型进行描述．由于这类模型含有不确定的随机因素，分析起来通常比确定性的模型困难．有的模型难以作定量分析，得不到解析的结果，或者是虽有解析结果，但计算代价太大以至不能使用．在这种情况下，可以考虑采用MonteCarlo方法。MonteCarlo方法是计算机模拟的基础，它的名字来源于世界著名的赌城——摩纳哥的蒙特卡洛，其历史起源于1777年法国科学家蒲丰提出的一种计算圆周π的方法——随机投针法，即著名的蒲丰投针问题。第41页,共61页，星期六，2024年，5月五：蒙特卡洛方法MonteCarlo方法的基本思想是首先建立一个概率模型，使所求问题的解正好是该模型的参数或其他有关的特征量．然后通过模拟一统计试验，即多次随机抽样试验（确定m和n），统计出某事件发生的百分比．只要试验次数很大，该百分比便近似于事件发生的概率．这实际上就是概率的统计定义．利用建立的概率模型，求出要估计的参数．蒙特卡洛方法属于试验数学的一个分支．第42页,共61页，星期六，2024年，5月五：蒙特卡洛方法蒙特卡洛方法适用范围很广泛，它既能求解确定性的问题，也能求解随机性的问题以及科学研究中的理论问题．例如利用蒙特卡洛方法可以近似地计算定积分，即产生数值积分问题．第43页,共61页，星期六，2024年，5月蒙特卡洛法求圆周率clearn=50000X=rand(n,1);Y=rand(n,1);k=0;fori=1:n;ifX(i)^2+Y(i)^2=1k=k+1;endend4*k/n第44页,共61页，星期六，2024年，5月普丰投针求圆周率的M程序functionpi_value=pinpi(k,a,l)%k投针次数%a线距%l针长（必须小于等于d）%pi_value返回pi值m=0;ifla;fprintf(error：针长必须小于等于%d\n,a);fprintf(请重新调用函数pinpi(k,d,l)\n);pi_value=0;elsefori=1:kifa*rand(1)=l*sin(pi*rand(1))m=m+1;endendp=m/k;pi_value=2*l/(a*p);fprintf(投针法求得pi=%d\n,pi_value);end第45页,共61页，星期六，2024年，5月formatlonggpinpi(100000,4,3)投针法求得pi=3.143797e+000ans=3.14379728795087第46页,共61页，星期六，2024年，5月任意曲边梯形面积的近似计算

一个古老的问题：用一堆石头测量一个水塘的面积．应该怎样做呢？测量方法如下：假定水塘位于一块面积已知的矩形农田之中．如图8．2所示．随机地向这块农田扔石头使得它们都落在农田内．被扔到农田中的石头可能溅上了水，也可能没有溅上水，估计被“溅上水的”石头量占总的石头量的百分比．试想如何利用这估计的百分比去近似计算该水塘面积？第47页,共61页，星期六，2024年，5月任意曲边梯形面积的近似计算

第48页,共61页，星期六，2024年，5月任意曲边梯形面积的近似计算

结合图8．2中的图形(1)分析，只要已知各种参数及函数（a，b，H，f(x)），有以下两种方法可近似计算水塘面积．1．随机投点法1）赋初值：试验次数n=0，成功次数m=0；规定投点试验的总次数N；2）随机选择m个数对xi,yi,1im，其中yaxib,0yiH，置n＝n＋l；3）判断n≤N，若是，转4，否则停止计算；