Discussion :

[LeTube]

yop yop


dans un programme, je cherche à prendre les valeurs non nulles de 2 matrices définies par bloc (X et Z) et à les stocker dans une autre matrice (de taille inférieure) où une colonne correspond à un bloc

mon problème est que mon programme ne garde que les valeurs négatives de la matrice Z et je ne comprends pas pourquoi

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
 
for k = 1:Nc                     %Nc= nombre de blocs dans la grosse matrice
    p = 0;                      
 
    for i = 1:Nl                   %Nl*Nr = taille des blocs des grosses matrices
 
        for j = 1:Nr
 
            if (X(i+(k-1)*Nl,j+(k-1)*Nr)~=0)
                p = p+1;                                  %p sert de repère dans la colonne k de xx et zz
                xx(p,k) = X(i,j);
                zz(p,k) = Z(i,j);                         
            end
 
        end
    end
end

à noter, seules les points de coordonnée x non nulle m'intéressent, je prends le z aussi mais c'est parce que j'ai besoin de leur ordonnée dans le plan xOz

Si quelqu'un a une idée je suis preneur.
Sinon, pas de bol
Cordialement
Tube

[Jerome Briot]

Pourrais-tu nous montrer ce que contient X ou nous fournir ce tableau dans un fichier .mat (mis dans une archive zip ou rar) ?

[LeTube]

Voilà pour les données

[Jerome Briot]

Le code suivant semble donner le résultat attendu, non ?

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
clc
clear all

load X
load Z

Nc = 10;
Nl = 10;
Nr = 5;

for k = 1:Nc                     %Nc= nombre de blocs dans la grosse matrice
    p = 0;                      

    for i = 1:Nl                   %Nl*Nr = taille des blocs des grosses matrices

        for j = 1:Nr

            if (X(i+(k-1)*Nl,j+(k-1)*Nr)~=0)
                p = p+1;                                  %p sert de repère dans la colonne k de xx et zz
                xx(p,k) = X(i+(k-1)*Nl,j+(k-1)*Nr);
                zz(p,k) = Z(i+(k-1)*Nl,j+(k-1)*Nr);                         
            end

        end
    end
end

[LeTube]

Un grand merci Dut (encore une fois et c'est surement pas la dernière )

en effet c'était bien un problème d'indice

Merci Merci et encore Merci

Tube, imbécile mais heureux