Discussion :

[genius4evers]

Salut.

je suis entraine de calculer l'information mutuelle entre 2 variables ,
apres un recherche j'ai trouve des functions dans le site du Matlab
en faisant :

Code :

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

z = mutualInformation(e',(original_message));

Error :

??? Error using ==> sparse
Sparse matrix sizes must be non-negative integers less than MAXSIZE as
defined by COMPUTER.  Use HELP COMPUTER for more details.

Error in ==> mutualInformation at 16
    Mx = sparse(idx,x,1,n,k,n);

qu'est ce qu'il faut faire ?

Merci d'avance

[ol9245]

k ou n doit être nul ou négatif.

menu debug / stop if error/warning

cocher : always stop ...

puis relance ton programme. Quand il se plantera, tu pourra inspecter la valeur des variables et savoir ainsi celle qui est pas dans les clous.

[genius4evers]

Salut.
l'erreur est dans l'instructin suivante :

Code :

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

Mx = sparse(idx,x,1,n,k,n);

voici les resultats dans l'attache , je ne comprends pas ou se fait l'erreur, il n'ya pas des nombres negatives,

[Invité]

Salut,

??? Error using ==> sparse
Sparse matrix sizes must be non-negative integers less than MAXSIZE as
defined by COMPUTER. Use HELP COMPUTER for more details.

[genius4evers]

Salut

1.comment je peux connaitre le size du sparse matrice ?
2.je ne sais pas comment utiliser la variable computer,
error dans l'attache.

[Invité]

Tu as aussi le droit de copier ces lignes et de les mettre dans des balises [code]
Sinon computer en minuscules :

Code :

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

[C MAXSIZE] = computer

[genius4evers]

Bonsoir
je ne comprends pas pourquoi il ya un erreur
j'ai vu le document du sparse :

Code :

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

1
2
3
S = sparse(i,j,s,m,n,nzmax) uses vectors i, j, and s to generate an m-by-n sparse matrix such that S(i(k),j(k)) = s(k), with space allocated for nzmax nonzeros. Vectors i, j, and s are all the same length. Any elements of s that are zero are ignored, along with the corresponding values of i and j. Any elements of s that have duplicate values of i and j are added together.
 
Note   If any value in i or j is larger than the maximum integer size, 2^31-1, then the sparse matrix cannot be constructed.

pour moi

Code :

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

1
2
3
4
5
 max(idx)
 
ans =
 
     1048576

j'ai pas depaase 2^32-1,
je pense qu'il ya un erreur dans l'alogorithme,

[Jean Dumoncel]

Bonjour,

c'est toujours pour l'instruction :

Code :

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

Mx = sparse(idx,x,1,n,k,n);

avec les valeurs de ta capture d'écran? Si oui, alors les agruments 3 et 4 sont les dimensions de ta matrice (et donc doivent être des entiers) or tu donnes k=3.9557 ...

[genius4evers]

Bonjour.
oui c'est vrai .

Mais si je prend

Code :

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

k=round(k) --> k=4

l'erreur est le meme.

moi je suis entraine d'utilise une fonctionne predefini , j'ai la trouve sur le site du Mathworks.

le code est le suivant :

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
 
function z = mutualInformation(x, y)
% Compute mutual information I(x,y) of two discrete variables x and y.
% Written by Mo Chen (mochen80@gmail.com).
 
    assert(numel(x) == numel(y));
    n = numel(x);
    x = reshape(x,1,n);
    y = reshape(y,1,n);
 
    l = min(min(x),min(y));
    x = x-l+1;
    y = y-l+1;
    k = max(max(x),max(y));
 
    idx = 1:n;
    Mx = sparse(idx,x,1,n,k,n);
    My = sparse(idx,y,1,n,k,n);
    Pxy = nonzeros(Mx'*My/n); %joint distribution of x and y
    Hxy = -dot(Pxy,log2(Pxy+eps));
 
    Px = mean(Mx,1);
    Py = mean(My,1);
 
    % entropy of Py and Px
    Hx = -dot(Px,log2(Px+eps));
    Hy = -dot(Py,log2(Py+eps));
 
    % mutual information
    z = Hx + Hy - Hxy;
endfunction