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.