Brevet FR2648251A1 — publié le 19901214

DESCRIPTION "METAODE D'APPRENTISSAGE ET STRUCTURE DE RESEAU DE NEURONES'. L'invention concerne une méthode d'apprentissage aise en oeuvre dans un réseau de neurones qui effectue des phases d'apprentissage pour adapter ses coefficients synaptiques'à partir d'exemples à l'aide de l'algorithme de rétropropagation du gradient de l'erreur. Elle concerne également une structure de réseau de neurones et un calculateur programmé pour mettre en oeuvre une telle méthode. Les réseaux de neurones trouvent leurs applications dans le traitement d'image, de la parole... Les réseaux de neurones sont constitués d'automates connectés entre eux par des synapses auxquelles sont affectés des coefficients synaptiques. Ils permettent de traiter des problèmes qui le sont difficilement par les ordinateurs séquentiels classiques. Pour effectuer un traitement donné, les réseaux de neurones doivent au préalable apprendre à l'effectuer. Cette phase appelée apprentissage fait appel à des exemples pour lesquels, à partir de données d'entrée, les résultats que l'on doit obtenir en sortie sont connus à l'avance. Dans un premier temps le réseau de neurones qui n'est pas encore adapté à.la tache prévue, va délivrer des résultats erronés. On détermine alors une erreur EP entre les résultats obtenus et ceux que l'on aurait dû obtenir et, à partir d'un critère d'adaptation, on modifie les coefficients synaptiques pour permettre au réseau de neurones d'apprendre l'exemple choisi. Cette étape est réitérée sur le lot d'exemples considérés comme nécessaires à un apprentissage satisfaisant du réseau de neurones. Une méthode très répandue pour opérer cette adaptation est celle de la rétropropagation du gradient. Pour cela on détermine les composantes du gradient gj,L de l'erreur Ep précédente (calculée sur la dernière couche L) par rapport à chaque état de neurone xjL. Ces composantes sont ensuite rétropropagées dans le réseau de neurones à partir de ses sorties afin de déterminer d'abord les composantes internes gjl (1L) puis les corrections à apporter aux coefficients synaptiques Wij, 1 des neurones os05 concernés. Cette méthode est par exemple décrite dans les documents: - D.E. Rumelhart, D.E. Hinton, and R.J. Williams, 'Learning Internal Representation by Error Propagation', in D.E. Rumelhart, and J.L. McClelland (Eds), 'Parallel Distributed Processing: Exploration in the Microstructure of Cognition', Vol.1, Foundations, MIT Press (1986). - 'Experiments on neural net recognition of spoken and written text.' D.J. Burr, IEEE Trans. on Acoustic, speech and signal processing Vol.36, n 7, July 1988 p.1162. Mais lorsqu'une telle méthode est mise en oeuvre dans un réseau de neurones, il peut arriver que pour certaines applications les temps d'apprentissage deviennent très longs. Ainsi par exemple dans le cas d'un problème de parité cette difficulté a été observée. Le problème de parité se présente par exemple dans le cas d'un réseau de neurones dont les entrées sont réunies à des signaux binaires 1/0 et dont la sortie doit délivrer un état 1 lorsque le nombre d'entrées à 1 est impair et un état O dans le cas inverse. La difficulté de l'apprentissage est due ici au fait que l'état de sortie doit changer lorsqu'une seule des entrées change d'état alors que lorsque c'est un nombre pair de changements d'état des entrées qui se produit, la sortie doit rester inchangée. De même, par exemple, lorsque le réseau de neurones est utilisé dans des problèmes de classification il peut être très difficile d'opérer la séparation des classes entre lesquelles la distance euclidienne minimale est petite car le réseau de neurones a besoin d'un temps très long pour apprendre à différencier ces différentes classes. Ceci gêne la séparation des données d'entrées codées en continu et en particulier lorsque certains des exemples, qui appartiennent & des classes différentes, ont des entrées qui diffèrent très peu les unes des autres. Le problème posé est donc de réduire les temps d'apprentissage du réseau de neurones tout en réduisant au maximum le matériel supplémentaire nécessaire. La solution consiste en une méthode d'apprentissage mise en oeuvre dans un réseau de neurones, formé de L couches, qui comprend les étapes de: détermination des états xj,1 des neurones d'une couche 1 à partir des potentiels de sortie Yi,l-1 délivrés par les neurones de la couche précédente qui leur sont connectés par des coefficients synaptiques Wij,1 ou à partir de données d'entrée yio pour la couche 1=1, tel que: Xj,1 = i Wij,l.Yi,1-1 détermination des potentiels Yj,1 des neurones de sortie par application d'une fonction non-linéaire F tel que: Yj,l = F(xj,l) avec: 1 indice de la couche considérée iel4L j indice du neurone de la couche de sortie 1 i: indice du neurone de la couche d'entrée 1-1 la méthode comprenant des phases d'apprentissage par itération à l'aide de P exemples qui sont introduits successivement sur les entrées du réseau de neurones avec: initialisation de la matrice de coefficients synaptiques Wij,1 du réseau de neurones, introduction des données d'entrée y de chaque exemple p destiné à l'apprentissage, comparaison des résultats obtenus Yj,L sur la couche de sortie L avec la sortie yjP prévue pour cet exemple p P présenté en entrée, pour définir une erreur partielle E, p détermination de la somme E de toutes les erreurs P partielles E observées pour chaque neurone de sortie J et pour chaque exemple p, détermination des différentes composantes du gradient g aE /3x de l'erreur E par rapport aux états ]nj,$ xjL pour la couche de sortie L, mise en oeuvre de la méthode de rétropropagation des composantes gj,L du gradient pour que le réseau de neurones détermine les composantes gjl1 du gradient pour les autres couches à partir de la matrice de coefficients synaptiques transposée, détermination des variations &xj,1 subséquentes qui ont un signe opposé à la composante correspondante gj,1, afin d'adapter le réseau de neurones mise à jour des coefficients synaptiques à partir de ces variations Axj1l, caractérisée en ce que pour déterminer les variations Axj,l subséquentes des états de neurones, la méthode comprend une étape de multiplication des composantes gj,l1 du gradient par des paramètres ej,l pour calculer des variations Axj,l proportionnelles à -$j,l.gj,1, ou ej,l1 dépend de l'état du neurone j de la couche 1, avec 6j,1=1 lorsque -gj,l1 et xj,l sont de signes différents, et 8j,l=$l+ lorsque -gj,1 et xj,l sont de même signe, avec 0601+41. Au cours de l'apprentissage on présente un exemple p donné. Les données propres à l'exemple p introduites sur les entrées du réseau de neurones vont délivrer sur la dernière couche L du réseau, un résultat Yj,L pour un neurone de sortie j donné. Or on connait au départ le résultat p y que l'on aurait d obtenir. On peut ainsi calculer pour un exemple une erreur tel que: p I{L)p - EP = % j= (y - y) j=1 j jL C'est l'expression du calcul d'une erreur quadratique moyenne. Un autre critère de comparaison peut étre utilisé. Selon la méthode connue de rétropropagation du gradient, on détermine les composantes gj,l du gradient de l'erreur pour chaque contribution d'un état xj,l1 d'un neurone tel que: g = úaE /ax ],l],1 o Xj,l représente l'état du neurone avant l'action de la fonction non linéaire. Pour cela on calcule les composantes g = 8E /ax se rapportant à la couche de sortie L tel que j,L jL g = (y - y) . F' o F' est la dérivée des jL jL) j,L j,L fonctions non linéaires de sortie. Le réseau de neurones est alors chargé avec la matrice de coefficients synaptiques transposée Wji,l1 et les composantes gj,L sont rétropropagées dans le réseau à partir de sa sortie. Le réseau détermine ainsi les autres composantes gj,l du gradient avec lsL. Ces composantes gj,l sont utilisées pour déterminer les variations Axj,l qui servent à corriger les coefficients synaptiques Wij,1 pour que le réseau s'adapte à l'exemple concerné. Généralement cette correction selon la méthode connue est effectuée tel que Wij,l(nouveau)=Wij,l(ancien)+k.axjl.Yi,l-1. Selon l'invention les composantes gj,l ne sont pas utilisées tel quel mais sont préalablement multipliées chacune par un paramètre propre ej,l1 qui dépend pour un neurone j donné du signe: - de l'état xj,l de ce neurone - et du signe de la composante du gradient gjl. Ces paramètres sont $j,l = 1 lorsque -gj,l et xj,l sont de signes différents, et 8j,l = el+ lorsque -gj,1 et xj,j sont de même signe avec 0 Mais pour accélérer l'apprentissage, préférentiellement au cours des premières itérations de l'apprentissage, pour chaque exemple donné, 31+ est choisi soit voisin de zéro soit égal à zéro. Il est possible complémentairement au cours des itérations ultérieures de l'apprentissage, pour chaque exemple donné, de faire croître 81+ vers la valeur 1. Avantageusement, la stratégie du signe selon l'invention permet au début de l'apprentissage de faire que les corrections opérées tiennent compte du signe de l'erreur observée, puis au fur et à mesure que l'apprentissage progresse, cette stratégie permet de provoquer graduellement des corrections moins brutales ayant une plus grande précision. Les fonctions non linéaires qui opèrent pour déterminer les potentiels de sortie peuvent être choisies faiblement ou fortement non linéaires. Pour accroître l'efficacité de la stratégie du signe selon l'invention il est possible de modifier en cours de l'apprentissage le choix de ces fonctions non linéaires. Mais les variations Axj,l qui se déduisent de la méthode de la rétropropagation du gradient ne doivent pas conduire à des variations trop importantes des coefficients synaptiques. Aussi selon un mode complémentaire de l'invention on effectue une normalisation tel que la somme des carrés des coefficients synaptiques reste quasiment constante. Ainsi les fonctions non linéaires sont choisies faiblement non linéaires en début d'apprentissage puis s'approchent de fonctions de type signe en fin d'apprentissage, et pour permettre ce choix, les coefficients synaptiques qui convergent vers un neurone donné j gardent une norme E (W)z2 quasiment constante. i 1j,1 Ces fonctions non linéaires F sont par exemple du type Yj, l=tanh(xj,l/T1) o T1 est un paramètre relatif à la couche 1 appelé température de la couche 1. Les variations apportées au cours de l'apprentissage au degré de non linéarité des fonctions non linéaires sont obtenues en faisant varier, pour chaque couche, les paramètres T1. La stratégie du signe selon l'invention, qui consiste d'abord à privilégier des corrections brutales (8+ petit et positif) basées sur le signe de l'erreur puis à l'adoucir par les paramètres e+ voisins de l'unité afin d'opérer des corrections ayant une plus grande précision, peut donner lieu à une action similaire au niveau de la structure P globale. Pour cela pour chaque exemple chaque composante E J p de l'erreur E est elle-mème multipliee par un paramètre Oj,L. Ainsi des corrections (paramètre nj,L) qui concernent simultanément la totalité des états neuronaux peuvent être superposées à chacune des corrections individuelles effectuées pour chaque état (paramètre e+). Pour favoriser la stratégie du signe présentée précédemment, on introduit un facteur de correction qj,L qui dépend de chaque neurone de sortie j de la dernière couche L. L'erreur Ep est alors déterminée telle que: p p E = En..E jL j p p avec E = (y - y)2 j j j,L Cette erreur est une fonction quadratique. Dans le cas général, pour chaque neurone de sortie j de la couche L (pour un exemple p donné) cette erreur p E est: J p p E =H (y - y) j] j,L o R est une fonction de l'écart entre le résultat obtenu Yj,L et le résultat prévisible yjP. Cette erreur Ep ainsi déterminée est utilisée pour déterminer les composantes gj,L et gj,l (avec l*L) du gradient comme cela a été développé précédemment. Ainsi pour déterminer ces composantes du gradient la méthode comprend une étape de détermination de l'erreur EP en appliquant préalablement un facteur de correction qj,L qui dépend du neurone j de la dernière couche L tel que E úP E IL jL. E. j=1 j afin de favoriser le début de l'apprentissage, p avec nj,L=1 si y et Yj,L sont de signes différents p et rj,L=O+si Y. et Yj,L sont de même signe avec Oen+(1. Dans le cas courant qj,L = ej,L Préférentiellement chaque erreur partielle p p 2 E est l'erreur quadratique (y - y) J 3 j,L Selon un mode subsidiaire la stratégie peut être développée au niveau de chaque couche du réseau de neurones.. Compte-tenu du rôle prépondérant dévolu aux couches d'entrée, il y a lieu d'accélérer l'apprentissage des couches d'entrée et de ralentir l'apprentissage des couches de sortie. Dans l'application classique de la méthode de rétropropagation du gradient l'état d'un neurone xj,l est modifié d'une quantité -Axj,1 qui tient compte de la composante du gradient gj,1 correspondante. Ceci s'effectue en multipliant chaque composante du gradient gj,1 par une constante de proportionnalité qui est la même pour toutes les couche.s du réseau de neurones. Selon le mode subsidiaire, l'invention propose d'opérer la correction en affectant à chaque neurone de chaque couche un coefficient de proportionnalité Bj,l de sorte que chaque correction -Axj,l soit proportionnelle à 1j,l.gj,1. Le paramètre Dj,l est relié à la stratégie du signe précédemment décrite en le prenant proportionnel au paramètre Oj, 1 servant à déterminer les corrections xj,. l Ainsi Bj,l est proportionnel à 8l.ej,j o B1 est un paramètre identique pour toute une couche 1 donnée. Selon ce mode subsidiaire, on affecte à chaque couche 1 un paramètre 01 qui lui est propre et qui permet de contrôler la vitesse d'apprentissage sur les couches d'entrée par rapport à celle sur les couches de sortie. Pour cela le paramètre B1 décroit quand 1 croit de la couche d'entrée vers la couche de sortie. Ainsi l'invention comprend une étape de multiplication des composantes ej, l.gj,l par des constantes 1 qui dépendent de chaque couche, pour que 8xjl soit alors proportionnel à l.Oj,l.gj,1 avec B1 strictement décroissante selon le numéro de la couche, en allant des couches d'entrée vers les couches de sortie, afin que les corrections appliquées aux états de neurones permettent d'accélérer l'apprentissage sur les couches d'entrée et de ralentir l'apprentissage sur les couches de sortie. L'invention sera mieux comprise à l'aide des figures suivantes données à titre d'exemple non limitatif qui représentent: Figure 1: un schéma indiquant les mécanismes du traitement opéré par une structure ayant une couche de neurones d'entrée et un neurone de sortie. Figure 2: un schéma représentant une structure à plusieurs couches: couche d'entrée, couches cachées, couche de sortie. Figure 3: un schéma représentant une structure de réseau de neurones mettant en oeuvre la méthode selon l'invention. La figure 1 représente le schéma classique d'un traitement opéré sur un réseau de neurones élémentaire formé d'une couche d'entrée à plusieurs neurones 10i...10I(l_) qui délivrent respectivement les signaux d'entrée Y,1-1, Y2,l-1...YI(l-1),l-1 à un seul neurone de sortie dont l'état est xj,. Cet état est déterminé par des moyens de calcul 11 tel que: xj,l = [ Wijl Yi,l-1 Cet état xj,l est soumis à l'action d'une fonction non linéaire (bloc 12) pour délivrer le potentiel de sortie Yj,l après application de cette fonction F. Yj,l = F(Xj,1) Ce potentiel de sortie Yj, l peut alors servir d'état d'entrée pour une couche suivante. On réalise ainsi un empilement en couches représenté sur la figure 2 avec une couche d'entrée 1=1, des couches cachées 1=2, 3, et une couche de sortie l=L. Les neurones d'une couche sont uniquement reliés à ceux de la couche suivante selon des coefficients synaptiques Wijl. Chaque état de neurones est déterminé cv q t L J l selon les expressions précédentes en commençant à la couche 1=1. Pour effectuer l'apprentissage, c'est-à-dire adapter les coefficients synaptiques Wij,l à une tâche donnée, on présente en entrée des exemples pour lesquels les résultats yjP sur la couche de sortie sont connus à l'avance. L'erreur Ep est calculée pour chaque exemple sur tous les états de sortie, puis on détermine ses variations relatives par rapport à des faibles variations bxj,1 de chaque état intermédiaire. Les composantes gj,1 du gradient sont telles que: j,1 E/a,1 Pour cela les composantes gjL sur la couche de sortie sont calculées puis sont rétropropagées dans le réseau de neurones qui restitue les autres composantes 9j,1 du gradient d'erreur. Ces composantes permettent de déterminer les variations axjl,1 qui s'en déduisent pour les états xj,1l pour que le réseau de neurones s'adapte à la tAche prévue. Ceci précède la mise à jour des coefficients synaptiques Wij,1 comme cela a été indiqué précédemment. Ces étapes de la méthode sont réalisées dans une structure de réseau de neurones dédiée tel que cela est représenté sur la figure 3 ou dans un calculateur programmé pour la mise en oeuvre de la méthode. La mémoire 30 stocke la matrice de coefficients synaptiques Wij,1, et la matrice transposée Wjil issues initialement des moyens d'entrée 29. Les coefficients synaptiques sont transmis à des moyens de calcul 31 qui reçoivent les potentiels d'entrée Yi,l-1 issus de la couche précédente. Ces moyens 31 déterminent: xj,l = E Wijl Yi,l-1 Au départ les états de neurones d'entrée Yi,l-1 sont des exemples YiO appliqués en entrée du réseau. Ces exemples sont issus d'une mémoire d'exemples 32. Un sélecteur 33 permet d'effectuer cette sélection. La mémoire d'exemples 32 stocke également les résultats yjP que l'on doit obtenir pour chaque exemple p et pour chaque potentiel de sortie j. Les états xj,1 des neurones de sortie sont soumis à l'action d'une fonction non linéaire dans l'organe 34 qui délivre pour chaque exemple les potentiels de sortie Yj,L de la dernière couche L tels qu'ils sont délivrés par le système. Pour opérer les étapes intermédiaires du calcul de couches en couches, les potentiels de sortie Yj,l1 d'une couche 1 sont stockés temporairement dans une mémoire d'état 37 et utilisés comme états d'entrée pour la couche suivante. Chaque potentiel YjL est comparé avec l'état prévu yjP dans le comparateur 35 qui, de plus, P stocke toutes les erreurs détectées E et les additionne pour ], délivrer une erreur E relative à chaque exemple. Les composantes du gradient gjL sont déterminées par le calculateur hôte 36. Pour cela il reçoit l'erreur EP, les potentiels de sortie Yj,L et les états p prévus y. Le calculateur hôte 36 détermine les composantes ] gj,L tel que: g = e.(y - y) . F' j,L j,L j,L j j,L avec lj la couche de sortie. Ces composantes gj,L sont introduites dans les moyens de calcul 31 qui permettent la mise en oeuvre de la méthode de la rétropropagation du gradient, c'est-à-dire que ces composantes gjL sont introduites sur les couches de sortie et leur action est rétropropagée sur les couches d'entrée. Les composantes du gradient gj,l=aEP/axj,1 (avec lsL) sont ainsi déterminées par les moyens de calcul 31 par rétropropagation du gradient de l'erreur EP. Les composantes gj,l sont délivrées au calculateur hôte 36 qui détermine les variations Axjl subséquentes pour chaque état de neurones. Pour cela il multiplie chaque composante gjl par son paramètre ejl selon l'invention. Toutes les variations àxj,l sont transmises à l'organe de mise à jour 38 qui détermine les nouveaux coefficients synaptiques Wij,1 et les délivre à la mémoire 30. Ce processus est réitéré pour effectuer la totalité de la phase d'apprentissage. Au cours de celle-ci le calculateur hôte 36 va pouvoir appliquer un paramètre de correction 81+ voisin ou égal à zéro pour les premières itérations puis le faire croître pour approcher la valeur 1 au cours des itérations ultérieures. De même le calculateur hôte p 36 va opérer la multiplication de E par les paramètres J Oj,L avant de calculer les composantes gj,L pour effectuer le rétropropagation du gradient dans les moyens de calcul 31. Lorsque les constantes Al propres à chaque couche sont appliquées aux corrections 8j l.gj, pour déterminer des variations -Axj,l proportionnelles à Bl.$j,l.gj,l, ceci est effectué par le calculateur hôte 36 préalablement à la remise à jour des coefficients synaptiques Wijl par l'organe de mise à jour 38. Ainsi la structure de réseau de neurones en couches selon l'invention comprend des moyens pour mettre en oeuvre la méthode d'apprentissage décrite et pour cela est munie: - de moyens pour stocker les coefficients synoptiques, - de moyens pour stocker des exemples à apprendre puis les introduire dans le réseau de neurones, - de moyens pour comparer pour chaque exemple les potentiels de neurones obtenus en sortie avec les résultats prévus pour chaque exemple et pour délivrer une erreur selon les écarts observés, - de moyens pour calculer des états de neurones de sortie à partir des potentiels des neurones d'entrée et pour effectuer la rétropropagation du gradient de ladite erreur et délivrer les composantes gj,l dudit gradient, - de moyens pour appliquer en sortie des fonctions non linéaires, - de moyens pour calculer des nouveaux coefficients synaptiques en tenant compte des composantes gjl du gradient et de paramètres multiplicateurs propres à la méthode qui permettent de contrôler l'importance attribuée à certaines itérations du cycle des itérations ou l'importance attribuée à certaines couches ou à certains neurones du réseau de neurones. Le schéma de la figure 3 est présenté sous la forme d'une structure de réseaux de neurones formée de blocs fonctionnels gérés par un calculateur hôte. Il est possible que l'ensemble des fonctions à réaliser soit intégré au calculateur lui-même. Dans ce cas l'invention concerne également un calculateur qui est programmé pour mettre en oeuvre les étapes de la méthode décrite. Le tableau 1 représente l'organigramme des principales étapes d'un exemple de programme selon l'invention. - L'étape 1 initialise n+ et e+ à de faibles valeurs positives et fixe les températures T1. Pour la couche 1=1 la valeur T1 est voisine de la moyenne des valeurs absolues des entrées relatives à l'exemple p et pour 1,1, la valeur T1 est de l'ordre de 1 (boucle sur 1). Les coefficients synaptiques Wij,l sont initialisés par un choix aléatoire ou à des valeurs connues (boucle sur i et j). - L'étape 2 introduit dans le réseau de neurones les valeurs des entrées Yio pour un exemple p. - L'étape 3 calcule des états xj,1 et des potentiels de sortie Yj,l. Le calcul des états xj,l peut faire intervenir un seuil sj,l qui peut également être introduit dans la fonction non linéaire F. - L'étape 4 applique la stratégie du signe sur l'erreur de p sortie. Pour cela on effectue le produit y. y et on ]],L considère son signe. Si le produit est négatif ou nul: nj,L prend la valeur 1. Dans le cas inverse qjL prend la valeur r+ L'erreur EP sur la couche de sortie est déterminée et les composantes du gradient gj,L sont calculées. - L'étape 5 Les dérivées F'j,l1 des fonctions non linéaires sont calculées. Les composantes du gradient gj,j-1 sont ensuite calculées par rétropropagation du gradient. On teste le produit -gj,l.xjj. Lorsque ce produit est négatif ou nul, ejL est mis égal à 1. Lorsque ce produit est positif, ej,L est égal à 8+ avec O On calcule ensuite Bj,l. - L'étape 6 Les composantes du gradient gj,l sont utilisées pour déterminer les variations bxj,l subséquentes. Cette étape donne un choix d'un exemple de fonction auto-adaptative qui permet de répercuter les composantes gj,1 sur les variations axjl. Cette fonction fait intervenir le module GZ du gradient gjl1, des facteurs v,Z qui règlent l'amplitude de la correction et la valeur moyenne B des termes 0jl associés aux différents neurones. - L'étae 7 Elle permet de répartir les variations axj,j calculées pour l'exemple p, à la fois sur les coefficients synaptiques Wij,l1 et sur les seuils sjl. Le taux de répartition est contrôlé par un paramètre ol qui fait intervenir la norme i Y i1,1-1 Cette étape 7 représente un exemple de taux de répartition qui permet de. garder quasiment constante la norme des coefficients synaptiques pour un neurone de sortie donné. On impose que les changements soient effectués avec des modifications des poids et des seuils aussi faibles que possibles. - L'étape 8 Lorsque l'erreur cumulée = E p=1 pour tous les exemples est inférieure ou égale à une valeur prédéterminée z, l'apprentissage est achevé. Lorsque cette erreur est supérieure à a, on continue selon les étapes suivantes. - L'étae 9 Les températures T1 sont légèrement abaissées et pour cela la valeur initiale est multipliée par un paramètre Qg1 compris entre O et 1. - L'tape 10 Les valeurs de n+ et 8+ sont réajustées. - L'étape 11 Un autre exemple p' est sélectionné et le traitement recommence à l'étape 2. TABLEAU 1 1. Init O( i,p 1= 2,.... L Tii oi=,.. I(1-1); j=1,.....I(l) L-Cwij,1 2. Ex(p) - Yio; i=1...,I(O) = l 3. I1=1. x =(11 W y j,l i=1 ij,l-l i,l-1 j,1 -4 = F(x) = tanh(x)/T 1 jL ]l/T 4. -> l=L: =....I(L) _ nj,L=q+; On(+ .I(L) g.,n 8E /x =n. (Yj -y).aF(x)/8x ],L,L ],L,L j j]/L j,L l=L,.....; j=l,...,I(1-1) F' j, l=aF(xj, l)/axj, l l 1, IWl)F 1 j,11ji l,l-1 -gj,lXj, lI 1= i+j=+1 j1j111 j= r1j, 6. OGv1; v l; t:1 T: (2-v) T l=LI(l) 2 Gz-l I] =1 l 1 'gjil l=L I(l) 22l=L I(l) 2 (1E-- El ej= l ( 'gj 1)/(lEl É i 'gj 1 1=1 joljl,' =p j=l j,1, 1=1,.....L; j=l,.....I(l) ax E 3J'l'gJ'l..DTD: Ax= - t. úP. G2 +. (SG2) 7.1: l....L f -2 i=j!-l) 2 tY-11'i,; ^1-1 O ) 1-1 i,1-1 i j: i,.. .I(1) o.A sj,]=.-2 2 2 IIYIl-(x +s a + o ll 1-1 3,1 j,1 1-1 ls = s + as l S,1j,1j,+ y.,Ax ti-1 j,1 Wij,1- 2 22 1-1 I j, j, 1-1 W.(x+s).8x ij,1j,1 j,1j,1..DTD: - 2 2 2 IIYIl- (x +s j +o 1-1 i,1j,1 1-1 ij,1 iij,1l iij,1l P=pp ( EP úP) -E END p=> >O 9. 1=1,.....L Ti=pl.T1 (O REVENDICATIONS 1. Méthode d'apprentissage mise en oeuvre dans un réseau de neurones, formé de L couches, qui comprend les étapes de: 05. détermination des états xjl des neurones d'une couche 1 à partir des potentiels de sortie Yil-1 délivrés par les neurones de la couche précédente qui leur sont connectés par des coefficients synaptiques Wijl ou à partir de données d'entrée Yi,o pour la couche 1=1, tel que: xj, l= WijlYi,l-1 détermination des potentiels Yj.l des neurones de sortie par application d'une fonction non-linéaire F tel que: Yj,l = F(xj,l) avec: 1: indice de la couche considérée 11(L j indice du neurone de la couche de sortie 1 i: indice du neurone de la couche d'entrée 1-1 la méthode comprenant des phases d'apprentissage par itération à l'aide de P exemples qui sont introduits successivement sur les entrées du réseau de neurones avec: initialisation de la matrice de coefficients synaptiques Wijl du réseau de neurones, introduction des données d'entrée y de chaque exemple j,o p destiné.à l'apprentissage, comparaison des résultats obtenus Yj,L sur la couche de sortie L avec la sortie yjP prévue pour cet exemple p p présenté en entrée, pour définir une erreur partielle E, ) P détermination de la somme E de toutes les erreurs p partielles E observées pour chaque neurone de sortie j et pour chaque exemple p, détermination des différentes composantes du gradient g = 8E /8x de l'erreur E par rapport aux états j,L j,L xj,L pour la couche de sortie L, L b ' 0 L v I mise en oeuvre de la méthode de rétropropagation des composantes gj,L du gradient pour que le réseau de neurones détermine les composantes gj,1 du gradient pour les autres couches a partir de la matrice de coefficients synaptiques transposée, détermination des variations &xj,1 subséquentes qui ont un signe opposé à la composante correspondante gj,1, afin d'adapter le réseau de neurones mise à jour des coefficients synaptiques à partir de ces variations Axj,l caractérisée en ce que pour déterminer les variations &xj,1 subséquentes des états de neurones, la méthode comprend une étape de multiplication des composantes gj,1 du gradient par des paramètres j,l1 pour calculer des variations àxj, l proportionnelles à -ejl.gjl, o $j,1 dépend de l'état du neurone j de la couche 1, avec Sj,1=1 lorsque -gj,l et xj,1 sont. de signes différents, et ej,1=$1+ lorsque -gjl1 et xj.l1 sont de même signe, avec Oe1+(1. 2. Méthode selon la revendication 1 caractérisée en ce que, au cours des premières itérations de l'apprentissage, pour chaque exemple donné, 31+ est choisi soit voisin de zéro soit égal à zéro. 3. Méthode selon la revendication 2 caractérisée en ce que, au cours des itérations ultérieures de l'apprentissage, pour chaque exemple donné, 81+ croit vers la valeur 1. 4. Méthode selon une des revendications 1 à 3 caractérisée en ce que les fonctions non linéaires sont choisies faiblement non linéaires en début d'apprentissage puis s'approchent de fonctions de type signe en fin d'apprentissage, et pour permettre ce choix, les coefficients synaptiques qui convergent vers un neurone donné j gardent une norme Z (W)2 quasiment constante. 1 1],1 5. Méthode selon la revendication 4 caractérisée en ce que les fonctions non linéaires F sont du type 21 2648251 YD, = tanh(xj,l/T1) o T1 est un paramètre relatif à la couche 1 appelé température de la couche 1. 6. Méthode selon la revendication 5 caractérisée en ce que les variations apportées au cours de l'apprentissage au degré de non linéarité des fonctions non linéaires sont obtenues en faisant varier, pour chaque couche, les paramètres T1. 7. Méthode selon une des revendications 1 à 6 caractérisée en ce que, pour déterminer des composantes du gradient gjL, la méthode comprend une étape de détermination de l'erreur Ep en appliquant préalablement un facteur de correction n jL qui dépend du neurone j de la p p dernière couche L tel que E X nj'L- E afin de favoriser le début de l'apprentissage, paecs p avec qj,L=l si y et YjL sont de signes différents et nj,L=n+si y. et YjL sont de même signe J avec 04n+41. 8. Méthode selon la revendication 7 caractérisée en ce que nj,L=ej,L. 9. Méthode selon une des revendications 1 à 8 caractérisée en ce que l'erreur partielle E est l'erreur p 2 quadratique % (y - y) 10. Méthode selon une des revendications 1 à 9 caractérisée en ce qu'elle comprend une étape de multiplication des composantes ejll.gj,1 par des constantes B1 qui dépendent de chaque couche, pour que -axj,1 soit alors proportionnel à lj,.,1 avec B1 strictement décroissante selon le numéro de la couche, en allant des couches d'entrée vers les couches de sortie, afin que les corrections appliquées aux états de neurones permettent d'accélérer l'apprentissage sur les couches d'entrée et de ralentir l'apprentissage sur les couches de sortie. 11. Structure de réseau de neurones caractérisée en ce qu'elle comprend des moyens pour mettre en oeuvre la métho- de d'apprentissage selon une des revendications 1 à 10 et pour cela est munie: - de moyens pour stocker les coefficients synaptiques, de moyens pour stocker des exemples à apprendre puis les introduire dans le réseau de neurones, - de moyens pour comparer pour chaque exemple les potentiels de neurones obtenus en sortie avec les résultats prévus pour chaque exemple et pour délivrer une erreur selon les écarts observés, - de moyens pour calculer des états de neurones de sortie à partir des potentiels des neurones d'entrée et pour effectuer la rétropropagation du gradient de ladite erreur et délivrer les composantes qjl dudit gradient, - de moyens pour appliquer en sortie des fonctions non linéaires, - de moyens pour calculer des nouveaux coefficients synaptiques en tenant compte des composantes gj,l1 du gradient et de paramètres multiplicateurs propres à la méthode qui permettent de contrôler l'importance attribuée à certaines itérations du cycle des itérations ou l'importance attribuée à certaines couches ou à certains neurones du réseau de neurones. 12. Calculateur caractérisé en ce qu'il est programmé pour mettre en oeuvre la méthode d'apprentissage selon une des revendications 1 à 10.