Brevet FR3030963A1 — publié le 20160624

La présente invention concerne le domaine technique de la transmission radiofréquence et plus particulièrement de l'égalisation aveugle de canal. Dans une transmission radiofréquence, un signal M est transmis entre un émetteur et un récepteur. Le canal radiofréquence par lequel il transite, du fait notamment des problèmes causés par le milieu où se propage les ondes hertziennes, déforme le signal émis M . Aussi c'est un signal X qui est effectivement reçu par le récepteur. Il est habituel de modéliser le signal X par une fonction H figurant les modifications apportées par le canal, appliquée par convolution au signal d'origine M, auquel s'ajoute un bruit B de mesure, Soit X = H(M)+ B. /0 Un exemple de perturbation est causé par les réflexions multiples ou multi- chemins. Suite à une réflexion sur un obstacle un signal voit son trajet se rallonger, relativement à un signal se propageant plus directement. Il s'ensuit qu'un récepteur peut recevoir, superposées à un signal propagé selon un chemin plus court, direct ou comportant moins de réflexions, une ou plusieurs occurrences du même signal, retardées 15 car propagées selon des chemins plus longs. Le but d'un égalisateur de canal, encore dénommé estimateur ou identificateur de canal, est de corriger a posteriori les modifications causées par un canal à un signal. Pour cela un égalisateur traite le signal d'entrée X reçu, et produit un signal de sortie Y, aussi proche que possible du signal M émis. Un égalisateur est ainsi un 20 dispositif qui applique au signal d'entrée X une fonction corrective F, afin de produire un signal de sortie Y, avec Y =F(X),--M . La qualité d'un égalisateur est d'autant meilleure que le signal de sortie Y reproduit, le cas échéant avec un retard T correspondant au délai de transmission, le signal d'origine M, soit que la quantité *(t)- M ) est la plus réduite possible. 25 Afin d'aider à résoudre la problématique d'égalisation, il peut être avantageux d'inclure dans le signal émis M , des éléments de référence connus, qui lorsqu'ils sont reconnus, le cas échéant modifiés, dans le signal reçu X, aident le récepteur à identifier le canal et facilitent ainsi la détermination de la fonction corrective F. Cependant il n'est pas toujours possible ou souhaitable d'inclure de tels éléments de référence. Un 30 égalisateur est dit aveugle ou autodidacte lorsqu'il n'utilise pas de tels éléments de référence. Pour réaliser un égalisateur aveugle il est connu d'employer un filtre adaptable. Un tel filtre produit le signal de sortie Y à partir du signal d'entrée X. Le filtre est adapté, par asservissement, en fonction d'une fonction de coût J. Cette fonction de coût J est typiquement construite autour d'un ou plusieurs critères, indicatifs de la grandeur Efr(t)-M(t ), que l'on cherche à minimiser et/ou annuler, alors même que M est inconnu. La problématique de l'égalisation aveugle a été théorisée par GODARD sous 5 le nom de « forçage à module constant » (ou en anglais : « constant modulus algorithm », CMA). Un des résultats, issu d'hypothèses statistiques sur le signal émis M, \ 2 \ 2 \ permet de remplacer *(t)-M(t ) ( , par E (d 2 -r) , où r est un rayon de convergence que l'on sait déterminer ou estimer à partir du signal d'entrée X. Ceci fournit un moyen d'estimer un critère pour la fonction de coût J. 10 D'autres critères ont encore été proposés pour réaliser une égalisation aveugle : SATO, KURTOSIS, BENVENISTE. Un autre résultat de GODARD est que ce seul critère n'est pas suffisant en ce qu'il ne garantit pas la convergence vers la solution optimale. Aussi nécessite-t-il des améliorations. 15 Muni d'une fonction de coût J, il est connu, en appliquant une méthode dite du gradient stochastique, d'adapter le filtre de manière convergente. Dans la pratique, l'implémentation d'une telle méthode au moyen d'un filtre numérique adaptable, nécessite une puissance de calcul importante. Dans le cas d'application à un autoradio, installé dans un véhicule automobile, il convient d'utiliser un 20 processeur durci « automotive » (c'est-à-dire apte à satisfaire les conditions imposées par les cahiers des charges du secteur automobile), adapté pour résister à un tel environnement d'utilisation. De tels processeurs durcis « automotive » présentent classiquement des performances en retrait par rapport aux processeurs « grand public » que l'on peut rencontrer dans des applications dont l'environnement est moins 25 contraignant : ordinateurs de bureau, téléphone cellulaire, etc. Dans tous les cas, une augmentation de performance pour un processeur induit une augmentation sensiblement exponentielle du coût et il est toujours intéressant de réduire les besoins en termes de quantité de calcul. Dans l'état de développement technologique actuel, les processeurs durcis 30 « automotive » sont capables d'implémenter une méthode GODARD de base. Cependant ils restent encore incapables d'implémenter une méthode comprenant des compléments pourtant nécessaires. Dans le cadre d'une utilisation automobile, du fait de la mobilité, le canal est modifié, et il apparait d'importantes variations de la puissance du signal reçu. Cette 35 variation est susceptible de causer une divergence de l'asservissement adaptant le filtre. Une réponse classique pourrait être d'appliquer préalablement une normalisation de la puissance au moyen d'un asservissement à gain automatique (en anglais : automatic gain control, AGC). Cependant une telle approche n'est pas souhaitable dans le cas d'utilisation présent en ce qu'elle modifie le signal d'entrée X , perturbant ainsi divers traitements et mesures réalisés en aval du dispositif égalisateur. Ceci est renforcé en ce que le dispositif égalisateur est optionnel et que son ajout ou retrait doit être transparent pour le récepteur radiofréquence. L'objet de la présente invention vise à répondre à ces différents inconvénients en proposant un égalisateur aveugle de canal, apte à compenser l'excursion en puissance du signal d'entrée X. Pour cela l'invention propose un dispositif égalisateur aveugle de canal, pour un récepteur radiofréquence adapté à une modulation à signal d'enveloppe constante à l'émission, comprenant : - un filtre numérique linéaire adaptable, défini à un instant par ses coefficients, apte à filtrer un signal d'entrée pour produire un signal de sortie, - un estimateur apte à estimer une puissance du signal d'entrée, - un adaptateur apte à adapter le filtre en calculant les coefficients du filtre à un instant, en soustrayant aux coefficients du filtre à un instant précédent le gradient d'une fonction de coût affecté d'un coefficient de correction, où la fonction de coût comprend un premier critère de distance entre le carré du signal de 20 sortie et la puissance, où le coefficient de correction est un produit comprenant un coefficient de convergence constant, et un coefficient de mise à l'échelle inversement proportionnel au carré de la puissance. Une telle caractéristique permet avantageusement, d'adapter le coefficient de 25 correction du gradient stochastique afin de tenir compte de la puissance effectivement reçue pour le signal d'entrée. Elle permet de plus d'éviter l'emploi d'un asservissement à gain automatique, préjudiciable aux autres modules disposés en aval. Selon une autre caractéristique, la fonction de coût est défini(e) par la formule J 2 - avec J la fonction de coût, y (n) l'échantillon du signal de sortie à 30 l'instant n, et P la puissance. Selon une autre caractéristique, le dispositif comprend encore un moyen de dispersion, remplaçant le gradient de la fonction de coût, défini par la formule W = y(n)I2 P) y (n) X (n) , avec V l'opérateur de gradient, J la fonction de coût, y (n) l'échantillon du signal de sortie à l'instant n, P la puissance, X (n) le vecteur des 35 échantillons x(k) du signal d'entrée, k e - K étant la dimension du filtre, et z désignant le conjugué d'un nombre complexe z, par la formule VJ = (n12 P) y (n + q) X (n) , avec q entier aléatoire, tel que -Q L'invention propose ainsi nombre de caractéristiques inédites, pouvant être mise en oeuvre seules ou de manière complémentaire, afin d'offrir un dispositif égalisateur aveugle de canal amélioré. Parmi ces améliorations, plusieurs visent à simplifier le dispositif, en vue de sensiblement réduire la puissance de calcul nécessaire, afin notamment de pouvoir l'implanter sur un processeur durci « automotive », le cas échéant à virgule fixe. D'autres caractéristiques, détails et avantages de l'invention ressortiront plus clairement de la description détaillée donnée ci-après à titre indicatif en relation avec des dessins sur lesquels : - la figure 1, présente un synoptique d'une chaîne de réception radiofréquence, comprenant un dispositif égalisateur aveugle, - la figure 2 détaille un tel dispositif égalisateur, - la figure 3 illustre le découpage de l'intervalle de puissance en zones afin de déterminer une approximation du coefficient de mise à l'échelle, - la figure 4 est un tableau contenant des valeurs numériques illustratives de l'approximation, - la figure 5 illustre l'erreur d'approximation du coefficient de mise à l'échelle. L'invention est applicable à un récepteur radiofréquence adapté à une 20 modulation à signal d'enveloppe constante à l'émission. Une telle modulation comprend la modulation de fréquence et la modulation de phase. Telle qu'illustrée à la figure 1, une chaîne de réception radiofréquence, du type adaptée à une modulation à signal d'enveloppe constante à l'émission, comprend classiquement, de gauche à droite, une antenne 2, un sélecteur 3, un démodulateur 4 et 25 un moyen de reproduction sonore 6. Le sélecteur 3 est en charge d'extraire, du signal hertzien complet reçu par l'antenne 2, un signal sensiblement centré sur la fréquence centrale d'un canal. Un sélecteur 3 peut être schématisé par un passe bande centré sur la fréquence centrale de modulation du canal, réalisant une extraction d'une partie du spectre. Le sélecteur 3 ne modifie pas la fréquence, qui reste la fréquence de base. Le 30 démodulateur 4 est en charge de réaliser une démodulation du signal sélectionné afin de le transformer en un signal audio apte à être reproduit. Le moyen de reproduction sonore 6 est en charge de cette reproduction. Il est figuré par un haut-parleur, mais peut comprendre tout système de reproduction comprenant un casque, un ou plusieurs haut-parleurs et le cas échéant un dispositif d'amplification. 35 Le dispositif égalisateur 1 s'intercale entre le sélecteur 3, duquel il reçoit un signal d'entrée X, en fréquence de base, et le démodulateur 4, auquel il fournit, après traitement d'égalisation, un signal de sortie Y de caractéristiques similaires à celles du signal d'entrée X. Une chaîne de réception radiofréquence peut encore comprendre, au moins un dispositif 5 de traitement du signal audio. Un tel dispositif 5 utilise typiquement des informations issues du signal (signal de sortie Y ou d'entrée X en l'absence de dispositif égalisateur 1) pour améliorer le rendu audio. Ces informations issues du signal, avant ou après démodulation, sont mesurées par des capteurs 7, tels un capteur de mesure de niveau de champ ou un capteur de mesure de niveau de modulation. Le dispositif égalisateur 1 étant optionnel, il est avantageux que le signal de 10 sortie Y présente une dynamique comparable à celle du signal d'entrée X, afin de ne pas perturber ces capteurs 7 et le dispositif 5 de traitement et que l'ajout éventuel d'un dispositif égalisateur 1 reste transparente pour ces modules. La figure 2 détaille le dispositif égalisateur 1 et son contenu. Le principe d'un tel dispositif égalisateur 1 est de filtrer le signal d'entrée X pour obtenir le signal de 15 sortie Y. Pour cela le dispositif égalisateur 1 comprend un filtre 10 adaptable. Le filtre 10 est adaptable afin que le signal de sortie Y reproduise, au plus près, le signal émis M. Le principe de l'invention peut être appliqué à tout type de filtre. Cependant la description qui suit est plus particulièrement illustrée dans le cas où le filtre est réalisé au moyen d'un filtre 10 numérique. 20 Une hypothèse simplificatrice suppose que la perturbation apportée par la transmission du signal est linéaire et que l'on peut la corriger au moyen d'un filtre linéaire. Un tel filtre 10 linéaire est entièrement défini par ses coefficients a(k) , k E 10, ..., où K est la dimension du filtre 10. Ces coefficients a(k) sont regroupés dans un vecteur A qui permet de les désigner globalement. Ces coefficients étant modifiés à chaque instant 25 d'échantillonnage, par l'adaptateur, on note A(n) le vecteur des coefficients à un instant n. Afin que le signal de sortie Y recopie au mieux le signal d'origine M, le dispositif égalisateur 1 comprend un adaptateur 12. Cet adaptateur 12 est apte à adapter le filtre 10 en calculant, les coefficients A(n) du filtre 10 à l'instant n. Selon la méthode 30 connue du gradient stochastique, une convergence du filtre 10 peut être obtenue en calculant, les coefficients A(n) du filtre 10 à l'instant n, de manière récurrente, en fonction des coefficients A(n -1) du filtre 10 à l'instant précédent n-1 et d'une fonction de coût J. Ce calcul est réalisé selon la formule A(n) = A(n -1)- /1VJ , où A(n) est le vecteur des coefficients a(k) , k E 10, ..., K}, du filtre 10, K étant la dimension du filtre 10, 35 à l'instant n, A(n -1) est le vecteur des coefficients du filtre 10 à l'instant précédent n- 1, 2 est un coefficient de correction, V est l'opérateur de gradient et J est la fonction de coût. La fonction de coût J à minimiser comprend un critère indicatif de la grandeur Efr(t)-M(t \12 ). Le critère est ici une distance entre le carré du signal de sortie Y et la puissance P, indicative d'un rayon de convergence r, avec une relation P = r2, issu de la méthode de GODARD. Puisque le critère dépend de la puissance P, le dispositif 1 comprend encore un estimateur 11 apte à estimer la puissance P du signal d'entrée X. Le coefficient de correction 2 est un élément important en ce qu'il détermine 10 la vitesse de convergence de l'asservissement d'adaptation du filtre 10. S'il est trop petit la convergence n'est pas obtenue ou est trop lente. Au contraire, s'il est trop grand, l'adaptation est instable ou divergente. De manière connue, le coefficient de correction 2 est pris égal à un coefficient de convergence y constant, caractéristique du dispositif 1 et classiquement 15 réglé empiriquement au cours de tests. Cependant, et cette remarque est au coeur de l'invention, même si par hypothèse, la puissance du signal émis M est constante, du fait des perturbations occasionnées par le canal, la puissance du signal reçu en entrée X peut être variable, et présenter un intervalle de variation étendu. La transmission ne s'effectue pas à gain 20 unitaire. Si le signal d'entrée X est affecté d'un gain g, X g.X , le gradient de la fonction de coût J se voit affecté d'un gain g4, Vif g4 .VJ . Aussi le coefficient de correction 2 se trouve multiplié par g4. Afin d'éviter que cet effet préjudiciable menace la convergence de l'adaptation, selon une caractéristique importante, le coefficient de correction 2 est corrigé par multiplication par un coefficient de mise à l'échelle G, 25 inversement proportionnel à la puissance 4 du rayon de convergence r ou, ce qui équivalent, inversement proportionnel au carré de la puissance P. On peut alors écrire r4 p2 La méthode de GODARD préconise un critère et une fonction de coût J n f définie par la formule J = Ily(k)12 P)2 , avec J la fonction de coût, y(n) l'échantillon k=n-M 30 du signal de sortie Y à l'instant n, et P la puissance. Le nombre M+1 de termes de cette somme peut être important (de quelques dizaines à plusieurs centaines) et conduit à un nombre important de calculs. Aussi selon un mode de réalisation préférentiel, cette formule est simplifiée et remplacée par sa version instantanée J=),(012 - P)2 , limitée à l'instant n courant, avec J la fonction de coût, y(n) l'échantillon du signal de sortie Y à l'instant n, et P la puissance. Avec une telle fonction de coût J, le gradient de la fonction de coût J, s'écrit VJ = y(n)I2 P) y(n) X (n) , avec V l'opérateur de gradient, J la fonction de coût, y(n) l'échantillon du signal de sortie Y à l'instant n, P la puissance, X (n) le vecteur des échantillons x(k) du signal d'entrée X, k e pour les K+1 derniers instants, K étant la dimension du filtre 10, et z désignant le conjugué d'un nombre complexe z. On vérifie au passage que si le signal d'entrée X est affecté d'un gain g, le signal de sortie Y est affecté d'un facteur g, et la puissance est affectée d'un facteur g 2 . Aussi le gradient de la fonction de coût J, devient : VJ = y(n12 - g2 P) g y(n) g X (n) = gly(n12 - P)y(n) X (n) et donc VJ g4 VJ Une autre caractéristique contributive est l'introduction dans l'adaptation d'une 15 dispersion. Pour cela le dispositif 1 comprend encore, typiquement intégré dans l'adaptateur 12, un moyen de dispersion. Ce moyen de dispersion remplace dans la formule précédente du gradient le deuxième terme y(n) , échantillon du signal de sortie Y à l'instant n, par un terme y(n + q), échantillon du même signal de sortie Y, mais à l'instant n + q, avec q entier, positif ou négatif. q est aléatoirement choisi dans un 20 intervalle [-Q ; +Q], avec Q entier. Ainsi, y(n) échantillon du signal de sortie Y à l'instant n, est aléatoirement remplacé par y(n + q), échantillon du même signal de sortie Y, mais pris à un instant n + q, aléatoirement varié de manière centré autour de l'instant n et selon une étendue de variation Q. La formule du gradient devient VJ = y(n)I2-P)y(n+q)X(n). 25 L'étendue de variation Q est préférentiellement égale à 10. Il a été décrit précédemment un premier critère, issu du module constant de GODARD sous la forme : ),(n)12 - P)2 . Afin d'améliorer encore la convergence du filtre 10, il est ajouté un autre critère. Ce critère formalise une autre propriété d'une modulation à signal d'enveloppe constante à l'émission, à savoir que la projection sur 30 l'axe réel du produit du signal et de sa dérivée est nulle. Ce second critère s'exprime ReG)(n) (y(n)- y(n -1))). La fonction de coût J, qui incorpore nécessairement le premier critère, devient alors J y(n)12 - PY + ke(y(n)(y(n)- y(n - 1))))2 , avec J la fonction de coût, y(n) , respectivement y(n -1), l'échantillon du signal de sortie Y à l'instant n, respectivement à l'instant précédent n -1 , P la puissance, Re(z) désignant la partie réelle d'un nombre complexe z, et z désignant le conjugué d'un nombre complexe z. Parmi les filtres numériques linéaires, il est connu les filtres à réponse impulsionnelle infinie, RII (en anglais IIR) et les filtres à réponses impulsionnelle fini, RIF, (en anglais FIR). Le filtre 10 peut être un filtre de type RII ou RIF. La formulation non récursive d'un RIF est cependant plus simple. Aussi selon un mode de réalisation, le filtre 10 est modélisé par un filtre à réponse impulsionnelle finie, RIF. Un tel filtre RIF est tel que y(n) =Ia(k) x(n - k)= A(n) X (n) , avec y(n) k=o l'échantillon du signal de sortie Y à l'instant n, x(k) la eme composante, correspondant à l'instant k, du vecteur X (n) représentant le signal d'entrée X à l'instant n, k e - K , n} pour les K+1 derniers instants, K étant la dimension du filtre 10, a(k) le eme coefficient du filtre 10, eme composante du vecteur A(n) regroupant ces coefficients à l'instant n, k E 10, ..., , K étant la dimension du filtre 10. La dimension K du filtre 10 est un paramètre dont le choix résulte d'un compromis. Une réduction de la dimension K permet de réduire le nombre de calculs à réaliser. Au contraire une augmentation de La dimension K permet de considérer un intervalle de temps plus long. Aussi, en relation avec la fréquence d'échantillonnage, relativement au problème de réflexion multiple, une augmentation de la dimension K permet d'augmenter le retard maximal, ou ce qui est équivalent, le chemin de réflexion maximal, qui peut être pris en compte et être égalisé par le dispositif égalisateur 1. Selon une autre caractéristique, le dispositif égalisateur 1 comprend encore un contrôleur 13. Ce contrôleur 13 a pour fonction de comparer les coefficients A(n) du filtre 10, avec deux seuils minimum et maximum. Ainsi pour chaque coefficient a(k) , k e 10, ..., Kl, le contrôleur 13 vérifie qu'il est compris entre un seuil minimum a min(k) et un seuil maximum a max(k), selon l'expression a min(k) Si tous les coefficients a(k) , k E 10, restent dans l'intervalle [a min(k) ; a max(k)1, le contrôleur 13 valide les coefficients a(k) et autorise le filtre 10 à les appliquer. Au contraire, si l'un au moins des coefficients a(k) sort de l'intervalle [a min(k); a max(k)1, le contrôleur 13 provoque une réinitialisation du filtre 10. Dans ce dernier cas, similaire à l'initialisation du filtre 10 lors de son démarrage, tous les coefficients a(k) sont annulés à l'exception du coefficient a(0) appliqué à x(n) qui est pris égal à 1, soit [a(0), a(k), a(K)]= 0, .., 0]. Le contrôleur 13 est typiquement disposé, entre l'adaptateur 12 et le filtre 10, de manière à contrôler les coefficients A(n) après leur calcul par l'adaptateur 12 et à ne transmettre au filtre 10 les coefficients A(n) calculés, que lorsqu'ils sont valides. Selon un mode de réalisation préféré, chaque seuil minimum a min(k) est 10 opposé et égal en valeur absolue au seuil maximum a max(k) correspondant, selon l'expression a mie= -a max(k), Vk E 10, ..., Kl. Il a été vu précédemment qu'une caractéristique importante consistait à corriger le coefficient de correction 2 en lui appliquant un coefficient de mise à l'échelle G s inversement proportionnel au carré de la puissance P. Une telle opération 15 d'inversion, si elle peut être effectuée en un cycle sur un processeur à virgule flottante, n'est en général pas disponible sur un processeur à virgule fixe. Afin de pouvoir approcher ce calcul, de manière réalisable par un processeur à virgule fixe, le dispositif 1 comprend encore un approximateur, apte à déterminer une approximation du coefficient de mise à l'échelle G s . 20 Le principe de cet approximateur est illustré en relation avec la figure 3. Il consiste à partager l'étendue ou intervalle de variation de la puissance P, représentée ici par l'axe des abscisses, en un nombre S de zones Zi, i = 1 S, avec S = 6 sur la figure 3. Il s'agit d'un partage ou partition, en ce que les zones Zi sont toutes disjointes, et que l'union de toutes les zones Zi couvre la totalité de l'intervalle de puissance P. Le 25 découpage des zones Zi est réalisé de telle manière que la puissance P d'une zone Zi varie d'un facteur 4 relativement à la zone précédente Z(i-1), ou, ce qui est équivalent, que le signal X ou encore le rayon de convergence r varie d'un facteur 2. Il s'ensuit que d'une zone Zi à la zone Z(i+1) immédiatement suivante, le coefficient de mise à l'échelle G s , proportionnel à -1, est divisé par 16, ou 24. Ainsi, si l'on note bi, i = 0 S, p2 30 les bornes entre les zones Zi, chaque zone Zi est encadrée par les bornes b(i-1) et bi, ou encore la zone Zi est définie par l'intervalle [b(i-1) ; bi]. Les bornes bi vérifient la relation bi = 4 x b(i-1), i = 1 ... S. Sur chacune des zones Zi, le coefficient de mise à l'échelle Gs est modélisé par une fonction approchante FAi. La fonction FAi associée à chaque zone Zi est avantageusement, en factorisant un facteur multiplicatif 24 cumulé, une fonction simple, telle une fonction polynomiale de degré d faible, ainsi définie par un nombre restreint d+1 5 de coefficients. La détermination de ces coefficients est obtenue en écrivant toutes les relations et les conditions aux limites, notamment que le coefficient de mise à l'échelle Gs découpé en fonction approchantes FAi par paliers est continu. Une fonction approchante FAi est ainsi associée à chaque zone Zi. La condition de continuité s'exprime par une égalité entre la valeur de la fonction approchante FAi de la zone Zi en bi et la 10 valeur de la fonction FA(i+1) de la zone suivante Z(i+1) au niveau de leur borne commune bi. Ainsi, sur chaque zone Zi, le calcul du coefficient de mise à l'échelle Gs , proportionnel à -1, est approximé par l'application de la fonction FAi(P) correspondante, pz ce qui donne lieu à des calculs plus simples et surtout réalisables par un processeur à 15 virgule constante. Il a été vu que le degré du polynôme des fonctions approchantes FAi est avantageusement réduit. Dans la pratique, il est pris au plus égal à 2, avec des polynômes des fonctions approchantes FAi quadratiques, définis par trois coefficients, ou encore préférentiellement égal à 1, avec des polynômes des fonctions approchantes FAi 20 affines de la puissance P, ainsi définis par deux coefficients a et p. Ainsi, dans ce dernier cas, une fonction approchante FAi est définie par une formule FAi(P)-(aP+13)28 , avec un facteur 28 introduisant, à chaque zone Zi, de gauche à droite, une division par 16, soit une division par 24, cumulée à chaque zone Zi, soit Ô - 4(i - S), avec i l'indice de la zone Zi et S le nombre total de zones, soit 25 FAi(P)-(aP+ )6)24(t-s) . On obtient alors le coefficient de mise à l'échelle Gs en appliquant pour chaque zone Zi, la fonction FAi associée à la zone Zi de la puissance P, soit Gs - FAi(P)-(aP+ )3)24(t-s) . Cette approximation, permet avantageusement de remplacer, par une fonction 30 affine, un calcul de 1- simple à évaluer par un processeur à virgule fixe et par une pz multiplication par une puissance de 2, simplement réalisée par un décalage binaire. Avantageusement encore, une mise en facteur par une puissance de 2, choisie pour chaque zone Zi, est appliquée, selon la formule aP+13=2:1 -P+-, de telle manière à ce que les nouveaux 217 2/7) 217 217 coefficients, a' et [3', soient tous les deux entiers. Afin de mieux détailler l'approximation, un exemple numérique est décrit ci-dessous, en relation avec le tableau de la figure 4. Il est supposé que la puissance P, varie dans un intervalle [4,77.10-7 ; 0,03125]. On se propose de partager cet intervalle comme décrit précédemment, ce qui résulte en S = 8 zones Z1 ... Z8, délimitées par 9 bornes b0 b8. Le tableau de la figure 4 donne pour chaque zone Zi, sa borne inférieure b(i-1) en 2ème colonne et sa borne supérieure bi en 3ème colonne. Pour chaque zone Zi, sont déterminés les deux coefficients a et [3 de la 10 fonction affine et Ô la puissance de 2, qui définissent la fonction FAi. a est égal à 15 et est noté dans la 4ème colonne. [3 est égal à -21 et est noté dans la 5ème 16 (b, -k_1) 16 colonne. Ô est égal à 4(i - S) et est noté dans la 6ème colonne. Le tableau comprend encore, les coefficients a' et [3' entiers, notés respectivement en 8ème et 9ème colonnes, obtenus par multiplication par 2, 1-1 étant noté 15 en 7ème colonne. Les puissances de 2, net Ô se simplifient en un unique décalage binaire. La figure 5 illustre l'erreur introduite par l'approximation du coefficient de mise à l'échelle G. Sur un intervalle [1 ; 2] figurant le domaine de variation d'une variable x, indicatif d'une zone Zi, dont la borne supérieure est double de la borne inférieure, sont représentées cinq courbes 21-25. La courbe 21 représente la droite identité y = x, ici 20 homogène au rayon de convergence r ou encore au signal d'entrée X. La courbe 22 représente la parabole y = x2, homogène à la puissance P, et dont la borne supérieure, égale à 4, est dans un rapport 4 avec la borne inférieure, égale à 1. La courbe 23 représente la courbe y = x4, homogène au carré de la puissance P. La courbe 24 représente le facteur de mise à l'échelle Gs approximé par une fonction affine de la 25 puissance P ou de x2, entre la borne inférieure 1 où elle prend comme valeur 1 et la borne supérieure 4 où elle prend comme valeur 1/16, soit la courbe y = -5/16x2 + 21/16 homogène à -5/16P + 21/16. La courbe 25 représente le produit du facteur de mise à l'échelle Gs approximé et de x4, soit la courbe y = Gs x4, homogène à y= Gs P2. Puisque le facteur de mise à l'échelle Gs approché est une approximation de 30 1 , la qualité de l'approximation est d'autant meilleure que le produit Gs P2 est proche p2 de 1. Il apparait ainsi que l'approximation est parfaite, par construction, aux bornes bi des zones Zi, et présente une variation d'un facteur présentant un maximum de l'ordre de 3,5. En l'absence de correction par le facteur de mise à l'échelle Gs, le coefficient de correction 2 présenterait une erreur de l'ordre de l'excursion de puissance, soit b8 / b0, supérieure à 65000 dans l'application numérique. Aussi une erreur d'au maximum 3,5 introduite par l'approximation, est-elle tout à fait acceptable. La puissance P est utilisée par le dispositif égalisateur 1 pour calculer le coefficient de mise à l'échelle Gs (selon la méthode exacte ou selon l'approximation) et pour calculer le premier critère. Aussi un estimateur 11 est utilisé pour estimer la puissance P du signal d'entrée X. Selon un mode de réalisation, cette puissance P est déterminée, par une moyenne sur une fenêtre temporelle glissante et est affectée d'un /0 gain fixe y. Ce gain y est configurable afin de pouvoir être initialisé en usine. La n puissance P est déterminée par la formule P = y Ilx(k)I2 , avec P la puissance, y le k=n-N gain fixe configurable, x(k) l'échantillon du signal d'entrée X à l'instant k, k e In - N,...,n}, N étant la taille de la fenêtre temporelle glissante. L'invention concerne encore un récepteur radiofréquence comprenant un tel 15 dispositif 1 égalisateur de canal selon l'un quelconque des modes de réalisation décrits. REVENDICATIONS1. Dispositif (1) égalisateur aveugle de canal, pour un récepteur radiofréquence adapté à une modulation à signal d'enveloppe constante à l'émission, comprenant : - un filtre (10) numérique linéaire adaptable, défini à un instant ( n ) par ses coefficients ( A(n)), apte à filtrer un signal d'entrée (X) pour produire un signal de sortie (Y ), - un estimateur (11) apte à estimer une puissance ( P ) du signal d'entrée ( X ), - un adaptateur (12) apte à adapter le filtre (10) en calculant, les coefficients (A(n)) du filtre (10) à un instant (n), en soustrayant aux coefficients (A(n -1)) du filtre (10) à un instant précédent (n-1) le gradient d'une fonction de coût (J) affecté d'un coefficient de correction (2), où la fonction de coût ( J ) comprend un premier critère de distance entre le carré du signal de sortie (Y) et la puissance (P), où le coefficient de correction ( ) est un produit comprenant un coefficient de convergence (p) constant, caractérisé en ce que le coefficient de correction (2) /5 comprend encore un coefficient de mise à l'échelle (Gs) inversement proportionnel au carré de la puissance (P ). 2. Dispositif (1) selon la revendication 1, où la fonction de coût ( J ) est défini(e) par la formule J =),(n12 - , avec J la fonction de coût, y(n) l'échantillon du signal de sortie ( Y ) à l'instant n, et P la puissance. 20 3. Dispositif (1) selon la revendication 2, comprenant encore un moyen de dispersion, remplaçant le gradient de la fonction de coût ( J ), défini par la formule Vif y(n)12 P) y (n) X (n) , avec V l'opérateur de gradient, J la fonction de coût, y(n) l'échantillon du signal de sortie (Y) à l'instant n, P la puissance, X (n) le vecteur des échantillons x(k) du signal d'entrée (X), ke N,...,n}, K étant la dimension du 25 filtre (10), et z désignant le conjugué d'un nombre complexe z, par la formule VJ = y(n)I2 P) y(n + q) X (n) , avec q entier aléatoire, tel que - Q 4. Dispositif (1) selon l'une quelconque des revendications 1 à 3, où la fonction de coût ( J ) comprend encore un deuxième critère de nullité de la projection sur l'axe réeldu produit du signal de sortie (Y) et de sa dérivée, la fonction de coût étant alors définie par la formule =)(11)12 P)2 ±(Re(3n)(Y(n)- Y(n -i))))2 avec J la fonction de coût, y(n) , respectivement y(n -1), l'échantillon du signal de 5 sortie (Y) à l'instant n, respectivement à l'instant précédent n -1, P la puissance, Re(z) désignant la partie réelle d'un nombre complexe z et z désignant le conjugué d'un nombre complexe z. 5. Dispositif (1) selon l'une quelconque des revendications 1 à4, où le filtre (10) est un filtre à réponse impulsionnelle finie, RIF. /0 6. Dispositif (1) selon l'une quelconque des revendications 1 à 5, comprenant encore un contrôleur (13), apte à vérifier après un calcul des coefficients (A(n)) du filtre (10), que chaque coefficient est encadré entre un seuil minimum et un seuil maximum, avec préférentiellement chaque seuil minimum égal à l'opposé du seuil maximum correspondant, ledit contrôleur (13) étant encore apte à réinitialiser les /5 coefficients ( A(n)) du filtre (10) dans le cas contraire. 7. Dispositif (1) selon l'une quelconque des revendications 1 à 6, comprenant encore un approximateur apte à déterminer une approximation du coefficient de mise à l'échelle (G,), comprenant un moyen de partage apte à partager l'étendue de puissance (P) du signal d'entrée (X) en une pluralité de zones (Zi), telles que la 20 puissance (P) varie sensiblement d'un facteur 4 entre deux zones (Zi) contiguës, un moyen d'association apte à associer à chaque zone (Zi) une fonction approchante (FAi) de calcul du coefficient de mise à l'échelle (Gs) sur ladite zone (Zi) en fonction de la puissance (P), sous forme d'un polynôme divisé par une puissance de 2 quadruple d'un indice correspondant à la zone (Zi), en continuité avec la fonction approchante (FAi) 25 associée à une zone (Zi) contiguë, et un moyen d'application apte à appliquer la fonction approchante (FAi) associée à une zone (Zi) correspondant à la puissance (P) pour calculer une approximation du coefficient de mise à l'échelle (Gs). 8. Dispositif (1) selon la revendication 7, où un polynôme d'une fonction approchante (FAi) est une fonction affine de la puissance ( P ). 30 9. Dispositif (1) selon l'une quelconque des revendications 1 à 8, où l'estimateur (11) détermine la puissance (P) du signal d'entrée (X), par une moyenne sur une fenêtre temporelle glissante, affectée d'un gain fixe ( y) configurable. 10. Récepteur radiofréquence caractérisé en ce qu'il comprend un dispositif (1) égalisateur de canal selon l'une quelconque des revendications précédentes.