数学建模专题三MonteCarlo模拟教材课程.pptVIP

数学建模专题三;内容提纲;引言(Introduction);举例;Solution;例1. 蒲丰投针问题; 一些人进行了实验，其结果列于下表 ：;基本思想;例1　在我方某前沿防守地域，敌人以一个炮排（含两门火炮）为单位对我方进行干扰和破坏．为躲避我方打击，敌方对其阵地进行了伪装并经常变换射击地点．; 需要模拟出以下两件事： ;i：要模拟的打击次数； k1：没击中敌人火炮的射击总数； k2：击中敌人一门火炮的射击总数； k3：击中敌人两门火炮的射击总数； E：有效射击比率； E1：20次射击平均每次毁伤敌人的火炮数．;模拟框图;模拟结果;Date;理论计算;结果比较;注：rand(n)=rand(n,n); name 的取值可以是;fix(x) : 截尾取整，直接将小数部分舍去 floor(x) : 不超过 x 的最大整数 ceil(x) : 不小于 x 的最小整数 round(x) : 四舍五入取整; 随机投掷均匀硬币，验证国徽朝上与朝下的概率是否都是 1/2 ; 随机投掷骰子，验证各点出现的概率是否为 1/6 ; 用蒙特卡罗 ( Monte Carlo ) 投点法计算 ? 的值 ; 在画有许多间距为 d 的等距平行线的白纸上，随机投掷一根长为 l ( l ? d ) 的均匀直针，求针与平行线相交的概率，并计算的 ? 值;n=100000; l=0.5; d=1; m=0; for i=1:n alpha=rand(1)*pi; y=rand(1)*d/2; if y=l/2*sin(alpha) m=m+1; end end fprintf(针与平行线相交的频率为：%f\n,m/n); fprintf(计算出来的 pi 为：%f\n’,2*n*l/(m*d));; 设某班有 m 个学生，则该班至少有两人同一天生日的概率是多少？;n=1000; p=0; m=50; % 设该班的人数为 50 for t=1:n a=[]; q=0; for k=1:m b=randperm(365); a=[a,b(1)]; end c=unique(a); if length(a)~=length(c) p=p+1; end end fprintf(‘至少两人同一天生日的概率为：%f\n,p/n);;clear; m = 50; p1= 1:365; p2= [1:365-m, 365*ones(1,m)]; p = p1./p2; p = 1- prod(p); fprintf(至少两人同一天生日的概率为：%f\n,p);;排队问题随机模拟;[1] 系统的假设： （1） 顾客源是无穷的； （2） 排队的长度没有限制； （ 3）到达系统的顾客按先后顺序依次进入服务， “先到先服务”。;w：总等待时间； ci：第i个顾客的到达时刻； bi：第i个顾客开始服务时刻； ei：第i个顾客服务结束时刻； xi：第i-1个顾客与第i个顾客之间到达的间隔时间； yi：对第i个顾客的服务时间。;c1;;用蒙特卡洛法解非线性规划问题;基本假设;框 图; 在Matlab软件包中编程，共需三个Ｍ－文件:randlp.m, mylp.m, lpconst.m.主程序为randlp.m.;% randlp.m function [sol,r1,r2]=randlp(a,b,n)　　 %随机模拟解非线性规划 debug=1; a=0;　　　　　　　　　　　　　 %试验点下界 b=10;　　　　　　　　　　　　 　%试验点上界 n=1000;　　　　　　　　　　　　 %试验点个数 r1=unifrnd(a,b,n,1);　　　　　　%n?1阶的［a,b］均匀分布随机数矩阵 r2=unifrnd(a,b,n,1); sol=[r1(1) r2(1)]; z0=inf; for i=1:n x1=r1(i); x2=r2(i); lpc=lpconst([x1 x2])； if lpc==1 z=mylp([x1 x2])； if zz0 z0=z； sol=[x1 x2]； end end end;Thanks!