Discussion :

[linda05]

Bonjour;
J'ai une boucle dans mon programme qui me donne une matrice par exemple :

Code :

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

B=[3 2 1;2 1 3];

et je veux remplir une matrice P :

Code :

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

P=zeros(4*size(B,1),4);

par la matrice B et 4 tel que P soit :

Code :

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

P=[4 3 2 1;3 4 2 1;3 2 4 1;3 2 1 4;4 2 1 3;2 4 1 3;2 1 4 3;2 1 3 4];

je remplis P pour trouver une matrice C:

Code :

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

C= [sum( A( sub2ind(size(A),rows,cols) ) ,2) P]

tel que :

Code :

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

1
2
rows = P(:,1:n-1);
cols = P(:,2:n);

A une matrice d'ordre 4.
et puis je ne garde que:

Code :

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

P=[C((C(:,1)==max(C(:,1))),2:end)];

comment jeux programmer le remplissage de P ?
Merci beaucoup
À bientôt

[Invité]

Bonjour,

Une solution :

Code :

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

1
2
3
4
5
B = [3     2     1
     2     1     3];
P = 4*[eye(4)  eye(4)];
P(~P) = kron(B, ones(4,1)).';
P = P.';

[Jerome Briot]

Une alternative, moins performante si B possède un grand nombre de lignes :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
Pl = zeros(4,4);
Pu = zeros(4,4);
P = zeros(4*size(B,1),4);
 
for n = 1:size(B,1)
 
    BB = repmat(B(n,:),3,1);    
 
    Pl(2:4,1:3) = tril(BB);
    Pu(1:3,2:4) = triu(BB);
 
    P((n-1)*4+1:n*4,:) = Pl+Pu+diag([4 4 4 4]);
 
end

[linda05]

Merci bien;
juste une remarque pour le diag ds le cas generale si

Code :

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

1
2
size(B,2)==n
alors

Code :

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

diag(4*ones(1,n+1))

j'ai une autre question pour la matrice C,si j'ai par exemple:

Code :

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

C(:,1)=[1 3 15 4 9 0];

ce que je voudrai, c'est comment arriver a retrouver le max de c(:,1) avec la valeur qui est juste inférieur à le max(donc 15 et 9 ds l'exemple)
merci

[linda05]

Bonjour a tous;
j'ai une autre question:si on a maintenent un vecteur:

Code :

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

V=[3 2 1];

et P une matrice:

Code :

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

P=zeros(4*size(V,2),4)

,ma question c'est comment remplir P tel que P soit:

Code :

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

P=[4 3 2 1;3 4 2 1;3 2 4 1;3 2 1 4;4 2 1 3;2 4 1 3;2 1 4 3;2 1 3 4;4 1 3 2;1 4 3 2;1 3 4 2;1 3 2 4];

Merci

[Invité]

Dans le même esprit que précédemment :

Code :

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

1
2
3
4
5
6
7
V = [3 2 1];
vCirc = cell2mat(arrayfun(@(k) circshift(V,[1 k]), ...
    (0:-1:-numel(V)+1).',...
    'UniformOutput', false));
P = repmat(4*eye(4), 1, 3);
P(~P) = kron(vCirc, ones(4,1)).';
P = P.'

[linda05]

Merci beaucoup

[linda05]

Bonjour a tous ;
je rencontre des difficultés :
par exemple pour la matrice :

Code :

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

1
2
3
4
B=[
     3     2     4     1
     2     4     1     3
     3     1     4     2];

Je cherche une matrice P tel que P soit :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
P = [
     5     1     3     2     4
     1     5     3     2     4
     1     3     5     2     4
     1     3     2     5     4
     1     3     2     4     5
     5     4     2     3     1
     4     5     2     3     1
     4     2     5     3     1
     4     2     3     5     1
     4     2     3     1     5 ];

Donc mon code est :

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
B=[3 2 4 1;2 4 1 3;3 1 4 2];
BCirc=zeros(size(B,2)*size(B,1),size(B,2))
for i=1:size(B,1)
    BCirc((i-1)*size(B,2)+1:i*size(B,2),:)=cell2mat(arrayfun(@(k) circshift(B(i,:),[1 k]), ...
    (0:-1:-numel(B(i,:))+1).',...
    'UniformOutput', false));
end
BCirc
BCirc =unique(BCirc ,'rows')
n=5;
m = size(BCirc,1);
P = zeros(n*m,n);
Pl = zeros(n,n);
Pu = zeros(n,n);
 
for i = 1:size(BCirc,1)
    BB = repmat(BCirc(i,:),n-1,1)  ;  
 
    Pl(2:n,1:n-1) = tril(BB);
    Pu(1:n-1,2:n) = triu(BB);
 
    P((i-1)*n+1:i*n,:) = Pl+Pu+diag(n*ones(1,n));
 
end

mes questions:
1. dans la matrice B on remarque que la ligne 3 est en sens inverse de la ligne 2, comment je le détecter pour le suprimer et B soit :

Code :

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

B=[3 2 4 1;2 4 1 3] ;?

2. il ’ya des répétitions ds la matrice BCirc ?et pour ca j’ai ajouté :
BCirc =unique(BCirc ,'rows') ;
J'aimerai trouver la matrice P d'une façon propre et optimisée
SVP j'ai besoin d'un coup de main pour optimizer le code
Merci