Matlab

% Possible solutions implemented in Matlab for some of the challenge problems

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% MATLAB CODE

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% PROBLEM 1.0

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

c = [12.1, 6.45, 73, 24.6, 15.2, 44.3, 19.0]

m = [58.22, 57.25, 53.23, 43.91, 42.69, 46.23, 43.08, 46.06, 44.6, 47.55, 53.91, 63.05, 43.62, 48.26, 62.83, 41.3, 39.05, 38.66, 55.3, 45.19, 68.7]

imean = 2.5

istd =0.5

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Method of matching moments, assuming independence

C = lognormal(mean(c), std(c))

M = normal(mean(m), std(m))

I = normal(imean, istd)

H = C .* I ./ M

1 - prob(H,5) % probability that 5 < H

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Confidence boxes, assuming independence

C = CBlognormal(c)

M = CBnormal(m)

I = normal(imean, istd)

H = C .* I ./ M

1 - prob(H,5) % probability that 5 < H

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% PROBLEM 1.1

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

C = env(CBlognormal(c), 0)

M = CBnormal(m)

I = normal(imean, istd)

H = C .* I ./ M

1 - prob(H,5) % probability that 5 < H

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% PROBLEM 1.2

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

C = env(CBlognormal(c), 0)

for i = 1 to many(m) do m.i = censor(m.i, [0,20])

M = CBnormal(m)

I = normal(imean, istd)

H = C .* I ./ M

1 - prob(H,5) % probability that 5 < H

Page updated

Google Sites

Report abuse