A=[3,4;1,0] % Matrice 2x2
B=(A>0)     % B = matrice tq Bij=1 si Aij>0 et 0 sinon
V=[5;5]   % Vecteur colonne de dimension 2, s'écrit aussi [5,5]'
length(V) % Renvoie la longueur du vecteur V
size(A)   % La fonction size renvoie la taille
%
%% Opérations matricielles classiques
%
B=A' % B = transposée de A 
A*B  % Multiplication matricielle de A et de B
A+B  % Addition matricielle de A et de B
A^99 % Puissance matricielle
expm(A)  % Exponentielle matricielle
inv(A)   % Matrice inverse
sqrtm(A) % Racine carrée matricielle cf. help sqrtm
logm(A)  % Logarithme matriciel. cf. help logm
A*V  % Image du vecteur colonne V par la matrice carrée A
V'*A  % Devinez...
A\V  % Solution du système linéaire AX=V (par pivot, et pas par inv(A))
help slash  % Aide explicative sur la division matricielle précédente
%
%% Opérations entrée par entrée
%
C=A.*B    % C(i,j)=A(i,j)*B(i,j)
C=A./B    % C(i,j)=A(i,j)/B(i,j)
C=A.^3    % C(i,j)=A(i,j)^3
C=3+A     % C(i,j)=3+A(i,j)
C=3*A     % C(i,j)=3*A(i,j)
C=A./3    % C(i,j)=A(i,j)/3
C=cos(A)  % C(i,j)=cos(A(i,j))
C=log(M)  % C(i,j)=log(A(i,j)).  Ne pas confondre avec logm(A) !
C=sqrt(M) % C(i,j)=sqrt(A(i,j)). Ne pas confondre avec sqrtm(A) !
C=exp(A)  % C(i,j)=exp(A(i,j)).  Ne pas confondre avec expm(A) !
C=abs(A)  % C(i,j)=|A(i,j)|
% Notez qu'à chaque fois, la variable C est redéfinie,
% sans que cela pose problème.
% Listons les variables actuellement définies
whos
% Consultons l'aide sur les opérations élémentaires
help ops
% et sur les opérations élémentaires matricielles
help elmat
%
% Pour finir...
%
n=4;
A=ones(n,n)  % Matrice de taille 4x4 dont tous les éléments valent 1
B=zeros(n,n) % Matrice nulle de taille 4x4
C=eye(n,n)   % Matrice identité de taille 4x4. Que donne eye(3,7) ?
J=[A,B,C]    % J = juxtaposition horizontale des matrices A, B et C
K=[A;B;C]    % K = juxtaposition verticale des matrices A, B et C
N=0*J        % Astuce pour obtenir une matrice N nulle de même taille que J
             % sans connaitre la taille de J.
% Création d'une matrice M 4x3 «creuse» à partir d'un vecteur de données
% [0.5,0.1,0.2], et de vecteurs de position des données [1,2,4] et [2,1,3]
% Donc M(i,j)=0 sauf M(1,2)=0.5, M(2,1)=0.1 et M(4,3)=0.2
M = sparse( [1,2,4], [2,1,3], [0.5,0.1,0.2], 4, 3)