Discussion :

[momo032114]

bonjour
J'ai une matrice 327723 lignes par 3 colonnes
j'aimerai en fait cree une autre matrice carree ou je met la troisieme colonne de la premiere matrice selon la place de la premiere colonne et de la deuxieme colonne de la premiere matrice
a mon avis un exemple s'impose pour comprendre :
premiere matrice

Code :

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

1
2
3
4
5
6
7
8
9
445100.000000 6686925.000000 -14.292036 
445112.500000 6686925.000000 -101.512749 
445125.000000 6686925.000000 -267.692780 
445137.500000 6686925.000000 -89.545677 
445150.000000 6686925.000000 -93.365463 
445012.500000 6686937.500000 -899.299683 
445025.000000 6686937.500000 -1176.459839 
445037.500000 6686937.500000 -1206.381470 
444900.000000 6686962.500000 -23.203894

deuxieme matrice

Code :

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

1
2
3
-14.292036  -101.512749  -267.692780   -89.545677   93.365463 
      0         -899.299683  -1176.459839 -1206.381470     0
      0              0            -23.203894          0                0

tout depend de mon Y si il ne varie pas je ne change pas de ligne sinon je change #
svp

[Caro-Line]

Tu aurais pu essayé de simplifier un peu l'exemple, là ça fait un peu mal aux yeux

[EDIT]
En fait j'ai pris une mauvaise hypothèse ci-dessous par rapport au placement des zéros
Donc la question préliminaire est en fait : quel est le critère pour le placement des zéros dans ta matrice résultat ?
[/EDIT]

Il faut donc :
Dans un 1er temps trouver les indices de la 2ème colonne où on a les mêmes données : pour cela il faut utiliser la fonction UNIQUE.

Connaitre aussi les éléments uniques de ta 1ère colonne afin de pouvoir placer tes zéros, toujours avec la fonction UNIQUE.

Question : est-ce que ta matrice est déjà triée ?
Dans ton exemple il semblerait qu'elle soit triée sur la 2ème colonne, si ce n'est pas le cas il faudrait le faire au préalable avec la fonction SORT.

Ensuite une méthode serait de boucler sur les éléments uniques de la 2ème colonne, de regarder pour chacun quels éléments de la 1ère colonne existent et remplir la ligne de ton résultat en fonction.

Bon la fin de mon post n'est peut-être pas très claire mais essaye déjà de commencer et quand tu bloques, poste ton code, on verra pour la suite

[momo032114]

desole pour l'exemple je voulais tente d'expliquer le probleme

donc pour repondre a tes questions ma matrice est deja ordonnee dans le sens croissant par rapport a Y

ensuite pour les zeros en fait cela depend du plus grand nombre de recurrence de Y ; ensuite les zeros vont se placer symetriquement par rapport aux chiffres de la troisieme colonne ( enfin je voudrais )

voila le debut de 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
fid=fopen('02-04 RMS2.txt','rt') % ouverture du fichier 
%%%%%%%%%%%%%%%%% enleve l entete %%%%%%%%%%%%
for n=1:20 % enleve les 20 premieres lignes
    fgetl(fid);
end 
%%%%%%%%%%%%%%%%%% lit le fichier %%%%%%%%%%%%
A=fscanf(fid,'%lf',[5 327723])';
%%%%%%%%%%%%%%%% supprime colonne 4 et 5 %%%%%
A(:,5)=[];
A(:,4)=[];
 
X=A(:,1);
Y=A(:,2);
Z=A(:,3);
 
n=length(Y)
 
b=unique (Y);

ensuite moi je ferai

Code :

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

1
2
3
4
for i=1:n
    z=Z(i)
    y=Y(i)
           if y=b

ben la je coince il faudrait cree la matrice et mettre les elements de Z dedans si Y est le meme
merci de m'aider caro

[Caro-Line]

Ce n'est pas la peine de passer à chaque fois par des variables temporaires

Comme tu ne m'as pas dit comment on choisissait où on met les zéros, je te propose une solution, à toi après de l'adapter :

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
M=[445100.000000 6686925.000000 -14.292036
445112.500000 6686925.000000 -101.512749
445125.000000 6686925.000000 -267.692780
445137.500000 6686925.000000 -89.545677
445150.000000 6686925.000000 -93.365463
445012.500000 6686937.500000 -899.299683
445025.000000 6686937.500000 -1176.459839
445037.500000 6686937.500000 -1206.381470
444900.000000 6686962.500000 -23.203894 ];
res2=unique(M(:,2));%éléments uniques de la 2ème colonne
%Solution sans zéros supplémentaires
%Comme toutes les lignes n'auront pas le même nombre d'éléments
%on passe par un cell array
resultat=cell(length(res2),1);
for kR=1:length(res2)
    %on recherche les indices lignes tels que M(idx,2)==res2(kR)
    idx = find(M(:,2)==res2(kR));
    resultat{kR}=M(idx,3);
end

[momo032114]

merci caro

comment je fais maintenant pour recuperer ma matrice car j'aimerai faire dessus une fft2 mais la le probleme c'est cell ?
pour les zeros je dois regarder la ligne ou il ya le plus de chiffres dans la matrice et ensuite placer les zeros de part et d'autre dans les autres lignes

[Caro-Line]

Voilà un code qui fait ce que tu souhaites :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
%Remplir de zeros (et remettre en ligne)
%1. Taille maxi des vecteurs dans resultat :
tailleMax=max(cellfun(@length,resultat));
for kR=1:length(resultat)
    tailleResK=length(resultat{kR});
    nbZeros=tailleMax-tailleResK;%nb de zeros a rajouter
    gauche=round(nbZeros/2);%au cas où le nombre de zéros est impait
    droite=nbZeros-gauche;
    resultat{kR}=[zeros(1,gauche) resultat{kR}' zeros(1,droite)];
end
resFinal=cell2mat(resultat);%transfo cellule => matrice

Mais est-ce que ça a vraiment une signification physique de faire ça ?
Ça me parait bizarre...
Bon par contre il serait bon que tu essayes par toi-même (en cherchant dans l'aide, en testant) de comprendre comment je suis arrivée là.
Ce n'est que comme ça que tu pourras progresser.

Si tu as des questions spécifiques sur certaines lignes du code, n'hésites pas à demander.

D'autre part il y a surement un tas de façons différentes de faire et surement des meilleures, j'ai juste fait au plus vite.

[momo032114]

merci beaucoup caro
t'es super
si tu veut plus de details je suis en train de travailller sur des donnees sismiques pour tracer des graphs 4D
voila merci encore je vais travailler le code que tu m'as donne

[soft001]

D'autre part il y a surement un tas de façons différentes de faire et surement des meilleures, j'ai juste fait au plus vite.

j'espère bien ça, surtout si on optimise au maximum le temps d'exécution.

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
 
r=[445100.000000 6686925.000000 -14.292036
445112.500000 6686925.000000 -101.512749
445125.000000 6686925.000000 -267.692780
445137.500000 6686925.000000 -89.545677
445150.000000 6686925.000000 -93.365463
445012.500000 6686937.500000 -899.299683
445025.000000 6686937.500000 -1176.459839
445037.500000 6686937.500000 -1206.381470
444900.000000 6686962.500000 -23.203894 ];
index=find(diff(r(:,2))~=0);
s=[diff([0; index]);size(r,1)-index(end)];
l=1;
m=1;
res=zeros(length(index)+1,max(s));
for i=1:length(index)+1
    if s(l)~=max(s)
        for j= fix(((max(s)-s(l))/2))+mod(max(s)-s(l),2)+1 :s(l)+fix(((max(s)-s(l))/2))+mod(max(s)-s(l),2)
            res(i,j)=r(m,3);
            m=m+1;
        end
    else
        for j= 1:max(s)
            res(i,j)=r(m,3);
            m=m+1;
        end  
    end
    l=l+1;
end

[Jerome Briot]

j'espère bien ça, surtout si on optimise au maximum le temps d'exécution.

Si r devient très grand, il est préférable de faire ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
[u,a,b] = unique(r(:,2));
 
nu = numel(u);
 
v = histc(b,1:nu);
mv = max(v);
 
res = zeros(nu,mv);
 
for n = 1:nu
   idx = (1:v(n))+(mv-v(n))/2;
   res(n,idx) = r(b==n,3);
end

[soft001]

Salut

Si r devient très grand, il est préférable de faire ceci :

Je demande pardon, "l'élève commente le travail de son professeur", j'arrive pas à croître que je peux le faire.
En tous cas, je vois si je me trompe pas que vous n'avez pas tenu compte de l'influence de la parité de (mv-v(n)).
Votre code après rectification :

Code :

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

1
2
3
4
5
6
7
8
9
10
[u,a,b] = unique(r(:,2));
nu = numel(u);
v = histc(b,1:nu);
mv = max(v);
res1 = zeros(nu,mv);
for n = 1:nu
   idx = (1:v(n))+ceil((mv-v(n))/2);
   res1(n,idx) = r(b==n,3);
end

De plus, je doute, même si r est très grand, que l'utilisation de votre code sera favorisée. Pourtant, pour s'assurer mieux de mes paroles j'ai modifié mon code comme suit :

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
 
 
index=find(diff(r(:,2))~=0);
s=[diff([0; index]);size(r,1)-index(end)];
m=1;
ms=max(s);
ls=length(s);
res=zeros(ls,ms);
for i=1:ls
    if s(i)~=ms
        for j= fix(((ms-s(i))/2))+mod(ms-s(i),2)+1 :s(i)+fix(((ms-s(i))/2))+mod(ms-s(i),2)
            res(i,j)=r(m,3);
            m=m+1;
        end
    else
        for j= 1:ms
            res(i,j)=r(m,3);
            m=m+1;
        end  
    end
end

Deuxième méthode "Inspirer de votre solution"

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
index=find(diff(r(:,2))~=0);tt=[0 index'];
s=[diff([0; index]);size(r,1)-index(end)];
m=1;
ms=max(s);
ls=length(s);
res3=zeros(ls,ms);
for i=1:ls
            m=tt(i)+1:s(i)+tt(i);
            idx = (1:s(i))+ceil((ms-s(i))/2);
            res3(i,idx)=r(m,3);
end

[Jerome Briot]

En tous cas, je vois si je me trompe pas que vous n'avez pas tenu compte de l'influence de la parité de (mv-v(n)).
Votre code après rectification :

Effectivement, j'ai fais l'hypothèse que le nombre de zéro de part et d'autre était toujours pair.

De plus, je doute, même si r est très grand, que l'utilisation de votre code sera favorisée.

Par rapport à la première solution, si.

Et pourtant si... Pourtant, pour s'assurer mieux de mes paroles j'ai modifié mon code comme suit :

Effectivement la seconde version est bien mieux optimisée.

Je demande pardon, "l'élève commente le travail de son professeur", j'arrive pas à croître que je peux le faire.

Ben heureusement que si...