FR 2805945 A1 20010907 FR 0002656 A 20000301 Surveillance et simulation perfectionnées de systèmes complexes, notamment de mécanismes et de contrôles de flux et de congestions dans des réseaux de communication L'invention concerne la surveillance et la simulation de systèmes complexes. Dans le cadre du contrôle de flux et de congestions dans des réseaux de communication, notamment de type internet, une analyse fine du débit offert est souhaitée pour estimer les influences respectives des paramètres du réseau. Avec le développement des techniques de communications en ligne et les problèmes rencontrés, notamment de congestion, différents protocoles de contrôle de flux et de congestion ont vu le jour, notamment le contrôle TCP (de l' anglais "transmission control protocole"). On connaît des procédés d' analyse de ces protocoles. Parmi ces procédés connus, un procédé basé sur une expression mathématique adaptée du débit dans un protocole de type TCP a permis une approche analytique du contrôle. Le principe sur lequel repose ce procédé est décrit notamment dans : M. Mathis, J. Semske, J. Mahdavi and T. Ott, "The Macroscopic Behavior of the TCP Congestion Avoidance Algorithm", Computer Communication Review, 27 (3), July (1997). Ce procédé, quoique prometteur, a montré ses limites dans des applications pratiques, notamment le fait que le caractère aléatoire du trafic ne soit que partiellement pris en compte, ou encore qu' il nécessite une approximation de tous les noeuds du réseau à un unique noeud équivalent, virtuel. Un autre procédé connu, prenant mieux en compte l'approche stochastique, a permis d'appréhender le caractère aléatoire du trafic. Ce procédé, plus récent, découle du principe décrit dans les travaux suivants : Padhye, J., Firiou, V., Towsley, D. and Kurose, J. "Modeling TCP throughput : simple model and its empiri- cal validation", Proc. of ACM SIGCOMM (1998) Padhye J., Firiou V., Towsley D., "A Stochastic Model of TCP Reno Congestion Avoidance and Control", Technical Report, 99-02, CMPSCI, Univ. of Massachusetts, Amherst (1999). Cependant, ce procédé a aussi montré ses limites, notamment le fait qu'il nécessite encore une approximation de tous les noeuds du réseau à un noeud équivalent. La présente invention vient améliorer la situation. Selon une approche différente, l'invention propose d'utiliser une représentation dans l'algèbre dite "max-plus" de systèmes complexes, tels que des réseaux de communication et notamment du contrôle de flux et de congestions. Pour obtenir le détail des principes mathématiques sur lesquels repose une telle représentation, on pourra se référer à l'ouvrage suivant : F. Baccelli, G. Cohen, G.J. Olsder, and J. P. Quadrat, Synchronization and Linearity, Wiley (1992). Globalement, l'algèbre max-plus scalaire est un semi-anneau sur la droite réelle où l'addition devient la fonction "max" (valeur la plus grande parmi un ensemble de valeur) et la multiplication, la fonction "plus" (somme). L'utilisation de l'algèbre max-plus permet de ramener les calculs d'un système compliqué à une simple représentation matricielle. La Demanderesse a montré et vérifié en pratique que l' utilisa- tion de l'algèbre max-plus s'adapte de façon très satisfaisante à la surveillance et à la simulation de systèmes tels qu'un réseau de communications, contrôlé ou non. Elle permet notamment de surmonter le caractère aléatoire des paramètres du réseau, tout en considérant une pluralité de noeuds. De plus, la Demanderesse a montré que la représentation d'un réseau utilisant un protocole TCP était linéaire dans l'algèbre max-plus, ce qui permet, en pratique, d'appliquer des traitements de données simples. La présente invention porte alors sur un dispositif d'aide à la surveillance et/ou de simulation d'un système complexe, notamment d'un réseau de communication. Selon une première importante de l'invention, le dispositif comprend : - une mémoire pour stocker de premières données représentatives de paramètres du réseau, ainsi que pour recevoir au moins de secondes données représentatives d'événements dans le réseau, une portion de ladite mémoire étant réservée au stockage de données sous forme matricielle, - un module de calcul, apte à effectuer sur au moins deux matrices de structure dynamiquement variable, une opération formant produit selon l'algèbre dite MAX-PLUS, - un module de modélisation pour construire au moins une première matrice et une matrice courante respectivement en fonction des premières données et des secondes données, selon un modèle choisi, et - un module pilote pour appliquer répétitivement la première matrice et la matrice courante au module de calcul, la matrice produit obtenue devenant une nouvelle première matrice. Préférentiellement, les premières données comprennent des informations sur la topologie du réseau, telles que le nombre de routeurs traversés par la connexion à surveiller ou à simuler, les propriétés de ces routeurs (tailles mémoires dites "buffers", ou autre), les propriétés statistiques des trafics offerts dans le réseau, etc. Selon une autre caractéristique avantageuse de l'invention, le module de modélisation comprend : - un sous-module de modélisation statique, pour construire la première matrice en fonction des premières données, et - un sous-module de modélisation dynamique, pour construire au moins une matrice courante en fonction des secondes données. Avantageusement, le dispositif de l'invention est apte à traiter des matrices comportant des coefficients dynamiquement variables, dont au moins la matrice courante précitée. Dans le traitement qu'effectue le dispositif selon l' invention, les matrices construites sont avantageusement de même dimension. Ces secondes données comprennent préférentiellement des informations relatives à des pertes dans le réseau, à des flux transverses dans le réseau, par rapport à une connexion contrôlée à surveiller ou à simuler, à des congestions dans le réseau, ou encore à des dépassements de délais dans le réseau dits "time-out". Selon une autre caractéristique préférentielle avantageuse, la matrice produit obtenue est un vecteur représenté par une matrice à colonne unique, ce qui permet de limiter les traitements et leur durée. La première matrice, représentative des paramètres du réseau, est avantageusement structurée au départ comme un vecteur. Selon une seconde caractéristique importante de l'invention, la matrice produit obtenue est représentative d'un débit dans le réseau associé à la connexion à surveiller ou à simuler, d'un débit moyen dans le réseau, ou encore de fluctuations d'un débit instantané dans le réseau. Selon une troisième caractéristique importante de l'invention, le module de modélisation est agencé pour construire successivement une pluralité de matrices, en nombre correspondant sensiblement au nombre de paquets dans le réseau. Le moèle choisi comprend préférentiellement la considération ie la taille variable d'une fenêtre utilisée pour contrôler le nombre de paquets dans le réseau. Il peut s'agir d'un réseau contrôlé par un protocole de type TCP, comprenant typiquement des routeurs à discipline de type "premier arrivé premier servi", ou encore des routeurs à discipline de type WFQ (de l' anglais "weighted fair queuing" ) . Le protocole TCP contrôlant le réseau peut aussi bien être basé sur un modèle de Reno ou un modèle de Tahoe, comme on le verra plus loin. Le service du réseau peut être déterministe, ou encore aléatoire, comme on le verra en détail plus loin. La présente invention vise aussi un procédé d'aide à la surveillance d'un système complexe, notamment d'un réseau de communication. Un tel procédé comprend globalement les étapes suivantes : a) obtenir des premières données représentatives de paramètres du réseau, b) construire une première matrice, selon un modèle choisi, en fonction desdites premières données, c) recevoir, à un instant choisi, au moins de secondes données représentatives d'événements dans le réseau, d) construire au moins une secorae matrice de structure dynamiquement variable, selon le mocèle choisi, en fonction des secondes données, et e) effectuer sur lesdites matrices une opération formant produit selon l'algèbre dite MAX-PLUS, la matrice produit obtenue étant représentative de l'état @@ réseau audit instant choisi. S'il est souhaité de suivre une évolution temporelle de l'état du réseau à des instants choisis, le procédé comporte avantageusement l'étape supplémentaire suivante : f) répéter, à des instants choisis, les étapes c), d) et e), tandis que la matrice produit obtenue devient la première matrice après l'étape e). La présente invention vise aussi un procédé de simulation d'un système complexe, notamment de mécanismes et de contrôles de flux et de congestions dans un réseau de communication. Ce procédé comprend globalement les étapes suivantes : a) obtenir des premières données représentatives de paramètres propres au réseau, b) construire une première matrice, selon un modèle choisi, en fonction desdites premières données, c) simuler des événements dans le réseau et prévoir au moins de secondes données représentatives desdits événements, d) construire au moins une seconde matrice selon le modèle choisi, en fonction desdites secondes données, et e) effectuer sur lesdites matrices une opération formant produit selon l'algèbre dite MAX-PLUS, la matrice produit obtenue étant représentative d'un état du réseau subissant lesdits événements. Pour prévoir une évolution de l' état du réseau en fonction des événements qu'il subit, ce procédé comporte avantageusement l'étape supplémentaire suivante : f) répéter, pour des événements successifs, les étapes c), d) et e), tandis que la matrice produit obtenue devient la première matrice après l'étape e). D'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après, et des dessins annexés sur lesquels : - la figure 1A représente schématiquement un dispositif au sens de la présente invention, - la figure 1 représente schématiquement un nombre K de files en tandem dans un réseau, avec contrôle de flux, - la figure 2A représente un exemple d'évolution pas-à-pas de dateurs et de la taille d'une fenêtre utilisée pour contrôler le nombre de paquets dans le réseau, - la figure 2 représente des interactions entre plusieurs paquets dans le réseau, - la figure 3 représente une variation du débit (courbe en traits pleins), obtenue par simulation, dans un réseau de protocole TCP basé sur le modèle de Tahoe sans la phase exponentielle, - la figure 4 illustre une interprétation graphique des débits asymptotiques dans un réseau de protocole TCP basé sur le modèle de Reno avec délais déterministes, - la figure 5 représente une variation du débit, obtenue par simulation et montrant une décroissance du débit en cas de pertes aléatoires dans le réseau, - la figure 6 représente une variation du débit (courbe en traits pleins), obtenue par simulation, dans un réseau de protocole TCP basé sur un modèle de Reno markovien, - la figure 7 représente des variations de débit comparées, obtenue par simulation, dans des réseaux de protocoles TCP basés respectivement sur un modèle de Reno déterministe RD (traits pleins), de Reno markovien RM (traits pointillés longs), de Tahoe déterministe TD (traits pointillés moyens) et de Tahoe markovien TM (traits pointillés courts), - la figure 8 représente une variation du débit (courbe en traits pleins), obtenue par simulation, dans un réseau de protocole TCP basé sur un modèle de Tahoe avec phase exponen- tielle, - la figure 9 représente des variations comparées de débits, obtenues par simulation, dans un réseau de protocole TCP basé sur un modèle de Tahoe avec des services #3 et 08 respective- ment constants et égaux à 1, et la figure 10 représente schématiquement un réseau avec ses files et ses routeurs. L'annexe I comprend les références bibliographiques indexées entre crochets dans la description ci-après. Les dessins et la description ci-après contiennent, pour l'essentiel, des éléments de caractère certain. Ils pourront donc non seulement servir à mieux faire comprendre la présente invention, mais aussi contribuer à sa définition, le cas échéant. En se référant à la figure 1, le dispositif se présente sous la forme d'un ordinateur comprenant une unité centrale UC munie d'un microprocesseur P qui coopère avec une carte-mère CM. Cette carte-mère est reliée à divers équipements, tels qu'une interface de communication COM (de type Modem ou autre), une mémoire morte ROM et une mémoire de travail RAM (mémoire vive). La carte-mère CM est reliée en outre à une interface graphique IG, laquelle pilote l' affichage de données sur un écran ECR que comporte le dispositif. Il est prévu en outre des moyens de saisie, tels qu'un clavier CLA et/ou un organe de saisie dit "souris" SOU, reliés à l'unité centrale UC et permettant à un utilisateur une interactivité avec le dispositif. La mémoire ROM, ou encore la mémoire RAM stocke les premières données précitées, représentatives des paramètres du réseau (topologie, propriétés des routeurs, etc). Dans l'exemple, la mémoire RAM reçoit les secondes données précitées, représenta- tives d'événements dans le réseau (flux transverses, conges- tions, pertes, etc). Dans le cadre d'une aide à la surveil- lance du réseau ces secondes données peuvent être reçues par l'interface de communication COM. Dans le cadre d'une simulation, l'acquisition de ces secondes données peut être effectuée par un calcul basé sur un modèle de simulation, comme on le verra plus loin. La mémoire RAM au moins peut être adressable en fonction de rangées et de colonnes de matrices et permettre ainsi un stockage de données sous forme matricielle. La mémoire ROM comprend un module de modélisation MOD qui, en coopération avec le microprocesseur P, permet de construire la première matrice précitée et une matrice courante, respectivement en fonction des premières données et des secondes données, selon un modèle choisi que l'on verra plus loin. La mémoire ROM comprend un module CAL qui, en coopération avec le microprocesseur P, permet d'effectuer sur au moins deux matrices de structure dynamiquement variable, une opération formant produit selon l'algèbre MAX-PLUS. Dans l'exemple, on entend par "matrices de structure dynamiquement variable", des matrices dont les coefficients au moins sont dynamiquement variables. Avantageusement, les modèles qui seront décrits ci- après permettent de ramener les matrices construites (et plus particulièrement les matrices courantes) à des matrices dont seuls les coefficients sont dynamiquement variables. Le module de modélisation MOD comprend alors : - un sous-module de modélisation statique ST, pour construire la première matrice en fonction des premières données, et - un sous-module de modélisation dynamique DYN, pour cons- truire, en fonction des secondes données, au moins une matrice courante dont les coefficients sont dynamiquement variables. La mémoire ROM comprend en outre un module PIL qui, en coopération avec le microprocesseur P, permet d'appliquer répétitivement la première matrice précitée (comprenant les premières données) et la matrice courante (comprenant les secondes données) au module de calcul CAL. La matrice produit obtenue est stockée en mémoire et devient une nouvelle première matrice. Elle peut être ensuite multipliée (dans l'algèbre MAX-PLUS) à une autre matrice courante, comportant de nouvelles secondes données représentatives de nouveaux événements dans le réseau. De diverses approches ont été proposées pour caractériser les propriétés clés du mécanisme de contrôle de flux à fenêtre du type TCP à partir notamment de considérations heuristiques, de simulations, d'approximations fluides ou encore d'analyses markoviennes [10, 11, 1, 12, 13, 141. Tous les modèles analytiques sont basés sur la réduction du réseau à un seul noeud représentant le goulot d'étranglement [9]. Par ailleurs il a été récemment démontré que le contrôle de flux à fenêtre d'un réseau multidimentionnel admet une représentation max-plus linéaire lorsque la taille de la fenêtre est constante [5]. Ici , nous nous intéressons aux modèles qui combinent le mécanisme de contrôle adaptatif de TCP et un réseau multidimensionnel constitué de plusieurs routeurs en série. Nous montrons que la dynamique d'un tel réseau contrôlé peut être décrite au niveau paquets via des itérations de produits matriciels dans l'algèbre max-plus. Nous considérons à la fois le cas où les temps de transmission des paquets sont déterministes et les divers modèles stochastiques qui ont été utilisés dans la littérature, en particulier, les cas où il y a des pertes aléatoires en plus des pertes dues au dépassement de la capacité des buffers, et le cas où les temps de transmission des paquets sont aléatoirement perturbés par les autres trafics. Tous les aspects clés du protocole peuvent être représentés : pertes de congestion, time-outs, pertes aléatoires, délais de propagation, délais dûs aux attentes ou encore au mécanisme de contrôle de flux etc. Nous montrons comment cette approche permet d'obtenir des formules explicites pour le débit maximum alloué, lorsque les perturbations sont déterministes ou aléatoires. Ces formules sont asymptotiquement compatibles avec les formules classiques lorsque la taille maximale de la fenêtre tend vers l'infini. De plus, ce cadre permet d'analyser les fluctuations instantanées et aléatoires du débit, ce qui peut être utile pour estimer la qualité de service offerte à une connexion. Ce modèle générique est aussi particulièrement bien adapté pour des simulations efficaces de la dynamique d'une session TCP opérant son contrôle de bout-en-bout sur un réseau de grande dimension. Dans la Section 2, nous donnons une représentation maxplus générale de notre modèle de base. Dans les Sections 3 et 4, nous considérons les services déterministes avec une évolution déterministe puis markovienne de la taille de la fenêtre. Nous montrons que pour ces modèles déterministes, le débit ne dépend que du temps aller-retour "round trip time (RTTjP et de la bande passante "bandvridth". Des variantes et extensions du modèle déterministe sont considérées dans la Section 5. Nous considérons notamment 1) le cas des pertes aléatoires en plus des pertes de congestion ; une description détaillée de la taille du buffer et des pertes de congestion ; le cas de services aléatoires ; le cas avec time-outs. Dans la Section 6, nous donnons une liste des diverses extensions qui peuvent être traitées dans ce cadre mathématique. et qui sont compatibles avec notre représentation. Elles conduisent toutes à des formules analytiques ou à de nouveaux concepts de simulations basés sur le calcul du produit d'un grand nombre de matrices ; particulier, nous montrons que le coût de la simulation par cette approche de la transmission de n paquets sur K routeurs est en 2n(KWmax)2,où Wmax est la taille maximum de la fenêtre. 2 Représentation max-plus 2. 1 Algèbre max-plus L' "algèbre" max-plus scalaire est un semi-anneau sur la droite réelle où l'on remplace l'addition par max (notée #) et la multiplication par plus (notée ). La loi # est distributive par rapport à # ce qui permet d'étendre les concepts classiques de l'algèbre linéaire à ce cadre, et en particulier la théorie des matrices. Ce semi-anneau est noté (Rmax,#,@), où Nmax = R U {#00} est la droite réelle complétée avec #oo, l'élément neutre de . Par la suite, nous notons (I , ffi, (À ED B)v 4,,t ED Bkj = max(A;k, Bkj), ('4 B)ij = # AiA: 1 2. 2 Modèle de base Nous considérons un réseau PAPS (Premier Arrivé Premier Servi) avec K files en tandem tel que le n-ième client arrivant à la station i reçoit un service #i(n). Dans le contexte de TCP, ce réseau modélise une seule source envoyant des paquets vers une seule destination, à travers un chemin composé de K routeurs. La variable #i(n) est le retard aléatoire causé par le trafic transversal (les autres utilisateurs) au niveau du routeur i sur le n-ième paquet. Ce retard n'inclut pas les délais d'attente, mais seulement le ralentissement de la vitesse du serveur dû à la présence du trafic transverse. Le délai de propagation du n-ième paquet entre les routeurs i à j sera noté di,j (n). Le débit du flux d'entrée est contrôlé par une fenêtre dynamique, dont la taille est égale au nombre total de paquets envoyés par la source à un instant donné et n'ayant pas atteint la destination (ou plus précisément les paquets qui n'ont pas encore été "acquittés"). 2. 3 Evolution de la taille de la fenêtre La taille de la fenêtre a une évolution générale définie par les équations suivantes: W(0) = 1, W(n + 1) = f(W(n), ACK(n , où ACK(n) est le signal de contrôle de flux/congestion donnant l'information sur l'état du système au temps n et où f est une fonction que nous préciserons ultérieurement. Par exemple, ACK(n) = 1 si aucune congestion (ou aucune perte de paquet) n'est observée par le n-ième paquet, sinon ACK(n) = 0. Par la suite, nous supposerons, dans certains cas, que f dépend aussi de W,, qui est, pour les modèles TCP, le seuil qui sépare la phase de croissance exponentielle - slow-start phase - de la phase de croissance linéaire -congeition-avoidance phase. La relation de récurrence (1) définit la taille de référence de la fenêtre. La taille effective est alors définie comme la partie entière de la taille de référence: wn = (int) W(n) = [w(n)]. (1) Nous supposons que la taille maximale de la fenêtre est finie et nous notons: tory = m87C L7, = (int) Wmax = ## . n#N Nous supposons également que l'évolution de la taille de la fenêtre peut être décomposée en deux phases qui dépendent de ACK(n) de la façon suivante: phase de croissance: {n : ACK(n) = 1}, 0 Voici le modèle que nous allons étudier: Contrôle de flux TCP ACK(n) E {0,1} est le signal d'acquittement du n-ième paquet qui détecte un état de congestion ou une perte de paquet. Les exemples usuels de la politique idéale de l'évolution de la fenêtre sont : 1. TCP Tahoe : 0 W.(n), wu (phase linéaire) et W. (n + 1) : = W,{n); f (W (n), 0, W, {n)) - 1; et W,(n + 1) : LaW(n)j. 2. TCP Reno : 0 Dans les exemples suivants, nous nous limitons à la valeur a = 1/2. 2. 4 La dynamique du vecteur des dateurs A partir de maintenant, nous supposons que la file à l'entrée est saturée (même si le cas nonsaturé peut être intégré facilement cf. 6). Alors le réseau se comporte comme un réseau fermé et son débit donne le taux maximum auquel la source peut envoyer des paquets tout en gardant un buffer d'entrée stable [5]. Soit x.(n) la date à laquelle le n-ième paquet arrivant au routeur i commence son service sur ce routeur. On se réfère à la figure 1 montrant un nombre K de files en tandem, avec contrôle de flux. Soit y; (re) = xi (n) + oi (n) la date à laquelle le n-ième client quitte le routeur i. Si on pose #0(n) = 0 pour tout n > 0, alors , pose yo(n) = K(n- n-1) dx,o(n - n-1) Yi(n) [Yi-l(n)t8Id.-l,.(n)E9y.(n-l)]t8IO".(n), ,K. Soient yen) = (y(n), Y2(n), ..,YK(n)) E et Zen) = (y(n),Y(n - 1),..,'(n -'tv' + 1))t E -a::...1. Ce vecteur sera appelé vecteur des dateurs par la suite. Soient Mi, i E {1, ..,w*}, des matrices données de RmaxK,K, Soit E la matrice de RmaxK,K dont tous les éléments sont égaux à -oo. Ci-dessous, on note (M1M2 - # \MW*) la matrice de 1 xw définie par blocs de taille K x K: tous les blocs sont égaux à la matrice E de RmaxK,K, sauf pour la première ligne de blocs qui est égale à Ml , MZ, ... , M,,,.. Si initialement le système est vide, l'évolution du vecteur des dateurs Z (n) donnée par la récurrence max-plus linéaire suivants: Z(0) ~ (0, .., 0)t, Z(n) = AW..-t (n) Zen - 1), \ln 1, (2) où Al(n) ~ (M(n) 0 Mi(n)..E) D, A2(n) = (M(n)IM2(n)IEI..IC) D, ..., Aw.(n) ~ (M(n)\ |..||Af,.(n)) ID D. Dans ces formules, M(n) et Ml(n) sont données par: i i-1 (M(n iJ = L C1k{n) + L dk,k+l (n), ifi>j, -oo, si i On peut définir l'exposant de Lyapunov de la station k par limn#+# yk(n)/n. A cause de la monotonie, il est clair que cette limite est indépendante de k. Cette propriété est plus généralement vraie sous l'hypothèse d'irréductibilité définie dans [8]. # Au niveau de la représentation donnée ici, nous ne prenons pas en compte le contrôle d'erreur et nous ne faisons aucune différence entre les paquets originaux et les paquets retransmis. En particulier, nous ne faisons pas distinction entre le send rate, throughput ou goodput [13]. # L'équation (2) est la base du schéma de simulation algébrique auquel nous avons fait allusion dans l'introduction. Comme les matrices Awn-1(n) sont de dimension Kw*, et comme seuls les produits matrice-vecteur sont nécessaires, on peut simuler la transmission contrôlée de n paquets à travers le réseau en 2n(Kw*)2 opérations sur un seul processeur. Le* Figure 2A représente un exemple d'évolution pas-à-pas explicite des dateurs et de la taille de la fenêtre: nous avons pris K = 5 avec {o~\ , # # # Q5) ~ (1,1,2,1,1). Donc, w' = 4. Nous avons considéré une évolution déterministe de la taille de la fenêtre (TCP Tahoe sans la phase de slow start, cf. 3): (Wl, W2, WS, ..,) = (1,2,2,3,3,3,4,1,2,2, ..). Par exemple, le paquet # 5 correspond à la taille de la fenêtre 3, ce qu'on a noté 5 (3) la figure: ceci veut dire que juste après la transmission du paquet # 5, il y a exactement 3 paquets dans le réseau qui n'ont pas encore été "acquittés". Comme le montre cette figure, le retards imposés aux paquets (voir par exemple le paquet # 6) dans les routeurs internes peuvent être assez complexes. Pour la simplicité de la présentation, nous allons d'abord considérer le cas où tous les délais de propagation di,j(n) valent 0, et ensuite nous montrerons comment les formules doivent être modifiées pour couvrir le cas de délais de propagation non nuls. 3 Cas déterministe Supposons que les services sont tous déterministes et constants, i.e. #i(n) = #i, pour tout n, où #i est une constante. Soit or* = maxlÍK (1, et S = 1 a,. Nous supposons que la congestion est détectée par la politique suivante: Cette détection peut être interprétée de la façon suivante: 1/#* est le taux de service du goulot d'étranglement du réseau et S est le temps d'aller-retour RTT. On détecte donc une congestion quand le taux d'émission moyen wn/RTT atteint le taux du goulot d'étranglement. Pour le modèle à fenêtre statique, w* est la taille de la fenêtre optimale donnée par la formule bien connue: bande passante x délai [9]. Sous les hypothèses ci-dessus, 1. w* est donné par w* = min {n : #*n > S} [s/#*] + 1 2. wn devient périodique avec période T qui peut être décomposée en: T = tl + t2 +### + tw*, où t, est le nombre d'occurrences de wn = i pendant une période. 3. 1 Tahoe et Teno 3.1.1 TCP Tahoe Nous considérons d'abord le modèle TCP Tahoe sans la phase exponentielle. Sous la politique (Ci), nous avons: Pour tout n est donné par Si tous les délais sont déterministes, A >. = (w. - 1)8 + a.. Par conséquent, le débit dépend des délais (ci,... a K) seulement à travers S et a. et est donné par: Lorsque w* # +00, le débit asymptotique est: Sur la Figure 2 (on suppose que w* = 6), nous avons représenté l'évolution pas-à-pas de la date d'entrée yo(n) et de la date de départ yx(n) (du paquet # n); les propriétés algébriques ci-dessus ont l'interprétation suivante : la détection de congestion, les paquets envoyés se comportent comme s'il n'y avait aucune interaction entre eux, sauf pour le couple de paquets envoyés en même temps lorsque la taille de la fenêtre croît d'une unité ; ces derniers paquets, le second paquet quitte toujours la station K avec un délai de #* par rapport au premier. On peut donc lire directement la valeur propre 5S + or* sur l'évolution graphique. Le débit de saturation pour le TCP Tahoe avec la phase exponentielle est donné par: Lorsque w* # +00, le débit asymptotique est: La figure 3 représente l'allure du rapport n/y4(n) (traits pleins) et de wn (traits pointillés) dans les conditions suivantes: TCP Tahoe sans la Phase exponentielle: quatre files en tandem avec #1 = 3.2, 0'2 = 4.61, u3 = 2.7, Q4 = 4.61. w' = 4, Wn E {1,2,3,4}. Le débit de saturntion est égal à 0. 140084 3. 1.2 TCP Reno Une évolution périodique déterministe de TCP Reno a été considérée dans [11] pour obtenir une valeur heuristique du débit. La représentation max-plus ci-dessus conduit à une nouvelle formule qui raffine celle de 11 il. Si tous les délais sont déterministes, le débit de saturation ne dépend que de S et de #*: Lorsque w*# +00, le débit asymptotique est: n se- rr:ér Q.: '(1 La formule bien connue du débit d'une connexion TCP en fonction du taux de perte pperte et du temps aller-retour RTT est de la forme [11]: où co est une constante réelle. Pour notre modèle déterministe, nous avons: RTT = S et Lorsque w* -3 oc, nous avons ,fp,,-,t, - -',w2: (case 1) and "Pperte '" J-8g -.!.- (cases 2 and 3). D'où: cas 1: Co = 0.71 ; cas 2 et 3: Co = fi 1.22. Par conséquent, pour les grandes valeurs de w* (ou les petites valeurs de pperte), la formule asymptotique du Corollaire 2 coîncide bien avec celle de [11]. Tous les résultats ci-dessus sont vrais pour des délais de propagation d;,j constants à condition de remplacer la valeur de S par S = Pour les modèles déterministes, le débit que nous obtenons peut être comparé à celui simulé par le simulateur NS en choisissant une taille de paquet arbitraire et en prenant une vitesse de service du routeur i correspondant à #i. Les écarts entre les débits obtenus par la simulation NS et par nos formules ne peuvent provenir ici que de différences sur les mécanismes de détection de la perte/congestion, ou du fait que nous avons pris ici la partie entière de W (n) dans f (comme dans [9]) tandis que NS utilise W(n). En effet, pour tous les modèles déterministes avec la même évolution périodique de W(n), les évolutions sont exactement les mêmes. Voici un exemple : sur NS avec source ftp. K = 10, taille d'un paquet 1250 (40 pour les accusés de réception), taille des buffers 2 ; tous les délais di,j sont égaux à O.lms sauf dK0 qui est égal à 1ms, vitesse de service: (10, 5, 4, 2, 5, 4, 5, 5, 4, 5, 5)Mb/s pour les liens 0 -1, .., 9 -10,10 - 0. A t = 100s, NS donne 152. 27 paquets/s. Pour cet exemple, S = 25. 5ms et #* = 5ms ; nous avons par (5) : 134 paquets/s. Cependant, nous remarquons que w* est en fait égal à 7 dans la simulation NS (un RTT est nécessaire pour détecter le triple-acks sur ce modèle NS), au lieu de 6 dans notre modèle ; évaluant (5) avec w* = 7, on a 152. 55 packets/s. Enfin, pour tous les modèles déterministes, la suite {wn} est déterministe, et on peut donc calculer l'évolution du débit d'émission de manière exacte à partir des produits de matrices maxplus correspondants. 4 Modèle de Markov avec services déterministes Nous supposons toujours que les services sont tous déterministes. Maintenant nous supposons que l'évolution de W (n) peut être décrite indépendamment des autres éléments du réseau par la matrice de transition markovienne P telle que: Vn6N. si W(n) W(n), g(W(n)) Si tous les services sont déterministes, pour toute évolution markovienne de W (n) avec la structure ci-dessus, le débit ne dépend des délais qu'à travers S et or*. De plus, l'exposant de Lyapunov (l'inverse du débit) est de la forme: [gamma] = #[gamma]## [gamma]#([gamma]), où # est la probabilité stationnaire de {yK(n + 1) - YK(n)}nEN et où A = {S - (k -l)a*,k = 1,.., w* - 1} U iu*l. 13 4. 1 Modèle avec deux transitions possibles Ici nous supposons que Wmax = w*, où w* est défini par la politique (CI), et que Ce modèle est à comparer à celui de [12], où une probabilité de perte globale est utilisée pour capturer en même temps les time-outs (TO) dûs aux pertes des paquets et les triple-duplications des ACKs (TD) dûs à la congestion. Dans notre modèle, ces deux mécanismes sont au contraire décrits séparément; les pertes de paquets qui génèrent les TO constituent une suite i.i.d. (indépendante et identiquement distribuée), independante de tous les autres éléments du réseau et sont capturées par le paramètre p-; les pertes dues à la congestion sont capturées par le paramètre w*. Pour TCP Tahoe avec une évolution markovienne comme ci-dessus (i.e. sous (Ml) et g(W(n)) = 1), le débit est donné par: Lorsque w' # oo, si p- # 2/(w*)2,le débit asymptotique est: On peut facilement vérifier que p- = 0 correspond au cas de TCP Tahoe déterministe (sans la phase exponentielle). Par ailleurs, si p- # 0 et p- #2/(w*)2 (ce qui veut dire que les pertes de paquets sont liées à la congestion d'une façon qui est en moyenne similaire à celles du cas du TCP Tahoe déterministe), nous avons: co # 0. 60. Donc, la dégradation des performances peut être significative (approximativement 15%) lorsque qu'on passe du modèle déterministe (où co # 0. 71) au modèle markovien. L'impact des TO sur les performances a déjà été remarqué par exemple dans [13, 10]. L'influence prépondérante de p- sur la probabilité de perte globale peut être quantifiée à partir de notre modèle analytique comme suit : probabilité de perte globale de ce modèle est: où h{vj*) est la perte due à la congestion et p~(1 - p(w')) est la perte due aux TO. Pour (1) fixé, la Figure 5 montre la décroissance du débit en p- (nous avons aussi fixé le taux du goulot d'étranglement : #* = 1 et posé S = w* - 1). On voit là dans quelle mesure l'impact des pertes par TO est prépondérant par rapport aux pertes dues à la congestion. Des résultats similaires peuvent être obtenus pour des modèles du type TCP Reno et aussi pour des modèles basés sur une évolution markovienne de la taille de la fenêtre avec trois transitions (ou plus). Voici quelques exemples. Soit K = 4 avec ci = 3.2, a2 = 4.61, ors = 2.7, U4 = 4.61. w* = 4, w" E {1,2,3,4}. La Figure 6 montre l'évolution de 4"n et w" pour TCP Reno markovien avec (p+,po,P-) = (0.8,0.1,0.1). La Figure 7 montre les débits de TCP Reno et TCP Tahoe dans les cas déterministe et markovien ((p+,po,P~) # (0.8,0.1,0.1). Dans le cas déterministe, les débits sont égaux à 0.140084 (Tahoe) et 0.172166 (Reno). Les débits du modèle déterministe semblent supérieurs à ceux des modèles markomens. CRI 5 Buffers finis et services aléatoires Buffers finis Nous considérons d'abord des files en tandem avec des services déterministes et avec une taille du buffer fixe b pour toutes les stations. Dans ce cas, d'après la formule de Little, il est naturel de définir la détection de la congestion par: Ce modèle conduit encore à un comportement déterministe et périodique de Wn et à une analyse similaire à celle de 3.1. On remarque que le modèle de 3 est un cas particulier de (C2) avec la taille du buffer b = 1. Nous considérons maintenant des files en tandem avec des services aléatoires et avec une taille du buffer fixe b pour toutes les stations. Soient u*(n) = max19K ui(n) et S (n) = K1 of(n). Nous supposons que les supports des temps de service sont bornés et nous posons: S = sup{S(n)} et (1. = suplor*(n)l. Soit r(n) = â n . Par la suite, nous prenons w* = 6 + sup{r(n)}. Par analogie avec ce que nous avons proposé ci-dessus (C3) : ACK(n) = 1, si wn ~ b + r(n). Dans ce modèle, W(n) a une évolution déterministe lorsque W(n) E { 1, .., b} et une évolution aléatoire lorsque W (n) E {b + 1, .., b + K}. L'interprétation de ce modèle est la suivante: pour une évolution de la taille de la fenêtre du type TCP, nous retenons comme date de début de la congestion le moment où W (n) = S " Ensuite, une perte de paquet a lieu lorsque le nombre de paquets dans le système dépasse la taille du buffer, ce que nous approximons par l'événement W (n) = b + r(n). Voici une autre interprétation: b est une taille arbitraire (par exemple, la moitié de la taille réelle du buffer) telle que si le nombre des paquets en attente dépasse ce seuil, la congestion est détectée. où uy , w2, ... , est la suite caractéristique de Tahoe (1,2,2,3,3,3, ... , w. - 1, w. - 1, w.). Sous la politique (C3), si les vecteurs de temps de service (a1(n),... ,00K(n))} sont i.z.d. en n, le débit de TCP Tahoe sans phase expo- nentielle est donné par: où 0 est le vecteur (0, 0)t E lfxl et Les mécanismes de time-out peuvent être pris en compte, par exemple, par la condi- tion S(n) > TO ou par la condition yx(n) - yx(ra - w,,-,) > TO. Par exemple, dès que Z(n) est markovien et que le support .6. des v.a. f yK (n) - yK (n - mn~1 ), n, E N} est fini, nous avons la même structure régénérative comme ci-dessus. La loi d'une cycle de régénération Tl peut être explicitement calculée par la récurrence suivante: Cette formule de récurrence est valide pour n Une formule similaire peut être obtenue pour TCP Tahoe avec phase exponentielle ou TCP Reno, et aussi pour diverses extensions du modèle ci-dessus en incluant des pertes de paquets comme dans 4. Ici, nous considérons quatre files en tandem avec une taille du buffer b = 50 et avec une distribution multinomiale indépendante pour les services. La Figure 8 montre l'évolution de n/y4(n) et W(n) pour le modèle TCP Tahoe avec la phase exponentielle et avec b = 50. Les #i(n) sont i.i.d. et mutuellement indépendantes, à valeurs dans {1,5,10} avec des probabilités respectives de 0.3,0.4,0.4. La Figure 9 montre une comparaison des débits de deux TCP Tahoe avec la phase exponentielle: b = 10 et K = 10. Dans le premier cas, #3(n) est constant égal à 1 et {#j(n),j # 3} une suite i.i.d., à valeurs dans il, 10, 20} avec les probabilités respectives 0.3,0.2,0.5. Dans le second cas, u8(n) est constant égal à 1 et les autres services comme ci-dessus. Le calcul des débits au moyen de la formule du Théorème 5, montre que la permutation des caractéristiques de deux routeurs peut changer la valeur du débit. Donc, déjà dans ce cas, nous ne pouvons pas réduire le modèle à un ,modèle avec un seul serveur, celui du goulot d'étranglement. Noter aussi que la connaissance des valeurs moyennes est insuffisante pour prédire la valeur du débit moyen. # 6 Exploitation de l'approche Dans les Sections antérieures, nous nous sommes limités aux valeurs moyennes du débit de saturation. En fait, à partir de notre analyse, on peut obtenir des résultats plus complets, soit analytiquement, soit par simulation algébrique ; concerne par exemple : 1. la loi des débits instantanés lininoo P (YK(n + 1) - yx(n))-1 t) = P(YK(F) > t), mais une formule plus précise peut être donnée en prenant en compte des retransimissions de paquets perdus. Dans ce cadre, nous pouvons aussi traiter les problèmes suivants au moins dans le cadre de la simulation algébrique: 4 les modèles ouverts (où la sourée n'est pas saturée : au lieu de FTP) où le processus des arrivées est décrit par ses caractéristiques statistiques; 5 les connexions multiples et les interactions entre plusieurs utilisateurs; 6 les connexions multipoints à travers un réseau avec une structure arborescente au lieu d'une structure linéaire de routeurs en série (voir [7] pour le cas à fenêtre constante). Ces dernières questions feront l'object de recherches futures. 7 Simulation Dans cette Section, nous donnons quelques précisions sur la structure d'un simulateur algébrique qui permettrait de prévoir les performances de TCP ou d'autres mécanismes de contrôle de flux issus de TCP dans des réseaux IP existant déjà ou en cours de conception. Le principe de base de ce simulateur est l'analyse d'une connexion contrôlée. La simulation de cette connexion commence par l'acquisition des données concernant le réseau: - nombre des routeurs traversés par la connexion, caractéristiques de chacun des routeurs: capacité en Mb/s, taille des buffers bi, scheduling etc.; - délais de propagation sur chacun des liens (di,i+1); - caractéristiques statistiques des flux transverses Fi lorsque ces dernières sont connues (comme par exemple dans le cas des routeurs du backbone Internet); - caractéristiques statistiques du trafic offert par la connexion contrôlée; - caractéristiques statistiques des trafics offerts par les autres flux en compétition avec la connexion contrôlée sur le routeur d'accès. Ces caractéristiques sont résumées sur la Figure 10 On construit alors en fonction de ces données et des caractéristiques de la version de TCP choisie (ou de l'algorithme de contrôle de flux retenu) un simulateur algébrique qui consiste en la construction pas à pas des matrices Awn du Théorème 1 et le calcul par récurrence des dateurs du sytème comme indiqué dans ce théorème. On déduit de l'analyse mathématique qui a été présentée dans cet article que l'on peut cal- culer à partir de la simulation les principales caractéristiques de l'effet de TCP sur la connexion contrôlée (débit moyen, débit moyenné sur une période de temps donnée, fluctuations du débit instantané etc.). Pour compléter cette description du simulateur, voici quelques indications sur la manière de prendre en compte certains mécanismes fins: Pour un routeur i, soit bi la taille de son buffer par flux; la détection du dépassement de cette capacité se fait par la condition suivante sur les dateurs: xi (n+bi)#yi(n). Comme déjà indiqué, les mécanismes précis des time-outs peuvent 'Être pris en compte par la condition yK(n)- yK(n- wn-1) > TO. Cette formulation vaut aussi lorsqu'il y a évolution de TO lui-même (dans certaines versions, la variable TO est mise à jour à chaque timeout : TO est multiplié par 2 en cas de timeout, jusqu'à ce qu'on atteigne la 64 fois la valeur d'initialisation). Nous avons clairement montré que l'on pouvait représenter les variations et adap- tations de la fenêtre dans le cadre max-plus linéaire. Les variations et adaptations de la variable TO sont de nature semblable et sont donc aussi représentables dans ce cadre. Dans certaines versions de TCP, il est recommandé de procéder au regroupement des accusés de réception. On attend par exemple que trois paquets soient arrivés à la destination pour renvoyer un seul accusé pour les trois paquets. Ceci peut aussi être représenté au moyen d'un chaînage de deux mécanismes de fenêtre variable, et rentre assez naturellement dans le cadre d'une description max-plus linéaire. Un mécanisme de scheduling dans un routeur (par exemple FIFO, ou Weighed Fair Queueing avec éventuellement des différenciations de services) peut être représenté dans cette simulation par la prise en compte de l'influence de ce mécanisme sur les durées de service #i(n); une fois calculées en fonction du trafic transverse, du trafic de la connexion contrôlée et du mécanisme de scheduling, ces durées aléatoires peuvent être injectées dans le simulateur. Pour des infrastructures de réseaux d'accès en cours de conception, on ne peut évidemment pas mesurer les caractéristiques statistiques des trafics transverses sur les routeurs d'accès. Supposons qu'on dispose d'une description statistique du trafic offert prévu sur les routeurs d'accès de ce réseau (par exemple N trafics HTTP). On peut alors utiliser le simulateur pour calculer le trafic transverse inconnu comme un point fixe du système : soitt la loi du trafic de la connexion contrôlée observé au niveau du routeur d'accès; t doit être tel que si on prend comme trafic offert sur la connexion contrôlée celui d'une session HTTP, - le trafic de la connexion contrôlée observé sur le routeur d'accès a pour loi t; - le trafic transverse observé en ce point est la somme de N trafics de loi t. Ainsi, dans les réseaux saturés ou non, le mécanisme de contrôle en feedback adaptatif de TCP est un feedback linéaire dans l'algèbre max-plus. Ceci conduit à une repré- sentation simple de l'effet de ce protocole sur n'importe quel réseau qui admet lui même une représentation max-plus en l'absence du contrôle, comme c'est le cas par exemple pour les files en tandem ou les réseaux fork-join que l'on trouve dans les arbres multipoints. Nous en avons déduit des formules explicites pour de divers modèles avec services déterministes qui raffinent les résultats bien connus de la littérature. Ces formules confirment que dans ce cas, le débit ne dépend que du RTT et du taux du goulot d'étranglement. De nouvelles formules ont aussi été obtenues lorsque les services sont aléatoires. Les aléas représentent ici l'effet du reste du trafic sur la connexion contrôlée. Dans ce cas, nous avons montré que le débit ne peut pas être obtenu seulement à partir de considérations en moyenne, et que l'ordre et le comportement statistique fin des routeurs ne peuvent être ignorés. L'ensemble des modèles entrant dans ce cadre est très riche : peut en effet choisir des services déterministes ou aléatoires, un contrôle de flux basé sur la perte ou congestion; les pertes peuvent provenir de la congestion ou des time-outs, ou être aléatoire, ou encore être une combinaison de ces trois possibilités ; version de Reno ou Tahoe peut être choisie, avec ou sans la phase exponentielle etc. Nous avons montré comment notre approche pouvait être utilisée pour analyser quelques une de ces combinaisons ; est important de souligner que toutes les combinai- sons peuvent être en principe analysées dans ce cadre. Plus généralement, cette approche fournit un cadre générique pour la simulation des protocoles du type TCP sur des réseaux qui peuvent être grands. La simulation est basée sur un algorithme simple qui exploite la linéarité et qui a une complexité maîtrisée. Bien entendu, l'invention n'est pas limitée à la forme de réalisation décrite précédemment à titre d' exemple, elle s'é- tend à d'autres variantes. ANNEXE Références [1] Altman, E., Bolot, J., Nain, P., Elouadghiri, D., Erramdani, M., Brown, P. and Collange, D. (1997) Performance Modeling of TCP/IP in a Wide-Area Network. INRIA Report, March, N 3142, INRIA. [2] Baccelli, F., Cohen, G. Olsder, G. J. and Quadrat, J. P. (1992) Synchronization and Linearity, Wiley. [3] Baccelli, F. and Hong, D. (1998) Analytic Expansions of (max, +) Lyapunov Exponents. INRIA Report, May, N 3427, INRIA Sophia-Antipolis, To appear in Annals of Appl. Prob. [4] Baccelli,F., Gaubert, S., Hong, D. (1999) Représentation and Expansion of (Max,Plus)- Lyapunov Exponents. Proc. of 37th Annual Allerton Conf. on Communscation, Control and Computing, September, Allerton Park. [5] Baccelli,F., Bonald, T. (1999) Window flow control in FIFO networks with cross traffic. Queueing Systems, 32,195-231. [6] Chang, C.S. (1999) Performance guarantees in Communication Networks. Springer Verlag. [7] Chaintreau, A, Baccelli, F., Diot, C. (2000) Impact of Network Delay Variation on Multicast Sessions Performance With TCP-like Congestion Control, soumission SIGCOMM 2000. [8] Hong, D. (1999) Equality of the maximum and minimum throughputs of irreducible stochastic max-plus linear systems. In préparation. [9] Keshav, S. (1997) An engineering approach to computer networking. Addison-Wesley Profes- sional Computing Séries. [10] Lakshman, T. V., Madhow, U. (1997) The performance of TCP/IP for networks with high bandwidth-delay products and random loss. IEEE/ACM Trans. Networking, June. [Il) Mathis, M. Semske, J. Mahdavi, J and Ott, T. (1997) The Macroscopic Behavior of the TCP Congestion Avoidance Algorithm. Computer Communication Review, 27 (3), July. [12] Padhye, J., Firiou, V., Towsley, D. and Kurose, J. (1998) Modeling TCP throughput : simple model and its empirical validation. Proc. of ACM SIGCOMM. [13] Padhye, J., Firiou, V., Towsley, D. (1999) A Stochastic Model of TCP Reno Congestion Avoidance and Control. Technical Report, 99-02, CMPSCI, Univ. of Massachusetts, Amherst. [14] http: //www.psc.edu/networking/tcp~friendly.html Revendications 1. Dispositif d'aide à la surveillance et/ou de simulation d'un système complexe, notamment d'un réseau de communication, caractérisé en ce qu'il comprend: - une mémoire (ROM,RAM) pour stocker de premières données représentatives de paramètres du réseau, ainsi que pour recevoir au moins de secondes données représentatives d'événements dans le réseau, une portion de ladite mémoire étant réservée au stockage de données sous forme matricielle, - un module de calcul (CAL), apte à effectuer sur au moins deux matrices de structure dynamiquement variable, une opération formant produit selon l'algèbre dite MAX-PLUS, - un module de modélisation (MOD) pour construire au moins une première matrice (Z(n)) et une matrice courante (A(n)) respectivement en fonction des premières données et des secondes données, selon un modèle choisi, et - un module pilote (PIL) pour appliquer répétitivement la première matrice et la matrice courante au module de calcul, la matrice produit obtenue (Z(n+1)) devenant une nouvelle première matrice. 2. Dispositif selon la revendication 1, caractérisé en ce que les premières données comprennent des données relatives à la topologie du réseau. 3. Dispositif selon la revendication 2, caractérisé en ce que les premières données comprennent le nombre (K) de routeurs du réseau traversés par une connexion à surveiller ou à simuler et des propriétés desdits routeurs. 4. Dispositif selon la revendication 3, caractérisé en ce que les premières données comprennent des tailles mémoire des routeurs ( bK ) . 5. Dispositif selon l'une des revendications 2 à 4, caractérisé en ce que les premières données comprennent des propriétés statistiques des trafics offerts dans le réseau. 6. Dispositif selon l'une des revendications précédentes, caractérisé en ce que le module de modélisation comprend : - un sous-module de modélisation statique (ST), pour construire la première matrice en fonction des premières données, et - un sous-module de modélisation dynamique (DYN), pour construire au moins une matrice courante en fonction des secondes données. 7. Dispositif selon l'une des revendications précédentes, caractérisé en ce que les secondes données comprennent des informations relatives à des pertes (p~,p+) dans le réseau. 8. Dispositif selon l'une des revendications précédentes, caractérisé en ce que les secondes données comprennent des informations relatives à des flux transverses dans le réseau, par rapport à une connexion contrôlée à surveiller ou à simuler. 9. Dispositif selon l'une des revendications 7 et 8, caractérisé en ce que les secondes données comprennent des informations relatives à des congestions dans le réseau. 10. Dispositif selon l'une des revendications 7 à 9, caractérisé en ce que les secondes données comprennent des informations relatives à des dépassements de délais (TO) dans le réseau. 11. Dispositif selon l'une des revendications précédentes, caractérisé en ce que ladite matrice produit est représentative d'un débit dans le réseau associé à la connexion à surveiller ou à simuler. 12. Dispositif selon la revendication 11, caractérisé en ce que ladite matrice produit est représentative d'un débit moyen dans le réseau. 13. Dispositif selon la revendication 11, caractérisé en ce que ladite matrice produit est représentative de fluctuations d'un débit instantané. 14. Dispositif selon l'une des revendications précédentes, caractérisé en ce que le module de modélisation est agencé pour construire une pluralité de matrices, en nombre correspondant sensiblement au nombre de paquets dans le réseau. 15. Dispositif selon l'une des revendications précédentes, caractérisé en ce que ledit modèle choisi comprend la considération de la taille variable d'une fenêtre (wn) utilisée pour contrôler le nombre de paquets dans le réseau. 16. Dispositif selon la revendication 15, caractérisé en ce que ledit modèle choisi comprend la considération d'un réseau comprenant des routeurs à discipline de type "premier arrivé premier servi". 17. Dispositif selon l'une des revendications 15 et 16, caractérisé en ce que ledit modèle choisi comprend la considération d'un réseau comprenant des routeurs à discipline de type WFQ. 18. Dispositif selon l'une des revendications précédentes, caractérisé en ce que le modèle choisi comprend la considération d'un réseau contrôlé par un protocole de type TCP. 19. Dispositif selon l'une des revendications 15 à 18, caractérisé en ce que ledit modèle choisi comprend la considération d'un service (a) déterministe. 20. Dispositif selon l'une des revendications 15 à 18, caractérisé en ce que ledit modèle choisi comprend la considération d'un service (a) aléatoire. 21. Dispositif selon l'une des revendications précédentes, caractérisé en ce que la matrice courante au moins comporte des coefficients (Alk) dynamiquement variables. 22. Dispositif selon l'une des revendications précédentes, caractérisé en ce que la matrice produit obtenue est un vecteur représenté par une matrice à colonne unique (Z(n)). 23. Procédé d'aide à la surveillance d'un système complexe, notamment d'un réseau de communication, caractérisé en ce qu'il comprend les étapes suivantes : a) obtenir des premières données représentatives de paramètres du réseau, b) construire une première matrice (Z(n)), selon un modèle choisi, en fonction desdites premières données, c) recevoir, à un instant choisi, au moins de secondes données représentatives d'événements dans le réseau, d) construire au moins une seconde matrice (A(n)) de structure dynamiquement variable, selon le modèle choisi, en fonction des secondes données, et e) effectuer sur lesdites matrices une opération formant produit selon l'algèbre dite MAX-PLUS, la matrice produit obtenue (Z(n+1)) étant représentative de l'état du réseau audit instant choisi. 24. Procédé selon la revendication 23, caractérisé en ce qu'il comporte en outre l'étape suivante : f) répéter, à des instants choisis, les étapes c), d) et e), tandis que la matrice produit obtenue devient la première matrice après l'étape e), ce qui permet de suivre une évolution temporelle de l'état du réseau auxdits instants choisis. 25. Procédé de simulation d'un système complexe, notamment de mécanismes et de contrôles de flux et de congestions dans un réseau de communication, caractérisé en ce qu'il comprend les étapes suivantes : a) obtenir des premières données représentatives de paramètres propres au réseau, b) construire une première matrice (Z(n)), selon un modèle choisi, en fonction desdites premières données, c) simuler des événements dans le réseau et prévoir au moins de secondes données représentatives desdits événements, d) construire au moins une seconde matrice (A(n)) selon le modèle choisi, en fonction desdites secondes données, et e) effectuer sur lesdites matrices une opération formant produit selon l'algèbre dite MAX-PLUS, la matrice produit obtenue (Z(n+1)) étant représentative d'un état du réseau subissant lesdits événements. 26. Procédé selon la revendication 25, caractérisé en ce qu'il comporte en outre l'étape suivante : f) répéter, pour des événements successifs, les étapes c), d) et e), tandis que la matrice produit obtenue devient la première matrice après l'étape e), ce qui permet de prévoir une évolution de l'état du réseau en fonction desdits événements.