silico.biotoul.fr
 

M1 MABS BBS Math TD Calcul Matriciel

From silico.biotoul.fr

(Difference between revisions)
Jump to: navigation, search
m
m
Line 1: Line 1:
= Création d'une matrice =
= Création d'une matrice =
-
<code lang="shell">
+
<source lang="rsplus">
-
A = matrix( c(11, 12, 13, 21, 22, 23, 31, 32, 33), ncol=3, byrow=TRUE )
+
A = matrix( c(11, 12, 13, 21, 22, 23, 31, 32, 33), ncol=3, byrow=TRUE )
-
</code>
+
</source>
A partir d'un fichier [[File:gold.metadata.txt]] :
A partir d'un fichier [[File:gold.metadata.txt]] :
-
<code lang="shell">
+
<source lang="rsplus">
-
G=read.table("gold.metadata.txt", sep="\t", header=TRUE)
+
G=read.table("gold.metadata.txt", sep="\t", header=TRUE)
-
class(G)
+
class(G)
-
</code>
+
</source>
G est un data frame ; les deux premières colonnes contiennent l'identifiant et le nom de l'organisme. Pour extraire la matrice de données (colonnes numériques), on fait :
G est un data frame ; les deux premières colonnes contiennent l'identifiant et le nom de l'organisme. Pour extraire la matrice de données (colonnes numériques), on fait :
-
<code lang="shell">
+
<source lang="rsplus">
-
as.matrix( G[ , 3:12] )
+
as.matrix( G[ , 3:12] )
-
</code>
+
</source>
-
Cas de la matrice unité d'ordre n (notée <math>I_n</math>), exemple avec n = 5 :
+
Cas de la matrice unité d'ordre n (notée <math>I_n</math>), exemple avec n = 3 :
-
<code lang="shell">
+
<source lang="rsplus">
-
diag(5)
+
I3 = diag(3)
-
</code>
+
</source>
Cas des vecteurs colonne ou ligne :
Cas des vecteurs colonne ou ligne :
-
<code lang="perl">
+
<source lang="rsplus">
-
V = c(1, 0, 5)
+
V = c(1, 0, 5)
-
class(V)
+
class(V)
-
# obtention de la matrice vecteur colonne :
+
# obtention de la matrice vecteur colonne :
-
as.matrix(V)
+
as.matrix(V)
-
# pour obtenir la matrice vecteur ligne :
+
# pour obtenir la matrice vecteur ligne, on fait la transposée du résultat pécédent :
-
t( as.matrix(V) )
+
t( as.matrix(V) )
-
</code>
+
</source>
= Opérations sur les matrices =
= Opérations sur les matrices =
* addition <tt>+</tt>, soustraction <tt>-</tt>, multiplication par un nombre <tt>*</tt>,
* addition <tt>+</tt>, soustraction <tt>-</tt>, multiplication par un nombre <tt>*</tt>,
-
* produit : <tt>A %*% B</tt>, avec A une matrice et B une matrice ou un vecteur
+
* produit : entre 2 matrices A et I3 : <tt>A %*% I3</tt>, entre une matrice A et un vecteur V <tt>A %*% V</tt>
* transposition : <tt>t(M)</tt>
* transposition : <tt>t(M)</tt>
*
*

Revision as of 07:59, 27 September 2011

Création d'une matrice

A = matrix( c(11, 12, 13, 21, 22, 23, 31, 32, 33), ncol=3, byrow=TRUE )

A partir d'un fichier File:Gold.metadata.txt :

G=read.table("gold.metadata.txt", sep="\t", header=TRUE)
class(G)

G est un data frame ; les deux premières colonnes contiennent l'identifiant et le nom de l'organisme. Pour extraire la matrice de données (colonnes numériques), on fait :

as.matrix( G[ , 3:12] )

Cas de la matrice unité d'ordre n (notée In), exemple avec n = 3 :

I3 = diag(3)

Cas des vecteurs colonne ou ligne :

V = c(1, 0, 5)
class(V)
# obtention de la matrice vecteur colonne :
as.matrix(V)
# pour obtenir la matrice vecteur ligne, on fait la transposée du résultat pécédent :
t( as.matrix(V) )


Opérations sur les matrices

  • addition +, soustraction -, multiplication par un nombre *,
  • produit : entre 2 matrices A et I3 : A %*% I3, entre une matrice A et un vecteur V A %*% V
  • transposition : t(M)