MATLAB kommandoer til eksempel F i avsnitt 3.3

 

% Vi vil lage et plott av tettheten for den bivariate normalfordelingen i eksempel F i avsnitt 3.3 i Rice.

% Vi setter f?rst opp ?et grid? av x- og y-verdier:

?

x=-3:0.05:3;

y=-3:0.05:3;

[x,y]=meshgrid(x,y);

 

% Vi beregner s? den binormale tettheten for gitte verdier av parametrene:

 

mx=0; my=0; sx=1; sy=1;

u=(x-mx)/sx; v=(y-my)/sy;

r=0.9;

konst=1/(2*pi*sx*sy*sqrt(1-r^2));

f=konst*exp(-(u.^2+v.^2-2*r*u.*v)/(2*(1-r^2)));

 

% Vi plotter s? den simultane sannsynlighetstettheten og setter navn p? aksene

mesh(x,y,f)

xlabel('x-akse')

ylabel('y-akse')

zlabel('f(x,y)')

 

% Vi lager plott med ?h?ydekurver?:

contour(x,y,f)

xlabel('x-akse')

ylabel('y-akse')