Ensait - A3 - Aide à la décision (2)

Date: Évaluation du 2008-12-03 (14h00)

durée 2h, tous documents autorisés
login: azevedo ; passwd fleury

Rappel des consignes

1. La première qualité que l'on attend d'un ingénieur est de savoir présenter ses conclusions. Les listings, les calculs, les graphiques et autres "sorties d'ordinateur" ne peuvent en aucun cas remplacer un relevé de conclusions, rédigé de façon précise et scientifique.
2. Ces différents types de matériaux doivent être ordonnés de façon à produire un document organisé, fournissant une liste de réponses argumentées à la liste des questions de l'énoncé.
3. La méthode RECOMMANDÉE est de rédiger les réponses à la main et de faire des couper-coller avec des ciseaux et de la colle pour insérer les "sorties d'ordinateur" sur votre copie. L'expérience confirme que cette méthode est à la fois la plus rapide et la plus robuste. C'est aussi celle qui permet de se concentrer sur les compétences évaluées et non sur les compétences/incompétences en dactylographie.
4. Un minimum de rigueur orthographique et grammaticale est attendu. Il sera considéré qu'un étudiant incapable de comprendre la grammaire élémentaire n'est probablement pas capable de comprendre des notions d'un niveau ingénieur.
5. Les copies sont évaluées en proportion de leur contenu et pas en fonction de leur poids. Dans cet ordre d'idées, il est demandé d'imprimer les procédures que vous avez modifiées (i.e dont vous êtes co-auteur) et pas la totalité des bibliothèques (waits=22) mises à votre disposition.
6. Chaque simulation doivent être initialisée avec la valeur ddmmyy où dd/mm/19yy est votre date de naissance.
7. Enfin, l'attention des étudiants est attirée sur le fait que le trafic réseau de leur ordinateur est susceptible d'être enregistré pendant la durée de l'évaluation.

1 Petit essai

1. Modifier la fonction petitessai de façon à ce que les services soient uniformément répartis dans $\left[18,\,28\right]$ et les inter-arrivées uniformément réparties dans $\left[9,\,22\right]$. Ne pas oublier de changer la valeur de initseed.
2. Lancer le programme et en commenter le résultat. Ce commentaire s'appuiera sur le chronogramme des événements.

.../...

2 Simulation M/Ga/1

Utiliser le programme 'système M/Ga/1' avec les réglages $marv=65,\, msrv=55,\, a=2.4,\, N=5000$ et $seed=ddmmyy$. Le processus d'interarrivée $X$ est alors une loi exponentielle de paramètre $\lambda$ et le processus de service $Y$ est une loi Gamma de paramètres $a$ et $b$ tels que $\mathrm{E}\left(X\right)=marv$ et $\mathrm{E}\left(Y\right)=msrv$.

1. Déterminer le débit d'arrivée $\lambda$ à partir de $marv$ et donner la distribution (pdf) de probabilité associée. Quelle est la variance du processus d'inter-arrivée ?
2. Rappeler quelle est la distribution de probabilité associée à la loi Gamma de paramètres $a$ et $b$. Déterminer la valeur de $b$ en fonction de $a$ et de $msrv$. Que vaut la capacité de service $1/\mu$ ? Quelle est la variance des services ?
3. Rappeler ce qu'est le paramètre $\rho$ et donner sa valeur.
4. Rappeler ce que signifient les graphes "1-Durée des états" et "2-États lors de l'arrivée". Pourquoi sont-ils si ressemblants dans l'exemple étudié ?
5. Le programme affiche diverses moyennes, écart-types et auto-corrélations. Pour ce qui est des services, donner l'intervalle de confiance à $96\%$ concernant la moyenne. Même question pour la variance (on signale que le $M^{4}$ du cours vaut $M^{4}=3a\left(a+2\right)b^{4}$ pour une loi Gamma).
6. Pour ce qui est des services résiduels, considérer le fait que la loi en masse associée à la loi $Gamma\left(a,b\right)$ est à son tour une loi Gamma, de paramètres $a+1,\, b$ et obtenir les intervalles de confiance correspondants.
7. Quelle serait la durée de séjour moyenne dans une file M/M/1 ayant les mêmes valeurs de $\lambda$ et $\mu$ ? Expliquer pourquoi la valeur obtenue pour la file M/Ga/1 est inférieure.
8. Comparer le nombre de services résiduels constatés et le nombre de clients entrés dans le système. Expliquer la proportion constatée en la reliant à $\lambda$ et $\mu$.

3 Chi-square test

1. Utiliser la commande
   hndl=mopen('dat_services.txt', 'rb'); sss=mget(N, 'f', hndl); mclose(hndl);
   pour récupérer les durées de service engendrées pendant la simulation précédente.
2. Procéder à un test du $\chi^{2}$ pour vérifier s'il est raisonnable d'admettre que ces nombres ont effectivement été simulés selon la loi gamma caractérisée par $a$ et $msrv$. On utilisera pour cela une partition en dix classes bien choisies.

4 Simulation à trois serveurs

On considère maintenant un système M/Ga/3, avec les mêmes arrivées exponentielles que précédemment ($marv=65$) et trois serveurs de loi Gamma, assurant la même capacité totale de service $\mu$.

1. On conserve la valeur du paramètre $a$. Quelle est la nouvelle valeur du paramètre $b$ (pour chacun des serveurs).
2. Modifier les lignes 5 et 7 du programme menu_MGa3 pour tenir compte des nouvelles valeurs de $marv$, $a$, $b$. Exécuter le programme menu_MGa3, avec $nb\_eve=1000$, $nb\_bat=40$ et la politique 2 (tour). Ce programme donne une série d'encadrements. Comparer avec les résultats précédents (serveur unique).
3. Exécuter le programme menu_MGa3, avec $nb\_eve=1000$, $nb\_bat=40$ et la politique 4 (best). Ce programme donne une série d'encadrements. Comparer avec les deux séries de résultats précédents.
4. Reproduire et commenter le graphe 6 (états cumulés)