20138 随机模拟32.pptxVIP

模拟的概念;模拟的方法; 在实际问题中，面对一些带随机因素的复杂系统，用分析方法建模常常需要作许多简化假设，与面临的实际问题可能相差甚远，以致解答根本无法应用．这时，计算机模拟几乎成为唯一的选择．;例1　在我方某前沿防守地域，敌人以一个炮排（含两门火炮）为单位对我方进行干扰和破坏．为躲避我方打击，敌方对其阵地进行了伪装并经常变换射击地点．; 需要模拟出以下两件事： ;2. 符号假设;4. 模拟结果;5. 理论计算;6. 结果比较;clear k1=0;k2=0;k3=0; for i=1:200 i R1=rand if R1=0.5 R2=rand if R2=3/6 k1=k1+1; else if R25/6 k3=k3+1; else k2=k2+1; end end else k1=k1+1; end end E=(k2+k3)/200 E1=(0*k1+1*k2+2*k3)/200;[1] 系统的假设： （1） 顾客源是无穷的； （2） 排队的长度没有限制； （ 3） 到达系统的顾客按先后顺序依次进入服务， 即“先到先服务”．; [2] 符号说明 w：总等待时间；ci：第i个顾客的到达时刻； 　　 bi：第i个顾客开始服务时刻； ei：第i个顾客服务结束时刻． xi：第i-1个顾客与第i个顾客到达之间的时间间隔 yi：对第i个顾客的服务时间;[3] 模拟框图;clear i=2; w=0; x(i)=exprnd(10); c(i)=x(i); b(i)=x(i); while b(i)=480 y(i)=unifrnd(4,15); e(i)=b(i)+y(i); w=w+b(i)-c(i); i=i+1; x(i)=exprnd(10); c(i)=c(i-1)+x(i); b(i)=max(c(i),e(i-1)); end i=i-2; t=w/i m=i ;clearcs=100;for j=1:cs j w(j)=0; i=2;x(i)=exprnd(10);c(i)=x(i);b(i)=x(i);while b(i)=480 y(i)=unifrnd(4,15);  e(i)=b(i)+y(i);  w(j)=w(j)+b(i)-c(i);   i=i+1;  x(i)=exprnd(10);  c(i)=c(i-1)+x(i);  b(i)=max(c(i),e(i-1)); endi=i-2;t(j)=w(j)/i;m(j)=i;end pt=0;pm=0;for j=1:cs pt=pt+t(j); pm=pm+m(j);endpt=pt/cspm=pm/cs