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.