FR 3022081 A1 20151211 FR 1459469 A 20141003 La présente invention est relative aux systèmes audio pour véhicule et particulièrement aux procédés de filtrage pour le filtrage sélectif dynamique de canaux de diffusion en modulation de fréquence (FM). Plus précisément, on s'intéresse aux systèmes audio capables de recevoir et 5 démoduler des ondes électromagnétiques diffusant des programmes radiophoniques sur la modulation de fréquence. Sur la bande de la modulation de fréquence, on peut trouver couramment des canaux utiles séparés de 200 kHz ou 300 kHz en terme de fréquence porteuse, mais il arrive aussi que deux canaux adjacents aient des fréquences porteuses séparées 10 seulement par 100 kHz, ce qui pose un problème de discrimination et de filtrage sélectif pour focaliser la réception sur un des canaux et rejeter les signaux de l'autre canal. Dans les systèmes connus, on applique un filtre d'entrée plus ou moins étroit en fonction des conditions de réception qui sont estimées après démodulation, c'est-à-dire sur la qualité audio ; la largeur du filtre est donc pilotée par une boucle de rétroaction. 15 Mais il se trouve que cette boucle de rétroaction peut introduire des instabilités inutiles, voire néfastes, dans le système, notamment en cas de changements de contenu FM et/ou de présence de canaux adjacents ou divers perturbateurs électromagnétiques et/ou de changements de conditions de réception. Il est donc apparu un besoin de proposer une solution plus performante pour 20 filtrer sélectivement les signaux d'entrée, afin d'assurer une réception optimisée du canal d'intérêt même avec un canal adjacent proche. Selon l'invention, il est proposé un procédé de filtrage sélectif dynamique d'ondes électromagnétiques diffusant des programmes radiophoniques sur la modulation de fréquence, le procédé comprenant les étapes : 25 /a/ sélectionner un premier canal fréquentiel (F0) portant des signaux d'une station d'intérêt destinée à être écoutée, /b1/ appliquer aux signaux reçus un filtre d'entrée centré sur le premier canal fréquentiel, opérer une transposition en bande de base, et numériser le signal filtré, 30 /b2/ appliquer une transformée de Fourier et rechercher dans la courbe formant le résultat de la transformée de Fourier, une fréquence de croisement (FC), définie comme l'abscisse du minimum local de la courbe le plus proche de la porteuse centrale d'intérêt, /c/ appliquer aux signaux numérisés un filtre sélectif asymétrique de part et 35 d'autre du premier canal fréquentiel, le filtre ayant une de ses deux bornes (FW) déduite de la fréquence de croisement, à savoir FW = FC +/- FH, FH étant une valeur de correction comprise entre 0 kHz et 40 kHz, l'autre de ses bornes étant laissée à une valeur par défaut plus éloignée de F0, /d/ répéter continuellement les étapes /b1/ /b2/ et /ci, jusqu'à un changement de canal fréquentiel. La détection d'un tel minimum local permet d'identifier le passage du canal d'intérêt à un autre canal adjacent perturbateur, et permet de déterminer si ce canal adjacent est au-dessus ou au-dessous du canal d'intérêt. Moyennant quoi on peut appliquer un filtrage asymétrique pour éliminer d'un côté ce canal adjacent perturbateur et conserver une largeur suffisante de l'autre côté. On obtient ainsi des signaux audio démodulés qui reproduisent les signaux diffusés par le premier canal d'intérêt en annulant ou limitant fortement les perturbations induites par la présence d'un second canal adjacent perturbateur, ou la présence d'interférences en dehors du contenu audio. On note aussi que, grâce à la transposition vers la bande de base, tous les 15 traitements numériques se font avantageusement dans la bande de base, assez facilement. Dans divers modes de réalisation de la méthode selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : - Le filtre numérique asymétrique est avantageusement un filtre complexe, à savoir 20 à coefficients complexes, traitant simultanément une partie réelle et une partie imaginaire de l'enveloppe complexe des signaux numérisés. Moyennant quoi, en fonction des coefficients réel et imaginaire du filtre on peut rendre le filtre asymétrique par rapport à la fréquence centrale d'intérêt F0. - Le filtre à coefficients complexes permet de définir une borne supérieure du filtre 25 et une borne inférieure du filtre, lesdites bornes étant définies indépendamment l'une de l'autre et non pas nécessairement symétriques par rapport à la fréquence d'intérêt. On peut ainsi ajuster indépendamment la borne inférieure et la borne supérieure du filtre numérique de manière à réagir sélectivement à la présence d'un canal perturbateur d'un seul des deux côtés de la fréquence 30 d'intérêt. - La détection d'un tel minimum local peut être symptomatique de la présence d'un second canal perturbateur ayant une fréquence porteuse décalée de 100 kHz par rapport à la fréquence porteuse du canal d'intérêt ; moyennant quoi le procédé proposé permet de discriminer deux canaux très proches, de façon plus efficace 35 que les systèmes de l'art antérieur, notamment à rétroaction. - La numérisation et la transformée de Fourier sont de préférence répétées au moins toutes les 500 ps. On dispose ainsi de moyens de décision rapides de manière à ajuster la décision de façon très dynamique. - Le filtre d'entrée présente une largeur d'environ 300 kHz. Ceci permet d'éliminer avantageusement en amont tous les canaux dont la porteuse est séparée de 200 kHz ou plus de la porteuse de la fréquence d'intérêt. - La valeur de la borne contrainte du filtre numérique asymétrique est réajustée à une cadence identique à la fréquence de récurrence des étapes /b1/, /b2/ et /c/, de préférence supérieure à 1 kHz. De sorte que le filtre peut être resserré si le contenu audio présente peu d'excursion ou si un canal adjacent est détecté et le filtre peut être élargi si le contenu audio est plus fourni et sans canal adjacent potentiellement perturbateur. - Le minimum local le plus proche est de préférence déterminé en calculant la dérivée de la courbe formant la transformée de Fourier et en identifiant les points où cette dérivée s'annule et en retenant le point le plus proche de la porteuse du canal d'intérêt. Ce qui représente un traitement mathématique bien maitrisé et facile à mettre en oeuvre. - La valeur de correction FH peut être choisie en fonction de la valeur de la dérivée seconde au point du minimum local. On peut ainsi ajuster la fréquence de coupure du filtre numérique en fonction de la puissance du canal adjacent perturbateur, qui est d'autant plus grande que la dérivée seconde est élevée. - On peut prévoir un lissage temporel et/ou un filtre passe-bas numérique concernant le mécanisme d'ajustement de la borne contrainte du filtre ; moyennant quoi l'ajustement, bien que rapide, ne comporte pas de saut. L'invention vise aussi un système audio capable de recevoir et de démoduler des ondes électromagnétiques diffusant des programmes radiophoniques sur la modulation de fréquence, remarquable en ce qu'il est configuré pour mettre en oeuvre le procédé selon l'une quelconque des caractéristiques précédentes ; un tel système présente les avantages énoncés ci-dessus. D'autres aspects, buts et avantages de l'invention apparaîtront à la lecture de la description suivante d'un de ses modes de réalisation, donné à titre d'exemple non limitatif. L'invention sera également mieux comprise en regard des dessins joints sur lesquels : - la figure 1 représente schématiquement un spectre électromagnétique partiel dans la bande FM et son évolution au cours du temps, avec un canal d'intérêt et un canal adjacent potentiellement perturbateur, - la figure 2 représente un schéma bloc d'un système audio selon l'invention, - la figure 3 représente schématiquement un spectre électromagnétique, transposé de celui de la figure 1 vers la bande de base, - la figure 4 illustre différents cas de figure de la courbe spectrale issue de la transformée de Fourrier, notamment ses variations au cours du temps, - la figure 5 représente schématiquement les contenus audio respectifs du canal d'intérêt et des canaux adjacents potentiellement perturbateurs, au cours du temps, - la figure 6 illustre un cas de figure où le filtre à coefficients complexes opère un filtrage asymétrique. Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. Les figures 1 et 2 illustrent un système audio 1 pour véhicule, aussi appelé autoradio, capable de recevoir et démoduler des ondes électromagnétiques diffusant des programmes radiophoniques sur la bande FM (« modulation de fréquence »). S'agissant des programmes radiophoniques, comme connu en soi, le principe de la modulation de fréquence est basé sur l'utilisation d'une excursion de fréquence par rapport à une fréquence porteuse (F0) pour encoder les signaux audio à diffuser. L'excursion de fréquence est proportionnelle au signal audio encodeur. Une pluralité de canaux peut être diffusée sur la totalité de la bande FM (88 MHz - 108 MHz dans beaucoup de pays, 76 MHz - 90 MHz au Japon). Chaque canal présente typiquement une largeur de bande utile d'environ 240 kHz, c'est-à-dire +/-120 kHz autour de la porteuse ; ceci correspond à un spectre audio composite de départ (bande de base Mono + Stéréo + RDS) borné par 60 kHz, modulant la porteuse en fonction de son amplitude et de son contenu spectral pour donner un canal d'au plus 240 kHz de large. On comprend dès lors que, lorsque deux canaux ne sont séparés, en termes de fréquence porteuse F0, que par 100 kHz, il peut y avoir un léger chevauchement des 30 signaux sur le spectre, problème auquel la présente invention s'attache à apporter une solution. Dans la suite du présent document, la fréquence FO désignera la fréquence porteuse d'un premier canal, aussi appelé canal d'intérêt, que l'utilisateur du système audio a sélectionné, cet utilisateur souhaitant écouter le programme radiophonique diffusé 35 par ce premier canal à fréquence porteuse F0. On remarque que le contenu spectral de ce premier canal évolue avec le temps, principalement en raison du contenu audio modulant qui évolue lui aussi avec le temps. Le spectre n'est pas toujours centré sur F0, et notamment le lieu du maximum 71 se déplace autour de FO en permanence en fonction du contenu audio. La figure 1 illustre l'évolution du contenu spectral du canal d'intérêt de porteuse FO ainsi que le contenu spectral d'un second canal F2, que l'on appellera ici 5 canal perturbateur. Le second canal présente une fréquence porteuse F2 peu éloignée de la fréquence porteuse F0, par exemple l'écart peut être de 100 kHz. Le canal perturbateur peut être situé soit au-dessus (comme illustré), soit au-dessous de la fréquence du premier canal F0. Le contenu spectral du second canal évolue lui aussi avec le temps de manière indépendante du contenu spectral du premier canal. 10 Sur la figure 2, le système audio 1 comprend une partie frontale analogique 2, une partie centrale numérique 4 avec une unité de commande logique 3 (coeur de commande logique), un démodulateur 7 et une partie aval analogique avec un amplificateur 19 et des haut-parleurs 9. Par ailleurs, le système audio comprend une interface utilisateur 6 avec un écran d'affichage et une surface tactile ou des boutons 15 physiques (non représentés) ; ce qui permet à l'utilisateur de choisir la station/le canal d'intérêt. Les signaux électromagnétiques sont reçus par un dispositif d'antenne 10, les signaux sont ensuite filtrés par un filtre d'entrée 11. Les signaux issus du filtrage sont ensuite numérisés par un convertisseur analogique digital 12. 20 La poursuite du traitement du signal est réalisée par des opérations numériques, notamment une mise en mémoire de données dans une mémoire tampon 13. Avantageusement, l'opération de numérisation est effectuée en bande de base et comporte une opération de détermination de l'enveloppe complexe du signal 25 modulé ramené en bande de base par la transposition -F0. Cette enveloppe complexe comporte une partie de signaux réels notée I (composante en phase, « In phase » en anglais) et une partie de signaux imaginaires notée Q (composante en « Quadrature ») ; la partie de signaux imaginaires Q étant l'image de la portion de spectre située à gauche de 0 ; la partie de signaux réels I est l'image de la portion de spectre située à droite de 0. 30 On peut représenter l'enveloppe complexe du signal modulé ramené en bande de base par X = I + jQ, j étant le nombre complexe dont le carré vaut -1. Un filtre complexe 15, à savoir à coefficients complexes, peut s'écrire sous la forme F = C + jD. L'opération de filtrage complexe consiste à faire le produit du signal enveloppe 35 par le filtre à savoir : Y= Fx X= (I + jQ) x (C + jD) = IC - QD + j(ID + QC) À partir des données mises en mémoire tampon (I et Q), il est procédé à une analyse de transformée de Fourrier 14 qui donne en sortie une courbe de densité spectrale de puissance dont il sera question en détail plus loin. Par ailleurs, les données mises en mémoire tampon font l'objet d'un filtrage numérique par le filtre numérique 15 ayant des bornes programmables grâce au paramétrage du filtre numérique complexe. Les valeurs respectives des borne inférieure FB1 et borne supérieure FB2 sont déterminées par les valeurs C et D filtre 15 déterminées en temps réel par la détection ou non d'un canal perturbateur adjacent. Un bloc de contrôle 16 du filtrage numérique est en particulier en charge de 10 commander, en fonction de l'analyse du contenu de la courbe de densité spectrale de puissance, la valeur des coefficients C et D du filtre numérique complexe, desquelles résultent les bornes inférieure FB1 et supérieure FB2, comme expliqué ci-dessus. Tout ou partie des opérations de mise en mémoire tampon, transformée de Fourrier, filtrage numérique et démodulation peuvent être effectués par un circuit dédié de 15 type DSP 5 (« DSP » étant la contraction du terme anglais « Digital Signal Processor », c'est-à-dire un processeur de signal numérique). On remarque que le filtrage d'entrée 11 présente une largeur de l'ordre de 200 à 300 kHz centré sur F0, ce qui permet d'éliminer les signaux provenant des canaux éloignés de F0, mais ne permet pas d'éliminer correctement les signaux d'un canal 20 adjacent proche dont la densité spectrale peut déborder à l'intérieur du filtre d'entrée. Comme illustré à la figure 3, dans le cas de figure où le second canal F2 potentiellement perturbateur présente une porteuse décalée à 100 kHz de F0, il existe un point de croisement noté FC au niveau duquel la densité spectrale de puissance de chaque canal se trouve au même niveau, à l'endroit où les courbes se croisent. 25 Les signaux en aval du filtrage d'entrée, après numérisation et transposition spectrale, puis traitement par transformée de Fourrier sont illustrés à la figure 4. Il faut comprendre que concernant les figures 3 et 4, elles illustrent le traitement en bande de base c'est-à-dire après transposition de fréquence de -F0. L'effet de la présence du second canal perturbateur F2 se manifeste par une 30 remontée locale de la densité spectrale à droite du point noté « FC+ » référence 92b. L'analyse de la courbe de densité spectrale en sortie du bloc 14 permet d'identifier des minimums locaux, et en particulier d'identifier le minimum local le plus proche de la fréquence centrale d'intérêt. Grâce à la détermination de FC, le bloc de commande du filtre 16 applique un 35 filtre numérique sélectif avec des valeurs de coefficients C et D en conséquence. L'une des bornes inférieure FB1 et supérieure FB2 peut être contrainte FW (plus proche de F0), tandis que l'autre borne est laissée à une valeur par défaut plus éloignée de FO (cf. figure 6). Avantageusement le filtrage numérique est asymétrique, car dans la majorité des cas pratiques, il n'y a qu'un seul canal perturbateur potentiel à proximité du canal FO ; 5 ce perturbateur peut être à gauche (fréquence inférieure) à droite (fréquence supérieure) par rapport à F0. On peut choisir une borne du côté du signal perturbateur par exemple égale à FC +/- FH, FH étant une valeur de correction comprise entre 0 kHz et 40 kHz, FH pouvant être un décalage paramétré faisant l'objet d'un mécanisme d'apprentissage, et/ou de 10 corrélation par rapport à la valeur de la dérivée seconde de la courbe spectrale issue de la transformée de Fourier. Si aucun minimum n'est identifié, comme par exemple en cas de la courbe 93 sur la figure 4, les coefficients complexes C, D du filtre et par conséquent les bornes du filtre numérique pourront alors être déduite du contenu spectral du canal d'intérêt, comme 15 cela sera détaillé plus loin. Sur la figure 4, en trait mixte, dans le cas de la courbe 91, le minimum local le plus proche de la fréquence centrale noté FC- référence 91b, se situe à gauche de cette dernière. Dans ce cas, la courbe présente également un maximum local c'est-à-dire une bosse 91a symptomatique de la présence d'un canal adjacent à gauche (cf. figure 5). 20 De même, dans le cas de courbe 92, on trouve un maximum local 92a à droite du minimum local 92b. La détection d'un minimum local est obtenue en calculant la dérivée de la courbe spectrale et en recherchant les racines c'est-à-dire les lieux où cette dérivée s'annule. De plus, on calcule la dérivée seconde, un minimum étant caractérisé par une 25 valeur positive de la dérivée seconde à l'endroit où la dérivée primaire s'annule. On remarque que généralement la hauteur de la bosse latérale sera liée à la valeur de la dérivée seconde, et la valeur de correction FH peut être obtenue comme une fonction de la valeur de la dérivée seconde à l'endroit le plus proche de la fréquence. Le mécanisme décrit plus haut est sans cesse répété ce qui conduit à une 30 prise de décision dynamique. La fréquence de récurrence est au moins égale à 1 kHz de préférence. Typiquement, la numérisation et la transformée de Fourier sont répétées au moins toutes les 500 ps. Comme ceci est illustré à la figure 5, le contenu spectral 70 du canal 35 d'intérêt FO évolue avec le temps, ainsi que le contenu spectral 72 d'un éventuel second canal, perturbateur F2 représenté en traits pointillés en partie basse de la figure, et le contenu spectral 74 d'un éventuel troisième canal, perturbateur, F3 représenté en traits pointillés sur la partie haute de la figure. Les références 71, 73 et 75 désignent respectivement le lieu des fréquences de densité spectrale maximum qui évoluent elles aussi avec le temps. Pendant la première période de temps illustrée 81, la courbe de densité 5 spectrale issue de la transformée de Fourrier ne révèle pas de minimum et par conséquent le filtrage numérique est opéré avec une largeur et des bornes [FB1 - FB2] de valeur par défaut par exemple 100 kHz de large. Au début de la seconde période de temps illustrée 82, l'analyse de la courbe spectrale permet d'identifier un minimum local FC induit par la présence du canal 10 adjacent F2 à droite dont le contenu spectral s'est rapproché de la fréquence centrale d'intérêt F0. La borne supérieure FB2 du filtre numérique 15 est alors réduite progressivement vers une valeur contrainte FW plus proche de FO au fur et à mesure que l'interférence entre les deux canaux augmente. Pendant ce temps, la borne inférieure FB1 15 ne change pas ; c'est donc uniquement le filtre 151 qui voit sa borne ajustée par le bloc de commande 16. Au début de la troisième période de temps illustrée 83, la borne supérieure FB2 du filtre numérique 15 est à nouveau augmentée au fur et à mesure que l'interférence entre les deux canaux F0, F2 diminue. 20 A la quatrième période de temps illustrée 84, en l'absence de détection d'un minimum local dans la courbe spectrale, on retrouve progressivement les bornes inférieure FB1 et supérieure FB2 par défaut. Dans l'exemple illustré, le second canal F2 disparaît. Au cours d'une période de temps ultérieure 85, il apparaît un autre canal 25 perturbateur baptisé troisième canal, de fréquence porteuse F3 située juste en dessous par exemple à 100 kHz de la fréquence d'intérêt F0. Pendant la période de temps 85, la courbe de densité spectrale issue de la transformée de Fourrier ne révèle pas de minimum et par conséquent le filtrage numérique est opéré avec des bornes [FB1 - FB2] de valeur par défaut. 30 Au début de la période de temps 86, l'analyse de la courbe spectrale permet d'identifier un minimum local FC induit par la présence du canal adjacent F3 à gauche dont le contenu spectral s'est rapproché de la fréquence centrale d'intérêt F0. La borne inférieure FB1 du filtre numérique 15 est alors augmentée vers une valeur contrainte FW plus proche de FO au fur et à mesure que l'interférence entre les 35 deux canaux F0, F3 augmente. Pendant ce temps, la borne supérieure FB2 ne change pas ; c'est donc uniquement le filtre 15Q qui voit sa borne ajustée. Au début de la période de temps 86, le point FC s'éloigne de FO et par conséquent la borne inférieure FB1 diminue et tend à revenir progressivement vers sa valeur par défaut. Pendant les périodes de temps 88, 89, il est à nouveau détecté un point de 5 minimum FC à gauche, et la borne inférieure FB1 du filtre est de nouveau augmentée pour exclure les signaux du troisième canal F3. Au cours de la période de temps suivante 8A, on retrouve une situation sans détection de minimum local dans le spectre similaire la borne inférieure FB1 revient, de préférence avec lissage, vers sa valeur par défaut. 10 II faut noter que les valeurs par défaut des bornes inférieure FB1 et supérieure FB2 peuvent être calculées en dépendance du contenu spectral du canal d'intérêt F0, la largeur dépendant de l'étalement du spectre du canal d'intérêt. REVENDICATIONS1. Procédé de filtrage sélectif dynamique d'ondes électromagnétiques diffusant des programmes radiophoniques sur la modulation de fréquence (FM), le procédé comprenant les étapes : /a/ sélectionner un premier canal fréquentiel (F0) portant des signaux d'une station d'intérêt destinée à être écoutée, /b1/ appliquer aux signaux reçus un filtre d'entrée (11) centré sur le premier canal fréquentiel (F0), opérer une transposition en bande de base, et numériser (12) le signal filtré, /b2/ appliquer une transformée de Fourier (14), et rechercher dans la courbe formant le résultat de la transformée de Fourier, une fréquence de croisement (FC), définie comme l'abscisse du minimum local de la courbe le plus proche de la porteuse centrale d'intérêt, /c/ appliquer aux signaux numérisés un filtre sélectif asymétrique (15) de part et d'autre de la fréquence du canal fréquentiel (F0), le filtre ayant une de ses deux bornes (FW) déduite de la fréquence de croisement (FC), à savoir FW = FC +/- FH, FH étant une valeur de correction comprise entre 0 kHz et 40 kHz, l'autre de ses deux bornes étant laissée à une valeur par défaut plus éloignée de la fréquence du canal fréquentiel (F0), /d/ répéter continuellement les étapes /b1/ /b2/ et /c/ jusqu'à un changement de canal fréquentiel, ce qui permet d'identifier le passage du canal d'intérêt à un autre canal adjacent perturbateur ; de sorte que les signaux audio démodulés reproduisent les signaux diffusés par le premier canal en annulant ou limitant fortement les perturbations induites par la présence d'un second canal adjacent situé sur un des côtés, ou d'interférence en dehors du contenu audio. 2. Procédé selon la revendication 1, dans lequel le filtre numérique asymétrique est un filtre complexe, à savoir à coefficients complexes, traitant simultanément une partie réelle (I) et une partie imaginaire (Q) de l'enveloppe complexe des signaux numérisés. 3. Procédé selon la revendication 1, dans lequel le filtre numérique (15) est de type complexe à coefficients complexes, permettant de définir une borne supérieure du filtre (FB2), et une borne inférieure du filtre (FB1), lesdites bornes étant définies indépendamment l'une de l'autre et non pas nécessairement symétriques par rapport à la fréquence du canal fréquentiel (F0). 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel la détection d'un minimum local est symptomatique de la présence d'un second canal perturbateur ayant une fréquence porteuse (F2) décalée de 100 kHz par rapport à la fréquence porteuse (F0) du canal d'intérêt. 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel la numérisation et la transformée de Fourier sont répétées au moins toutes les 500 ps. 6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel le filtre d'entrée (11) présente une largeur d'environ 300 kHz. 7. Procédé selon l'une quelconque des revendications 1 à 6, dans la valeur de la 10 borne contrainte (FW) du filtre numérique asymétrique est réajusté à une cadence identique à la fréquence de récurrence des étapes /b1/, /b2/ et Ici, de préférence supérieure à 1 kHz. 8. Procédé selon l'une quelconque des revendications 1 à 7, dans lequel le minimum local le plus proche est déterminé en calculant la dérivée de la courbe formant 15 la transformée de Fourier et en identifiant les points où cette dérivée s'annule et en retenant le point le plus proche de la porteuse du canal d'intérêt. 9. Procédé selon la revendication 7 et la revendication 8, dans lequel la valeur de correction (FH) est choisie en fonction de la valeur de la dérivée seconde au point du minimum local. 20 10. Système audio capable de recevoir et de démoduler des ondes électromagnétiques diffusant des programmes radiophoniques sur la modulation de fréquence, caractérisé en ce qu'il est configuré pour mettre en oeuvre le procédé selon l'une quelconque des revendications précédentes.