DS de modélisation

IUT informatique de Bordeaux

janvier 2019

Le fichier Maple sera ouvert en worksheet mode mais des zones de texte seront insérés pour la
rédaction. Le travail sera rendu sur Moodle sous forme de fichier Maple intitulé DS − maple − nom,
après avoir averti l’enseignant de la remise. Les nom, prénom et groupe seront aussi indiqués en en-
tête du fichier. Le fichier sera agrémenté de zones de texte où toutes les explications nécessaires seront
fournies. Les documents de travail des TP Maple sont autorisés s’ils sont sur la session de l’étudiant.
L’accès internet n’est pas autorisé. (Logs conservés et sessions vérifiables par les administrateurs)

1 Chaı̂ne de Markov, simulation et probabilité

14 points
On lance n fois une pièce de monnaie bien équilibrée et on note la suite des résultats des n lancers
en indiquant P pour pile et F pour face.
Pour une suite de n lancers donnée on observe les séquences obtenues, une séquence est une succession de
Piles ou de Faces consécutifs, la longueur de la séquence étant le nombre de Piles ou de Faces consécutifs.
On donne les états suivants pour cette suite :
— A : elle ne contient pas de séquence de longueur au moins 6 et finit par une séquence de longueur
1.
— B : elle ne contient pas de séquence de longueur au moins 6 et finit par une séquence de longueur
2.
— C : elle ne contient pas de séquence de longueur au moins 6 et finit par une séquence de longueur
3.
— D : : elle ne contient pas de séquence de longueur au moins 6 et finit par une séquence de longueur
4.
— E : elle ne contient pas de séquence de longueur au moins 6 et finit par une séquence de longueur
5.
— F : elle contient une séquence de longueur au moins 6.
On appellera An l’évènement correspondant à l’état A pour n lancers et an = p(An ) etc.

Exemple 1 Liste de 6 lancers [P P F P F P ], liste des longueurs des séquences [2, 1, 1, 1, 1], état A.
Exemple 2 Liste de 7 lancers [P F P F F F F ], liste des longueurs des séquences [1, 1, 1, 4], état D.
Exemple 3 Liste de 12 lancers [P F F F F F F F P F P P ], liste des longueurs des séquences [1, 7, 1, 1, 2],
état F.
1. Simulation

(a) Ecrire une procédure ”Simulation” qui retourne une liste de n lancers aléatoires d’une pièce
de monnaie bien équilibrée. L’appliquer à n = 100, n = 300.
(Création aléatoire d’un booléen : with(ListTools) ;rand(2)). (1)
(b) Ecrire une procédure ”Séquence” qui retourne la liste des longueurs des séquences d’une liste
de n lancers. Tester la procédure sur les exemples donnés. (2)
1
 (c) En déduire une procédure ”Etat” qui retourne l’état d’une liste de n lancers. Tester la procédure
sur les exemples donnés. (2)
(d) Conclure par une procédure ”Fréquence” utilisant les procédures précédentes et qui simule 1000
suites de n lancers et retourne les fréquences des états A,B,C,D,E,F sur ces 1000 expériences.
L’appliquer à n = 100, n = 300. (2)
Pour simplifier les procédures, on pourra identifier l’état à la longueur de la séquence retenue, a=1,b=2,etc.
2. Probabilité
Les questions a, b, c, d, e sont à rendre dans le fichier Maple en mode texte.
an ; bn ; cn ; dn ; en ; fn sont 
donc les
 probabilités respectives qu’une suite de n lancers soit à l’état
an
 bn 
 
 cn 
A ;B ;C ;D ;E ;F et Xn =    , n ∈ N∗ .
 dn 

 en 
fn
(a) Donner a1 ; b1 ; c1 ; d1 ; e1 ; f1 soit X1 . (0.5)
(b) Donner a2 ; b2 ; c2 ; d2 ; e2 ; f2 soit X2 . (0.5)
(c) Déterminer P (An ) en fonction des probabilités des états pour n-1 lancers. (1)
(d) En déduire les relations de récurrence liant Xn et Xn−1 . (1)
(e) Montrer
 que la matrice de transition Tentre Xn et Xn−1 est :
1/2 1/2 1/2 1/2 1/2 0
 1/2 0 0 0 0 0 
 
 0 1/2 0 0 0 0 
T =
  (1)
 0 0 1/2 0 0 0  
 0 0 0 1/2 0 0 
0 0 0 0 1/2 1
(f ) En déduire une procédure permettant de calculer Xn en fonction de X1 . (1)
(g) Utiliser cette procédure pour calculer an ; bn ; cn ; dn ; en ; fn , n = 100, n = 300
Interpréter la valeur de fn obtenue pour ces deux valeurs. (1)
3. Comparer en norme le vecteur fréquence obtenue en 1.d et le vecteur loi de probabilité obtenue en
2.f pour n=300. (1)

2
2 Méthode des trapèzes
7 points
On découpe un intervalle [a; b] en n intervalles égaux. On approche l’intégrale de f sur chacun de ces
intervalles par l’aire du trapèze ayant pour sommets les extrémités de l’intervalle et les images par f de
ces extrémités.

B1 + B2
On rappelle que l’aire d’un trapèze de bases B1; B2 et de hauteur h est donnée par h × .
Rb 2
L’intégrale a f (t)dt, qui est l’aire du domaine sous la courbe, est donc approchée par la somme des
aires des trapèzes (aire totale) ayant chacun pour hauteur h = b−a
n si l’intervalle [a; b] est partagé en n
intervalles égaux.
1. Donner la formule en mode texte sur Maple de l’aire totale Jtrap(f, a, b, n) calculée ainsi entre
a et b, en fonction de f ; a ; b et n. (1)
Cette formule sera donnée sous forme d’une somme.
2. Dans Maple, écrire une procédure Jtrap(f,a,b,n), qui prend en entrée f ; a; b et n et calcule la
somme précédente. (1)
On utilisera la commande sum
Pour les questions suivantes, on prendra

f (x) = x2 + ln(x + 1), a = 0, b = 5

3. Calculer la valeur approchée de l’intégrale à l’aide de la procédure pour n=10. (0.5)

3
 4. En utilisant la bibliothèque Student[Calculus1]
et la commande ApproximateInt(· · · , · · · ,partition = 10, method = trapezoid), retrouver la valeur
approchée de l’intégrale. (0.5)
5. A l’aide de la commande int, calculer la valeur exacte de l’intégrale. (1)
6. Comparer pour n = 10, 100, 500, en écrivant une procédure retournant une liste, les valeurs
approchées de l’intégrale à celles de la commande Maple int qui calcule l’intégrale de f . (1)
7. Expliquer pourquoi il existe un entier p tel que :
Z 5
f (x) dx − Jtrap(f, a, b, p) < 0, 1
0
(1)
8. Écrire une procédure qui détermine le plus petit entier p vérifiant cette propriété. (1)

3 Listes
4 points
1. En faisant attention aux répétitions, écrire la liste des nombres qui sont la somme de deux carrés
d’entiers a2 + b2 avec 0 ≤ a ≤ 5, 0 ≤ b ≤ 5 (2)
2. Créer une nouvelle liste en extrayant de la liste précédente les nombres impairs. Quelle semble
être la forme générale de ces nombres ? (2)
REMISE SUR MOODLE A LA FIN DE LA SEANCE