clear; close all; clc; rand('twister',1000); %rand('twister',sum(100*clock)); m1=25; %Mittlere Zahl ankommender Kunden pro Stunde m2=30; %Mittlere Zahl bedienbarer Kunden pro Stunde lambda=m1/60; mu=m2/60; anzahl=1000000; gleich=rand(2,anzahl); ankunftsabstand=-1./lambda.*log(1-gleich(1,:)); bedienzeit=-1./mu.*log(1-gleich(2,:)); ankunft=cumsum(ankunftsabstand); abgang=zeros(1,anzahl); bedienbeginn=zeros(1,anzahl); bedienbeginn(1)=ankunft(1); abgang(1)=bedienbeginn(1)+bedienzeit(1); for k=2:anzahl if ankunft(k)>=abgang(k-1); bedienbeginn(k)=ankunft(k); else bedienbeginn(k)=abgang(k-1); end abgang(k)=bedienbeginn(k)+bedienzeit(k); end % Teilaufgabe a) % Teilaufgabe b) % Teilaufgabe c) lsg_c=mean(bedienbeginn-ankunft) % Teilaufgabe d) lsg_d=mean(abgang-ankunft) % Teilaufgabe e) % Teilaufgabe f) lsg_f=sum(bedienbeginn-ankunft>=15)/anzahl % Teilaufgabe g) lsg_g=sum(abgang-ankunft>=20)/anzahl