FICHIER Prof CoursLatexScilabR

Logiciel de traitement de texte scientifique : LATEX 2ε Logiciel Scilab Logiciel R
Logiciels Scientifiques (Latex, Scilab et R)
Mouhamadou DOSSO1
1 Université Félix Houphouët-Boigny de Cocody-Abidjan, Côte d’Ivoire 2
UFR de Mathématiques et Informatique

1 Logiciel de traitement de texte scientifique : LATEX 2ε

Overview
Formules Mathématiques
constructions mathématiques
Images
Bibliographie, Index
2 Logiciel Scilab
Fonctionnement général
Type de données
Programmation
Entrées/Sortie sous Scilab
Sorties graphiques
3 Logiciel R
Interface d’utilisation sous windows
Les objets : Vecteurs, Matrices, Matrices à plus de deux dimensions, liste et structures de données
Quelques fonctions usuelles
Quelques fonctions de statistique exploratoire
Construction d’une nouvelle fonction
Quelques éléments de programmation
Les entrées/Sorties et gestion des objets crées
Les graphiques avec R
Overview
Édition d’un document
On édite un document Latex comme dans les logiciels de

programation. Un exemple d’édition de document en Latex est
ci-dessous :
\documentclass[french,12pt]{article}
\usepackage[T1]{fontenc}
\usepackage{babel}
\begin{document}
bonjour, A tous ! Ici, vous devez écrire le texte du document.
\end{document}
Le texte du document doit-être écrire entre \begin{document}

et \end{document}. Puis, on le sauve sous un nom terminant
par .tex. Exemple: premier.tex
Overview
Compilation et visualisation d’un document tex.
Compilation
Pour la compilation ( avec WinEdit), on clique sur le menu
LATEX ou en forme de tête de chat en couleur verte situé
dans la barre en haut de l’éditeur WinEdit. On peut également
faire la compilation à l’aide du clavier en tapant Shift+Ctrl+L
pour LATEX et Shift+Ctrl+X pour la tête de chat.
N.B.:
En cliquant sur le menu LATEX, on a une simple compilation
qui vérifie les éventuelles erreurs ; mais il ne permet pas la
visualisation du document. Cependant, avec le menu en forme
de tête de chat en couleur verte, on a en plus de la
compilation et la vérification, la visualisation du document s’il
n’y a pas d’erreur.
Overview
Pendant la compilations
Créations d’autres fichiers

Au cour de la compilation d’un fichier LATEX, plusieurs fichiers
seront créés. Pour un fichier qui porte le nom premier.tex, on a
la création de :
premier.log :C’est la transcription détaillée de tout ce qui
s’est passé à la compilation. En particulier des éventuelles
erreurs.
premier.aux :C’est un fichier auxiliaire (parmi plusieurs
autre pour certrains documents plus complexes) très
précieux à LATEX.
premier.dvi : C’est le résultat de la compilation. Celui que
tu visualiseras.
Overview
Après la compilation
Visualisation
Pour visualiser votre document Latex, après compilation sans
visualisation, il vous suffit de cliquer sur le menu en forme de
loupe juste en bas de LATEX. Votre document aura
l’extention .dvi
Exemple : premier.dvi
Vous pouvez également visualiser votre document en fichier

PDF ou en fichier PS ou EPS en cliquant sur les autres icônes
qui suivent les icônes en forme de chat et LATEX.
Overview
Structuration du document
Les classes et les options

Les classes:Il s’agit des termes report, article, book et
letter qui définissent les types de document que que nous
voulons avoir.
Les options : Les options sont placées entre crochets
juste après le nom de la commande ; Exemple :
⋆ \documentclass[12pt]{report} permet de charger la classe
de document report en 12 points.
⋆ \documentclass[12pt]{report}[2012/09/03] permet de
demander à LATEX 2ǫ d’utiliser le style report avec l’option
12pt dans une version datant d’àprès le 03 septembre
2012. On aurra un message lors de compilation si une
version plus ancienne est trouvée.
Overview
Structure du document
Les packages
Les packages: Ce sont les macros de fichiers que nous
souhaitons utiliser ; comme exemple : pour epsfig.sty ou
fancyhdr.sty, on place \usepackage{epsfig} et
\usepackage{fancyhdr} entre \documentclass et
\begin{document}.
N.B.
Pour éviter la répetition de la commande \usepackage, il est
possible de donner une liste de packages à charger, en
séparant leurs noms par des virgules.
Exemple: \usepackage[german,english]{babel} ou
\usepackage{babel,indentfirst}
Overview
Structure de documents types - Rapport type
En-tête du document
\ documentclass[french,twoside,openright]{report}
\ usepackage[T1]{fontenc}
\ usepackage{babel,indentfirst}
\ begin{document}
fontenc :permettre la césure convenable de toute langue

ayant des caractères accentués.
babel : se charge de faire des traductions utiles (exemple
"chapter" en "chapitre")
indentfirst : pour indenter le premier paragraphe suivant
un titre
Overview
Structure de documents type - Rapport type
french : Pour pouvoir passer à tous les packages qui en

auront besoin. C’est-à-dire les packages qui sont
susceptible de produire du texte.
twoside : demande à ce que l’on travaille en recto-verso.
openright : indique que les chapitres doivent commencer
sur des pages impaires ( de doite )
Overview
Le titre
\ title{rapport bidon}
\ author{Benjamin \ textsc{Bayart}
\ and \ textsc{M \∧ eme}
\ and Moi \ TEXTSC {Aussi}
\ and Personne \ TEXTSC {d’Aussi}
\ date{Le \today}
\ maketitle
Overview
Le debut
\ tableofcontents
\ chapter∗ {Introduction}
\ addcontentsline{toc}{chapter}{Introduction}
\ markboth{\uppercase{Introduction}}
{\uppercase{Introduction}}
Overview
Le corps du document
\ part{Étude préliminaire}
\ chapter{on a commencé}
\ section{Le commencement}
\ subsection{Une sous section}
\ paragraph{Un paragraphe}
La fin du rapport
\ chapter∗ {Concusion}
\ end{document}
Les environnements.
Pour écrire une formule (symbole) ou une équation

mathématique, nous avons besoin d’un environnements
mathématiques.
Les environnements
Il y a 6 méthodes pour passer en mode maths :
L’environnement displaymath.
L’environnement math.
\[. . . \], équivalent à displaymath.
$. . . $, équivalent à math.
$$ . . . $$, équivalent à displaymath.
$ . . . $, équivalent à math.
Les environnements.
N.B.:
"displaymath" est prévu pour faire apparaître une
équation seule centrée sur une ligne.
"math" est prévu pour mettre une petite équation dans le
texte.
Package latexsym.
• Lettres Grecques:
\alpha α \beta β \gamma γ

\epsilon ǫ \varepsilon ε \zeta ζ
\theta θ \vartheta ϑ \iota ι
\lambda λ \mu µ \nu ν
\pi π \varphi ϕ \rho ρ
etc · · ·
\varsigma ς \tau τ \upsilon υ
\varphi ϕ \chi χ \psi ψ
\Gamma Γ \Delta ∆ \Theta Θ
\Xi Ξ \Pi Π \Sigma Σ
\Phi Φ \Psi Ψ \Omega Ω
Package latexsym
•Opérateurs binaires.
\pm ± \mp ∓ \times ×

\ast ∗ \star ⋆ \circ ◦
\cdot · \cap ∩ \cup ∪
\oplus ⊕ \ominus ⊖ \otimes ⊗
\bircirc \dagger † \ddagger ‡ etc · · ·
\sqcap ⊓ \sqcup ⊔ \bigtriangleup △
\vee ∨ \wedge ∧ \bigtriangledown ▽
\setminus \ \wr ≀ \triangleleft ⊳
\diamond ⋄ \odot ⊙ \triangleright ⊲
Package latexsym
• Symboles de relation.
\leq ≤ \le ≤ \geq ≥

\equiv ≡ \models |= \prec ≺
\sim ∼ \perp ⊥ \preceq
\simeq ≃ \mid | \ll ≪
\subset ⊂ \subset ⊂ \subseteq ⊆
etc · · ·
\sqsubset ⊏ \sqsupset ⊐ \sqsubseteq ⊑
\parallel k \approx ≈ \bowtie ⊲⊳
. ∼
\neq 6 = \doteq = \cong =
\in ∈ \ni ∋ \smile ⌣
\vdash ⊢ \dashv ⊣
Package latexsym
• Frèches
\leftrightarrow ↔ \leftharpoondown ↽
\longleftrighttarrow ←→ \hooklefttarrow ←֓
\Leftrightarrow ⇔ \leftarrow ←
\Longleftrightarrow ⇐⇒ \uparrow ↑
\longleftarrow ←− \Leftarrow ⇐
etc · · ·
\longrightarrow −→ \Uparrow ⇑
\Longleftarrow ⇐= \mapsto 7→
\Longrightarrow =⇒ \nearrow ր
\Updownarrow m \nwarrow տ
\leftharpoonup ↼ \updownarrow l
Package latexsym
• Constructions mathématiques
\widetilde{abc} g
abc \widehat{abc} d
abc
←−− −−→
\overleftarrow{abc} abc \overrightarrow{abc} abc
\overline{abc} abc \underline{abc} abc
z}|{ \underbrace{abc} abc
\overbrace{abc} |{z}
√abc √
n
\sqrt{abc} abc \sqrt[n]{abc} abc
abc
f′ f′ \frac{abc}{xyc} xyc
Quelques exemples de constructions mathématiques
Voyons des exemples de constructions mathématiques.

Sommes
Pour écrire
n
X
ui
i=0
V
on tape : \[\sum− {i = 0} {n}u− i\]
n
X
D’autre part, si on souhaite taper ui en plein dans le texte,
i=0
on fait appel à la commande \displaystyle en tapant :
V
$\displaystyle \ sum− {i = 0} nu− i$
Une autre commande

P est donnée par \textstyle : équation
dans le texte ni=0 ui
Opérations et fonctions
Cas «Particulier» : la limite
lim un = ℓ
n→+∞
s’obtient en tapant
\[\lim− {n \ rightarrow + \infty}u− n = \ell\]

Fraction, racines et accolades

X
+∞ n
x p
= ex = e2x
n!
n=0
n
X
ui = n0 + u1 + · · · + un
| {z }
i=0 n+1 termes
s’obtient à partir de
\[ \ sum− {n = 0}∧ {+ \ infty} \ frac{x ∧ n}{n!}

= e∧ x = \sqrt{e∧ {2x}}\]
\[ \ sum− {i = 0}∧ nu− i =
\ underbrace{n− 0 + u− 1 + \cdots + u− n}− {n + 1
\ hbox{\scriptsize termes}}\]
On aime avoir de grandes parenthèses en Mathématiques. On

procède de la manière suivante :
Délimiteurs
Pour avoir !
n
X
ui
i=0
on tape
\[\left(\sum− {i = 0}∧ n u− i \ right)\]
N.B.: Le délimiteur est utilisé pour les systèmes d’équations

(une accolade à gauche et rien à droite), comme dans
l’exemple suivant :
\[\left \ {system \ right.\]

Les matrices
Il est très facile d’écrire une matrice. Exemple
\ [A = \left(
\ begin{array}{ccccc}
1 & 2 & 3 & \ cdots & n \ \
2 & 3 & 4 & \ cdots & n + 1 \ \
3 & 4 & 5 & \ cdots & n + 2 \ \
\ vdots & \ vdots & \ vdots &
\ vdots & \ vdots \ \
n & n + 1 & n + 2 & \ cdots & 2n + 1 \ \
\ end{array}
\ right)∧ {2}\]
Les matrices
 2
1 2 3 ··· n
 2 3 4 ··· n+1 
 
 3 4 5 ··· n+2 
A= 
 .. .. .. .. .. 
 . . . . . 
n n + 1 n + 2 ··· 2n + 1
Le ’c’ dans l’argument du ’array’ indique que l’on souhaite une

colonne centrée. On peut préférer indiquer ’l’ pour une colonne
gauche ou ’r’ pour une colonne droite.
Images
Dessins à inclure
Pour inclure un fichier produit par certains logiciels comme Xfig

et gnuplot, on utilisera les commandes prévues par LATEX en
standard :
Inclusion de code LATEX
\ begin{figure}
\ input nom− du− fichier
\ caption{Titre de la figure}
\ end{figure}
Images
Dessins à inclure
Inclure du PostScript (graphics)

Pour l’inclusion des graphics, on prendra soin d’inclure le
package
\usepackage[dvips]{graphics}
Puis la commande
\includegraphics{fichier.ps}
viendra placer à l’endroit courant le dessin contenu dans le

fichier placé en paramétre.
Approche manuelle
Elle repose sur la commande \bibitem et l’environnement

thebibliography
Code Latex
Exemple :
\ begin{thebibliography}{99}
\ bibitem{Jai-Kha} Jairus M. KHALAGAI,
\ textit{Analyse 2}, African Virtual university.
\ bibitem{Fe-Rong} Felice Ronga
\ textit{Analyse réelle néo-post-élémentaire},
Genève, 2004 ap. J.-C..
\ end{thebibliography}
Approche manuelle
donnera :
Visualisation
1 Jairus M. KHALAGAI, Analyse 2, African Virtual university.
2 Felice Ronga, Analyse réelle néo-post-élémentaire,
Approche manuelle
En générale, on utilise la bibliographie pour les références.

Pour ce faire, dans le cadre de ce exemple, on a les
commandes suivantes
\ cite{Jai-Kha}
\ cite{Fe-Rong}
\ cite{Jai-Kha, Fe-Rong}
D’autre part, on peut mettre les trois premières lettres et du

nom de l’auteur et l’année de parution comme référence, on
procède comme suite :
Approche manuelle
Code Latex
Exemple :
\ begin{thebibliography}{WWW 99}
\ bibitem[Kha06]{Jai-Kha} Jairus M. KHALAGAI,
\ textit{Analyse 2}, African Virtual university.
\ bibitem[Ron04]{Fe-Rong} Felice Ronga,
\ textit{Analyse réelle néo-post-élémentaire,
\ end{thebibliography}
Approche manuelle
donnera :
Visualisation
Kha06 Jairus M. KHALAGAI, Analyse 2,African Virtual university.
Ron04 Felice Ronga, Analyse réelle néo-post-élémentaire,
Genève, 2004 ap. J.-C.
Les commandes LATEX
Retenons la commandes suivante permettant de faire une

citation :
\cite[commentaire]{clef1, clef1, · · · }
faire référence au(x) livre(s) ayant comme clef d’accès
clef1,clef2,... avec le commentaire donné, par exemple,
\cite[page 30]{texbook} donnera [58, page 30].
Puis cette autre commande qui indique les fichiers de la base
de données bibliographies à exploiter :
\bibliography{fichier1, fichier2, ...}
indique que les données utilisées sont dans les fichiers

fichier1.bib, fichier2.bib,...
Chapitre 2
LOGICIEL SCILAB
Introduction
Scilab est un logiciel de calcul numérique développé par

l’Institut National de Recherche en Informatique et en
Automatique (INRIA) et distribué gratuitémént sur presque tout
type d’ordinateur (PC Windows, Linux, Unix, MacIntosh). Pour
plus d’information et pour télécharger ce logiciel, vous pouvez
consulter le site Internet de INRIA :
http ://www.scielab.org.
Depuis la version 2.7, Scilab propose un éditeur intégré

permettant d’écrire les scripts sans sortir de l’environnement.
N.B.: On peut dans un premier temps utiliser Scilab comme
une calculette matrielle.
Fonctionnement général
Utilisation de l’aide en ligne
Une aide en ligne sur toutes les commandes SCILAB, les

fonction preprogrammées, les types, les type de données, etc,
est disponible en cliquant sur le bouton help puis sur les
différents sous-menus. On peut aussi taper les demandes
d’aides au clavier, par exemple :
help() : donne des informations générales sur l’aide en
ligne
apropos motclef : donne toutes les fonctions scilab qui
ont un un rapport avec le mot-clef spécifié
help nomdefonction : ouvre la fenêtre d’aide pour utiliser
la fonction Scilab nomdefonction
Type de données
Constantes spéciales
Elles sont précédées du caractères %. Voici la listes des plus

utilisées. La liste complète est obtenue en tapant la commande
who
%F %T %s %nan %inf %t %f %eps %io %i %e %pi

Type de données
Vecteurs
Pour définir un vecteur, la syntaxe est comme suite :

−− > v = [2, −3 + %i, 7] // vecteur ligne
2. − 3. + i 7.
−− > v ′ // vecteur transposé conjugué
ans=
2.
−3. − i
7.
−− > v.′ // vecteur transposé
ans=
2.
−3. + i
7.
Type de données
vecteurs
−− > w = [−3; −3 − %i; 2] //vecteur colonne

w=
2.
−3. − i
7.
′
−− > v + w // somme de deux vecteurs
ans =
−1.
−6. − 2.i
9.
−− > v + w //produit scalaire euclidien
ans =
18.
−− > w ′ . ∗ v //produit des composantes
ans=
−6. 8. + 6.i 14.
Type de données
Matrices
Les matrices suivent la même syntaxe que les vecteurs. Les

composantes des lignes sont séparées par des virgules et
chaque ligne est séparée de l’autre par un point virgule.
−− > // une manière de définir une matrice 3 × 3:
−− > A = [1, 2, 3; 0, 0, atan(1); 5, 9, −1];
−− > // une autre syntaxe pour faire la même chose
−− > A = [1 2 3
0 0 atan(1)
5 9 −1]
Type de données
Principales opérations sur les matrices.
Fonction Description
ones(i, j) matrice remplies de 1
zeros matrice nulle
eye(i, j) matrice identité
toeplitz(u) crée une matrices de Toeplitz
diag(u) matrice diagonale
diag(U) extrait la diagonale de U
triu(A) matrice triangulaire supérière
tril(A) matrice triangulaire inférière
linspace(a, b, n) vecteur de n composante
linsolve(A, b) resoudre Au + b = 0
A\b résoudre Au = b
cond (A) conditionnement d’une matrice
det(A) déterminant d’une matrice
Type de données
rank(A) rang d’une matrice

inv(A) inverse d’une matrice
pinv(A) pseudo inverse d’une matrice
svd (A) valeurs singulière d’une matrice
norm(A) norme matricielle ou vectorielle
u′ transposé conjugué de u
u.′ transposé conjugué de u
u∗v multiplication matricielle
u+v adition matricielle
u−v soustraction matricielle
u. ∗ v multiplication point par point
u./v division point par point
Type de données
Autres objets mathématiques
Un grand nombre d’ojets mathématiques ou logiques

interviennent dans un programme de simulation (entier, réel,
fonction, polynôme, matrice, booléen, ets...)
Instruction Description
a=2 affectation de 2 à a
a = %t affectation du booléen vrai à a
p = poly([%i − %i],′ z ′ ) affectation du polynôme
z = poly(0,′ z ′ );
q = 1 + 3 ∗ z + 4.5 ∗ z ∧ 2 affectation du polynôme
r = p/q affectation d’une fraction ractionnelle
A = [a + 1 2; atan(1) − 3] affectation d’une matrice réelle
rand (3, 3,′ u ′ ) matrice aléatoire 3 × 3
( loi uniforme sur [0, 1]
Type de données
Autres objets mathématiques
rand (3, 3,′ n′ ) matrice aléatoire 3 × 3

(loi gaussienne centrée)
eye(4, 4) matrice identité
ones(1, 2) ; 1 : 2 : 7 vecteurs lignes (1, 1) et (1, 3, 5, 7)
B = toeplitz matrice trangulaire symétrique
C = B′ ∗ B calcul de
D = C. ∗ B multiplication terme à terme
f = C(1 : 4, 1) ou f = C(:, 1) extraction des termes de
la première colonne de C
g =D\f division à gauche de f par D
deff (′ y = f (a, b)′ ,′ y = a + b ′ ) définition de la fonction
f (a, b) = a + b
Programmation
Ediction de scripts
Scilab possède un editeur intégré qu’on obtient en cliquant sur

le bouton Editor( ou Editeur pour la version française). Il gère
les fichiers que vous créez :
- S’il s’agit d’une fonction qui se termine par Endfunction,
Le fichier prend automatiquement le suffixe sci.
- Si c’est un script, il prend le suffixe sce.
Le bouton Execute dans la fenêtre d’édiction permet de lancer
l’exécution du script courant dans la fenêtre de commandes.
Toutes les fonctions, doivent d’abord être chargées ( avec la
commande getf) et utilisées avec la syntaxe d’appel.
Programmation
Fonctions ou macros ("function" or "macros’...)
De manière générale, la syntaxe de définition d’une fonction

externe est :
Visualisation
function [y1 , ..., ym ] = toto(x1 , ..., xn )

..
..
..
endfunction
où toto est le nom de la fonction, x1 , ..., xn , les n arguments

d’entrée et [y1 , ..., ym ] les m argurments de sortie.
Programmation
On peut crée le fichier "angle.sci" contenant les lignes :

Exemple de fonction
function [s] = angle(x, y) //calcul de l’angle en dégrés

s = 180 ∗ atan(y/x)/%pi; //visualise le trangle
xpoly([x, 0, x], [0, 0, y], ”lines”)
endfunction
puis dans la fenêtre de commandes on tape

getf("angle.sci","c");
Programmation
Si une des variables de la procédure n’est pas définie à

l’intérieur de celle-ci et non plus en argument d’entrée, elle
prend la valeur définie dans le programme appelant. Ceci
permet d’appeler les fonctions avec moins de paramètres
d’entrée que prévu. Ainsi
Exemple de fonction
− − > clear
− − > getf("angle.sci")
− − > angle(4)
! − −error
undefined variable : y
at line 2 of function angle called by :
angle(4)
Programmation
donne une erreur car le paramètre y n’est pas affecté. En

revanche
Exemple de fonction
− − > clear
− − > getf("angle.sci")
−−>y =2
y =
2.
− − > angle(4)
ans =
1.1071487
−−>
Programmation
Boucles
Il y a deux type de boucle en Scilab : Les boucles while et les

boucles for. La boucle for parcout un vecteur d’indices et
effectue à chaque pas toutes les instructions délimités par
l’instruction end.
Exemple
−− > x = 1; for k = 1 : 4, x = x ∗ k, end
x= x=
1. 6.
x= x=
2. 24.
Programmation
Boucles
La boucle for peut parcourir un vecteur ( ou une matrice) en

prenant comme valeur à chaque pas les éléments ( ou
colonnes ) successifs.
Exemple
−− > V = [−1 3 0]
−− > x = 1; k = v, x = x + k, end
x=
0
x=
3.
x=
3.
Programmation
Boucles ( while)
Exemple
La boucle while effectue une suite de commandes tant qu’une
condition est satisfaite.
−− > x = 1; while x < 14, x = x + 5, end

x=
6.
x=
11.
x=
16.
Ces deux types de boucle peuvent être interrompus par

l’instruction break.
Programmation
Tests
Dans Scilab on a la structure conditionnelle if-then-else-end

agrémenté du elseif parfois bien utile. La syntaxe est par
exemple
Exemple
−− >x = 16
x=
16.
−− >if x > 0 then, y = −x, else y = x, end
y=
− 16.
Dans le cas où le test doit départager un grand nombre (>2) de

possibilité,...
Programmation
Test (suite)
on peut utiliser le select-case, équivalent du switch-case en

C.
Exemple
−− >n = round (10 ∗ rand (1, 1))

−− >select n
−− >case 0 then
−− > printf (′ cas numero 0’)
−− >case 1 then
−− > printf (′ cas numero 1′ )
−− >else
−− > printf (′ autre cas′ )
−− >end
Programmation
Tests - Les opérations logiques dans le tests
fraçais anglais test Scilab

et and &
ou or |
non not ∼
égal equal ==
différent diferent <> ou ∼=
plus petit que lower than <
plus grand que greater than >
plus petit ou égal à lower than ou equal <=
plus grand ou égal à larger than or equal >=
Programmation
Utilisation de fonctions Scilab
On a les commandes fsolve et ode qui servent à résoudre des

système d’équations non linéaires et équations différentielles
ordinaires respectivement.
- Pour fsolve : La syntaxe la plus simple est :
x=fsolve(x0,fct)
et peut être précisée avec des arguments optionnels :
[x [, v [, inf 0]]] = fsolve(x0, fct [, fjac] [, tol]).
x0 : valeur initiale de l’argument de la fonction
fct et fjac : des fonctions externes
tol :un scalaire réel pour la tolérence sur l’erreur
x : vecteur contenant la solution du système
v: vecteur contenant la valeur de la fonction en x
inf0: indicateur de fin : ( 0 : mauvais paramètre ;1 : calcul
correct avec une erreur inférieure à tol ; 2 :nombre
maximum d’appels à la fonction atteinte ; etc ...)
.
Programmation
Utilisation de fonctions Scilab
- Pour ode : la syntaxe d’appel la plus simple est :

y=ode(y0,t0,t,f) où y0 est le vecteur des conditions
initiales, t0 : le temps initial, t :le vecteur de temps et y : la
matrice contenant la solution
y = [y(t(1)), y(t(2)), ...]
L’argument f de ode est une fonction externe ( syntaxe
d’appel : ydot = f (t, y)) où t est un scalaire réel ( temps )
et y un vecteur de réels.
Syntaxe d’appel avec des paramètres optionnels :
[y, w , iw ] = ode([type], y0, t0, t [, rto1 [, atol]), f [, jac])
[y, rd , w , iw ] =
ode(”root”, y0, t0, t, [, rto1 [, ato1]], f [, jac])
où type (≈ "adams","stiff","rk", "discrete","roots’,etc)
rtol,atol : vecteurs de constantes réelles de même taille
que y.
.
Interation avec l’utiliseur
Pour communiquer avec l’utilisateur du programme Scilab, on

utlise la commande input dont la syntaxe est :
[x] = input(”message”, [”string”])

"message" : est une chaîne de caractère qui s’affiche dans la
fenêtre de commande.
x est la vrariable où le programme stocke la valeur rentrée par
l’utilisateur ;
Pour afficher des résultats avec un format, on utlise la
commande
printf("format",variable)
Interation avec l’utiliseur
Exemple
−− >a = 1/3;
−− >printf(”valeur scalire = %f ”, a);
valeur scalaire = 0.333333
−− >a = 2;
−− >printf(”valeur entière = %d ”, a);
valeur entière = 2
−− >a = ”une chaîne de caractères”;
−− >printf(”pour afficher %s”, a);
pour afficher une chaîne de caractères
−− >b = 1;
Interaction avec l’utilisation
Exemple
−− >a = ”chaîne de caractères”;

−− >printf(”pour afficher %d %s”, b, a);
pour afficher 1 chaîne de caractères
Des "interfaces utilisations" plus élaborées peuvent être

programmées avec les menus. (Voir l’aide en ligne sur les
commandes uimenu, addmenu, etc...)
Sauvegarde de résultats en binaire

bf who :
Pour avoir la liste

save : Pour en avoir la liste
load : Pour sauver des variables dans un fichier.
clear : Pour effacer toutes ou une partie des variable.
Exemple
−− >a = [1 : 2 : 10] //définition du tableau a

a =
1. 3. 5. 7. 9.
−− >//sauvegarde dans le fichier binaire ′ tab− a′
−− >save(′ tab− a′ , a)
−− >clear a //effacement du tableau a
Sauvegarde de résultats en binaire
Exemple
−− >a //a n’existe plus

! − −error 4
undefined variable :a
−− >//récuperation du tableau a sauvé dans ’tab− a’

−− >load(′ tab− a′ )
−− >a
a =
1. 3. 5. 7. 9.
Sorties graphiques
Les tracés en 2 dimensions
.plot2d est la commande de base pour tracer une courbe.

plot2d(x,y,[style,chaîne,legende,rectangulaire]):
- x,y : sont des matrices de tailles [npt, nc] contenant les
abscisses et les ordonnées. nc est le nombre de courbes
et npt le nombre de points ; le même,pour toutes les
courbes.
- style est un vecteur de dimension (1, nc).
- chaîne est une chaîne de trois caractères "xyz".
Sorties graphiques
Ce premier exemple calcule automatiquement le repère

Exemple 1
x = [0 : 1 : 100]
y = sin(2 ∗ %pi ∗ x./100);
z = 3 ∗ y;
plot2d (x, y, [1],′ 121′ ,′ y(x)′ );
on a la sortie graphique suivante :

Sorties graphiques
Dans un autre exemple, le repère est imposé par les valeurs de

cadre:
Exemple 2
xbasc()
cadre = [min(x), min(y), max(x), max(y)];
plot2d (x, y, [1],′ 111′ ,′ y(x)′ , cadre);
on a la sortie graphique suivante :

Sorties graphiques
Graphiques dans l’espace
La commande plot3d représente graphiquement une matrice à

trois dimensions dont l’appel la plus simple est plot3d(x,y,z) où
x, y et z sont trois matrices et z contient les valeurs des points
de coordonnées (x, y).
- param3d pour représenter des courbes paramétrées en
trois dimensions.
- contour pour les courbes d’une fonction 3d donnée par
une matrice.
- fec pour les courbes de niveaux d’une fonction donnée par
ses valeurs nodales sur un maillage trangulaire.
Sorties graphiques
Instruction Desciption
plot2d(x,y) tracé de la courbe passant
par les points (x, y)
plot2d1(’ol1’,x,y) idem avec echelle logarithmique
sur les deux axes
fplot2d(x,f) tracé de la courbe (x, f (x))
champ(x,y,fx,fy) tracé du champ de vecteurs
(fx(x, y), fy(x, y))
plot3d(x,y,z) tracé de la surface passant
par les points (x, y, z)
param3d(x,y,z) tracé de la courbe paramétré
passant en (x, y, z)
contour(x,y,z,n) tracé de n courbes de niveau d’2 surface
Sorties graphiques
histplot(n,data) histogramme de l’échantillon data

divisé en n classes
xbasc() effacement des fenêtres graphiques
xset() modification des options graphiques
xsetech([x1,y1,x2,y2]) découpage d’une fenêtre graphiques
xstring(x,y,’coucou’) inscription de caractères sur une figure
Chapitre 3
LOGICIEL R
Fonction des menus
Ce logiciel R forme un interface utilisateur simple. Elle est

structurée autour d’une barre de menu et de diverses fenêtres.
- Le menu file ( ou fichier ) : contient les outils nécessaires à
la gestion de l’espace de travail ;
- Le menu Edit ( ou Edition ): contient les habituelles
commandes de copier-coller et la boîte de dialogue pour la
personnalisation de l’apparence de l’interface ;
- Le menu Misc traite de la gestion des objets en mémoire
et permet d’arrêter une procédure en cours de traitement.
- Le menu Package : automatise la gestion et le suivi des
librairies de fonctions ;
- Les menus windows ( ou fenêtre ) et Help ( ou Aide ) :
assurent des fonctions similaires à celles qu’ils occupent
dans les autres applications Windows ( les fénêtres, l’aide
en ligne et manuels de références de R)
D’autres fenêtres
Parmi les fenêtres, on a la console qui est la fenêtre principale

où on réalise par défaut les entrées de commandes et sorties
de résultats en mode texte. On y ajoute également un certain
nombre de fenêtres telles que les fenêtre graphiques et les
fenêtres d’informations ( historique des commandes, aide,
visualisation de fichier, etc...).
Ces dernières sont toutes appelées à partir de la console par
des commandes spécifiques.
Objets
Les éléments de base du language R sont des objets qui

peuvent être des données ( vecteurs, matrices,...), des
fonctions, des graphiques,...
Les objets se différencient en mode (qui sont :null ( objet vide
), logical, numeric, complex, character) qui décrivent leur
contenu et leur classe.
Les principales classes d’objets sont : vector, matrix, array,
factor, time-series, data.frame, list
N.B.: On peut avoir des vecteurs, matrices, tableaux,
variables,... de mode null (objet vide), logical, numeric,
complex, character.
Les vecteurs
C’est l’objet de base dans R. La syntaxe de création d’un

vecteur est :
c(a1 , a2 , a3 , ..., an )
où les (ai )1≤i≤n sont les coefficients du vecteur.
Exemple:
> a=c(5,5.6,1,4,-5) Création de l’objet a recevant un
vecteur numérique de dimension 5
et de cordonnées 5, 5.6, 1, 4, −5.
> a Affichage du vecteur a
> a[1] Affichage de la première coordonnée
du vecteur a
> b=a[2:4] création du vecteur numérique b de
dimension 3 et de coordonnées 5.6, 1, 4
Les vecteurs
> d=a[c(1,3,5)] Création d’un vecteur numérique d

de dimension 3
et de coordonnées 5, 1, −5
> sum(d) Calcul de la somme du vecteur d
> length(d) Affichage de la dimension de d
> t(d) Transposition du vecteur d
> e=c(1,1,1) ; t(d)%*%e Produit scalaire entre les vecteurs
d et e
> b=sample(1:10,10) Création d’un vecteur numérique
b de dimension 10 dont les coefficients
sont pris aléatoirement entre les entier
Les matrices
Les matrices comme les vecteurs, sont de modes quelconque,

mais elles ne peuvent pas contenir des éléments de nature
différente. La syntaxe de création d’une matrice est :
matrix(vec,nrow=n,ncol=p,byrow=T)
où vec est le vecteur contenant les éléments de la matrice
rangés en colonne ( sauf si "byrow=T" est choisie).
Exemple
> a=1:20 Création d’un vecteur d’entier
de 1 à 20
> x1=matrix(a,nrow=5) Création d’une matrice numérique
x1 de dimension 5 × 4 ayant pour
prémière ligne 1, 6, 11, 16
> x2=matrix(a,nrow,byrow=T) Création d’une matrice numérique
x2 de dimension 5 × 4 ayant pour
prémière ligne 1, 2, 3, 4
Les matrices
Exemple (Suite)
> x3=t(x2) Tranposition de la matrice x2
> b=x3%*%x2 produit de matriciel entre x3 et x2
> dim(x1) Affichage de la dimension de x1
> eigen(x1) Détermine les valeurs propres et
vecteurs propres de la matrices x1
Les matrices à plus de deux dimensions
Les matrices à plus de deux dimensions sont crées à l’aide de

la commande suivante :
array(vec,c(n,p,q,...))
où
- vec est le vecteur contenant les éléments de la matrice
- c(n,p,q,...) désigne les dimensions :
• n est le nombre de lignes,
• p le nombre de colonnes,
• q le nombre de cotes,
..
.
Les matrices à plus de deux dimensions
> x=array(1:50,c(2,5,5)) Construction de la matrice x

> x affichage de de x
> x[1,2,2] extrait le coef. placé aux indices
1,2 et 2
> dim(x) détermine les dimensions de x
> aperm(x) Transposition généralisée de x,
x[i,j,k] devient x[k,j,i]
Les listes
La construction d’une liste passe par la fonction
list(nom1=el1,nom2=el2,...) ,
l’utilisation des noms étant étant facultative.
>li = list(num = 1 : 5, y = ”couleur ”, a = T )

>li
>li$num ou li[[1]]
>li$a ou li[[3]]
>a = matrix(c(6, 2, 0, 2, 6, 0, 0, 0, 36), nrow = 3)
>res = eigen(a, symmetric = T ) Diagonalisation de a
>res$values
>res$vectors
Les structures de données
Les tableaux de données (data.frame) constituent une classe

partitulière de listes consacrée au stockage des données
destinées à l’analyse. Pour créer un tableau de données, on
peut régrouper les tableaux de même longueur à l’aide de la
commande
data.frame(nom1=var1,nom2=var2,...).
On peut ainsi transformer une matrice en tableau de données

en utilisant la commande as.data.frame(mat).
Les structures de données
Exemple
>v1 = sample(1 : 12, 30, rep = T ) Echatillonnage avec

remise dans les entiers de 1 et 12
>v2 = sample(LETTERS[1 : 10], 30, rep = T )
>v3 = runif (30) 30 réalisations indépendantes d’une loi
uniforme sur [0, 1]
>v4 = rnorm(30) 30 réalisations indépendentes d’une loi
normale de moyenne 0 et variance 1
>v1; v2; v3; v4
>xx = data.frame(v1, v2, v3, v4)
>xx
Les distributions usuelles
Loi Nom Paramètres Val. par défaut

Beta beta shape1,shape2
Binomiale binom size,prob
Cauchy cauchy location,scale 0,1
Khi-Deux chisq df
Explonentielle exp 1/mean 1
Fisher f df1,df2
Gamma gamma shape,1/scale -1
Géométrique geom prob
Hypergéométrique hyper m,n,k
Log-Normale lnom mean,sd 0,1
Logistique logis location,scale 0,1
Normale norm mean,sd 0,1
Poisson pois lambda
Les distributions usuelles
Student t dt
Uniforme unif min,max 0,1
Weibull weibull shape
Pour chacune de ces distributions, on dispose de quatre
commande préfixées par une des lettres d,p,q,r et suivi du
nom de distribution :
- dnomdist : fonction de densité pour une distribution de
probabilité et de la fonct. de probabilité (P(X = k));
- pnomdist :fonction de répartition (P(X ≤ x));
- qnomdist : fonction de quantité ;
- rnomdist :génère des réalisations aléatoires
indépendantes de la distribution nomdist.
> data(women)
> names(women)
> attach(women)
> mean(height) Calcul de la moyenne empirique de
variable quantitative height
> var(height) Calcul de la varance empirique de height
estimateur non biaisé (diviseur n − 1)
> sd(height) Calcul de l’écart-type de height
> median(height) Calcul de la médiane empirique de height
> quantile(height) Calcul des quantiles empiriques de height
> summary(weight) Résumé de weight
> summary(women) Résumé de women
> hist(weight,nclass=15) Histogramme de weight constitué

de 15 classes
> cor(height,weight) Calcul du coefficient de corrélation
linéaire
> empirique entre weight et height
> v1=rnorm(100)
> hist(v1)
> v2=factor(sample(letters[1:4],100,rep=T))
> table(v2) Résumé de la variable qualitative v2

> barplot(table(v2)) Diagramme en barre de v2
> sd(height) Diagramme en secteur v2
On peut définir une nouvelle fonction soit directement à partir

de la console, soit via un éditeur de texte externe grâce à la
commande fic(nom− fonction). La seconde possibilité permet la
correction du code en cours d’édition, tandis que la prémière
s’éffectue ligne par ligne, sans retour en arrière.
La syntaxe générale de la définition d’une nouvelle fonction par
par l’expression
nom− fonction=function(arg1[−expr 1], arg2[−expr 2], ...){

bloc d’intructions
}
Les accolades définissent le début et la fin du code source de

la fonction ; les crochets ne font pas partie de l’expression mais
indiquent le caractère facultatif des valeurs par défaut des
arguments.
On peut créer également une fonction personnalisée à partir

d’une fonction existente grâce
nom− fonction2=edit(nom− function1); fix(nom− fontion2)
Exemple
>x = 2
>carre = function(x) { x = x ∗ x; x }
>carre(2)
>x
>fix(carre)
On peut ajouter des commentaire au code en les faisant

précéder du symbole #
Cette partie concerne les commandes if, while, for

Exemple
>bool = T
>i = 0 ; while(bool == T ) {i = i + 1; if (i > 10) {bool = F }}
>i
>s = 0 ; for (i in 1 : 1000) {s = s + x[i]}
>s
>un = rep(1, 10000)
>t(un)% ∗ %x
>s = 0
Exemple (suite)
>system.time(for (i in 1 : 10000) {s = s + x[i]})[3]

>system.time(t(un)% ∗ %x)[3]
N.B.: Le logiciel R peut avoir des problème de mémoire de ce

fait, il est préférable de les remplacer par les outils de calcul
matriciel.
Les entrées/Sorties
Les données que l’on souhaite analyser proviennet de source

externes sous forme de fichiers. Les objets crées doivent
pouvoir être sauve gauvegardés dans des fichiers afin d’être
transportables.
- Formats propriétaire: save() autorise la sauvegarde de
n’importe quelle liste d’objet en mémoire ;
- Fichiers textes ASCII: ils sont prisent en charge par la
commande scan() ;
- Logiciels statistiques : La librairie foreign offre ces outils
pour une sélection des logiciels statistiques plus courants,
à savoir MINITAB, S-PLUS,SAS, SPSS et STATS.
gestion des objets crées
Lors de l’exécution de R, les fichiers .RData et .Rhistory sont

automatiquement crées dans le répertoire de travail. Losqu’on
quitte R à l’aide de la commande q(), les nouveaux objets
crées peuvent être sauvegardé dans les fichier : .Rdta ou
.Rhistory en mode texte.
- history(): permet de visualiser la suite de commandes que
l’on a tapée
- getwd(): permet de conaitre le répertoire de travail courant
au cours d’une session ;
- ls() : permet de visualiser la liste des objets crées
- rm(): permet de détruire des objets
Du graphiques avec R
Les principaux graphiques sont donnés dans le tableaux

suivant:
plot(x) Tracer le graphe des valeurs de x ordronnées
sur l’axe des abscisses.
plot(x,y) Tracer le graphe de y en fonction de x
sunflowerplot(x,y) Idem mais les points superposés sont
dessinés sous forme de fleurs dont le nombre
de pétales correspond au nombre de points.
pie(x) Tracer un graphe en camembert
boxplot(x) Tracer le graphe en boîtes et moustaches
de x.
stripplot(x) Tracer le graphe des valeurs de x sur une ligne.
pairs(x) Tracer le graphe des moyennes de x en
fonction des valeurs des facteurs f 1 (sur l’axe
des abscisses) et f 2 ( plusieurs graphes).
plot.ts(x),ts.plot(x) Tracer le graphe bivarié de x et y pour

chaque valeur de z.
hist(x,fred=T) Tracer un histogramme des fréquences de x.
barplot(x) Tracer un histogramme des valeurs de x.
qqnorm(x) Tracer un histogramme des valeurs de x.
qqplot(x,y) Tracer les quantiles de x en fonction de ceux
attendus d’une loi normale.
contour(x,y,z) Tracer des courbes de niveau.
filled.contour(x,y,z) Idem mais les aires entre les contours
sont colorées.
image(x,y,z) Idem mais en couleur.
persp(x,y,z) item mais en 3D.
symbols(x,y,...) Dessiner aux coordonnées données par x et
y des symboles (étoiles, cercles, boxplots...).
Pour ces fonctions, les options principales sont les

suivantes :
axes=TRUE ou FALSE Si TRUE, les axes et le cadre sont
tracés.
type="n","p","l","b,"h","s"... Précise le type de graphe dessiné.
type="n" supprime le graphe.
col="blue", col.axis, col.main... Précise la couleur du graphe,
des axes, du titre...
bg="yellow" Précise la couleur du fond.
xlim=c(0,10),ylim=c(0,20) Précise les limites des axes.
xlab="abscisse", Précise les annotations des axes
ylab="ordonnée" Précise les annotations des axes.
main="titre" Précise le titre du graphe.
sub="sstitre" Précise le sous-titre du graphe.
bty="n","o"... Contrôle comment le cadre est
tracé. bty="n" supprime le cadre.
cex=1.5, cex.axis, cex.main... Contrôle la taille des caractères.

font=1, font.axis... Précise la police du texte.
las="0" Contrôle comment sont orientées
les annotations des axes.
lty="1" Contrôle le type de lignes tracées.
lwd=1.5 Contrôle la largeur des lignes.
pch="+","o"... Contrôle le type de symbole utilisé
pour le tracé des points.
ps=1.5 Contrôle la taille en points du texte
et des symboles.
tck, tcl Précise la longueur des graduations
sur les axes.
mfcol=c(3,2), mfrow=c(3,2) Partitionne le graphe en 3 lignes
et 2 colonnes. Les figures sont remp
colonnes par colonnes ou lignes par

FICHIER Prof CoursLatexScilabR

Transféré par

Droits d'auteur :

Formats disponibles

FICHIER Prof CoursLatexScilabR

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

FICHIER Prof CoursLatexScilabR

Transféré par

Droits d'auteur :

Formats disponibles

Logiciel de traitement de texte scientifique : LATEX 2ε Logiciel Scilab Logiciel R

Logiciels Scientifiques (Latex, Scilab et R)

UFR de Mathématiques et Informatique

1 Logiciel de traitement de texte scientifique : LATEX 2ε

Édition d’un document

On édite un document Latex comme dans les logiciels de

Le texte du document doit-être écrire entre \begin{document}

Compilation et visualisation d’un document tex.

Créations d’autres fichiers

Vous pouvez également visualiser votre document en fichier

Les classes et les options

Structure de documents types - Rapport type

fontenc :permettre la césure convenable de toute langue

Structure de documents type - Rapport type

french : Pour pouvoir passer à tous les packages qui en

Structure de documents types - Rapport type

Structure de documents types - Rapport type

Structure de documents types - Rapport type

Pour écrire une formule (symbole) ou une équation

\alpha α \beta β \gamma γ

\pm ± \mp ∓ \times ×

\leq ≤ \le ≤ \geq ≥

Quelques exemples de constructions mathématiques

Voyons des exemples de constructions mathématiques.

Une autre commande

Quelques exemples de constructions mathématiques

\[\lim− {n \ rightarrow + \infty}u− n = \ell\]

Quelques exemples de constructions mathématiques

Fraction, racines et accolades

\[ \ sum− {n = 0}∧ {+ \ infty} \ frac{x ∧ n}{n!}

Quelques exemples de constructions mathématiques

On aime avoir de grandes parenthèses en Mathématiques. On

N.B.: Le délimiteur est utilisé pour les systèmes d’équations

\[\left \ {system \ right.\]

Quelques exemples de constructions mathématiques

Quelques exemples de constructions mathématiques

Le ’c’ dans l’argument du ’array’ indique que l’on souhaite une

Pour inclure un fichier produit par certains logiciels comme Xfig

Inclure du PostScript (graphics)

viendra placer à l’endroit courant le dessin contenu dans le

Elle repose sur la commande \bibitem et l’environnement

En générale, on utilise la bibliographie pour les références.

D’autre part, on peut mettre les trois premières lettres et du

Les commandes LATEX

Retenons la commandes suivante permettant de faire une

\bibliography{fichier1, fichier2, ...}

indique que les données utilisées sont dans les fichiers

Scilab est un logiciel de calcul numérique développé par

Depuis la version 2.7, Scilab propose un éditeur intégré

Utilisation de l’aide en ligne

Une aide en ligne sur toutes les commandes SCILAB, les

Elles sont précédées du caractères %. Voici la listes des plus

%F %T %s %nan %inf %t %f %eps %io %i %e %pi

Pour définir un vecteur, la syntaxe est comme suite :

−− > w = [−3; −3 − %i; 2] //vecteur colonne

Les matrices suivent la même syntaxe que les vecteurs. Les

Principales opérations sur les matrices.

rank(A) rang d’une matrice

Autres objets mathématiques

Un grand nombre d’ojets mathématiques ou logiques