随机模拟与系统仿真课件.pptVIP

随机模拟 一、随机模拟介绍 随机模拟是一种数量技术，它利用计算机 一、随机模拟介绍 （1）复杂的实际问题难于用解析理论处理， 一、随机模拟介绍 一、随机模拟介绍 （1）问题本身就是随机的； 一、随机模拟介绍 一、随机模拟介绍 二、随机数的产生 用随机模拟方法模拟某过程的时候，需要 产生各种概率分布的随机变量。 二、随机数的产生 设连续型随机变量 具有概率密度 则称 在区间[a，b]上服从均匀分布。 二、随机数的产生 二、随机数的产生 例1: 随机变量X表示随机到达的人数,分布列如下： xk 0 1 2 pk 0.4 0.3 0.3 模拟十分钟内到达的状况。 二、随机数的产生 用MATLAB模拟随机事件 rand(20, 1) %生成20个均匀随机数向量。 r=rand(1,10); for i=1:10; if r(i)0.4 n(i)=0; elseif r(i)0.7 n(i)=1; else n(i)=2; end; end 二、随机数的产生 r r =0.5678 0.7942 0.0592 0.6029 0.0503 0.4565 0.0185 0.8214 0.4447 0.6154 n n = 1 2 0 1 0 1 0 2 1 1 r r =0.2311 0.6068 0.4860 0.8913 0.7621 0.7919 0.9218 0.7382 0.1763 0.4057 n n = 0 1 1 2 2 2 2 2 0 1 二、随机数的产生 例2：利用[0，1]区间均匀分布的随机数r1，r2……表示服从负指数分布的随机数。 二、随机数的产生 二、随机数的产生 三、模拟法分类 三、模拟法分类 三、模拟法分类 L=zeros(1,31);% L 等待的顾客人数， T1=zeros(1,31);%T1总合等待时间， T2=zeros(1,31); %T2总合服务时间, L1=zeros(1,31); %L1总合到达顾客人数。 t=1;tau=1;x=0:30; r=rand(1,30);%随机数 %随机模拟 for i=1:30; t=t+1; if 0=r(i) r(i)0.4 n=0; elseif 0.4=r(i) r(i)0.7 n=1; else n=2; end; %排队分析 if L(t-1)==0 n==0 L(t)=L(t-1);T1(t)=T1(t-1);%模型 T2(t)=T2(t-1);L1(t)=L1(t-1); else L(t)=L(t-1)+n-1;T1(t)=T1(t-1)+L(t); T2(t)=T2(t-1)+tau; L1(t)=L1(t-1)+n; end; end; r=[0 r]; a=[x,r,L,L1,T1,T2] eL=T2(end)/tau %已被服务的人数 L2=(find(L1eL)) L3=sum(L(L2))%未被服务的顾客等待时间总和 g1=(T1(end)-L3)/eL %平均等待时间 g2=g1+tau %平均逗留时间 g3=eL/30 %平均每分钟服务的顾客人数 t RND n L T1 T2 1 .81 2 1 1 1 2 .02 0 0 0 1 3 .60 1 0 0 1 4 .04 0 0 0 0 5 .46 1 0 0 1 6 .31 0 0 0 0 7 .67 1 0 0 1 8 .25 0 0 0 0 9 .24 0 0 0 0 10 .10 0 0 0 0 11 .40 1 0 0 1 12 .02