Discussion :

[Choupi]

Bonjour à tous,

Je suis nouveau dans le monde de Matlab. J'utilise davantage R pour le calcul matriciel et je souhaite faire simplement :

matrice - vecteur ... sous R c'est possible, le vecteur est soustrait à chaque ligne ou à chaque colonne selon.

Comment faire ça sous matlab ? je ne vois pas l'interet de dupliquer mon vecteur pour avoir les mêmes dimenstions...

Merci bien !

Choupi

[Jerome Briot]

Quatre possibilités :

Code :

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

1
2
3
4
5
6
7
8
9
M=[1 2
    3 4
    5 6]
 
X=[4 6]
 
for n=1:size(M,1)
   M(n,:)=M(n,:)-X;
end

Ou

Code :

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

1
2
3
4
5
6
7
M=[1 2
    3 4
    5 6]
 
X=[4 6]
 
M=M-repmat(X,size(M,1),1);

Ou

Code :

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

1
2
3
4
5
6
7
M=[1 2
    3 4
    5 6]
 
X=[4 6]
 
M=reshape(cell2mat(arrayfun(@(n){M(n,:)-X},1:size(M,1))),size(M,2),[]).'

Ou

Code :

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

1
2
3
4
5
6
7
M=[1 2
     3 4
     5 6]
 
 X=[4 6]
 
M = bsxfun(@minus, M, X)

Pour le choix, tout dépend de la taille de la matrice à traiter et si on veut conserver un code lisible ou pas

[Choupi]

Merci c'est exactement ça.

Je vais faire des tests pour voir ce qui est le plus rapide. Pour le moment j'utilise une boucle.

Encore merci,

Choupi

[Jerome Briot]

Il est également possible d'utiliser un fichier MEX, voici l'idée :

Le code C, test.c:

Code C :

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
#include "mex.h"
 
void mexFunction( int nlhs, mxArray *plhs[],
                  int nrhs, const mxArray *prhs[])
{
int i, j, m ,n;
double *M, *X, *P;
 
m = mxGetM(prhs[0]);
n = mxGetN(prhs[0]);
 
M = mxGetPr(prhs[0]);
X = mxGetPr(prhs[1]);
 
plhs[0] = mxCreateDoubleMatrix(m, n, mxREAL);
 
P = mxGetPr(plhs[0]);
 
j=0;
for (i = 0; i < m*n; i++) 
    {
    j = (i / m);     
    P[i]=M[i]-X[j];
    }
}

A compiler avec MEX :

Code :

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

mex test.c

Et à utiliser comme ceci :

Code :

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

1
2
3
4
5
6
7
M=[1 2
    3 4
    5 6]
 
X=[4 6]
 
P=test(M,X)

Note : c'est l'idée qui est présenté ici.

Après quelques tests rapides sur les différentes méthodes, je dirais que les deux meilleures sont BSXFUN et MEX

[navhpf]

La solution suivante marche aussi. Je soupçonne que ce soit la plus rapide.

M - ones(3,1) * X


On peut bien sûr remplacer le 3 par size(M,1).

[Jerome Briot]

Je soupçonne que ce soit la plus rapide.

As-tu fais des tests avec des matrices de grandes dimensions ?