M=ones(3,3)+eye(3,3)
V=[1:10]
det(M)   % Déterminant de M
rank(M)  % Rang de la matrice 
trace(M) % Trace de M
D=eig(M) % Renvoie le vecteur colonne des valeurs propres de M
[P,D]=eig(M) % Diagonalise M = P*D*P^-1 avec D diagonale
null(M) % Renvoie une BON du noyau («null space» en anglais)
tril(M)    % Partie triangulaire inférieure (lower) de M
tril(M,-1) % Idem sans la diagonale, cf. help tril
triu(M)    % Partie triangulaire supérieure (upper) de M
triu(M,+1) % Idem sans la diagonale, cf. help triu
diag(M)    % Diagonale de M, dans un vecteur colonne
diag(V)    % Matrice diagonale dont la diagonale est le vecteur V
blkdiag(eye(3,3),M) % Matrice block diagonale de blocks eye(3,3) et M
%
sum(V)   % Somme des éléments du vecteur V
sum(M)   % Somme des colonnes de la matrice M (renvoie un vecteur ligne)
sum(M,2) % Somme les éléments de la matrice M selon la dimension 2
sum(sum(M)) % Somme totale des éléments de la matrice M
cumsum(V)   % Sommes cumulatives des entrées du vecteur V
cumsum(M)   % Matrice des sommes cumulatives des colonnes de M
cumprod(V)  % Produits cumulatifs des entrées du vecteur V
cumprod(M)  % Matrice des produits cumulatifs des colonnes de M
max(M)      % Renvoie un vect. ligne = au max sur chaque col. de M
max(max(M)) % Maximum des entrée de la matrice M. Que donne min ?
range(V)    % Renvoie l'étendue de V, c'est-à-dire max(V)-min(V)
%% Réplication de matrice par blocs
repmat(M,3,4) % Matrice par block 3x4 où chaque block = M
              % M n'a pas besoin d'être carrée.
%% Redimentionnement de matrice, colonne par colonne
reshape([1,2,3,4],2,2) % Donne la matrice 2x2 [1,3;2,4]
reshape([1,2,3;4,5,6],3,2) % Donne la matrice 3x2 [1,5;4,3;2,6]
%% Aide sur les opérations matricielles élémentaires 
%% et sur les fonctions matricielles
helpwin elmat
helpwin matfun
