TMA4240 Statistikk, høsten 2016
MTBYGG, MTING, MTMT, MTTEKGEO
Matlabscript: Estimerer forventingsverdi ved gjennomsnitt
% % Simulerer trekning av k kuler fra urne med totalt n kuler, % nR er røde og nG er gule. Trekker kuler gjentatte ganger og % estimerer forventningsverdier til X (antall røde kuler trukket ut) % og Y (antall gule kuler trukket ut) ved å regne ut gjennomsnitt. % Plotter resultatet som funksjon av antall forsøk utført. % n = 12; nR = 3; nG = 4; k = 3; m = 2000; meanx = []; meany = []; sumx = 0; sumy = 0; for i=1:m draw = randperm(n,k); xx = sum(draw <= nR); yy = sum(draw > nR & draw <= nR + nG); sumx = sumx + xx; sumy = sumy + yy; meanx = [meanx sumx/i]; meany = [meany sumy/i]; figure(1); plot(meanx); figure(2); plot(meany); if (i == 1) pause end pause(0.05); end [meanx(length(meanx)) meany(length(meany))]