Moyenne sur les lignes d'une matrice en négligeant certaines valeurs

Version imprimable

13/08/2008, 12h37
momo032114

Moyenne sur les lignes d'une matrice en négligeant certaines valeurs

bonjour
voila j'ai une matrice 1024*1024 je voudrait prendre le mean de cette matrice mais sans prendre les Zeros qui sont a l'interieur
et ensuite faire la soustraction que sur les nombres differents de zeros

ex :

0 0 12 15 18 0 0 0 ... me donne 15
0 11 6 2 3 8 0 0 .. me donne 7
0 0 0 0 3 6 0 0 me donne 4.5

et ensuite

0 0 -3 0 3 0 0 0 ..
0 4 -1 -5 -4 1 0 0 ..
0 0 0 0 -1.5 1.5 0 0...
13/08/2008, 14h09
Jerome Briot
Tu fais une boucle FOR-END qui scanne chaque ligne de ta matrice, tu utilises l'indexage logique pour trouver les valeurs non nulles sur chaque ligne et la fonction MEAN pour faire la moyenne

Sinon :
Code:

1 2 3 4 5 X = [0 0 12 15 18 0 0 0 0 11 6 2 3 8 0 0 0 0 0 0 3 6 0 0] bsxfun(@minus,X,sum(X,2)./sum(X~=0,2)).*(X~=0)
13/08/2008, 15h51
momo032114
j'ai tente la premiere methode mais ensuite j'aimerai enlever ce que le mean a chaque chiffre pour les lignes
Code:

1 2 3 4 5 6 for i=1:948 x=resfinal5(:,i); if x~=0 B=mean(x) end end
ensuite je voudrai que mean(x) de chaque ligne soit soustrait a chaque x de la meme ligne et recuper ensuite ma matrice avec mes zeros et x-mean(x)
13/08/2008, 16h24
Jerome Briot

Citation:

Envoyé par Dut

tu utilises l'indexage logique pour trouver les valeurs non nulles sur chaque ligne

=> :faq: : Qu'est-ce que l'indexage logique ?
13/08/2008, 16h26
Caro-Line
Avant de faire quoi que ce soit d'autre, revoit ton code : il est faux...
Pour plein de raisons :
- Tu boucles sur les colonnes et pas sur les lignes
- x~=0 sera vrai ssi tous les éléments de x sont différent de 0
  - soit B est la moyenne de toute la colonne qui n'a pas de zéros
  - soit il y a un zéro et la moyenne n'est pas calculée
Et fait un tour par la :faq:
Qu'est-ce que l'indexage logique ?
13/08/2008, 17h09
momo032114
alors j'ai regarder l'aide mais le probleme est que je n'y arrive pas
a mon avis il faut faire
Code:

1 2 3 4 5 for i=1:422 for j=1:948 x=resfinal(i,j); idx=find(x>0.1);
mais je vois pas comment on fait pour appliquer mean(car il est ecrit que maintenant la matrice idx cest comme du binaire ) et le soustraire pour ensuite recuperer la matrice du debut avec les x differents de zeros changes

Allez, puisqu'apparemment tu n'arrives pas à appréhender la faq (qui précise que l'indexage logique permet justement de ne pas utiliser FIND) :

Code:

1
2
3
4
5
6
7
8
9
10
11
X = [0 0 12 15 18 0 0 0
0 11 6 2 3 8 0 0
0 0 0 0 3 6 0 0]
Xmean=zeros(size(X,1),1);%vecteur colonne qui a le même nombre de lignes que X
%et qui va contenir les moyennes de chaque ligne
%Boucle sur les lignes :
for i=1:size(X,1)
   maLigne=X(i,:);%ième ligne
   Xmean(i)=mean(maLigne(maLigne~=0));%moyenne de maLigne pour maLigne~=0
end
Xmean %affichage

et si je fais cela ?

Code:

1
2
3
4
5
6
7
8
9
10
 
for i=1:422
    for j=1:948
        x=resfinal(i,j);
        if x~=0
           B=mean(x);
           C(i,1)=B ;
        end  
    end
end

et ensuite

Code:

1
2
3
4
5
6
7
8
9
10
11
for i=1:422
    for j=1:948
        x=resfinal(i,j);
 
        if x~=0
            y=x-C(i,1);
            A=zeros(422,948);
            B=A+y;
        end
    end
end

13/08/2008, 17h45
Caro-Line

1. Je suppose que tu n'as pas vu mon post précédent.

2. Essaye de comprendre pas à pas ce que tu fais.
Dans ton premier code : pour i et j données x est juste une valeur.
Si x est différent de 0 (par exemple 12 pour i=1 et j=3 dans l'exemple) tu as donc B qui est la moyenne d'une valeur (donc 12) => ça ne te choque pas ?
13/08/2008, 17h51
momo032114

si en fait j'ai reflechi ensuite a ce que j'avais poste mais c'etait trop tard
13/08/2008, 18h04
momo032114
merci
ce qui m'interresser au debut c'etait en fait de recuperer la matrice X mais avec les element (differents de zeros) changes de x-mean(x)

bon vais retenter un code mais ca foire car j'ai plus de zeros
Code:

1 2 3 4 5 6 7 8 9 10 for i=1:422 for j=1:948 x=resfinal(i,j); y=resfinalmean(i,1); if x~=0 M(i,j)=x-y; end end end
13/08/2008, 18h18
Caro-Line

Citation:

Envoyé par momo032114

bon vais retenter un code mais ca foire car j'ai plus de zeros

Pourquoi tu n'as plus de zéros ? Où ça ?
C'est quoi resfinal (ça n'est plus la matrice de départ), c'est quoi resfinalmean ?

Pour avoir une réponse il faut :
- donner le maximum d'informations
- essayer d'expliquer le plus clairement possible

As-tu essayer avant de coder de "pseudo-coder" ce que tu veux faire, c'est-à-dire d'écrire "à la main" l'algorithme que tu cas utiliser pour faire ce que tu souhaites ?
Peut-être que ce serait plus compréhensible :?
13/08/2008, 18h39
momo032114

en fait c bon j'ai trouve tout seul
merci de ton aide
13/08/2008, 19h16
Jerome Briot

Citation:

Envoyé par momo032114

en fait c bon j'ai trouve tout seul

On pourrait en savoir un peu plus ?

ben ouais
voici mon code
resfinal est ma matrice et resfinal mean mon mean de chaque ligne
Code:

1 2 3 4 5 6 7 8 9 for i=1:422 for j=1:948 x=resfinal(i,j); y=resfinalmean(i,1); if x~=0 resfinal(i,j)=x-y; end end end

Et pour éviter la double boucle, en 1 seule fois (sauf si garder la moyenne t'intéresse), et en utilisant l'indexage logique :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
X = [0 0 12 15 18 0 0 0
0 11 6 2 3 8 0 0
0 0 0 0 3 6 0 0]
%objectif : obtenir Y tel que Y(i,j)=X(i,j)-mean(X(i,:)) sauf pour X(i,j)=0
%Initialisation
Y=zeros(size(X));
%Boucle sur les lignes :
for i=1:size(X,1)
   maLigne=X(i,:);%ième ligne
   %soustraction pour les x ~=0
   Y(i,maLigne~=0)=maLigne(maLigne~=0)-mean(maLigne(maLigne~=0));
end
Y

Voici une autre solution sans recours aux structures itératives imbriquées.:D
Code:

1 2 3 4 5 6 7 8 r=[0 0 12 15 18 0 0 0;0 11 6 2 3 8 0 0;0 0 0 0 3 6 0 0]; [m n]=size(r); moy=sum(r,2)./sum(r~=0,2) u=(r~=0); b=repmat(moy,1,n); res=r-b.*u

Fuseau horaire GMT +2. Il est actuellement 12h28.