Discussion :

[Toma-HC]

Bonjour à tous,

Mon problème est le suivant:
J'ai une matrice carrée composée majoritairement de valeurs nulles. J'ai quelques valeurs non nulles par colonnes et je chercher à récupérer la première valeur non nulle de chaque colonne.

En gros j'ai cela:

Code :

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

1
2
3
4
5
 
[0 0 0 0 0; 
0 1 0 2 0;
0 0 0 0 0;
3 4 5 6 7;]

et je veux obtenir cela:

Code :

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

1
2
3
4
5
 
[0 0 0 0 0; 
0 1 0 2 0;
0 0 0 0 0;
3 0 5 0 7;]

voir cela:

Code :

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

1
2
 
[3 1 5 2 7]

J'ai testé cela (repris sur un autre fil de discussion):

Code :

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

1
2
3
4
 
idx=V1~=0;
OK=V1(idx)
toc

Mais cela permet de récupérer toutes les valeurs non-nulles. Peut-on adapter pour récupérer la première valeur non-nulle de chaque colonnes?

Merci d'avance!
Toma

[Invité]

Bonjour,

Une solution :

Code :

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

1
2
3
4
5
6
7
8
9
10
M = [0 0 0 0 0;
0 1 0 2 0;
0 0 0 0 0;
3 4 5 6 7];
 
V = zeros(1,size(M,2));
for i = 1:size(M,2)
   V(i) = M( find(M(:,i)~=0,1,'first') , i); 
end
V

[Toma-HC]

Bonjour,

Merci pour ta réponse.

Pour le moment, j’obtiens une message d'erreur:

Code :

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

Improper assignment with rectangular empty matrix

Je vais essayer de trouver ce qui cloche...

[Invité]

Alors c'est que tu as des colonnes qui ne contiennent que des 0.

Code :

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

1
2
3
4
5
6
7
V = zeros(1,size(M,2));
for i = 1:size(M,2)
    idx = find(M(:,i)~=0,1,'first');
    if ~isempty(idx)
        V(i) = M( idx , i); 
    end
end

[ol9245]

En 3 lignes, et sans boucle for :

Code :

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

1
2
3
4
5
6
7
8
9
10
M = [0 0 0 0 0;
0 1 0 2 0;
0 0 0 0 0;
3 4 5 6 7];
 
[a, b, c] = find(M) ;
[x, y, t] = unique (b, 'first') ;
I = sub2ind(size(M), a(y), b(y)) ;
 
disp(M(I))

[Toma-HC]

Merci!

ça fonctionne parfaitement!