Procédé de mesure d’une cavitation érosive d’une machine hydraulique (1) en fonctionnement, la machine hydraulique (1) comportant des parois (13, 16, 17) délimitant un conduit (12) de passage d’eau liquide et au moins un organe mobile (11, 3, 9), qui est situé dans le conduit (12) de passage d’eau liquide et qui est mis en mouvement lors du fonctionnement de la machine hydraulique (1), caractérisé en ce que le procédé comporte les étapes suivantes : au cours d’une première étape (E1), émission d’au moins un premier signal ultrasonore (e) impulsionnel à large bande par au moins un transducteur ultrasonore (20) émetteur, qui est positionné sur un premier côté (131, 171) d’au moins une première (13, 17) des parois (13, 16, 17) et qui est orienté au travers du conduit (12) de passage d’eau liquide vers un autre côté (172) de la première paroi (17) ou vers le au moins un organe mobile (11, 3, 9), au cours d’une deuxième étape (E2), réception, par au moins un transducteur ultrasonore (30) récepteur, qui est positionné sur le premier côté (131, 171) de la première paroi (13, 17), d’au moins un deuxième signal ultrasonore (x) reçu en réponse au premier signal ultrasonore (e) ayant été émis, au cours d’une troisième étape (E3), calcul, par une unité (40) de calcul reliée au au moins un transducteur ultrasonore (20) émetteur et au au moins un transducteur ultrasonore (30) récepteur, à partir du premier signal ultrasonore (e) ayant été émis et à partir du deuxième signal ultrasonore (x) ayant été reçu, d’au moins un indicateur de cavitation érosive de la machine hydraulique (1), à partir du calcul par l’unité (40) de calcul d’au moins un temps (TOA n ) d’arrivée du deuxième signal ultrasonore (x) ayant été reçu par rapport au premier signal ultrasonore (e) ayant été émis et/ou à partir du calcul par l’unité (40) de calcul d’une bande spectrale (BW k ) instantanée définie par une plage de fréquence, pour laquelle une valeur représentative (DSP (k ) (f)) de la densité spectrale de puissance (DSP i (k ) (f)) du deuxième signal ultrasonore (x) ayant été reçu est supérieure à un premier pourcentage (A) non nul donné de la valeur moyenne (DSP mean ) calculée de la valeur représentative (DSP (k ) (f)) de la densité spectrale de puissance (DSP i (k ) (f)) au cours du temps. Procédé suivant la revendication 1, caractérisé en ce que le premier signal ultrasonore (e) impulsionnel à large bande comporte une pluralité d’impulsions (n) émises avec une période prescrite (T r ) de répétition des impulsions, non nulle. Procédé suivant la revendication 2, caractérisé en ce que chaque impulsion (n) a, pendant une durée (L p ) de modulation, inférieure à la période prescrite (T r ) de répétition des impulsions, une modulation de fréquence supérieure ou égale à une valeur basse (f 1 ) de modulation de fréquence, non nulle, et inférieure ou égale à une valeur haute (f 2 ) de modulation de fréquence, non nulle. Procédé suivant la revendication 3, caractérisé en ce que la modulation de fréquence est linéaire ou cubique et va de la valeur basse (f 1 ) de modulation de fréquence à la valeur haute (f 2 ) de modulation de fréquence sur la durée (L p ) de modulation ou de la valeur haute (f 2 ) de modulation de fréquence à la valeur basse (f 1 ) de modulation de fréquence sur la durée (L p ) de modulation. Procédé suivant l’une quelconque des revendications 2 à 4, caractérisé en ce que lors de la troisième étape (E3), l’unité (40) de calcul calcule pour chaque impulsion (n) du premier signal ultrasonore e(t) ayant été émis le temps TOA n d’arrivée du deuxième signal ultrasonore x(t) ayant été reçu par rapport au premier signal ultrasonore e(t) ayant été émis suivant les équations suivantes : où n est l’indice de chaque impulsion. Procédé suivant l’une quelconque des revendications 2 à 4, caractérisé en ce que lors de la troisième étape (E3), l’unité (40) de calcul calcule pour chaque impulsion (n) du premier signal ultrasonore e(t) ayant été émis le temps TOA n d’arrivée du deuxième signal ultrasonore x(t) ayant été reçu par rapport au premier signal ultrasonore e(t) ayant été émis suivant les équations suivantes : où n est l’indice de chaque impulsion, FFT{x} est la transformée de Fourier rapide du deuxième signal ultrasonore x(t) ayant été reçu, est le conjugué complexe de la transformée de Fourier rapide du premier signal ultrasonore e(t), est la transformée de Fourier rapide inverse. Procédé suivant la revendication 5 ou 6, caractérisé en ce que lors de la troisième étape (E3), l’unité (40) de calcul calcule le au moins un indicateur de cavitation érosive de la machine hydraulique (1) comme comprenant la moyenne ( ) des temps TOA n d’arrivée pour les impulsions (n) du premier signal ultrasonore e(t) ayant été émis et/ou la variance ( ) des temps TOA n d’arrivée pour les impulsions (n) du premier signal ultrasonore e(t) ayant été émis et/ou l’écart-type ( ) des temps TOA n d’arrivée pour les impulsions (n) du premier signal ultrasonore e(t) ayant été émis. Procédé suivant la revendication 7, caractérisé en ce que lors de la troisième étape (E3), l’unité (40) de calcul détermine qu’une cavitation (C ) est présente, lorsque la moyenne (m TOA ) des temps TOA n d’arrivée est inférieure à un premier seuil (S m ) déterminé. Procédé suivant la revendication 7 ou 8, caractérisé en ce que lors de la troisième étape (E3), l’unité (40) de calcul détermine qu’une cavitation (C ) est présente, lorsque la variance ( ) des temps TOA n d’arrivée est supérieure à un seuil déterminé (V, V C ). Procédé suivant l’une quelconque des revendications 7 à 9, caractérisé en ce que lors de la troisième étape (E3), lorsque l’unité (40) de calcul a déterminé qu’une cavitation (C ) est présente, l’unité (40) de calcul détermine une position (P C ) de la cavitation (C ) dans la région définie par le produit (c.m TOA /2) entre la vitesse (c ) de propagation des signaux ultrasonores (e, x) et de la moyenne ( ) des temps TOA n d’arrivée, divisée par 2. Procédé suivant l’une quelconque des revendications 7 à 10, caractérisé en ce que lors de la troisième étape (E3), lorsque l’unité (40) de calcul a déterminé qu’une cavitation (C ) est présente, l’unité (40) de calcul détermine une durée de vie (DV C ) de la cavitation (C ) à partir de la variance ( ) des temps TOA n d’arrivée et/ou à partir de la moyenne (m TOA ) des temps TOA n d’arrivée. Procédé suivant l’une quelconque des revendications précédentes, lorsqu’elles dépendent au moins de la revendication 3, caractérisé en ce que lors de la troisième étape (E3), l’unité (40) de calcul calcule: la densité spectrale de puissance DSP i (k ) (f) du deuxième signal ultrasonore x k (t) ayant été reçu pour chacune de plusieurs fenêtres temporelles F i de réception ayant une durée W prescrite couvrant l’impulsion n, les fenêtres temporelles F i de réception étant décalées temporellement l’une à la suite de l’autre pour i allant de 1 à M, selon l’équation suivante : où M est un nombre prescrit de fenêtres temporelles F i de réception, la valeur représentative DSP (k ) (f) de la densité spectrale de puissance DSP i (k ) (f) du deuxième signal ultrasonore x k (t) ayant été reçu étant égale à la densité spectrale moyenne DSP (k ) (f) de puissance pour la densité spectrale de puissance DSP i (k ) (f) sur les M fenêtres temporelles F i de réception, selon l’équation suivante : la valeur moyenne DSP mean calculée de la valeur représentative DSP (k ) (f) de la densité spectrale de puissance DSP i (k) au cours du temps selon l’équation suivante : la bande spectrale BW k instantanée définie par la plage de fréquence supp f , pour laquelle la valeur représentative DSP (k ) (f) de la densité spectrale de puissance DSP i (k ) (f) du deuxième signal ultrasonore x k (t) ayant été reçu est supérieure au premier pourcentage A non nul donné de la valeur moyenne DSP mean calculée de la valeur représentative DSP (k ) (f) de la densité spectrale de puissance DSP i (k ) (f) au cours du temps selon l’équation suivante : où f 1 est la valeur basse de modulation de fréquence, f 2 est la valeur haute de modulation de fréquence, f est la fréquence, t est le temps. Procédé suivant la revendication 12, caractérisé en ce que les fenêtres temporelles F i de réception se chevauchent l’une à la suite de l’autre. Procédé suivant la revendication 12 ou 13, caractérisé en ce que la durée W prescrite de la fenêtre temporelle F i de réception est supérieure ou égale à 1 fois ou 2 fois ou 3 fois la période prescrite (T r ) de répétition des impulsions. Procédé suivant l’une quelconque des revendications 12 à 14, caractérisé en ce que lors de la troisième étape (E3), lorsque l’unité (40) de calcul a déterminé qu’une cavitation (C ) est présente, l’unité (40) de calcul détermine une direction (D C ) de déplacement de la cavitation (C ) à partir de la bande spectrale instantanée (BW k ). Procédé suivant l’une quelconque des revendications précédentes, caractérisé en ce que le premier pourcentage A non nul donné est supérieur ou égal à 50 % ou à 70 % ou à 90 % ou à 95%. Procédé suivant l’une quelconque des revendications précédentes, caractérisé en ce que le au moins un transducteur ultrasonore (20) émetteur comporte un mode actif, pouvant être commandé par l’unité (40) de calcul et dans lequel la première étape (E1), la deuxième étape (E2) et la troisième étape (E3) sont effectuées, et un mode passif, pouvant être commandé par l’unité (40) de calcul et dans lequel la première étape (E1), la deuxième étape (E2) et la troisième étape (E3) ne sont pas effectuées et le au moins un transducteur ultrasonore (20) émetteur n’émet pas, l’unité (40) de calcul commandant par défaut lors de la troisième étape (E3) le au moins un transducteur ultrasonore (20) émetteur dans le mode passif, dans lequel le au moins un transducteur ultrasonore (30) récepteur reçoit un troisième signal (x y ), l’unité (40) de calcul calculant lors de la troisième étape (E3) la densité spectrale de puissance (DSP y (k ) (f), DSP’ y (k ) (f)) du troisième signal (x y ) ayant été reçu dans chacun de plusieurs créneaux temporels (k-1, k) se suivant et ayant chacun une durée identique (D) de créneau temporel, l’unité (40) de calcul commandant lors de la troisième étape (E3) le au moins un transducteur ultrasonore (20) émetteur dans le mode actif, lorsque la densité spectrale de puissance (DSP y (k ) (f), DSP’ y (k ) (f)) du troisième signal (x y ) ayant été reçu dans un créneau temporel actuel (k) diffère de plus d’un deuxième pourcentage (B) non nul donné par rapport à la densité spectrale de puissance (DSP y (k ) (f), DSP’ y (k ) (f)) du troisième signal (x y ) ayant été reçu dans un créneau temporel (k-1) antérieur au créneau temporel actuel (k). Procédé suivant la revendication 17, caractérisé en ce que le créneau temporel (k-1) antérieur au créneau temporel actuel (k) est le créneau temporel (k-1) précédant le créneau temporel actuel (k). Dispositif (100) de mesure d’une cavitation érosive d’une machine hydraulique (1) en fonctionnement, la machine hydraulique (1) comportant des parois (13, 16, 17) délimitant un conduit (12) de passage d’eau liquide et au moins un organe mobile (11, 3, 9), qui est situé dans le conduit (12) de passage d’eau liquide et qui est mis en mouvement lors du fonctionnement de la machine hydraulique (1), caractérisé en ce que le dispositif (100) de mesure comporte : au moins un transducteur ultrasonore (20) émetteur, qui est positionné sur un premier côté (131, 171) d’au moins une première (13, 17) des parois (13, 16, 17) et qui est orienté au travers du conduit (12) de passage d’eau liquide vers un autre côté (172) de la première paroi (17) ou vers le au moins un organe mobile (11, 3, 9), le au moins un transducteur ultrasonore (20) émetteur étant apte à émettre au moins un premier signal ultrasonore (e) impulsionnel à large bande, au moins un transducteur ultrasonore (30) récepteur, qui est positionné sur le premier côté (131, 171) de la première paroi (13, 17), et qui est apte à recevoir au moins un deuxième signal ultrasonore (x) reçu en réponse au premier signal ultrasonore (e) ayant été émis, une unité (40) de calcul reliée au au moins un transducteur ultrasonore (20) émetteur et au au moins un transducteur ultrasonore (30) récepteur, l’unité (40) de calcul étant apte à calculer à partir du premier signal ultrasonore (e) ayant été émis et à partir du deuxième signal ultrasonore (x) ayant été reçu, au moins un indicateur de cavitation érosive de la machine hydraulique (1), à partir du calcul par l’unité (40) de calcul d’au moins un temps (TOA n ) d’arrivée du deuxième signal ultrasonore (x) ayant été reçu par rapport au premier signal ultrasonore (e) ayant été émis et/ou à partir du calcul par l’unité (40) de calcul d’une bande spectrale (BW k ) instantanée définie par une plage de fréquence, pour laquelle une valeur représentative (DSP (k ) (f)) de la densité spectrale de puissance (DSP i (k ) (f)) du deuxième signal ultrasonore (x) ayant été reçu est supérieure à un premier pourcentage (A) non nul donné de la valeur moyenne (DSP mean ) calculée de la valeur représentative (DSP (k ) (f)) de la densité spectrale de puissance (DSP i (k ) (f)) au cours du temps. Programme d’ordinateur, comportant des instructions de code pour la mise en œuvre de la troisième étape (E3) du procédé de mesure d’une cavitation érosive d’une machine hydraulique (1) suivant l’une quelconque des revendications 1 à 18, lorsqu’il est mis en œuvre sur un calculateur (40).