Méthode (100) de surveillance du régime moteur d’une machine tournante, la méthode (100) comprenant les étapes suivantes : Construction (110) d’un spectrogramme (1100) brut à partir d’un signal vibratoire mesuré par un capteur de vibrations ; Débruitage (120) du spectrogramme (1100) brut pour obtenir un spectrogramme (1200) égalisé, l’étape de débruitage (120) comprenant : Une première sous-étape (121) de détermination d’un pied du spectrogramme (1100) brut, la première sous-étape (121) comprenant : Une première sous-sous-étape (121a) de construction d’une fonction de régression robuste aux pics du spectrogramme (1100) brut ; Une deuxième sous-sous-étape (121b) de détermination du pied du spectre à partir de la fonction de régression ; Une deuxième sous-étape (122) de détermination du spectrogramme (1200) égalisé à partir du pied du spectrogramme ; Séparation de sources (130) dans le spectrogramme (1200) égalisé par détermination d’un estimateur des fréquences d’intérêt ; Détermination (140) d’une vitesse de rotation instantanée du régime moteur à partir de l’estimateur des fréquences d’intérêt. Méthode (100) de surveillance du régime moteur d’une machine tournante selon la revendication précédente, la méthode (100) comprenant en outre une étape (170) de détection d’une anomalie de fonctionnement du moteur de la machine tournante en fonction de la vitesse de rotation instantanée. Méthode (100) de surveillance du régime moteur d’une machine tournante selon l’une des revendications précédentes dans laquelle le spectrogramme (1100) brut est construit par une Transformée de Fourier appliquée au signal vibratoire sur une pluralité de fenêtres temporelles successives et de courte durée pour pouvoir considérer le signal vibratoire stationnaire dans chaque fenêtre temporelle de la pluralité de fenêtres temporelles. Méthode (100) selon l’une des revendications précédentes, dans laquelle l’étape de séparation de sources (130) est répétée jusqu’à satisfaction d’un critère de convergence, ladite étape de séparation de sources (130) comprenant une première sous-étape (131) de détermination de fréquences d’intérêt en fonction d’amplitudes des fréquences d’intérêt, et une deuxième sous-étape (132) d’estimation des amplitudes des fréquences d’intérêt en fonction des fréquences d’intérêt, lesdites première et deuxième sous-étapes (131,132) étant alternativement exécutées à chaque itération de ladite étape de séparation de sources (130), et l’estimateur des fréquences d’intérêt étant en outre déterminé en fonction des fréquences d’intérêt et des amplitudes des fréquences d’intérêt. Méthode (100) selon la revendication précédente, dans laquelle l’étape de séparation de sources (130) est mise en œuvre par un algorithme Espérance-Maximisation, où la première sous-étape (131) de détermination de fréquences d’intérêt est la E-étape de l’algorithme Espérance-Maximisation, et où la deuxième sous-étape (132) d’estimation des amplitudes des fréquences d’intérêt est la M-étape de l’algorithme Espérance-Maximisation. Méthode (100) selon l’une des revendications 4 et 5, dans laquelle la sous-étape (131) de détermination de fréquences d’intérêt en fonction d’amplitudes des fréquences d’intérêt comprend une première sous-sous-étape (131a) de détermination d’une probabilité a posteriori sur les fréquences d’intérêt, et une deuxième sous-étape (131b) de détermination d’une espérance d’une log-vraisemblance des amplitudes des fréquences d’intérêt. Méthode (100) selon l’une des revendications 4 à6, dans laquelle la sous-étape (132) d’estimation des amplitudes des fréquences d’intérêt en fonction des fréquences d’intérêt comprend une première sous-sous-étape (132a) de détermination d’une probabilité a posteriori sur les amplitudes des fréquences d’intérêt, une deuxième sous-sous-étape (132b) de détermination d’une probabilité a posteriori sur les amplitudes des fréquences d’intérêt, et une troisième sous-sous-étape (132c) de détermination d’une variance d’une distribution Gaussienne généralisée lissée du spectrogramme (1200) égalisé. Méthode (100) selon la revendication 6 ou selon les revendications 6 et 7, dans laquelle l’estimateur des fréquences d’intérêt est déterminé, à la satisfaction d’un critère d’arrêt, en fonction de la probabilité a posteriori sur les fréquences d’intérêt. Méthode (100) selon l’une des revendications précédentes, la méthode (100) comprenant en outre une étape (150) d’interpolation de l’estimateur des fréquences d’intérêt. Méthode (100) selon la revendication précédente, la méthode (100) comprenant en outre une étape (160) de lissage de l’interpolation de l’estimateur des fréquences d’intérêt.