Discussion :

[NELLLY]

Voilà un exp

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
format long e
p1=0.07;
TCmin=30000;
for n1=1:5
    for n2=4:7
        for L1=3:5
            for L2=4:5
                n1
                n2
                L1
                L2
                q11=binocdf(L1,n1,p1)
                q12=binocdf(L2,n1,p1)-binocdf(L1,n1,p1)
                q21=binocdf(L1,n2,p1)
                q22=binocdf(L2,n2,p1)-binocdf(L1,n2,p1)
                 Q=[q11 q12;q21 q22];
                 QQ=eye(2)-Q
                 if det(QQ)==0
                     TC=NaN;
                 else
                     A=[0 1]*inv(QQ)*ones(2,1)
                     B=A-2
                     C=(3*A)+B
                     TC=A+C+B
                   end
                 end
             end
         end
     end

si det(QQ)==0 on ne fait pas le calcul de TC.
Lors de l'execution
une erreur s'affiche

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
 
n1 =
 
     5
 
 
n2 =
 
     5
 
 
L1 =
 
     3
 
 
L2 =
 
     5
 
 
q11 =
 
    9.998866727999994e-001
 
 
q12 =
 
    1.133272000005681e-004
 
 
q21 =
 
    9.998866727999994e-001
 
 
q22 =
 
    1.133272000005681e-004
 
 
QQ =
 
    1.133272000005681e-004   -1.133272000005681e-004
   -9.998866727999994e-001    9.998866727999994e-001
 
Warning: Matrix is close to singular or badly scaled.
         Results may be inaccurate. RCOND = 6.776264e-021.
> In C:\MATLAB6p5\work\vssmc\exp55.m at line 22

le determinant dans ce cas =0 et pourtant il fait le calcul

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
 
A =
 
   -7.379533931401575e+019
 
 
B =
 
   -7.379533931401575e+019
 
 
C =
 
   -2.951813572560630e+020
 
 
TC =
 
   -4.427720358840945e+020

où est le problème?

[christophe_halgand]

Pour moi, il a un problème dans le calcul de "det(QQ)"

Peux tu nous donner le résultat de matlab stp

(au cas où tu ne serais pas, tu écris "det(QQ)" juste avant le if et sans le ";")

Merci

[Merel]

Bonjour Nelly,

Pour moi, il n'y a aucun problème avec ton programme. Le déterminant n'est pas nul. Il est juste très petit et matlab affiche 0...

Au premier warning,

Code :

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

1
2
 
det(QQ)=2.140692026037594e-021

Donc ton programme va bien calculer tes valeures A, B, C et TC.

Maintenant il faut voir si ce warning va te gêner ou non...

[mr_samurai]

utilise plutot :

Code :

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

abs(det(QQ))<eps

ou, le pseudo-inverse :

Code :

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

pinv(QQ);

[paradize3]

A la place de

Code :

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

1
2
 
abs(det(QQ))<eps

je conseillerais plutot la fonction cond(...)

Salutations,

Gregoire