MATLAB kommandoer til eksempel C i avsnitt 3.3

 

% Vi vil lage et plott av den simultane sannsynlighetstettheten i eksempel C i avsnitt 3.3 i Rice for ulike verdier av a (alfa). [Vi bruker U og V for � betegne variablene, ikke X og Y som i Rice.]

 

% Vi� setter f�rst opp �et grid� av verdier for u og v:

u=0:0.01:1;

v=0:0.01:1;

[u,v]=meshgrid(u,v);

 

% Vi beregner s� den simultane sannsynlighetstettheten c(u,v)� [dvs. copula-tettheten]

a=-1;

c=(1+a)+a*(4*u.*v-2*u-2*v);

 

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

mesh(u,v,c)

xlabel('u-akse')

ylabel('v-akse')

zlabel('c(u,v)')

 

 

% Plott selv copula-tettheten for a=1 og a=0 og legg merke til hva verdien til a betyr for �fasongen� til tettheten.