MATLAB kommandoer som illustrerer de teoretiske beregningene i eksempel C og D p? sidene 61-62 i Rice?

 

 

Eksempel C side 61

 

% Generere f?rst n=1000 standardnormalfordelte observasjoner og

% tegner histogram av dem:

 

n=1000;

Z=normrnd(0,1,1,n);

hist(Z)

 

% Kvadrerer s? de normalfordelte observasjonene og tegner histogram av

% de transformerte observasjonene:

 

X=Z.^2;

hist(X)

 

% Merk at histogrammet av de transformerte observasjonene ser ut som en

% gammafordeling med formparameter alfa < 1 (jf. figur 2.11 p? side 53 i Rice).

% Det er si samsvar det vi finner teoretisk i eksempel C.

 

 

 

Eksempel D side 62

 

% Generere f?rst n=1000 standardnormalfordelte observasjoner og

% tegner histogram av dem:

 

n=1000;

U=unifrnd(0,1,1,n);

hist(U)

 

% Regner s? ut de inverse av observasjonene og tegner histogram av

% de transformerte observasjonene (vi m? fikse opp aksene ved

% axis-kommandoen for ? f? et rimelig utseende av histogrammet):

 

V=1./U;

hist(V,0.5:1:21)

axis([0,20,0,600])

 

% Merk at histogrammet av de transformerte observasjonene f?r et utseende som

% er i samsvar med den tettheten vi finner i eksempel D.