Discussion :

[kwatz]

Bonjour,

Je souhaiterais extraire d'une matrice les lignes non-nulles.
J'ai vu que la fonction sparse correspond, mais elle ne me renvoie pas une matrice, mais des adresses et leur valeurs associées.

Voici mon code :

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
 
%Calcul de N
clear
clc
Tfusion=input ('Entrer la Température de fusion? \n')
 
L=load ('C:\Documents and Settings\Sartorius\Mes documents\Sartorius_Stedim\sortie_matlab\Tmilieu.txt','-ASCII'); %charge le tableau des températures
 
m=round(length(L(1,:))/2);
 
for k=1:length(L(:,1)) %met à zeros les lignes où T milieu<Tfusion
    if L(k,m) < Tfusion
        L(k,:)=0;
    end
end
L = sparse(L)
pause
temps=L(:,1);
tf=length(temps);
T=L(:,2:length(L(1,:)));
 
A=7.3e-6;
Ea=6e4;
 
[Derive]=calcul_expo(A,Ea,T);
 
Simpson(Derive,temps(1),temps(tf));

En fait, je charge un tableau qui comprend sur la 1ère colonne des temps, et sur les colonnes suivantes l'évolution de la température à différentes coordonnées.
Je souhaite avoir un tableau qui ne comporte que les cellules dont la température est supérieure à une valeur fixe ( Tfusion dans mon cas). C'est pourquoi, j'ai choisi de mettre à 0 toutes les cellules ne satisfaisant pas cette condition, et que je veux enlever les lignes avec les 0.

[duf42]

Bonjour,

En utilisant FIND tu devrais y arriver:

Code :

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

1
2
3
4
5
6
7
8
9
 
% Trouve les indices des lignes où il y a des 0 (seulement pour les colonnes de 2 à la fin)
[line_idx, col_idx] = find(L(:, [2:end])<Tfusion);
 
% On s'assure de ne garder qu'une fois chaque ligne
line_idx = unique(line_idx);
 
% On supprime les lignes ainsi définies
L(line_idx, :) = [];

Voilà j'espère que ca pourra t'aider

Bonne journée,
Duf

P.S. Il doit y avoir une solution plus courte mais j'ai pas le courage de la chercher

[nahouto]

Essaie un truc comme ça :

Code :

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

L(L(:,2)==0,:) =[]

[kwatz]

Merci, c'est top!

[nahouto]

je me suis trompé, je n'avais pas compris que L contenait plus de 2 colonnes.

Regarde plutot ça :

Code :

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

L(~all(L(:,2:end),2),:) = [] ;

[kwatz]

En fait, j'ai l'impression que la première solution fonctionne bien, car elle elimine les ligne où la deuxième valeur est nulle, non?

[nahouto]

En fait, j'ai l'impression que la première solution fonctionne bien, car elle elimine les ligne où la deuxième valeur est nulle, non?

oui, la première solution que j'ai proposée efface la ligne si et seulement si la deuxième colonne (pour cette ligne) est nulle.

La deuxième solution que j'ai proposée efface la ligne si et seulement si une des colonnes (sauf la première) (pour cette ligne) est nulle.

[Jerome Briot]

La deuxième solution que j'ai proposée efface la ligne si et seulement si une des colonnes (sauf la première) (pour cette ligne) est nulle.

Dans ce cas, tu aurais dû utiliser ANY, non ?

[kwatz]

oui, la première solution que j'ai proposée efface la ligne si et seulement si la deuxième colonne (pour cette ligne) est nulle.

Oui, donc ca fonctionne bien car je met à zero toute la ligne quand la condition n'est pas remplie.

[nahouto]

Dans ce cas, tu aurais dû utiliser ANY, non ?

c'est équivalent dans ce cas, puisqu'il met toute la ligne à zéro (je n'avais pas saisi ce point là)

Oui, donc ca fonctionne bien car je met à zero toute la ligne quand la condition n'est pas remplie.

ok, je comprends mieux.

[kwatz]

D'ailleurs, avec cette astuce

Code :

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

L(L(:,2)==0,:) =[]

Je peux même enlever ma boucle for, mon code devient du coup :

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
 
%Calcul de N
clear
clc
Tfusion=input ('Température de fusion? \n')
L=load ('C:\Documents and Settings\Sartorius\Mes documents\Sartorius_Stedim\sortie_matlab\Tmilieu.txt','-ASCII');
 
m=round(length(L(1,:))/2);
 
L(L(:,2)<Tfusion,:) =[];
 
temps=L(:,1);
tf=length(temps);
T=L(:,2:length(L(1,:)));
 
A=0.0000073;
Ea=20000;
 
[Derive]=calcul_expo(A,Ea,T);
 
S=Simpson(Derive,temps(1),temps(tf));

Merci pour votre aide!

[Jerome Briot]

D'ailleurs, avec cette astuce

Et cette astuce s'appelle l'indexage logique :

• Qu'est-ce que l'indexage logique ?
• III-C. Indexage logique

[Jerome Briot]

Tu peux même remplacer le bloc suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
temps=L(:,1);
tf=length(temps);
T=L(:,2:length(L(1,:)));
 
A=0.0000073;
Ea=20000;
 
[Derive]=calcul_expo(A,Ea,T);
 
S=Simpson(Derive,temps(1),temps(tf));

par :

Code :

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

1
2
3
4
5
6
A=0.0000073;
Ea=20000;
 
[Derive]=calcul_expo(A,Ea,L(:,2:end));
 
S=Simpson(Derive,L(1,1),L(end,1));

A vérifier

[kwatz]

Tu peux même remplacer le bloc suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
temps=L(:,1);
tf=length(temps);
T=L(:,2:length(L(1,:)));
 
A=0.0000073;
Ea=20000;
 
[Derive]=calcul_expo(A,Ea,T);
 
S=Simpson(Derive,temps(1),temps(tf));

par :

Code :

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

1
2
3
4
5
6
A=0.0000073;
Ea=20000;
 
[Derive]=calcul_expo(A,Ea,L(:,2:end));
 
S=Simpson(Derive,L(1,1),L(end,1));

A vérifier

Merci, je viens de tester, ca a l'air de fonctionner !
A ce rythme, je vais finir par poster tout mon programme pour que tu l'optimise !

[Jerome Briot]

A ce rythme, je vais finir par poster tout mon programme pour que tu l'optimise !

Je dirais plutôt que c'est à toi de bien comprendre chacune de nos "optimisations" et d'essayer de les appliquer toi-même

Dans ce sens, pense à conserver les anciennes versions de tes programmes "non-optimisés" pour pouvoir comparer leur évolution

[kwatz]

Héhé, si je comprends bien, je vais devoir faire le boulot

Merci beaucoup en tous cas!