Discussion :

[InfoMatlab]

Alors je suis toujours en train de bosser sur mes problèmes d'info. J'en ai une dizaine, il y en a que je comprends bien mais parfois je bloque notamment sur un qui consiste à programmer la fonction reshape de Matlab.

Déjà j'ai bien compris que la matrice générait par reshape comporte le même nombre d'élément que la matrice source. Je dois donc créer une condition pour déclencher le programme :

Code :

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

if RxC=MxN

Où [R,C] sont des variables d'entrées de ma fonction et [M,N] sont les lignes et les colonnes de ma matrice source que je trouve par la commande :

Code :

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

[M,N]=size(A)

Maintenant là où je bug c'est que je veux créer une nouvelle matrice de taille [R,C]. J'ai beau lire mon livre de syntaxe en long en large et en travers je ne vois pas comment créer cette matrice à part en créant une matrice zeros ou ones pour ensuite faire une opération avec ma matrice d'origine. Mais les opérations ne sont possibles que dans certains cas.

Je ne code rien pour le moment, je n'ai pas Matlab sur mon ordi perso, je réfléchis juste à la structure de la fonction. Mais je me demandais si en faisant ça, ça marcherait :

Code :

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

1
2
3
 
B = ones(R,C)
C = B.*A

Après je ne sais pas si Matlab va prendre tous les éléments de A et les multiplier avec tous les éléments de B alors que leur taille est différente ?

[tachmou]

Salut,

C'est une idée intéressante, mais je ne pense pas que ce soit possible. Personnellement, je ferais la résolution en travaillant sur les éléments.

Cordialement,

[Jean Dumoncel]

Bonjour,

Code :

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

1
2
3
 
B = ones(R,C)
C = B.*A

Après je ne sais pas si Matlab va prendre tous les éléments de A et les multiplier avec tous les éléments de B alors que leur taille est différente ?

Non, cela provoquera une erreur, car dans la documentation de .*:

c = a.*b
...
Inputs a and b must have the same size unless one is a scalar.

Le plus simple est de commencer par des boucles imbriquées pour remplir ta nouvelle matrice.

[InfoMatlab]

Bon je viens de sortir un code. Ce n'est pas le definitif mais c'est juste pour comprendre comment y arriver :

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
30
31
%%On transforme ma matrice A en vecteur C de n elements
A = [ 8 2 3; 4 5 6];
[r,c]=size(A);
%creation de mon vecteur B d'element egal aux elements de A
B = zeros(1,r*c);
%On balaie les lignes de A
for i=1:r
%On balaie les colonnes de A
for j=1:c
%Un element i,j se place dans mon vecteur B 
B(c*(i-1)+j)=A(i,j);
end
end
% On affiche B
C = B
 
%On transforme les elements de A en une matrice D de taille 3x2
%creation de la matrice D
D = zeros(3,2);
n=0;
%On balaie les lignes de D
for i=1:3
%On balaie les colonnes de D
for j=1:2
%A chaque boucle on passe a un element n + 1 de C 
n=n+1;
D(i,j)=C(n);
end
end
%On affiche D
F = D

Je me suis dit qu'il était plus facile de sortir tous les éléments de la matrice origine en un vecteur puis de travailler dessus. C'est peut être pas le plus efficace mais c'est la seule manière que je comprends.

[tachmou]

Bonjour,

Moi je trouve ça excellent.

Cordialement,

[le fab]

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
%%On transforme ma matrice A en vecteur C de n elements
A = [ 8 2 3; 4 5 6];
[r,c]=size(A);
%creation de mon vecteur B d'element egal aux elements de A
B = zeros(1,r*c);
%On balaie les lignes de A
for i=1:r
    %On balaie les colonnes de A
    for j=1:c
        %Un element i,j se place dans mon vecteur B 
        B(c*(i-1)+j)=A(i,j);
    end
end
% On affiche B
C = B

peut s'écrire plus simplement :

Code :

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

1
2
3
A = [ 8 2 3; 4 5 6];
B = A.';
B = B(:).'

[shaiHulud]

Bonjour,

L'assignation linéaire me semble plus simple pour régler le problème. Cette version marche aussi pour des ND arrays

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
function y= my_reshape(x,new_size)
assert(prod(size(x)) == prod(new_size))
y= nan(new_size);
y(:)= x(:);
 
function DEMO
x= rand(4,4)
y= my_reshape(x,[8 2])
assert(all(x(:)==y(:)))
yy= my_reshape(y,[4 4])
isequal(x,yy)

[InfoMatlab]

Ok merci pour les astuces, et j'ai corrigé le code. Il y avait une erreur fondamentale la-dedans meme si ca fonctionnait.

[Jerome Briot]

Pour éviter tout confusion, et bien que ce ne soit pas le cœur de l'exercice demandé, je rappelle que la fonction reshape ne réordonne pas les éléments d'un tableau.

Pour plus d'informations : [Défis][5] RESHAPE plus rapide que son ombre ?

[shaiHulud]

je rappelle que la fonction reshape ne réordonne pas les éléments d'un tableau

Si cette remarque porte sur le code que j'ai posté, je ne la comprends pas

[Jerome Briot]

Si cette remarque porte sur le code que j'ai posté, je ne la comprends pas

Non, pas du tout, la solution que tu proposes est très bonne
On peut même se passer d'un des opérateurs (:) :

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
>> a = [1 3 5 ; 2 4 6]
 
a =
 
     1     3     5
     2     4     6
 
>> b = zeros(3,2);
>> b(:) = a
 
b =
 
     1     4
     2     5
     3     6

Je faisais juste une remarque globale pour dire que la meilleure façon de redimensionner un tableau est d'utiliser la fonction reshape.
Car son temps d'exécution est petit quelque soit la taille du tableau à redimensionner.

L'énoncé de cet exercice (pédagogique au niveau des techniques d'indexage et/ou des structures itératives) ne doit pas conduire à une méconnaissance des mécanismes de MATLAB.

[shaiHulud]

Merci:-) Tout ca m'inspire quelques réflexions :

Dans une optique POO, si on voulait séparer les données de l'interface (ici l'indicage), on pourrait créer une classe qui surcharge subsref et subsasgn pour indicer les données suivant la forme crée, sans reshaper les données elles mêmes.

On aurait ainsi des données reshapées tres lentes en accès, mais dont la création serait très rapide, et ne nécessitant pas de recopier les données.