Discussion :

[AC88mm]

Bonsoir à tous,

hé bien voici mon problème, je pars des densités de proba de deux distributions Gamma ci-dessous, et je voudrais créer la densité de proba jointe, mais comme je construis le vecteur densité en deux morceaux (xin négatif, xin positif), je ne sais pas comment construire ma matrice de densité jointe à l'aide du meshgrid. Je suis parvenu à le faire pour des distributions gaussiennes jointes, mais ici je n'y arrive pas à cause de ce "xin négatif, xin positif", et "yin négatif, yin positif". Pourriez-vous m'aider s'il vous plaît ?

Merci d'avance et bonne soirée

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
N = 40;
nN = (N-1)/2;
lamb1 = 1;
lamb2 = 3;
alph1 = 2;
alph2 = 2;
 
figure
 
grid
hold on
xin=-nN:nN/1000:nN;
h(1:1000)=0;
for i=1001:length(xin)
h(i)=(lamb1^alph1/gamma(alph1))*xin(i)^(alph1-1)*exp(-lamb1*alph1);
end
plot(xin,h,'LineWidth',2)
xlabel('x1')
ylabel('fX1(x1)')
title('densité de probabilité de la distribution Gamma #1')
 
 
figure
 
grid
hold on
yin=-nN:nN/1000:nN;
h(1:1000)=0;
for i=1001:length(yin)
h(i)=(lamb2^alph2/gamma(alph2))*yin(i)^(alph2-1)*exp(-lamb2*alph2);
end
plot(yin,h,'LineWidth',2)
xlabel('x2')
ylabel('fX2(x2)')
title('densité de probabilité de la distribution Gamma #2')
 
%figure
% 
% grid
% hold on
% [xin,yin]=meshgrid(-nN:nN/50:nN); % et là je ne sais plus quoi faire

[AC88mm]

J'ai résolu mon problème en passant par une bête boucle pour faire des produits terme à terme
Il doit y avoir moyen d'utiliser reshape aussi, je pense; si jamais vous avez des idées...

Voici le code et bonne soirée ^^

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
N = 15;
nN = (N-1)/2;
lamb1 = 1;
lamb2 = 3;
alph1 = 2;
alph2 = 2;
 
figure
 
grid
hold on
xin=-nN:nN/100:nN;
hx(1:100)=0;
for i=101:length(xin)
hx(i)=(lamb1^alph1/gamma(alph1))*xin(i)^(alph1-1)*exp(-lamb1*xin(i));
end
plot(xin,hx,'LineWidth',2)
xlabel('x1')
ylabel('fX1(x1)')
title('densité de probabilité de la distribution Gamma #1')
 
 
figure
 
grid
hold on
yin=-nN:nN/100:nN;
hy(1:100)=0;
for i=101:length(yin)
hy(i)=(lamb2^alph2/gamma(alph2))*yin(i)^(alph2-1)*exp(-lamb2*yin(i));
end
plot(yin,hy,'LineWidth',2)
xlabel('x2')
ylabel('fX2(x2)')
title('densité de probabilité de la distribution Gamma #2')
 
figure 
 
grid
hold on
for i=1:length(xin)
    hxy(i,:)=hx(i)*hy(:);
end
subplot(2,1,1)
[xin,yin]=meshgrid(-nN:nN/100:nN);
surfl(xin,yin,hxy)
xlabel('x1')
ylabel('x2')
zlabel('fX1X2(x1,x2)')
subplot(2,1,2)
contourf(xin,yin,hxy)
xlabel('x1')
ylabel('x2')
zlabel('fX1X2(x1,x2)')
title('densité de probabilité jointe')

[Jerome Briot]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
for i=1:length(xin)
    hxy(i,:)=hx(i)*hy(:);
end

Pourquoi pas simplement :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

hxy = hx.'*hy;

[AC88mm]

En effet pourquoi pas! Merci beaucoup Dut ^^