Brevet FR2799321A1 — publié le 20010406

PROCÉDÉ DE CONTROLE EN CONTINU DE LA QUALITÉ DES SONS NUMÉRIQUES EN DISTRIBUTION L'invention concerne un procédé contrôle en continu de la qualité des sons numériques distribution. Les processus de codage audionumériques utilisés les services de radio- ou de télédiffusion ont permis réduire la quantité de données à transmettre. En contrepartie, cette réduction est susceptible d'entraîner perte irrémédiable de la qualité du son. par rapport au signal source d'origine. L'importance des défauts engendrés dépend à la fois du débit alloué au codeur, de la complexité du contenu du signal sonore, ainsi que des problèmes liés à la transmission du signal. Pour des raisons techniques ou de responsabilité de télédiffusion, il est nécessaire d'évaluer de façon continue le niveau de qualité du signal audio. Des méthodes d'évaluation subjectives d'équipements, par surveillance et appréciation humaine, sont d'une mise en #uvre lourde, et peu fiables. En particulier, parmi les inconvénients plus spécifiques des processus ou méthodes l'art antérieur, on peut citer la mise en ceuvre d'évaluations subjectives longes et coûteuses ; # l'absence de complétude des informations nécessaires pour effectuer le contrôle de la qualité sonore perçue, lorsque ces informations sont fournies par des analyseurs de flux binaire ; # l'absence d'analyse objective du contenu sonore, seule capable de refléter la qual-*té finale des signaux sonores perçus ; les défauts inhérents à l'analyse différentielle tels que - mise à disposition de la source non codée, comme source de référence ; - séquences analysées de courte durée, 20 secondes au plus, lesquelles ne sont pas représentatives du service analysé ; - transparence de certains défauts à ce type d'analyse ; - analyse généralement discontinue et non totalement significative. particulier, les processus d'analyse différentielle, basés sur le système d'audition humain, entre une source sonore de référence et la source sonore à évaluer, peuvent permettre une mise en ceuvre automatique. Toutefois cette solution apparait peu pratique car il est nécessaire de disposer de la source sonore de référence. La présente invention a pour objet de remédier aux inconvénients précités des processus ou méthodes l'art antérieur par la mise en #uvre d'un procédé basé une étude fine du signal numérique et du comportement en continu de ce dernier, afin de permettre, à partir des méthodes classiques, d'apprécier le niveau de qualité globale du signal. Le procédé de contrôle en continu de la qualité sons en distribution, objet de la présente invention, sons numériques étant disponibles en mode stéréophonique selon un signal numérique représentant au :-.oins une voie drcite et une voie gauche, est remarquable en ce qu'il consiste à effectuer une analyse statistique du contenu de ce signai numérique sur chacune de ces voies. L'analyse statistique consiste découper temporellement le signal nurr.érique selon. des suites successives d'échantillons, comportant un nombre déterminé d'échantillons, et à effectuer, lors de l'existence d'un programme de sons numériques, une analyse spectrale de chacune des suites d'échantillons pour observer les variations d'énergie et d'enveloppe du signal numérique dans le domaine temporel et fréquentiel et calculer un. indice qualité global. Un indice de quali " final est calculé ' partir des variations d'énergie et d'envelcppe et de 1 indice de qualité global, sous -`orme d'une valeur bornée continue dans le temps, cet indice de qualité final étant représentatif de la qualité des sors numériques perçus. procédé, objet de la présente invention trouve application à la surveillance opérationnelle et continue des composantes sonores des services audio- et audiovisuels, avant et après distribution secondaire notamment aux prestations d'expertise d'équipement, codeurs et multiplexeurs, d'expertise de la qualité de service, plates-formes d'expérimentation. Ce procédé, objet de la présente invention, sera mieux compris à la lecture de la description et à l'observation des dessins ci-après dans desquels - la figure la représente, sous forme de schéma synoptique, un organigramme général du procédé de contrôle en continu de la qualité des sors numériques en distribution, objet de - présente invention ; - la figure lb représente, à purement illustratif, un processus de création de suites d'échantillons du signal numérique, permettant la mise en oeuvre du procédé objet de la présente invention ; - la figure 2 représente, sous fo-T.e d'un organigramme, un détail d'un mode de réalisation préférentiel de l'étape de calcul à partir des variations d'énergie et d'enveloppe de indice de qualité final ; - la figure 3a représente un organigra=e 'atif à un mode de calcul préférentiel non lirr,itatif d'une valeur Cb(t) liée à la bande passante u signal numérique et permettant la mise en ceuvre d,.i mode de réalisation préférentiel du procédé objet de la présente invention représenté en figure 2 ; - la figure 3b représente un organigramme re;atif à un mode de calcul préférentiel non limitatif d'une valeur Cs(t) liée aux propriétés stéréophoniques du signal numérique temporel et permettant _a mise en oeuvre du mode de réalisation préférentiel procédé objet de la présente invention représenté figure - figure 3c représente un organigramme relatif à un mode de calcul préférentiel non limitatif d'une valeur Cw(t) liée au blanchiment du signal numérique temporel pour chaque voie du signal numérique temporel et permettant la mise en oeuvre du procédé objet de la présente invention représenté en figure - figures 4a et 4b représentent un processus détection d'un signal de brève coupure ; - les figures 5a et 5b représentent un processus détection d'un signal parasite de sifflement ; - figures 6a et 6b représentent un processus détection d'un signal parasite de bourdonnement ; - figure 7 représente un processus de détection déphasage intervoies entre les signaux numériques véhiculés par les voies d'un signa stéréophonique. Une description plus détaillée du procédé contrôle en continu de la qualité des sons numériques distribution, objet de la présente invention, sera maintenant donnée en liaison avec les figures la, lb et les figures suivantes. D'une manière générale, on indique que le procédé objet de la présente invention permet d'obtenir une valeur d'indice de qualité bornée, s'échelonnant par exemple entre deux limites supérieures de qualité, excellent à mauvais, cette valeur bornée étant continue dans le temps et indicative de la qualité du signal sonore. Par valeur continue dans le temps, on comprend bien entendu que cette valeur consiste en fait en des valeurs discrètes successives calculées sur des intervalles de temps suffisamment courts pour que ces valeurs successives soient représentatives d'une valeur de qualité considérée comme continue dans le temps. Ainsi qu'on l'a représenté en figure la, le procédé objet de l'invention s'applique à des sons numériques, lesquels sont disponibles en mode stéréophonique selon un signal numérique, noté ADS représentant au moins une voie droite et une voie gauche, le procédé objet de la présente invention pouvant le cas échéant être appliqué à des signaux sonores quadriphoniques ou autres. D'une manière générale, le procédé objet de présente invention consiste à effectuer une analyse statistique du contenu du signal numérique précité chacune des voies. En référence à la figure la, l'analyse statistique peut consister, en une étape 1, à découper temporellement le signal numérique selon des suites successives d'échantillons, SE" comportant un nombre aéterminé d'échantillons puis, à une étape 2, à effectuer analyse spectrale de chacune des suites d'échantillons pour observer les variations d'énergie, notées 4W, enveloppe, notées AE, du signal numérique dans domaine temporel et fréquentiel et calculer un indice qualité global I(t)=f(AW,AE) à partir des variations énergie et d'enveloppe. Les étapes précitées sont suivies d'une étape consistant à calculer, à partir des variations d'énergie et d'enveloppe et de l'indice de qualité global I(t), indice de qualité final, noté If(t), lequel est constitué par une valeur bornée et continue dans le temps. Cet indice est représentatif de la qualité des signaux numériques précités. En ce qui concerne l'étape 1 de découpage temporel, on indique que les suites d'échantillons peuvent etre constituées par des suites d'échantillons présentant taux de recouvrement rapport du nombre d'échantillons communs à deux suites consécutives S.,__, S, au nombre d'échantillons constitutif de chaque suite d'échantill'On s, ce taux pouvant être compris entre 0 et 75%. On indique en particulier que le découpage temporel précité peut être effectué par une mémorisation séqÜentielie de ces suites d'échantillons puis une relecture échantillon par échantillon mémorisés, le processus de relecture étant réalisé par adressage en recouvrement échantillons successifs pour réaliser le taux recouvrement considéré. Sur la figure lb, on a représenté à titre illustratif les suites successives d'échantillons, les suites successives Sn_1, Sn et Sn_1 étant en recouvrement de deux échantillons sur cent par exemple. Une description plus détaillée étapes 2 d'analyse spectrale de variations d'énergie d'enveloppe et de calcul d'un indice de qualité global et étape 3 de calcul d'indice de qualité final à partir des variations d'énergie et d'enveloppe AW et DE et du facteur de qualité global I (t) sera maintenant donnée en liaison avec la figure 2. D'une manière générale, on indique que l'étape 2 précitée, selon la figure la, consiste à calculer un indice de qualité global I(t) à partir d'au moins un critère fréquentiel et d'un critère temporel de variation d'énergie et d'enveloppe. En référence à la figure 2 précitée, l'étape 2 peut comprendre une étape 20 de détection de l'existence d'un programme radio- ou télédiffusé dans le signal numérique. Sur réponse négative à l'étape 20 précitée, une valeur arbitraire est allouée à l'indice de qualité final If(t)=1 à l'étape 21, la qualité en l'absence de programme étant réputée excellente. Au contraire, sur réponse positive à 'L'étape 20 précitée, l'étape 2 précédemment mentionnée consiste à prendre en compte les critères de qualité liés aux variations d'énergie OW et d'enveloppe AE, ces critères pouvant consister en le calcul de valeurs telles que valeurs Cb(t) liées à la bande passante du signal numérique, valeurs cs(t) liées aux propriétés stéréophoniques du signal numériques enfin, valeurs Cw(t) basées sur le blanchiment du signal temporel. étape 22 précitée est alors suivie d'une étape 23 consistant à calculer la valeur de 1 indice de qualité global, lequel est défini par une co7.binaison linéaire des valeurs (t), Cs(t) et Cw(t). A titre d'exemple non limitatif, l'indice de qualité global vérifie la relation (1) . La valeur de l'indice de qual' global ainsi obtenue pour une suite d'échantillons considérée est comprise entre 0, pour une qualité globale mauvaise, et 1 pour une qualité globale excellente. Suite à l'étape 2 précitée, l'étape 3 de calcul d'indice de qualité final peut alors être misse en #uvre ainsi que représenté dans le :rode réalisation préférentiel non limitatif de la figure 2. D'une manière générale, l'étape 3 consiste à pondérer la valeur de l'indice de qualité global I(t) en fonction de l'apparition de signaux de défaut susceptibles de gêner l'audition des signaux sonores, ces défauts constituant alarmes susceptibles d'inciter l'opérateur à prendre mesures pour assurer la qualité de la radio ou de la télédiffusion. D' manière générale, on indique que les signaux de défaut ou alarmes retenus sort les suivants - le sifflement ou la saturation, - le phénomène micro-coupure, - le bourdonnement, - le déphasage inter-voles. En ce qui concerne l'absence de programme, on rappelle que cette situation est réglée par l'étape 20 de l'étape 2 précédemment mentionnée dans la description. Ainsi sur la figure 2, on a représenté, dans un mode de réalisation préférentiel non limitatif, l'étape 3 comme consistant à détecter l'existence sur le signal numérique ADS d'au moins une perturbation de transmission 'u signal numérique, cette pert,,:rbation de transmission étant détectée à l'étape 3C pour l'existence d'un sifflement d'une saturation, à l'étape 31 pour l'existence d'un phénomène de à l'étape 32 pour l'existence d'un bourdonnement. Outre la détection de l'existence d'au moins une perturbation transmission du signal numérique aux étapes 30, 31 32 précitées, le procédé objet de la présente invention peut consister, pour la mise en ceuvre de l'étape à détecter la présence d'un déphasage inter voies à étape 33, la présence d'un tel déphasage n'étant toutefois pas considérée comme une perturbation de transmission en raison de déphasages relatifs introduits, dans certains cas, par les opérateurs sur la voie de gauche, respectivement la voie de droite des signaux audionumériques. Suite à la détection d'au moins une perturbation de transmission du signal numérique aux étapes 30, 31 32 précitées, le procédé objet de la présente invention consiste à affecter à l'existence de cette perturbation coefficient de pondération spécifique représentatif de contribution de cette perturbation à la dégradation de qualité des signaux numériques. Ainsi, en référence à la figure 2, pour la mise oeuvre de l'étape 3, on indique que sur réponse positive l'étape 30 de détection d'un sifflement ou d' saturation, un coefficient ps supérieur à 1 est affecté au phénomène de sifflement ou de saturation à l'étape , alors que sur réponse négative à l'étape 30, coefficient de pondération ps=l est alloué à l'étape 30b à ce meure phénomène de sifflement ou saturation. I1 en est de même pour le phénomène de micro coupure à l'étape 31 pour lequel, sur réponse positive, c'est-à-dire lors de l'existence d'une micro-coupure, un coefficient de pondération pur supérieur à 1 est alloué au phénomène précité à l'étape 31a, alors que sur réponse négative en l'absence de micro-coupure, un coefficient de pondération pm=l est affecté à ce même phénomène à l'étape 31b. De la même manière, pour le phénomène de bourdonnement à l'étape 32, sur réponse positive à l'étape de détection du bourdonnement précité, un coefficient pondération Pb supérieur à 1 est alloué au bourdonnement et un coefficient de pondération pb=l est alloué bourdonnement sur réponse négative à l'existence de ce phénomène à l'étape 32b. Compte tenu de la valeur des coefficients de pondération ps, pm et Pb affectés aux signaux de perturbation ou d'alarme de sifflement ou saturation, micro-coupure ou de bourdonnement, un coefficient pondération global, produit des coefficients pondération affectés à chacun des signaux de perturbation précités, est calculé à l'étape 34, lequel vérifie la relation (2) . P = Ps x Pm x Pb Ainsi que représenté en outre en figure 2, suite à 1a détection sur le signal numérique ADS d'un déphasage de valeur d à l'étape 33, ce déphasage correspondant à un déphasage inter-voies, le procédé objet de la présente invention consiste à affecter une valeur de critère déphasage D à cette valeur de déphasage lorsque cette valeur de déphasage est supérieure à 0, c'est-à-dire réponse positive au test 33, et une valeur de critère de déphasage D égale à 0 sinon à l'étape 33b, c'est-à-dire sur réponse négative au test 33. A titre d'exemple non limitatif, on indique que pour l'existence d'un déphasage détecté à l'étape 33, valeur de critère de déphasage peut avoir la valeur D = d/170 et D = 0 sinon, la valeur de d étant exprimée millisecondes par exemple. L'étape 34 est alors suivie d'une étape 35 consistant à calculer et déterminer l'indice de qualité final If(t) par comparaison de la différence entre l'indice de qualité pondéré, cet indice de qualité pondéré prenant la valeur de l'indice de qualité global divisé _e coefficient de pondération. p obtenu à l'étape 34, et va.eur du critère de déphasage D attribuée à ''étape 33a ou cette différence étant alors comparée à valeur Ainsi, pour attribuer l'indice de qualité final l'étape 35, celui-ci, en présence d'un programme radio télédiffusé, vérifie la relation (3) . I f (t) = sup(I(t) / p - D,0) . La relation (3) indique qu'à l'indice de qualité final attribuée la valeur la plus grande entre valeurs constituées par la différence précitée et la valeur 0. En ce qui concerne la valeur des coefficients de pondération, des essais ont montré que - s il y a détection de sifflement ou de saturation = 1,75 et p, =1 sinon ; - s'il y a détection d'une micro-coupure . pr, = 1,5 et p = 1 sinon ; - s il y a détection d'un bourdonnement Pb =l,25 et = 1 sinon ; - s'il y a un déphasage de valeur d en ms, alors D = d/170 et D = 0 sinon. On indique que la relation (3) réalisée à l'étape 35 est utilisée, puisque par hypothèse l'indice de qualité final ne peut pas avoir de valeur négative. Une description plus détaillée des processus de calcul des valeurs Cb(t) liées à la bande passante, Cs(t) liées aux propriétés stéréophoniques du signal numérique temporel et Cw(t) liées au blanchiment du signal numérique temporel, processus mis en #uvre à l'étape 22 représentée en figure 2, sera maintenant donnée en liaison avec les figures 3a, 3b, 3c. En référence à la figure 3a, l'étape de calcul de la valeur Cb(t) liée à la bande passante du signal numérique temporel est mise en oeuvre à partir une analyse statistique de la largeur de bande passante du signal audionumérique. En effet, en codage audionumérique à bas débit il existe certaine corrélation entre le débit alloué et la largeur de bande passante du signal codé. En fait, plus le débit alloué est faible et moins bonne est la qualité de ce dernier. processus permettant de détecter strictement la bande passante du signal ne se révèle pas suffisant afin d'estimer la qualité perçue puisque un signal dont le contenu de bande passante faible, signal codé ou non codé, risque d'être considéré à tort comme dégradé. Compte tenu de 'observation précédente, il est donc nécessaire d'évaluer la fréquence critique de ce signal au-delà de laquelle un codeur ne peut plus procéder au processus de codage et non pas la bande passante du signal numérique en tant que telle. Selon un aspect particulièrement remarquable du procédé objet de la présente invention, cette approche est rendue possible en constatant que le spectre d'un signal codé possède, généralement, pour caractéristique une forte décroissance d'énergie à l'endroit de la coupure la fréquence critique précitée. Parallèlement, les spectres des signaux à faible contenu en haute fréquence se caractérisent en général pas par une telle cassure, mais au contraire par une décroissance lente de l'énergie, laquelle ne permet pas de discerner une séquence de référence d'une séquence codée. Le procédé objet de l'invention, en particulier le processus de calcul de la valeur Cb(t) liée à la bande passante du signal numérique, permet de vérifier la cassure précédemment mentionnée existe bien avant de considérer l'estimation du facteur de qualité comme valable. Une telle contrainte améliore considérablement la pertinence du procédé, objet de l'invention, dans le cadre de la définition d'un critère d'acceptabilité lié au défaut de codage. une manière générale, on indique que le procédé objet la présente invention n'est valable que pour les zones de signal contenant de l'information, c'est ' dire hors des zones de silence. En effet, le but est d'estimer en moyenne la dernière fréquence codée et non la bande passante instantanée du signal. Dans ce but, le signal temporel, ainsi que représenté en figure 3a, est soumis à une décomposition fréquentielle, transformation temps/fréquence, par transformée de @'ourier discrète par exemple sur N points du signal temporel pondéré par une fenêtre, telle qu'une fenêtre Hamming. La décomposition en fréquence est indiquée l'étape 220 à la figure 3a. Le spectre de puissance résultant de cette transformation comprend +1 points. L'étape 220 précitée peut alors etre suivie avantageusement d'une étape 22_ ccnsistant ' déterminer l'existence d'une zone de silence. Le test réalisé à l'étape 221 peut consister à comparer l'énergie du spectre obtenu à une valeur de seuil. Sur réponse négative au test 221, ce dernier est suivi d'une étape 222 consistant à découper en P sous- bandes de K raies spectrales d'énergie déterminée la décomposition fréquentielle du signal numérique temporel obtenu à l'étape 220. Chaque sous-bande de la décomposition contient K raies d'énergie ek. Les raies et les sous-bandes vérifient la relation : KxP = N/2. L'étape 222 précitée est alors suivie, pour les voies de gauche et de droite véhiculant le signal numérique ADS, d'une étape de calcul 223 l'énergie moyenne E1 contenue dans chaque sous-bande de rang i. énergie moyenne contenue dans chaque sous-bande de rang vérifie la relation (41, Dans la relation précédente, on indique que _1 désigne l'énergie de chaque raie spectrale considérée, constitutive de la sous-bande de rang i correspondante. étape 223 précitée est alors suivie d'un processus consistant à déterminer le rang i, spécifique de la sous bande de rang i correspondante, pour laquelle se produit la fréquence de coupure, ou cassure precédemment mentionnée, par au moins une comparaison du rapport de l'énergie contenue dans la dernière sous-bande prise comme niveau de référence de bruit de fond à l'énergie contenue dans les P-1 autres sous-bandes à une première valeur de seuil. A titre d'exemple non limitatif, pour la mise oeuvre du processus de détermination du rang i, spécifique de la sous-bande de rang i pour laquelle se produit fréquence de coupure, ce processus peut être mis en ceuvre partir d'une étape 224 consistant à lire la valeur rang i de la sous-bande considérée, valeur arbitraire i=P, à vérifier si la sous-bande de rang correspondant correspond à la sous-bande de fréquence de coupure et à comparer, en une étape de test 225, le niveau d'énergie contenu dans la sous-bande de rang i correspondant, niveau énergie noté El, à celui, noté Ep, contenu dans les P-1 autres sous-bandes à une valeur de seuil notée Seuill. opération de comparaison s'écrit Sur réponse négative au test 225, le rang de la sous-bande i est décrémenté à la valeur i-1 à l'étape La valeur de l'indice de sous-bande i est alors soumise à l'étape 229, à une comparaison à la valeur 1 permettant de vérifier si toutes les sous-bandes ont été prises considération. Sur réponse négative au test 229, le processus repris, l'énergie de la sous-bande de rang i correspondant, différent de 1, étant soumise à nouveau au test 225. Selon un premier mode de réalisation du processus représenté en figure 2a, on indique que l'étape 225 peut alors être suivie, sur réponse positive au test à l'étape 225 précitée, d'une étape 228 consistant à mémoriser le rang i,=i de la sous-bande de fréquence pour laquelle la fréquence de coupure est détectée. Cette mémorisation intervient de manière particulièrement avantageuse dans un tableau de valeurs de rang à une étape notée 230. L'étape 230 précitée est alors suivie d'une étape consistant à rechercher dans le tableau de valeurs mémorisées, par un programme de tri, la valeur du rang i,. dont l'occurrence est la plus grande. L'étape 231 est alors suivie d'une étape 232 permettant en fait de déterminer la fréquence de coupure Fc la plus probable pour les voies de droite et de gauche. comprend en particulier que la détermination de fréquence de coupure F, la plus probable, F,gauche, ,droite, est réalisée par conversion du rang i, en valeur de la sous-bande de fréquence correspondante. L'étape précitée est alors suivie d'une étape 2 consistant à calculer la valeur moyenne Q des fréquences de coupure gauche et droite normalisée par la fréquence coupure théorique maximale P, la valeur moyenne Q précitée vérifiant la relation (5) . Dans la même étape 233, la valeur moyenne fréquences Q peut alors être soumise à une normalisation sur critère psycho-acoustique défini par au moins une valeur de seuil de bonne qualité de codage audionumérique, notée Seui13, et une valeur de seuil de mauvaise qualité de codage audionumérique, notée Seuil4. A l'étape 233 précitée, la valeur moyenne Q peut alors être comparée par comparaison de supériorité à la valeur Seui14 et d'infériorité à la valeur Seuil3 suivant la relation Seui14 Q Seui13 ?. A titre d'exemple non limitatif, indique qu'une fréquence de coupure de l'ordre de 17 kHz implique une bonne qualité de codage audionumérique, alors qu'une fréquence de coupure de l'ordre de 10 implique un codage avec énormément de dégradations. valeurs pour Seui14 et Seui13 peuvent par exemple correspondre à des fréquences de 10 kHz et 17 kHz respectivement. L'étape 233 précitée peut alors être suivie d'une étape 234 consistant en fait à calculer une valeur réduite constituant la valeur Cb(t) liée à la bande passante, la valeur précitée vérifiant la relation (6) . La valeur réduite est ainsi obtenue par une translation et une mise à l'échelle pour obtenir la valeur Cb(t) liée à la bande passante et dont la valeur est comprise entre 0 et Ainsi qu'on l'a en outre représenté en figure 3a, et de manière particulièrement avantageuse, le processus de calcul de la valeur liée à la bande passante peut en outre comporter, dans un deuxième mode de réalisation, une étape supplémentaire permettant de s'assurer que la coupure détectée correspond bien à une cassure au niveau de l'énergie spectrale. Cette é-ape supplémentaire consiste en une deuxième condition introduite à l'étape 226, insérée entre les étapes 225 et 228 précédemment citées. Ainsi, outre la première comparaison de l'étape 225, le procédé et le processus de calcul représentés en figure 3a comportent, sur réponse positive à la première comparaison de l'étape 225, une deuxième étape de comparaison du rapport E1/E._1 l'énerg@le de la sous bande de rang i à l'énergie la sous-bande de rang suivant i+1 à une deuxième valeur de seuil, désignée par Seuil2. Ainsi, l'étape suivante de mémorisation du rang i_=i référencée 228, mémorisation de la sous-bande de fréquence pour laquelle la fréquence de coupure est détectée, est alors conditionnée à la réponse positive à la première et à la deuxième comparaison réalisée à l'étape 225 et 226. La réponse négative au premier et deuxième test de comparaison 226 est suivie, si i#1, d'un retour au premier test comparaison et d'un appel de l'étape de recherche du rang i, dont l'occurrence est la plus grande sinon, à l'étape A la suite d'essais réalisés, on indique que pour N, nombre de points de la décomposition fréquentielle égal à 2048, le nombre N pouvant toutefois être compris dans une plage de valeurs comprises entre [256,4096], le processus de calcul de la valeur Cb(t) liée à la bande passante est optimum pour les valeurs ci-après P = 32 ([2;N/2]) K = 32 ([1;N/4]) Seuil = 100 ([10;1000]) Seui12 = 17 ([5;50 1 ) Seui13 = 0,7 ;[C,51;1):; Seui14 = 0,4[0;C,49)). Dans les valeurs numériques précitées, on indique que les valeurs entre parenthèses et crochets indiquent des plages de valeurs possibles susceptibles de convenir pour les différents paramètres précités. Une description plus détaillée d'un processus de calcul de la valeur Cs(t) liée aux proprîétés stéréophoniques du signal numérique temporel sera maintenant donnée en liaison avec la figure 3b. Le processus de calcul de la valeur Cs(t) précitée est basé sur le principe selon lequel les voies de gauc-:e et de droite véhiculant les signaux sonores sont codées indépendamment. Ceci implique que les erreurs de codage sont décorrélées entre les deux voies, tandis que le contenu sonore des deux voies reste, sauf exception, relativement similaire. Le processus de calcul mis en oeuvre repose donc sur le fait que le signal résiduel différence des énergies des voies gauche et droite est proportionnel à l'erreur de codage s'il y a eu codage. L'intérêt d'une telle approche réside dans le passage d'une analyse sans référence à une analyse oseudo différentielle dans laquelle le signal d'erreur est déduit par comparaison des signaux numériques véhiculés par les deux voies. Toutefois, un tel processus ne permet pas d'évaluer la qualité du codage pour un signal fortement stéréophonique ou, au contraire, strictement monophonique. Pour cette raison, le processus de calcu représenté en figure 3b relatif au calcul de 1a valeur @s(t) liée aux propriétés stéréophoniq-aes du signal numérique temporel est basé sur le spectre d'énergie du signal numérique obtenu après décomposition fréquent,ielle par une transformée de Fourier sur N points du signal temporel, pondéré par une fenêtre de Hamming par exemple. spectre fréquentiel ainsi obtenu comprend N +1 raies. En conséquence, le signal temporel ainsi que représenté en figure 3b, est soumis à la transformée de Tourier sur N points à l'étape 220 telle que décrite précédemment en liaison avec la figure 3a. L'étape 220 précitée est alors suivie d'une étape consistant à calculer, pour chaque raie spectrale de rang k obtenue suite à la décomposition, fréquentielle, un facteur Qk représentatif de la qualité stéréophonique du signal à partir des spectres de fréquence de 1a voie gauche et SD de la voie droite. Le facteur Qk constitue en fait une différence normée des énergies des voies droite et gauche vérifiant la relation (7) . D'une manière plus spécifique, on indique que la valeur Qk=0 correspond à une raie de rang k et une fréquence strictement monophonique, alors la valeur Qk=1 correspond à une raie de rang k et à une fréquence fortement stéréophonique. Le processus de calcul de la valeur Cs(t) liée aux propriétés stéréophoniques du signai r@mérique consiste ensuite à déterminer _e pourcentage K(t) des raies appartenant à une bande de fréquence J_` do7:née dont le facteur Qk dépasse une valeur de seuil déterminée, notée SI, le pourcentage R(t) vérifiant la relation R(t) = n/K où n désigne le nombre de fois où le facteur Qk représentatif de la qualité stéréophonique du signal est supérieur à une valeur de seuil S, pour toute valeur de K appartenant à Of, la bande de fréquence précitée. A titre d'exemple non limitatif, afin de déterminer le pourcentage R(t), ainsi que représenté en figure 3b, ce processus peut consister à initialiser, à une étape , suite à l'étape 235 précitée, valeur de k indice raies de fréquences à la vale-r 0 la valeur de n à la valeur 0. L'étape 236 est suivie d' étape 237 consistant à comparer la valeur de l'indice de raies courant k à la valeur K nombre de raies issues de la décomposition spectrale. Sur réponse négative au test 237, ce test est suivi d'une étape 241 consistant à affecter à la valeur du pourcentage R(t) la valeur r./K pour la valeur de n. Au contraire, sur réponse positive au test 237, ce test est suivi d'un test 238 consistait à comparer la valeur du facteur Qk représentatif de la qualité stéréophonique du signal à la vale-,#r de seuil SI précédemment citée dans la description.. La comparaison s'écrit Qk > S1 ? . Sur réponse négative au test 238 de comparaison précité, la valeur de k désignart le rang de la raie spectrale est incrémenté d'une un--té à l'étape 240 et le processus de calcul est ramené à l'étape 23^ de vérification de comparaison d'infériorité du rang . à la valeur . Au contraire, sur réponse positive au ` 238, ce test est suivi d'une étape 239 d'incrémentat de la valeur d'une unité, cette étape d'@lncrémentation 239 étant elle-méme suivie de l'étape incrément at ion de l'indice k de la raie spectrale cc-rsidérée. L'étape 241 est alors s-_ivie d'une étape 242 consistant à corriger la valeur du pourcentage R(t; par une fonction spécifique A telle que la valeur de cette fonction du pourcentage R(t) soit comprise entre G et 1. - .a fonction A de la forme A(R:t)) est une fonction monotone croissante de la valeur du pourcentage R(t). A titre d'exemple non limitatif, la fonction A(R(t)) peut vérifier la relation L'étape 242 permet d'engendrer une valeur de pourcentage M(t), moyenne d'un nombre P déterminé de valeurs de pourcentage corrigées vérifiant la relation (8) Le processus de calcul de la valeur Cs(t) liée aux propriétés stéréophoniques du signal numérique temporel comporte égaiement une étape consistant à .é terminer, da-s une fenêtre temporelle de durée déterminée, fenêtre temporelle de s secondes, le nombre de fois F où une valeur de seuil d'alarme S2 a été franchie par la valeur de pourcentage corrigée A(R(t) . L'étape peut consister en une étape 245 de définition de la fenêtre et d'initialisation du nombre de fois F à la valeur 0, suivie d'une étape 246 de comparaison de supériorité de la valeur de la fonction A(R(t)) à la valeur S2 constituant un seuil d'alarme. La relation de comparaison s'écrit A(R(i)) S,- ? i désignant des instants successifs pendant la fenêtre de durée s. L'étape 246 est suivie d'une étape 24î consistant, sur réponse positive au test 246, à incrémenter la valeur du nombre de fois F d'une unité à l'étape 247, la réponse négative au test 246 ramenant à l'étape 245 pour passage à l'instant suivant appartenant à la fenêtre de durée s secondes. Les étapes 243 et 247 sont alors suivies d'une étape 244 consistant à calculer la valeur Cs(t) liée aux propriétés stéréophoniques du signal numérique temporel à partir d'une fonction de la valeur moyenne M (t) donnée à la relation (8), cette fonction vérifiant la relation (9) . En définitive, à un instant t, la valeur Cs(t) d'acceptabilité stéréophonique est donnée par la relation (9) précédemment mentionnée. Dans un exemple de :-lise en #uvre du processus de calcul représenté en figure 3b, on indique que pour N=2 , N pouvant être compris entre 256;9096], alors, le procédé est optimum pour les valeurs ci-après = [0;14,4 kHz] pour K nombre de raies spectrales obtenues = 614 ; = 0,99 ( [0,51;1]) = 1 seconde ([0,1;100]) P =100 ([ 1 ;1000]) S2 ( [0,01;l] ) . Dans les valeurs numériques précitées, on i 'que que les valeurs entre parenthèses et crochets désignent des plages de valeurs susceptibles d'être utilisées. Une description plus détaillée du processus de calcul de la valeur Cw(t) liée au blanchiment du signal numérique sera maintenant donnée en liaison avec la figure 3c. L'introduction du blanchiment du signal numérique permet d'effectuer une comparaison du signal numérique avant et après blanchiment. Le processus blanchiment est réalisé au moyen d'un filtre de blanchiment. Les propriétés d'un tel filtre sont les suivantes Pour un vecteur X constitué par les Ne échantillons temporels d'entrée du signal et pour le vecteur Y constitué par les Ne échantillons temporels de sortie du filtre de blanchiment, on désigne par W la matrice contenant les coefficients du filtre de blanchiment précité. L'expression du vecteur de sortie à partir du vecteur d'entrée est obtenue par la relation Y = WHX, -e symbole H indiquant les opérations da transposition, et conjugaison. Pour un signal numérique codé de qualité, le signal numérique soumis au blanchiment obtenu après passage dans le filtre de blanchiment correspond sensiblement à un bruit blanc dont la matrice de covariance RYY vérifie la relation désigne la puissance de ce bruit blanc et I la matrice identité. Toutefcis, Ryy est la valeur moyenne de la matri,-e YY.. notée . La matrice W contenant les coefficients du filtre étant considérée comme constante pendant la durée calcul de la valeur moyenne précitée, on obtient alors Relation (10) . Dans la relation précédente, RXx désigne la matrice covariance du signal temporel d'entrée. Cette matrice vérifie la relation (11) Etant admis que la matrice W possède une symétrie hermitienne, de la forme ',r; = W, la relaticn (11; précitée . s'écrit selon la relation #12) Des résultats expérimentaux ont montré qu'une approximation du type fournissait alors de bons résultats tout en simplif@lant de manière très importante les calculs. Globalement, le processus de calcul de la valeur Cw liée au blanchiment du signal numérique est réalisé de 1 manière suivante - Calcul de la matrice de covariance Rix_ du signa' numérique reçu ; - Filtrage passe-bas anti-repliement et décima-ion d' facteur 2 de ces signaux ; - Filtrage du signal décimé par la matrice de covariance inverse du signal initial. Le processus de filtrage ainsi mis en ceuvre correspond à un filtrage empirique pour lequel aucune justification théorique ne peut pour l'instant être établie. Ce processus n'est valablement mis en #uvre pour les zones de signal numérique reçu contenant l'information, c'est-à-dire hors des zones de silence. Dans ce but, suite à une étape de détection d' zone de silence 221, telle que décrite précédennent dans la description, le processus de calcul proprement dit est mis en oeuvre sur réponse négative à l'étape 221 précitée. Le processus est mis en couvre pour la fie de gauche, respectivement la voie de droite. Pour chacune des voies précitées le processus consiste alors à calculer la matrice de covariance Rg, Rd signal d'entrée et d'un signal a-'éatoire compris entre valeurs -1 et +1 aux étapes 250g, 250d. Cette opération peut être réalisée, ainsi que représenté de fière illustrative sur la figure 3c, par adjonction au signal numérique d'entrée de la voie gauche, respectivement de la voie droite, ::'un signal aléatoire engendré à une étape 248, ce signal aléatoire étant un signal de valeur comprise entre -1 et +1. Ce mode opératoire permet d'obtenir une matrice covariance toujours inversible. A partir des échantillons obtenus suite à la mise #uvre des étapes 249g et 249d, le calcul proprement dit de la matrice de covariance Rg et R.d aux étapes 250g et 250d peut être obtenu à partir du signal X, suite d'échantillons obtenue par la mise en #uvre des étapes 249g et 249d respectivement. La matrice X comprend 2xN2 échantillons et le calcul de la matrice de covariance Rg, Rd désignée sous la forme RXX est donnée la relation (13) Les éléments des matrices de covaria.^.ce Rg et Rd sont réels. Les étapes 250g et 250d sont alors suivies d'étapes de calcul des matrices de covariance inverses 251g et 251d respectivement. Les étapes précitées peuvent alors être suivies d'étapes de filtrage passe-bas anti-repliement 252g, 252d appliquées au signal numérique d'entrée sur les voies gauche et droite respectivement. Les étapes 252.: et 252d sont alors suivies d'une étape de décimation 253g, 253d, d'un facteur 2 pour engendrer une matrice d'entrée gauche et droite Eg, Ed respectivement. Ces opérations sont référencées aux étapes 254g et 254d respectivement. Les matrices Eg et Ed, matrices d'entrée, sont ebten-._es par rangement dans les matrices correspondantes des coefficients obtenus suite à l'opération de décimation 253g, 253d précitée. Suite à la création des matrices d'entrée -_- g et les étapes de filtrage permettant d'engendrer une matrice de sortie Sg à l'opération 255g et une matrice de sortie Sd à l'opération 255d est a@-ors réalisée à partir matrices d'entrée gauche Eg, respectivement droite Ed. Le signal de sortie pour la voie gauche respectivement droite est alors obtenu par 'opération vérifiant la relation (14) . Dans la relation précA-dente, S, R et E doivent être compris comme désignant Sg, Sd Rg, Rd et Eg, Ed respectivement. En référence à la figure 3c, on indique que le processus de calcul consiste alors, suite aux étapes 255g et 255d, à calculer à l'étape 256, à partir des matrices d'entrée et de sortie gauche et droite précitées, un rapport entre l'énergie du signal de sortie et l'énergie du signal d'entrée. Ce rapport, désigné par r, vérifie la relation (15) . relation précédente exprime le rapport en entre l'énergie du signal de sortile et l'énergie d',: signal 12 d'entrée, S9J'2 , Sd)2 désignant l'énergie du signal sortie les voies gauche, respectivement droite, 12 12 Ce k J e k I désignant l'énergie du signal d'entrée après décimation sur la voie gauche, respectivement droite, N désignant le nombre de lignes des matrices traitées lié au nombre d'échantillons par la relation Ne = 2xNxN. L opération 256 est alors suivie d'une opération 257 consistant, à partir des L dernières valeurs de rapport, rapport moyen entre l'énergie du signal de sortie et l'énergie du signal d'entrée, ce rapport moyen vérifiant la relation (16) . ce rapport moyen étant calculé dans une fenêtre glissante contenant les L derniers résultats. désignent l'énergie du signal d'entrée sur la voie gauche et droite, et N2 N2 Eg = E S2 (gauche) et = @ S2 (droite) s k=1 k k=1 k désignent l'énergie du signal de sortie sur la voie gauche et droite. L'étape 257 est alors suivie d'une étape consistant à soumettre la valeur ce rapport moyen à une comparaison de supériorité une première valeur de seuil S': et d'infériorité à une deuxième valeur de seuil S'2. Sur critère de comparaison précité satisfait, une étape de calcul de la valeur Cw(t) liée au blanchiment du signal numérique d'entrée est effectuée, cette valeur étant définie comme le rapport augmenté d'une unité de la différence du rapport moyen et de la deuxième valeur de seuil S'2 à la différence entre la deuxième S'-# et la première valeur de seuil S'Î La valeur Cw(t) liée au blanchiment du signal numérique d'entrée vérifie alors relation (17; . Sur la figure 3c, on a représenté les étapes consistant à soumettre la valeur rapport moyen à une comparaison de supériorité à la première et à la deuxième valeur de seuil S'i et S'., et de ca';cu-' de la valeur Cw(t) liée au blanchiment en une seule et même étape en raison. du fait que le calcul la valeur Cw;t? conditionné à la réussite la double comparaison de la valeur du rapport moyen valeurs de seuil S'1 et S'2 précitées. On obtient ainsi une valeur Cw(t) liée au blanchiment du signal d'entrée comprise entre la valeur 0 e` 1. Au contraire, en présence d'une zone de silence sur réponse positive au test 221, le rapport moyen n'est pas actualisé et la valeur Cw(t) liée au blanchiment du signal numérique d'entrée garde la valeur à l'instant précédent t-1. La valeur à l'instant précédent est donc utilisée comme valeur à l'instant courant. Des résultats expérimentaux ont permis de montrer que pour N=16, la matrice d'entrée contient 512 échantillons et le procédé est optimum pour les valeurs suivantes du filtre passe-bas anti-repliement utilisé pour réaliser les opérations aux étapes 252g et 252d. Ces valeurs sont données dans le tableau ci-après, pour un filtre anti-repliement comprenant K=93 coefficients. -0.0006 -0.0017 -0.0022 0.0010 0 .OIC6 0.0253 0.0376 0.U372 C.0193 -0.0082 -0.0268 -0.0203 0.0 0 87 C.0358 0.0323 -0.0086 -.C572 -0.0626 0.0089 0.14_3 @0.2@.7 0.32y4 C.2707 J.1913 0. 0089 -0.0626 -0.0572 -0.0086 0.r323 @. 0358 C . C037 -0.0268 -0.0082 0.0 1 93 0.0372 0.0376 0.0253 0.C106 0.0010 -0.0022 -0.0017 -0.0006 I La fenêtre glissante contenant les 1_ derniers résultats est L=100, la valeur L pouvant tc,_:=efois être comprise entre ([10;1000]). La valeur de seuil S'1 est égale à -60 dB et S'., = dB. Une description plus détaillée des opérations de détection de micro-coupure, de sifflement ou saturation, de bourdonnement et d'existence d'un déphasage entre voies mises en ceuvre à l'étape 3 par les étapes 31, 30, 32 et de figure 2 sera maintenant décrite en liaison avec figures 4a, 4b, 5a, 5b, 6a, 6b et 7. En ce qui concerne l'étape 31 de détection d' micro-coupure, encore désignée par brève coupure, on indique que celle-ci peut avantageusement consister à détecter, sur une suite d'échantillons successifs du signal numérique ADS, une décroissance rapide du niveau d'énergie de ce signal audionumérique vers une énergie nulle révélant une absence de réverbération du signal audionumérique précité. Sur la figure 4a, l'axe des abscisses est gradué en millisecondes et l'axe des ordonnées en amplitude, la brève coupure, également désignée sous le nom de mute étant représentée comme la décroissance rapide du niveau d'énergie du signal audionumérique vers une énergie nulle. En référence à la figure 4b, on indique de manière non limitative que l'étape de détection d'un signal parasite telle qu'une brève coupure peut comprendre étape 401 consistant à déterminer séparément sur chaque voie stéréophonique, pour une pluralité de suites de M échantillons successifs, l'énergie moyenne E., du signal transporté par cette voie, n désignant le rang de chaque suite d'échantillons Sr. L'étape 401 est suivie d'une étape consistant à comparer l'évolution de l'énergie moyenne pour les suites de M échantillons successifs. L'étape précitée peut être réalisée par comparaison de l'énergie moyenne En du signal transporté à _ valeur 0 à l'étape , puis d'une comparaison 403 d' ou plusieurs des énergies moyennes précitées à une val de seuil AdB. Ainsi, l'existence d'un signal parasite de brève coupure est révélée si l'une au moins des énergies moyennes est nulle et si une ou plusieurs énergies moyennes voisines de cette énergie moyenne nulle sont supérieures à une valeur de seuil donnée, la valeur 0. En ce qui concerne l'étape 3C de détection de sifflement ou de saturation, on indique cette étape sera décrite dans le cas de la détection un sifflement, une saturation étant le plus souvent accompagnée d'un sifflement. référence à la figure 5a, on indique que la détect' d'un signal parasite tel qu'un sifflement dans le signal audionumérique ACS peut consister avantageusement à détecter dans ce signal une augmentation subite transitoire de l'énergie spectrale de ce dernier dans bande de fréquences dont la fréquence basse est comprise entre 4,5 kHz et 6,5 kHz et dont la fréquence haute peut atteindre jusqu'à 20 kHz. Sur la figure 5a, l'axe des abscisses est gradué en fréquences et l'axe des ordonnées en énergies correspondantes pour les bandes de fréquences considérées. En référence à la figure 5b, on indique que le processus de détection d'un signal parasite tel qu'un sifflement peut comprendre une étape 501,502 consistant à -alculer sur une su,;te d'échantillons signal audionumérique ADS la composition. spectrale ce signal définie comme la valeur SI(i) de composantes fréquentielles en sous-bandes de fréquence centrale f; et de largeur de bande Af, n désignant le rang de la suite d'échantillons. Les étapes 501 et 502 sort alors suivies d'une étape 503,504 consistant à calculer la valeur moyenne de l'énergie En(sb) d'une plage des sous-bandes précitées pour la suite d'échantillons de rang r. considéré. Une étape de calcul 506 d'une valeur de contraste auditif alors réalisée, Cn,s-b à partir de valeur du rapport . Ce rapport calculé à l'étape 505 désigne le rapport entre l'énergie En(sb) de cette plage pour la suite courante et pour une pluralité de suites précédentes En_S(sb) d'échantillons. La valeur de contraste auditif vérifie la relation (18) . Dans cette relation, Rn(sb+i, désigne, pour i=-v, la valeur du rapport pour les sous-bandes voisines de la même suite d'échantillons de rang n et du même spectre Sr. En outre, à l' étame une ----.paraisor. de la valeur de contraste audit--- sb à une premiére valeur de seuil de sifflement, - " S51, est réalisée, la comparaison étant notée C--,,@ > L'étape 506 précitée suivie d'une étape 507 de calcul d'un paramètre de proximité, noté P,,Sb, vérifiant la relation (19) En outre, à l'étape une comparaison. du paramètre de proximité P..,Sb deuxième valeur de sifflement Ssz est réalisée, la comparaison étant notée Pr,sb > Ss2. La présence d'un signal parasite de sifflement est révélée si les comparaisons de supériorité de la valeur de contraste auditif et du paramètre de proximité sont toues deux vérifiées. En ce qui concerne l'étape de détection d'un signal parasite de bourdonnement réalisée à l'étape 32, on indique que cette étape, en référence à 'a figure 6a, peut consister à détecter un signal parasite constitué par un bruit rose dans une bande de fréquences comprise entre 0 et 1100 Hz et de niveau sensiblement constant dans 1a bande de fréquences précitée. la figure 6a, l'axe des abscisses est gradué en fréquences et "axe des ordonnées en niveau d'énergie du signal exprimé en décibels. On constate que dans la bande fréquences précitée, un niveau sensiblement constant, voisin de 40 dB, peut être mis en évidence en présence d' bourdonnement. En référence à la Figure cc, le processus de mise en évidence d'un signa: parasite de bourdonnement peut comprendre, sur au :oins une vole anche e'.: droite de ce signal une étape 701 consistant à calculer, la suite d'échantillons du signal numérique ADS, la composition spectrale de ce signal définie comme la valeur Sn(i) de composantes fréquentielles en sous-bande, fréquences centrales fi où n désigne le rang la suite d'échantillons considérée. L'étape 701 est suivie d'une étape 702 pour un nombre déterminé k fréquences centrales fi du domaine des basses fréquences 'étape 702 consistant à calculer un premier et un deuxième rapport des valeurs de composantes fréquentielles en sous-bande pour suite d'échantillons courante la suite d'échantillons précédente, ce premier rapport étant S (i) désigné par oc. = n et le deuxième rapport pour la i,n Sn_lO) suite échantillons courante et la suite échantillons S (i) suivante étant désigné par (3 = n Sn + I (1) L'étape 702 consiste également à comparer valeur des premier et deuxième rapports précités à première valeur de seuil de bourdonnement, notée Sb:. Sur répcnse négative à la comparaison précitée, l'étape 702 est rebouclée, 703, par une incrémentation de 1 indice i en i=i+l. Sur réponse positive à l'étape 702, cette dernière est suivie d'une étape 709 consistant à soumettre la comparaison des premier et deuxième rapports à un critère de proportion du nombre p de comparaisons vérifiées par rapport ' la totalité des k comparaisons effec-uées pour les k fréquences centrales f1. L'étape 704 consiste à réaliser un test de vérification que Pô raies fréquentielles remplissent la condition précédente sur la suite courante Sn. Sur réponse négative au test 704, un bouclage 708 permet de passer à la suite d'échantillons suivante rang n+1. Sur réponse positive au test 704, une étape 705 est réalisée, consistant à discriminer parmi les valeurs Sn(i) de composantes fréquentielles en sous-bandes, la valeur maximale Sn(imax) des valeurs de composantes fréquentielles relatives à la suite d'échantillons courante. L'étape 705 est elle-même suivie d'une étape 706 consistant à calculer le rapport de la valeur maximale avec la valeur correspondant à l'index du spectre de la suite précédente Sn_1 (lmax) . Ce rapport est noté S (i ) - n, i n max . En outre, ce rapport est comparé à une Sn-1(lmax ) deuxième valeur de seuil de bourdonnement notée Sb2 par comparaison d'infériorité. Ainsi, on comprend que, sur au moins une voie de transmission en mode stéréophonique signal audionumérique ADS, la détection d'un signal parasite de bourdonnement consiste à détecter l'existence d'une comparaison de supériorité des premier et deuxième rapports ai,n et Pi,n à la première valeur seuil de bourdonnement Sbl et l'existence d'une comparaison d'infériorité du rapport des valeurs maximales Mn,i à la deuxième valeur de seuil de bourdonnement Sb2. Suite à 'étape 706 précitée, une analyse statistique est réalisée par répétition opérations précédentes et mémorisation périodique sur durée s' déterminée d'une variable binaire de prédétection de ''existence d'un signal parasite de bourdonnement. A la variable binaire de prédétection est attribuée la valeur ? lorsque les critères de comparaison de supériorité et d'infériorité sont satisfaits et la valeur 0 sinon. L'analyse statistique consiste à décompter, à l'étape 707, dans la durée s' déterminée, le nombre d'occurrences valeur 1 de la variable binaire de prédétection comparer ce nombre à une troisième valeur de seuil de bourdonnement, tee S,3. Ainsi, lorsque, sur une observation de s' secondes, un nombre d'occurrences est supérieur à Sb3, la présence d'un signal parasite de bourdonnement est révélée lorsque la comparaison précitée est vérifiée. En ce concerne la mise en oeuvre de l'étape 33 de calcul du déphasage d, on indique, en référence à la fïgure 7, que cette étape peut consister à calculer à l'étape A la valeur du déphasage entre voies du signal audionumérique ADS à partir de la fonction d'inter corrélation du signal audionumérique présent sur chacune des voies, puis à comparer à l'étape B la valeur de déphasage d a une valeur de seuil. Sur la figure 7, les valeurs de déphasage et de seuil sort notées d respectivement dmax. Pour ce qui concerne la mise en oeuvre des étapes de détection de sifflement ou de saturation 30, de micro coupure 31, de bourdonnement 32 et de déphasage inter voies 33, d'autres procédures peuvent être mises en #uvre. Toutefois, les procédures indiquées dans la présente demande de brevet apparaissent particulièrement satisfaisantes. Pour une description plus détaillée de la mise en #uvre de ces procédures, on pourra utilement se reporter à la demande de brevet français n 99 09179 déposée le 08/03/1999 au nom titulaires de la présente demande. REVENDICATIONS 1. Procédé de contrôle en continu de qualité des sons en distribution, les sons numériques étant disponibles en mode stéréophonique selon signal numérique représentant au moins une voie droite et une voie gauche, caractérisé en ce qu'il consiste à effectuer une analyse statistique du contenu de ce signal numérique sur chacune desdites voies, ladite analyse statistique consistant . - à découper temporellement ledit signal numérique selo:@ des suites successives d'échantillons, comportant un nombre déterminé d'échantillons, et - à effectuer, lors de l'existence d'un programme de sons numériques, une analyse spectrale de chacune des suites d'échantillons pour observer les variations d'énergie et enveloppe dudit signal numérique dans le domaine temporel et fréquentiel et calculer un indice de qualité global ; - à calculer à partir desdites variations d'énergie et d'enveloppe et de l'indice de qualité global un indice de qualité final, valeur bornée et continue dans le temps, représentatif de la qualité desdits signaux numériques. 2. Procédé selon la revendication 1, caractérisé en ce lesdites suites d'échantillons sont constituées par suites d'échantillons présentant un taux de recouvrement, rapport du nombre d'échantillons communs à deux suites consécutives au nombre d'échantillons constitutif de chaque suite d'échantillons, compris entre 0 et -75%. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que 'Ladite étape cons-stant à calculer lors de l'existence d'un programme de sons numériques distribués un indice de qualité global consiste au moins à calculer indice de qualité global à partir d'au moins un critère fréquentiel et d'un critère temporel de variation 'énergie et d'enveloppe. 4. Procédé selon l'une des revendications 1, 2 ou 3, caractérisé en ce que ladite étape consistant à calculer à partir desdites variations d'énergie et d'enveloppe et de l'indice de qualité global un indice de qualité final consiste au moins . - à détecter l'existence sur ledit signal numérique d'au moins une perturbation de transmission dudit signal numérique et à affecter à l'existence de cette perturbation un coefficient de pondération spécifique, représentatif de la contribution de cette perturbation à la dégradation de la qualité desdits signaux numériques, la valeur de ce coefficient de pondération étant égale à 1 sinon ; - à pondérer la valeur dudit indice de qualité global par la valeur du produit de l'ensemble coefficients de pondération, pour obtenir un indice qualité global pondéré; - à détecter la valeur d'un déphasage inter-voies et à affecter une valeur de critère de déphasage spécifique à cette valeur de déphasage lorsque cette valeur de déphasage est supérieure à zéro et une valeur de critère de déphasage égale à zéro sinon ; - à déterminer ledit coefficient de qualité final par comparaison de la différence entre ledit coefficient de qualité pondéré et ladite valeur de critère de déphasage à la valeur zéro et à attribuer valeur égale à 1 audit coefficient de qualité global en 'absence de programme de sons numériques distribués. 5. Procédé selon la revendication 3, caractérisé en que l'étape de calcul dudit indice de qualité global I(t) est effectuée à partir d'un critère de valeur Cb(t) liée à la bande passante, d'un critère de valeur Cs(t) liée aux propriétés stéréophoniques et d'un critère de valeurs Cw(t) lié au blanchiment du signal numérique temporel, lesdites valeurs Cb(t), Cs(t) et Cw( ) étant constituées par des valeurs réelles positives comprises entre 0 et 1, ledit indice de qualité global I étant défini par une combinaison linéaire desdites valeurs et étant constitué par une valeur réelle, comprise entre 0 et 1. 6. Procédé selon la revendication 4, caractérisé en ce que l'étape consistant à détecter l'existence sur ledit signal numérique d'au moins une perturbation de transmission consiste à détecter une perturbation choisie parmi les perturbations de sifflement ou saturation, de micro-coupure et de bourdonnement respectivement. 7. Procédé selon la revendication 5, caractérisé en ce que l'étape de calcul de la valeur Cb (t) liée à la bande passante consiste à partir d'une décomposition fréquentielle du signal numérique temporel, à - discriminer l'existence d'une zone de silence et, en 'absence de zone de silence, - découper en P sous-bandes de K raies spectrales énergie déterminée ladite décomposition fréquentielle du siqnal numérique temporel ; - calculer pour les voies gauche et de droite l'énergie moyenne Ei contenue dans chaque sous-bande de rang i ; - déterminer le rang i, spécifique de la sous-bande de rang i correspondant pour laquelle se produit la fréquence de coupure, par moins une comparaison du rapport de l'énergie contenue dans la dernière sous- bande, prise comme niveau référence de bruit de fond, à l'énergie contenue dans les P-1 autres sous- bandes à une première valeur de seuil ;et sur réponse positive à cette comparaison, - mémoriser le rang ic = i de la sous-bande de fréquences pour laquelle la fréquence de coupure est détectée, dans un tableau de valeurs de rang ; - rechercher dans ce tableau par un programme de tri de la valeur du rang i dont l'occurrence est 1a plus grande, puis déterminer la fréquence de coupure Fc la plus probable pour les voies de droite et de gauche ; - calculer la valeur moyenne fréquences de coupure gauche et droite normalisée fréquence de coupure théorique maximale, P, - normaliser ladite valeur moyenne des fréquences sur critère psycho-acoustique défini par au moins une valeur de seuil (Seui13) de bonne qualité de codage audionumérique et une valeur de seuil (Seui14) de mauvaise qualité de codage audionumérique par décalage et calcul d'une valeur réduite constituant ladite valeur Cb(t) liée à la bande passante et vérifiant la relation 8. Procédé selon la revendication 7, caractérisé en ce que l'étape consistant à déterminer 1e rang i, spécifique de la sous-bande de rang i correspondant pour laquelle fréquence de coupure se produit comporte, outre première comparaison du rapport Ei/Ep de l'énergie contenue dans la dernière sous-bande, à l'énergie contenue dans les P-1 autres sous-bandes à une première valeur de seuil, Seuill, sur réponse positive à cette première comparaison, une deuxième étape de comparaison du rapport Ei/Ei+l, de l'énergie de la sous bande rang i à l'énergie de la sous-bande de rang suivant i+1 à une deuxième valeur de seuil, Seui12 l'étape suivante de mémorisation du rang i`=i de la sous bande de fréquences pour laquelle la fréquence de coupure est détectée étant conditionnée à la réponse positive auxdites première et deuxième comparaisons, la réponse négative auxdits premier et deuxième tests de comparaison étant suivie, si i$1, d'un retour au premier test comparaison et d'un appel de l'étape de recherche du rang i, dont l'occurrence est la plus grande sinon. 9. Procédé selon la revendication 5, caractérisé en ce que l'étape de calcul de la valeur Cs(t) liée aux propriétés stéréophoniques du signal numérique temporel consiste, partir d'une décomposition fréquentielle en raies de rang k du signal numérique temporel, à: - calculer, pour chaque raie de rang '.K, un facte-r Qk représentatif de la qualité stéréophonique signal, à partir des spectres de fréquence Sk de la voie gauche et SD de la voie droite, différence normée des énergies des voies droite et gauche de la forme - déterminer le pourcentage R(t) des raies appartenant à une bande de fréquences AF donnée dont facteur Qk dépasse une valeur de seuil déterminée S:, (t) = n/K, :-: etant le nombre de fois où Qk > Si b'k E 4F; - corriger la valeur du pourcentage R(t) par une fonction spécifique A telle que 0 5 A(R(t)) S 1, pour engendrer valeur de pourcentage M (t) moyenne ci un nombre P determiné de valeurs de pourcentage corrigé - déterminer dans une fenêtre temporelle de durée déterminée le nombre de fois F où une valeur de seuil d'alarme Sz a été franchie par la valeur de pourcentage corrigé A(R(t)) ; - calculer la valeur Cs(t) à partir d'une fonction de ladite valeur moyenne, de la forme 10. Procédé selon la revendication 5, caractérisé en ce que l'étape de calcul de la valeur Cw(t) liée blanchiment du signal numérique temporel consiste, partir dudit signal temporel, pour chacune des voies, en l'absence de détection d'une zone de silence - à calculer la matrice de covariance (Rg,Rd) du signal entrée et d'un signal aléatoire, compris entre les valeurs -1 et +1 ; calculer la matrice inverse de la matrice covariance ; - ' soumettre le signal d'entrée à un filtrage passe- anti-repliement et à une décimation d'un facteur deux, pour engendrer une matrice d'entrée gauche et droite (Eg, Ed) ; - à calculer, à partir de la matrice d'entrée gauche et droite, une matrice de sortie gauche et droite (Sg, Sd) ; - à calculer, à partir des matrices d'entrée et de sortie gauche et droite, un rapport entre l'énergie du signal de sortie et l'énergie du signal d'entrée ; - à calculer, à partir des L dernières valeurs rapport, un rapport moyen (r) entre l'énergie du signal de sortie et l'énergie du signal d'entrée ; - à soumettre la valeur de ce rapport moyen à comparaison de supériorité à une première valeur de seuil S'1 et d'infériorité à une deuxième valeur de seuil S'2 ; - à calculer la valeur Cw(t) liée au blanchiment comme le rapport, augmenté d' unité, de la différence du rapport moyen r et de deuxième valeur de seuil S'2 à la différence entre la deuxième S'z et la première S': valeur de seuil.