FR 2991056 A1 20131129 FR 1254781 A 20120524 La présente invention concerne un système électronique à capteurs intégrés, un procédé d'estimation de valeur de grandeur physique de fonctionnement et un programme d'ordinateur correspondant. La présente invention s'applique en particulier au domaine des circuits intégrés. La publication intitulée « A fully integrated 32 nm Multiprobe for dynamic PVT measurements with complex digital SoC », de L. Vincent, E. Beigne, L. Alacoque, S. Lesecq, C. Bour, Ph. Maurine, publiée en 2011 dans « 2nd European Workshop on CMOS Process Tuning and Variability, Vari 2011 » décrit un système électronique à capteurs intégrés, comportant : un circuit électronique présentant, en fonctionnement, une valeur d'une grandeur physique de fonctionnement, et un dispositif de mesure comportant : - des capteurs intégrés dans le circuit électronique de manière à fournir des mesures brutes sensibles à ladite valeur de la grandeur physique de fonctionnement, et - une unité de sortie destinée à fournir des mesures utiles à partir des mesures brutes. Dans cette publication, l'unité de sortie ne réalise aucun traitement particulier de sorte que les mesures utiles sont égales aux mesures brutes, et sont sous la forme de fréquences, sensibles à la température des capteurs, à la tension d'alimentation des capteurs et au processus de fabrication des capteurs. Ces trois variables de fonctionnement forment donc la grandeur physique de fonctionnement mentionnée précédemment. Chaque fréquence dépend de ces trois variables de fonctionnement de sorte que la valeur précise de chacune des variables de fonctionnement n'est pas directement accessible à partir des fréquences. Cela signifie que les mesures brutes ne renseignent pas directement sur la valeur de la grandeur physique de fonctionnement. II peut ainsi être souhaité de prévoir un système électronique à capteurs intégrés, par exemple comme celui de la publication précédente ou de tout autre type, qui permette d'estimer la valeur de la grandeur physique de fonctionnement à partir des mesures utiles. L'invention a donc pour objet un système électronique à capteurs intégrés, comportant : un circuit électronique présentant, en fonctionnement, une valeur v d'une grandeur physique de fonctionnement, un dispositif de mesure comportant : - des capteurs intégrés dans le circuit électronique de manière à fournir des mesures brutes {X11,1 sensibles à la valeur v de la grandeur physique de fonctionnement, et - une unité de sortie destinée à fournir des mesures utiles IX'.} à partir des mesure brutes IXJ, comportant en outre : un dispositif de calcul de répartition destiné à déterminer une fonction de répartition courante F(x) associée aux mesures utiles une mémoire dans laquelle des associations {Ak} sont enregistrées, chaque association Ak associant une fonction de répartition modèle Fk (x) à une valeur modèle v*k de la grandeur physique de fonctionnement, un dispositif de calcul d'estimation destiné à déterminer une estimation V de la valeur 7 de la grandeur physique de fonctionnement à partir de la fonction de répartition courante F(x) et d'au moins une partie des associations {Ak } enregistrées. Ainsi, l'estimation V peut être déterminée de manière fiable, c'est-à-dire que l'estimation V obtenue est proche de la valeur De façon optionnelle, l'unité de sortie est destinée à combiner les mesures brutes {X11,1 pour calculer les mesures utiles en sommant deux à deux les mesures brutes {X11,1 de sorte que chaque mesure utile IX'.} soit la somme de deux mesures brutes {X11,1. De façon optionnelle également, le dispositif d'estimation comporte : une unité de sélection préalable destinée à sélectionner des fonctions de répartition modèles {Fkl(x)} dans la mémoire, une unité de test d'hypothèse destinée à déterminer, pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), une valeur distinctive Cki , et à réaliser un test d'hypothèse utilisant les valeurs distinctives Kkil afin de sélectionner statistiquement, parmi les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)}, des fonctions de répartition modèles {Fk2(x)}, - une unité de calcul d'estimation destinée à déterminer l'estimation V à partir des valeurs modèles {v:21 associées aux fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)} et des valeurs distinctives {Ck2} correspondantes. De façon optionnelle également, l'unité de sélection préalable est destinée à sélectionner les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)} en fonction d'une estimation intermédiaire V' de la valeur 7 de la grandeur physique de fonctionnement, et de manière prédéterminée en cas d'absence d'estimation intermédiaire, et l'unité de calcul d'estimation est destinée, d'une part, à déterminer l'estimation intermédiaire V' à partir des valeurs {v:21 associées aux fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)} et des valeurs distinctives {Ck2} correspondantes, et, d'autre part, suivant le résultat d'un test d'arrêt, soit à fournir l'estimation intermédiaire V' à l'unité de sélection préalable de manière à former une boucle récursive de détermination d'estimations intermédiaires successives, soit à fournir la dernière estimation intermédiaire V' déterminée comme estimation V . De façon optionnelle également, l'unité de test d'hypothèse comporte : une unité de calcul de valeur distinctive destinée, pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), à déterminer une valeur discriminante tu d'une fonction discriminante Tki(F) prédéterminée dépendant de la fonction de répartition courante F(x) et de la fonction de répartition modèle préalablement sélectionnée Fk(x), et à déterminer la valeur distinctive Cki à partir de la valeur discriminante tu , et une unité de sélection statistique destinée à sélectionner, parmi les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)}, les fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)} à partir des valeurs distinctives {Cu}. De façon optionnelle également, l'unité de calcul de valeur distinctive est destinée, pour chaque valeur discriminante tu , à déterminer la probabilité critique pcki associée à cette valeur discriminante tu , et dans lequel la valeur distinctive Cki est la probabilité critique pcki . De façon optionnelle également, l'unité de sélection statistique est destinée à sélectionner, parmi les fonctions de répartition préalablement sélectionnées {Fu (d, celles pour lesquelles la valeur distinctive Cki est supérieure à un seuil a . De façon optionnelle également, le seuil a est égal à un pourcentage prédéterminé de la plus grande des valeurs distinctives {Cu}. L'invention a également pour objet un procédé d'estimation de valeur de grandeur physique de fonctionnement, comportant : recevoir des mesures utiles IX'.} d'un dispositif de mesure comportant des capteurs intégrés dans un circuit électronique de manière à fournir des mesures brutes {Xm} sensibles à une valeur v d'une grandeur physique de fonctionnement que le circuit électronique présente en fonctionnement, les mesures utiles IX'.} étant déterminées à partir des mesure brutes {Xm}, déterminer une fonction de répartition courante F(x) associée aux mesures utiles déterminer une estimation 'if- de la valeur 7 de la grandeur physique de fonctionnement à partir de la fonction de répartition courante F(x) et d'au moins une partie d'associations {Ak} enregistrées dans une mémoire, chaque association Ak associant une fonction de répartition modèle Fk (x) à une valeur modèle v*k de la grandeur physique de fonctionnement. L'invention a également pour objet un programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions pour l'exécution des étapes d'un procédé selon l'invention lorsque ledit programme est exécuté sur un ordinateur. L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels : la figure 1 représente schématiquement la structure générale d'un premier système électronique à capteurs intégrés selon l'invention, la figure 2 illustre les étapes successives d'un procédé de fonctionnement du système électronique à capteurs intégrés de la figure 1, - la figure 3 représente schématiquement la structure générale d'un second système électronique à capteurs intégrés selon l'invention, et - la figure 4 illustre les étapes successives d'un procédé de fonctionnement du système électronique à capteurs intégrés de la figure 3. En référence à la figure 1, un premier système électronique 100 à capteurs intégrés va être décrit. Le système électronique à capteurs intégrés 100 comporte tout d'abord une puce 102 sur laquelle des composants sont intégrés afin de former un circuit électronique 104 (indiqué en pointillés sur la figure 1) destiné à réaliser une ou plusieurs fonctions prédéterminées. La puce 102 est formée d'une tranche de semi- conducteur (« wafer » en anglais) sur laquelle sont gravés ou imprimés les composants du circuit électronique 104. Le système électronique à capteurs intégrés 100 comporte en outre un dispositif de mesure 106. Le dispositif de mesure 106 comporte tout d'abord des capteurs 108 intégrés à la puce 102 de manière à surveiller le fonctionnement d'une zone 110 du circuit électronique 104 s'étendant autour des capteurs 108. Les capteurs 108 sont ainsi eux aussi imprimés sur la tranche de semi-conducteur formant la puce 102. Cependant, les capteurs 108 ne participent pas à la réalisation de la ou des fonctions prédéterminées du circuit électronique 104. Plus précisément, lors du fonctionnement du circuit électronique 104, la zone surveillée 110 présente une valeur v d'une grandeur physique de fonctionnement, et les capteurs 108 sont destinés à surveiller cette valeur v . Une grandeur physique est toute propriété ou ensemble de propriétés de la nature qui peuvent être quantifiées par la mesure ou le calcul, et dont les différentes valeurs possibles peuvent par exemple s'exprimer à l'aide d'un nombre, par exemple réel ou complexe, généralement accompagné d'une unité de mesure. Les capteurs 108 sont destinés à être sensibles à la même grandeur physique de fonctionnement que la zone surveillée 110, c'est-à-dire à une grandeur physique de même nature. De manière idéale, lors de leur fonctionnement, les capteurs 108 présentent la même valeur v que la zone surveillée 110. Cependant, en pratique, il arrive qu'ils présentent une valeur y différente de la valeur 7, mais néanmoins proche de la valeur v , ou du moins présentant les mêmes variations que la valeur 7. Ainsi, les capteurs 108 sont destinés à fournir, à chaque fois qu'ils sont interrogés, M mesures, appelées mesures brutes {X11,1 ( m =1...M ), qui dépendent de la valeur y de la grandeur physique de fonctionnement. Par exemple, chaque capteur fournit une des mesures brutes IXJ, de sorte qu'il y a autant de capteurs 106 que de mesures brutes {Xm}. Comme expliqué précédemment, puisque la valeur y est sensiblement égale à la valeur v , ou du moins présente les mêmes sensibilités et/ou variations, les mesures brutes {X11,1 sont donc sensibles à la valeur v de la grandeur physique de fonctionnement de la zone surveillée 110. Dans l'exemple décrit, la grandeur physique de fonctionnement est bidimensionnelle : elle regroupe la température et la tension d'alimentation. Ainsi, la valeur 7 est composée de la valeur i de la température de la zone surveillée 110 et de la valeur V de tension d'alimentation de la zone surveillée, soit : 7 = (T, V), tandis que la valeur y est composée de la valeur T de la température des capteurs 108 et de la valeur V de la tension d'alimentation des capteurs 108, soit y = (T, V). Dans l'exemple décrit, l'intégration des capteurs 108 signifie que les capteurs 108 sont imprimés sur la même puce 102 que le circuit électronique 104, en particulier selon le même processus d'impression. Cela signifie en outre que les capteurs 108 sont alimentés par la même alimentation, fournissant la tension d'alimentation, que le circuit électronique 104 et en particulier que la zone surveillée 110. Dans l'exemple décrit, les capteurs 108 sont tous de même type et produisent donc le même type de mesure brute. Cependant, le comportement de chaque capteur (c'est-à-dire la mesure brute qu'il fournit) soumis aux mêmes variations de la grandeur physique de fonctionnement que les autres capteurs est de préférence le plus différent possible de celui des autres capteurs. Ainsi, les informations apportées par les capteurs 108 ne sont pas redondantes. Par exemple, le dispositif appelé « Multiprobe » et décrit dans la publication intitulée « A Fully Integrated 32 nm MultiProbe for Dynamic PVT Measurements within Complex Digital SoC » de L. Vincent, E. Beigné, L. Alocoque, S. Lesecq, C. Bour et P. Maurine, et publiée dans « 2nd European Workshop on CMOS Variability », VARI 2011, 2011, est utilisé pour former les capteurs 108. Comme cela est décrit dans cette publication, chaque capteur 108 est alors un oscillateur en anneau dont la fréquence d'oscillation est sensible à sa valeur de température et à sa valeur de tension d'alimentation. Un tel oscillateur en anneau fournit, en tant que mesure brute, une mesure de fréquence d'oscillation. Le dispositif de mesure 106 comporte en outre une unité de sortie 112 destinée à déterminer N mesures, appelées mesures utiles IX'.} ( n =1...N ), à partir des mesure brutes {Xm}. De préférence, les mesures utiles IX'.} sont plus nombreuses que les mesures brutes {Xm}. En effet, cela permet d'augmenter la puissance du test statistique qui sera réalisé par les éléments décrits plus loin, c'est-à-dire que cela permet d'augmenter la probabilité de rejeter l'hypothèse nulle de ce test statistique. Dans l'exemple décrit, l'unité de sortie 112 est destinée à combiner les mesures brutes {X.,} pour calculer les mesures utiles Par exemple, l'unité de sortie 112 est destinée à sommer deux à deux des mesures brutes {Xm} de sorte que chaque mesure utile IX'.} soit la somme de deux mesures brutes En variante, l'unité de sortie 112 peut être destinée à faire la différence deux à deux des mesures brutes {X.,} de sorte que chaque mesure utile IX'.} soit la différence de deux des mesures brutes IX.J. Lors des études réalisées, il a été trouvé que la somme deux à deux des mesures brutes donnait de meilleurs résultats que la différence deux à deux. Encore en variante, l'unité de sortie 112 peut être destinée à faire la différence trois à trois des mesures brutes {X11,1, ou tout autre type de combinaison possible. Le système électronique à capteurs intégrés 100 comporte en outre un dispositif de calcul de répartition 114 destiné à déterminer une fonction de répartition, appelée fonction de répartition courante F(x), associée aux mesures utiles Dans l'exemple décrit, la fonction de répartition courante F(x) représente la proportion des mesures utiles IX'.} dont la valeur est inférieure ou égale à x, et est donnée par : où 11- \ %--IN sinon. F(X) 1 -22100 1(K.), N est une fonction valant 1 lorsque X'n est compris dans 1- C.°, x] et 0 Ainsi, la fonction de répartition courante F(x) est une fonction à paliers croissante, et le nombre de paliers est égal au nombre N de mesures utiles Le système électronique à capteurs intégrés 100 comporte en outre une mémoire 116 dans laquelle K associations, appelées associations disponibles {Ak} (avec k =1...K ), sont enregistrées. Chaque association disponible Ak associe une fonction de répartition modèle disponible Fk(x) à une valeur modèle disponible v*k de la grandeur physique de fonctionnement. Les associations disponibles {Ak} sont par exemple enregistrées sous la forme d'une table de correspondance (ou LUT de l'anglais « Look-Up Table »). Chaque fonction de répartition modèle disponible Fk(x) et valeur modèle disponible vk* associée est par exemple obtenue au moyen d'une simulation du circuit électronique 104 et des capteurs 108, et/ou par une caractérisation des capteurs 108, et/ou par calibration des capteurs 108. Les valeurs modèles disponibles Ivkl forment ainsi un pavage de l'espace des valeurs possibles. Dans un mode de réalisation, le pavage est régulier. Dans l'exemple décrit, cela signifie que l'écart entre deux valeurs modèles de température ou de tension est toujours le même. Dans un autre mode de réalisation, le pavage est irrégulier. Dans l'exemple décrit, cela signifie que l'écart entre deux valeurs modèles de température ou de tension varie suivant les deux valeurs modèles considérées. En particulier cet écart peut être plus resserré pour les zones de l'espace des valeurs possibles les plus indécises. Dans encore un autre mode de réalisation, le pavage est en quinconce, ce qui permet de limiter le nombre de valeurs modèles disponibles tout en minimisant la distance entre ces valeurs. Dans l'exemple décrit, chaque fonction de répartition modèle disponible Fk(x) est obtenue en considérant des valeurs de mesures utiles modèles, ces dernières étant obtenues de manière empirique sur un système électronique à capteurs intégrés de test et/ou par modélisation et calculs. Ainsi, la fonction de répartition modèle disponible Fk(x) est également une fonction à paliers croissante, le nombre de paliers, noté Mk , étant égal au nombre de mesures utiles modèles. Dans un mode de réalisation particulier, toutes les fonctions de répartition modèles disponibles Fk(x) présentent le même nombre de paliers, par exemple égal au nombre N de mesures utiles Le système électronique à capteurs intégrés 100 comporte en outre un dispositif d'estimation 118 destiné à déterminer une estimation V de la valeur v de la grandeur physique de fonctionnement, et donc également de la valeur v de la grandeur physique de fonctionnement, à partir de la fonction de répartition courante F(x) et d'au moins une partie des associations disponibles {Ak } enregistrées dans la mémoire 116. Dans l'exemple décrit, le dispositif d'estimation 118 comporte tout d'abord une unité de sélection préalable 120 destinée à sélectionner dans la mémoire 116, parmi les fonctions de répartition modèles disponibles {Fk(x)}, des fonctions de répartition modèles, appelées fonctions de répartition modèles préalablement sélectionnées {Fkl(x)} (avec kle [1,1(]). Dans l'exemple de la figure 1, l'unité de sélection préalable 120 est destinée à sélectionner toutes les fonctions de répartition modèles disponibles {Fk (x )}. Dans l'exemple décrit, le dispositif d'estimation 118 comporte en outre une unité de test d'hypothèse 122 destinée à réaliser un test d'hypothèse afin de sélectionner statistiquement, parmi les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)}, des fonctions de répartition modèles, appelées fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)} (avec k2E [1,1(]). Pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), le test d'hypothèse présente comme hypothèse nulle : « la fonction de répartition courante F(x) est égale à la fonction de répartition Fkl(x) », et comme hypothèse contraire : « la fonction de répartition courante F(x) n'est pas égale à la fonction de répartition Fkl(x) ». Le test d'hypothèse est par exemple le test d'hypothèse de Smirnov (variante moderne ou bien variante historique) ou bien le test d'hypothèse de Cramér-von Mises. Dans l'exemple décrit, l'unité de test d'hypothèse 122 comporte tout d'abord une unité de calcul de valeur distinctive 124 destinée à déterminer, pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), une valeur distinctive Cki . La valeur distinctive Cki représente un écart entre la fonction de répartition courante F(x) et la fonction de répartition modèle préalablement sélectionnée Fu (X) . Dans l'exemple décrit, l'unité de calcul de valeur distinctive 124 est destinée, pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), à déterminer une valeur, appelée valeur discriminante tu , d'une fonction prédéterminée, appelée fonction discriminante Tki(F). La fonction discriminante Tkl(F) dépend de la fonction de répartition courante F(x) et de la fonction de répartition modèle préalablement sélectionnée Fu (X) . Dans l'exemple décrit, la fonction discriminante Tkl(F) est choisie de sorte que la valeur tkl de Tkl(F) a tendance à s'éloigner de 0 lorsque l'hypothèse nulle n'est pas vraie (test d'hypothèse bilatéral). Par exemple, dans le cas du test d'hypothèse de Smirnov, la fonction discriminante Tkl(F) est basée sur l'écart maximal entre les deux fonctions de répartition testées. Ainsi, le test d'hypothèse de Smirnov s'appuie sur l'observation d'un seul point des fonctions de répartition testées, à savoir celui correspondant à la distance maximale entre les deux fonctions de répartition. La fonction discriminante Tkl(F) est alors : Tkl (F) : tkl = suplFki(x)-F(xl . Dans le cas du test de Cramér-von Mises, la fonction discriminante Tkl(F) mesure la différence entre les deux fonctions de répartition testées sur l'intégralité de leurs domaines, ce qui en fait une bonne alternative au test d'hypothèse de Smirnov. La fonction discriminante Tkl(F) vaut alors : Tkl(F): tkl = NM" T[F" (X) - F(X)12dF(X) . N±M" _00 Dans l'exemple décrit, les fonctions de répartition étant des fonctions à paliers, appelées couramment fonctions de répartition empiriques, la fonction discriminante Tkl(F) est calculée par : 4NM" - 1 T" (F) : NM"Mkl)16(1\4" N) avec : N MkI (I; -i)2 + (si - j)2 , i=1 j=1 où, les valeurs utiles étant classées par ordre croissant dans un premier vecteur et les valeurs utiles modèles étant classées par ordre croissant dans un second vecteur, ces deux vecteurs étant combinés en un unique vecteur, appelé vecteur combiné, de sorte que les valeurs que le vecteur combiné contient sont également classées par ordre croissant, i est le rang des valeurs utiles dans le premier vecteur et r1 est le rang de ces valeurs utiles dans le vecteur combiné, et j est le rang des valeurs utiles modèles dans le second vecteur et s est le rang de ces valeurs utiles modèles dans le vecteur combiné. L'unité de calcul de valeur distinctive 120 est en outre destinée à déterminer la valeur distinctive Ckl à partir de la valeur discriminante tu . Dans un mode de réalisation, la valeur distinctive Ckl est la valeur discriminante tkl. C'est par exemple le cas pour le test d'hypothèse de Smirnov historique. Dans un autre mode de réalisation, la valeur distinctive Ckl est une probabilité, appelée probabilité critique pckp permettant d'obtenir une valeur de la fonction discriminante Tkl(F) égale à la valeur tkl ou plus extrême (c'est-à-dire dans le cas du test d'hypothèse bilatéral, une valeur t telle que t > C'est par exemple le cas pour le test d'hypothèse de Smirnov moderne (communément appelé test de Kolmogorov-Smirnov), où la probabilité critique pcki vaut : pcki = (-1)J+l e(-202) avec : Y NM" Tki (F.) - N + Mid En pratique, comme la somme sur j converge très rapidement vers sa valeur à l'infini, seul les J premiers termes de la somme sont déterminés, de sorte que pcki est déterminé par : pcki = (-1) 1 (-2k2J2) Par exemple, J est compris entre 3 et 100, ce qui suffit pour approcher de manière convenable la valeur à l'infini. Dans le cas du test d'hypothèse de Cramér-von Mises, la valeur distinctive Ckl peut également être la probabilité critique pckp Dans ce cas, la probabilité critique pcki vaut : Tki (F)- eTki (F) 1 pcki - .V45. Var(Tki (F)) + 6 avec : 1 eT = 1 + 6 6(Mki + N) et : 1 Mki N+1 4MkiN(Mki +N)-3(1VC +N2)-2MkiN Var(Tki(F))- 45 (Mu + N)2 4MkiN Dans l'exemple décrit, l'unité de test d'hypothèse 122 comporte en outre une unité de sélection statistique 126 destinée à sélectionner, parmi les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)}, à partir de leurs valeurs distinctives {Ckl}, des fonctions de répartition modèles, appelées fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)}. Plus précisément, l'unité de sélection statistique 126 est destinée à comparer la valeur distinctive Ckl à un seuil pour sélectionner, en fonction de la comparaison, les fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)}. Dans l'exemple décrit, l'unité de sélection statistique 126 est destinée à fournir des identifiants des fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)}, par exemples les indices k2, et les valeurs distinctives {Ck2} correspondantes. Dans le cas où la valeur distinctive Ckl est la probabilité critique pckp l'unité de sélection statistique 126 est destinée à sélectionner les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)} dont la valeur distinctive Ckl est supérieure à un seuil a . En effet, si la probabilité critique pcki est inférieure au seuil , cela signifie que l'hypothèse nulle est rejetée. Dans l'exemple décrit, le seuil a est égal à un pourcentage prédéterminé de la plus grande des valeurs distinctives {Ckl}. De préférence, le pourcentage prédéterminé est supérieur à 50%. Par exemple, le pourcentage prédéterminé est de 90% ou plus, ce qui permet de ne sélectionner que les fonctions de répartition modèles les plus probables. Il a en effet été trouvé que ce dernier cas donnait de meilleurs résultats. Dans le cas où la valeur distinctive Ckl est la valeur discriminante tu , afin de garantir qu'au moins une fonction de répartition modèle est sélectionnée et afin de diminuer la taille de mémoire nécessaire en évitant l'enregistrement de tables de fractiles, le seuil a est de préférence défini de la même manière que précédemment comme un pourcentage prédéterminé de la plus grande des valeurs distinctives {Ckl}. Dans l'exemple décrit, le dispositif d'estimation 118 comporte en outre une unité de calcul d'estimation 128 destinée à déterminer l'estimation V à partir, d'une part, des valeurs modèles {V*k2} enregistrées dans la mémoire 116 et associées aux fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)} et, d'autre part, des valeurs distinctives {Ck2} correspondantes. Par exemple, l'estimation V est égale à la moyenne des valeurs sélectionnées tyk*21 pondérée par les valeurs distinctives {Ck2} correspondantes. Dans le mode de réalisation décrit, le système électronique à capteurs intégrés 100 comporte en outre un dispositif de commande 130 destiné à commander le fonctionnement du circuit électronique 104 en fonction de l'estimation . Le dispositif de commande 130 a par exemple pour but d'ajuster au mieux les performances du circuit électronique 104 dans un contexte de faible consommation électrique. Par ailleurs, dans un mode de réalisation, tout ou partie du dispositif de calcul de répartition 114 et du dispositif d'estimation 118 est implémentée sous la forme d'un programme d'ordinateur exécuté par une unité de traitement (non représentée). On notera que, de manière générale, tous les éléments décrits précédemment sauf les capteurs 108, la mémoire 116 et le dispositif de commande 130 peuvent être implémentés sous la forme d'un programme d'ordinateur exécuté par une unité de traitement (non représentée). En référence à la figure 2, un procédé 200 de fonctionnement du système électronique à capteurs intégrés 100 va être décrit. Au cours d'une étape 202, les capteurs 108 sont interrogés et fournissent les mesures brutes {X11,1. Au cours d'une étape 204, l'unité de sortie 112 combine les mesures brutes {X11,1 pour calculer les mesures utiles Au cours d'une étape 206, le dispositif de calcul de répartition 114 reçoit les mesures utiles IX'.} et détermine la fonction de répartition courante F(x) associée aux mesures utiles Au cours d'une étape 208, le dispositif d'estimation 118 détermine une estimation V de la valeur y de la grandeur physique de fonctionnement, et donc également de la valeur v de la grandeur physique de fonctionnement, à partir de la fonction de répartition courante F(x) et d'au moins une partie des associations disponibles {Ak} enregistrées dans la mémoire 116. L'étape 208 comporte les étapes 210 à 216 suivantes. Au cours d'une étape 210, l'unité de sélection préalable 120 sélectionne dans la mémoire 116, parmi les fonctions de répartition modèles disponibles {Fk (x)}, des fonctions de répartition modèles préalablement sélectionnées {Fkl(x)}. Dans l'exemple décrit, la sélection est prédéterminée. Par exemple, toutes les fonctions de répartition modèles disponibles {Fk(x)} sont sélectionnées. Au cours d'une étape 212, l'unité de test d'hypothèse 122 réalise un test d'hypothèse afin de sélectionner, parmi les fonctions de répartition modèles préalablement sélectionnées IFki (d, des fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)}. L'étape 212 comporte les étapes 214 et 216 suivantes. Au cours d'une étape 214, l'unité de calcul de valeur distinctive 124 détermine, pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), la valeur distinctive Ckl correspondante. Comme expliqué précédemment, suivant le mode de réalisation, cette valeur distinctive Ckl est soit la valeur discriminante tu , soit la probabilité critique pcki associée à cette valeur discriminante tu . Au cours d'une étape 216, l'unité de sélection statistique 126 sélectionne, parmi les fonctions de répartition modèles préalablement sélectionnées IFki (d, à partir de leurs valeurs distinctives {Ckl}, les fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)}. L'unité de sélection statistique fournit en outre des identifiants des fonctions de répartition modèles statistiquement sélectionnées {Fk2(x )} et les valeurs distinctives {Ck2} correspondantes. Au cours d'une étape 218, l'unité de calcul d'estimation 128 détermine l'estimation V à partir, d'une part, des valeurs Ivk*21 obtenues dans la mémoire 116 grâce aux identifiants et, d'autre part, des valeurs distinctives {Ck2}. Au cours d'une étape 220, le dispositif de commande 130 commande le fonctionnement du circuit électronique 104 en fonction de l'estimation V . Le procédé reprend ensuite à l'étape 202 pour une nouvelle estimation. Les estimations sont par exemple réalisées à intervalles de temps réguliers. En référence à la figure 3, un second système électronique 300 à capteurs intégrés va être décrit. Le système électronique à capteurs intégrés 300 est identique à celui de la figure 1, si ce n'est que l'unité de sélection préalable, portant à présent la référence 302, est destinée à sélectionner, parmi les fonctions de répartition modèles disponibles {Fk (x)} de la mémoire 116, les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)} en fonction d'une estimation intermédiaire V' de la valeur y de la grandeur physique de fonctionnement, et donc également de la valeur 7 de la grandeur physique de fonctionnement, et de manière prédéterminée en cas d'absence d'estimation intermédiaire. Par exemple, en cas d'absence d'estimation intermédiaire, l'unité de sélection préalable 302 est destinée à sélectionner les fonctions de répartition modèles dont les valeurs modèles associées forment, dans l'exemple décrit, un pavage régulier de l'espace des valeurs possibles. Cependant, comme indiqué précédemment, le pavage pourrait être irrégulier. Lorsqu'une estimation intermédiaire V' est disponible, l'unité de sélection préalable est par exemple destinée à sélectionner des fonctions de répartition modèles dont les valeurs modèles associées se trouvent autour de l'estimation intermédiaire V' , par exemple à une distance inférieure ou égale à une distance prédéterminée. De préférence, les valeurs modèles associées aux fonctions de répartition modèles sélectionnées sont séparées par des distances décroissantes pour chaque nouvelle estimation intermédiaire V' . Cela permet de resserrer de manière récursive le pavage autour de la valeur la plus probable. L'unité de sélection préalable 302 permet ainsi de limiter les calculs réalisés, et donc le coût en temps de calcul, tout en gardant de bons résultats. Par exemple, il a été trouvé que la même précision que pour le système électronique à capteurs intégrés de la figure 1 est obtenue avec une économie de calcul d'environ 60%. En outre, l'unité de calcul d'estimation, portant à présent la référence 304, est destinée à déterminer l'estimation intermédiaire V' à partir des valeurs distinctives {Ck2} déterminées par l'unité de sélection statistique 126. L'estimation intermédiaire . 1 V' est par exemple égale à la moyenne des valeurs sélectionnées {Vk2} pondérée par les valeurs distinctives {Ck2} correspondantes. L'unité de calcul d'estimation 304 est en outre destinée à évaluer un test d'arrêt portant sur l'estimation intermédiaire et, suivant le résultat du test d'arrêt, soit à fournir l'estimation intermédiaire V' à l'unité de sélection préalable 302 de manière à former une boucle récursive de détermination d'estimations intermédiaires V' , soit à fournir la dernière estimation intermédiaire V' déterminée comme estimation V . En référence à la figure 4, un procédé 400 de fonctionnement du système électronique à capteurs intégrés 300 va être décrit. Le procédé de fonctionnement 400 est identique à celui de la figure 2, si ce n'est que l'étape 210 est remplacée par une étape 402 au cours de laquelle l'unité de sélection préalable 302 sélectionne, parmi les fonctions de répartition modèles disponibles {Fk (x)} de la mémoire 116, les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)} en fonction d'une estimation intermédiaire V' de la valeur y de la grandeur physique de fonctionnement, et donc également de la valeur 7 de la grandeur physique de fonctionnement, et de manière prédéterminée en cas d'absence d'estimation intermédiaire. En outre, l'étape 218 est remplacée par les étapes 404 à 410 suivantes. Au cours d'une étape 404, l'unité de calcul d'estimation 304 détermine l'estimation intermédiaire V' à partir des valeurs modèles tyk*21 et des valeurs distinctives {Ck2}. Au cours d'une étape 406, l'unité de calcul d'estimation 304 évalue un test d'arrêt portant sur l'estimation intermédiaire V' . Au cours d'une étape 408, si le résultat du test d'arrêt est négatif, l'unité de calcul d'estimation 304 fournit l'estimation intermédiaire V' à l'unité de sélection préalable 302 et le procédé 400 retourne à l'étape 402 de manière à former une boucle récursive de détermination d'estimations intermédiaires V' successives. Au cours d'une étape 410, si le résultat du test d'arrêt est positif, l'unité de calcul d'estimation 304 fournit la dernière estimation intermédiaire V' déterminée comme estimation V . Il apparaît clairement qu'un système électronique à capteurs intégrés tel que celui décrit précédemment permet d'estimer la valeur d'une grandeur physique de fonctionnement d'une zone d'un circuit électronique, et l'évolution de cette valeur au cours du temps. Le système électronique à capteurs intégrés permet en particulier de suivre une grandeur physique multidimensionnelle quand bien même les capteurs fournissent des mesures de nature identique. Par ailleurs, l'invention n'est pas limitée aux modes de réalisation décrits précédemment. Il apparaîtra en effet à l'homme de l'art que diverses modifications peuvent être apportées aux modes de réalisation décrits ci-dessus, à la lumière de l'enseignement qui vient de lui être divulgué. En particulier, la grandeur physique de fonctionnement n'est pas obligatoirement bidimensionnelle. Dans d'autres modes de réalisation, elle pourrait être monodimensionnelle ou multidimensionnelle à plus de deux dimensions. En outre, la grandeur physique de fonctionnement pourrait inclure, en plus ou à la place de la température et de la tension d'alimentation, le processus de fabrication du circuit électronique. En effet, des circuits électroniques réalisés en suivant un même mode opératoire peuvent différer les uns des autres, par exemple pour la raison que des séries différentes de pièces ont été utilisées. Dans ce cas, les différentes valeurs caractérisant le processus de fabrication peuvent par exemple s'exprimer par un mot qualifiant les performances du circuit par rapport aux performances attendues, et non par un chiffre. En outre, les mesures brutes ne sont pas obligatoirement combinées. Dans ce cas, l'unité de sortie peut être un simple fil et les mesures utiles sont alors égales aux mesures brutes. En outre, dans une variante, l'unité de sélection statistique est destinée à déterminer si la plus grande des valeurs distinctives est inférieure ou non à un seuil prédéterminé. Si la plus grande des valeurs distinctives est inférieure au seuil prédéterminé, cela signifie qu'aucune fonction de répartition modèle ne ressemble à la fonction de répartition courante, c'est-à-dire qu'il est peu probable que les mesures utiles soit issues de l'une des fonctions de répartition modèles. Dans ce cas, l'unité de sélection statistique est destinée à modifier le fonctionnement de l'une ou plusieurs parmi : l'unité de sortie, l'unité de sélection préalable, l'unité de calcul de valeur distinctive et l'unité de sélection statistique, et à faire recommencer la détermination de l'estimation (finale ou intermédiaire suivant le mode de réalisation) à partir des mesures brutes. Par exemple, l'unité de sélection statistique peut modifier le type de combinaison utilisé pour obtenir les mesures utiles à partir des mesures brutes. Ou bien l'unité de sélection statistique peut modifier la manière avec laquelle l'unité de sélection préalable sélectionne des fonctions de répartition modèles disponibles. Ou bien l'unité de sélection statistique peut modifier le test utilisé pour le calcul de la valeur distinctive. Ou bien l'unité de sélection statistique peut modifier le seuil qu'elle utilise pour sélectionner les fonctions de répartition modèles statistiquement sélectionnées. Dans les revendications qui suivent, les termes utilisés ne doivent pas être interprétés comme limitant les revendications aux modes de réalisation exposés dans la présente description, mais doivent être interprétés pour y inclure tous les équivalents que les revendications visent à couvrir du fait de leur formulation et dont la prévision est à la portée de l'homme de l'art en appliquant ses connaissances générales à la mise en oeuvre de l'enseignement qui vient de lui être divulgué. REVENDICATIONS1. Système électronique à capteurs intégrés (100 ; 300), comportant : un circuit électronique (104) présentant, en fonctionnement, une valeur v d'une grandeur physique de fonctionnement, un dispositif de mesure (106) comportant : - des capteurs (108) intégrés dans le circuit électronique (104) de manière à fournir des mesures brutes {Xm} sensibles à la valeur de la grandeur physique de fonctionnement, - une unité de sortie (112) destinée à fournir des mesures utiles IX'.} à partir des mesures brutes IXJ, caractérisé en ce qu'il comporte en outre : un dispositif de calcul de répartition (114) destiné à déterminer une fonction de répartition courante F(x) associée aux mesures utiles IX'. 1, une mémoire (116) dans laquelle des associations {Ak} sont enregistrées, chaque association Ak associant une fonction de répartition modèle Fk(x) à une valeur modèle vk* de la grandeur physique de fonctionnement, un dispositif de calcul d'estimation (118) destiné à déterminer une estimation V de la valeur 7 de la grandeur physique de fonctionnement à partir de la fonction de répartition courante F(x) et d'au moins une partie des associations {Ak} enregistrées. 2. Système électronique à capteurs intégrés (100 ; 300) selon la revendication 1, dans lequel l'unité de sortie (112) est destinée à combiner les mesures brutes {Xm} pour calculer les mesures utiles IX'.}, en sommant deux à deux les mesures brutes {Xm} de sorte que chaque mesure utile IX'.} soit la somme de deux mesures brutes {Xm}. 3. Système électronique à capteurs intégrés (100 ; 300) selon la revendication 1 ou 2, dans lequel le dispositif d'estimation (118) comporte : une unité de sélection préalable (120 ; 302) destinée à sélectionner des fonctions de répartition modèles {Fki(x)} dans la mémoire (116),une unité de test d'hypothèse (122) destinée à déterminer, pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), une valeur distinctive Cki , et à réaliser un test d'hypothèse utilisant les valeurs distinctives Kkil afin de sélectionner statistiquement, parmi les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)}, des fonctions de répartition modèles {F' (x)}, une unité de calcul d'estimation (128 ; 304) destinée à déterminer l'estimation V à partir des valeurs modèles Ivk*21 associées aux fonctions de répartition modèles statistiquement sélectionnées {Fk2(x )} et des valeurs distinctives {Ck2} correspondantes. 4. Système électronique à capteurs intégrés (300) selon la revendication 3, dans lequel : l'unité de sélection préalable (302) est destinée à sélectionner les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)} en fonction d'une estimation intermédiaire V' de la valeur 7 de la grandeur physique de fonctionnement, et de manière prédéterminée en cas d'absence d'estimation intermédiaire, et l'unité de calcul d'estimation (304) est destinée, d'une part, à déterminer l'estimation intermédiaire V' à partir des valeurs {v*k2} associées aux fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)} et des valeurs distinctives {Ck2} correspondantes, et, d'autre part, suivant le résultat d'un test d'arrêt, soit à fournir l'estimation intermédiaire V' à l'unité de sélection préalable (302) de manière à former une boucle récursive de détermination d'estimations intermédiaires successives, soit à fournir la dernière estimation intermédiaire V' déterminée comme estimation V. 5. Système électronique à capteurs intégrés (100 ; 300) selon la revendication 3 ou 4, dans lequel l'unité de test d'hypothèse (122) comporte : une unité de calcul de valeur distinctive (124) destinée, pour chaque fonction de répartition modèle préalablement sélectionnée Fkl(x), à déterminer une valeur discriminante tu d'une fonction discriminanteTkl(F) prédéterminée dépendant de la fonction de répartition courante F(x) et de la fonction de répartition modèle préalablement sélectionnée Fk(x), et à déterminer la valeur distinctive Ckl à partir de la valeur discriminante tu , - une unité de sélection statistique (126) destinée à sélectionner, parmi les fonctions de répartition modèles préalablement sélectionnées {Fkl(x)}, les fonctions de répartition modèles statistiquement sélectionnées {Fk2(x)} à partir des valeurs distinctives {Cu}. 6. Système électronique à capteurs intégrés (100 ; 300) selon la revendication 5, dans lequel l'unité de calcul de valeur distinctive (124) est destinée, pour chaque valeur discriminante tu , à déterminer la probabilité critique pcki associée à cette valeur discriminante tu , et dans lequel la valeur distinctive Ckl est la probabilité critique pcki . 7. Système électronique à capteurs intégrés (100 ; 200 ; 300) selon la revendication 5 ou 6, dans lequel l'unité de sélection statistique (126) est destinée à sélectionner, parmi les fonctions de répartition préalablement sélectionnées {Fkl(x)}, celles pour lesquelles la valeur distinctive Ckl est supérieure à un seuil a . 8. Système électronique à capteurs intégrés (100 ; 300) selon la revendication 7, dans lequel le seuil a est égal à un pourcentage prédéterminé de la plus grande des valeurs distinctives {Cu}. 9. Procédé d'estimation de valeur de grandeur physique de fonctionnement, comportant : - recevoir (206) des mesures utiles IX'.} d'un dispositif de mesure (106) comportant des capteurs (108) intégrés dans un circuit électronique (104) de manière à fournir des mesures brutes {X11,1 sensibles à une valeur v d'une grandeur physique de fonctionnement que le circuit électronique (104) présente en fonctionnement, les mesures utiles IX'.} étant déterminées à partir des mesure brutes IXJ, caractérisé en ce qu'il comporte : - déterminer (206) une fonction de répartition courante F(x) associée aux mesures utiles IX'.},déterminer (208) une estimation V de la valeur 7 de la grandeur physique de fonctionnement à partir de la fonction de répartition courante F(x) et d'au moins une partie d'associations {Ak} enregistrées dans une mémoire (116), chaque association Ak associant une fonction de répartition modèle Fk(x) à une valeur modèle vk* de la grandeur physique de fonctionnement. 10. Programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, caractérisé en ce qu'il comprend des instructions pour l'exécution des étapes d'un procédé selon la revendication 9 lorsque ledit programme est exécuté sur un ordinateur.