FR 3022083 A1 20151211 FR 1459472 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 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 numérique sélectif centré ayant 35 une largeur déduite de la fréquence de croisement, à savoir une largeur de 2 x (FO - FC) + FH, FH étant une valeur de correction comprise entre 0 kHz et 10 kHz, /d/ répéter continuellement les étapes /b1/ /b2/ et /c/ 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. Moyennant quoi on peut appliquer un filtrage de largeur adaptée pour éliminer ce canal adjacent perturbateur. 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 10 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 traitements numériques se font avantageusement dans la bande de base, assez facilement. 15 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 : - 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é 20 proposé permet de discriminer deux canaux très proches, de façon plus efficace 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. 25 - 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 d'intérêt. - La largeur du filtre numérique à savoir 2 x (FO - FC) + FH peut prendre des valeurs comprises entre 40 kHz (+/-20 kHz) et 300 kHz (+/-150 kHz), et peut être 30 réajusté à une cadence identique à la fréquence de récurrence des étapes /b/ 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. 35 - 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 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 largeur 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 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 et 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 du canal adjacent potentiellement perturbateur, au cours du temps. 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 de démoduler des ondes électromagnétiques diffusant des programmes radiophoniques sur la bande FM. 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 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é 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 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. 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 numérique), 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 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. Avantageusement, le filtre d'entrée peut être complété par une opération de changement de fréquence (décalage de F0) pour transposer les signaux en bande de base. Les signaux issus du filtrage sont ensuite numérisés par un convertisseur analogique digital 12. La poursuite du traitement du signal est réalisée par des opérations numériques, notamment une mise en mémoire dans une mémoire tampon de 5 données 13. À partir des données mises en mémoire tampon, 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 un filtre numérique 15 de largeur programmable. Un bloc de contrôle de largeur de filtre 16 commande la largeur 10 du filtre W en fonction de l'analyse du contenu de la courbe de densité spectrale de puissance. Tout ou partie des opérations de mise en mémoire tampon, transformée de Fourrier, filtrage numérique peuvent être effectuées par un circuit dédié de type DSP 5. Le terme « DSP » vient de l'anglais « Digital Signal Processor », soit un processeur de signal 15 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 adjacent proche dont la densité spectrale peut déborder à l'intérieur du filtre d'entrée. 20 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. Les signaux en aval du filtrage d'entrée, après numérisation et transposition 25 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 remontée locale de la densité spectrale à droite du point noté « FC+ » référence 92b. 30 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 de largeur de filtre 16 applique un filtre numérique sélectif centré sur FO avec une largeur paramétrée déduite 35 de la fréquence de croisement FC. On peut choisir par exemple une largeur égale à 2 x (FO - FC). De manière générale, on prendra comme valeur de largeur de filtre 16 la valeur W telle que : W = 2 x (FO - 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 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 5 sur la figure 4, la largeur du filtre numérique pourra alors être déduite du contenu spectral du canal d'intérêt, comme 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érencé 91b, se situe à gauche de cette dernière. Dans ce cas, la courbe présente également un maximum local c'est-à-dire une 10 bosse 91a symptomatique de la présence d'un canal adjacent à gauche. De même, dans le cas de courbe 92, on trouve un maximum local 92a à droite du minimum local 92b ; on peut même trouver un second minimum local 92c, mais qui ne joue pas de rôle car il est plus éloigné par rapport à la fréquence centrale 0/F0. La détection d'un minimum local est obtenu en calculant la dérivée de la 15 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 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 20 fonction de la valeur de la dérivée seconde à l'endroit le plus proche de la fréquence. Dans le cas où le canal d'intérêt est encadré par deux canaux potentiellement perturbateurs, l'analyse peut conduire à trouver 2 minimums, un de chaque côté, et on retiendra celui qui est le plus proche de la fréquence centrale (0 en bande de base). Le mécanisme décrit plus haut est sans cesse répété ce qui conduit à une 25 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 30 d'intérêt FO évolue avec le temps, ainsi que le contenu spectral 72 d'un éventuel canal perturbateur F2 représenté en traits pointillés. Les références 71 et 73 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é 35 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 W 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 induit par la présence du canal adjacent F2 dont le contenu spectral s'est rapproché de la fréquence centrale d'intérêt F0. La largeur W du filtre numérique 15 est alors réduite progressivement au fur et 5 à mesure que l'interférence entre les deux canaux augmente. Au début de la troisième période de temps illustrée 83, la largeur du filtre numérique 15 est à nouveau augmentée au fur et à mesure que l'interférence entre les deux canaux diminue. Au début de la quatrième période de temps illustrée 84, en l'absence de 10 détection d'un minimum local dans la courbe spectrale, la largeur du filtre numérique 15 est à nouveau augmentée vers la valeur par défaut. Au cours de la période de temps suivante 85, on retrouve une situation sans détection de minimum local dans le spectre similaire à la première période de temps 81. Au cours de la période de temps suivante 86, une nouvelle séquence de 15 chevauchement provoquée par une excursion négative sur le second canal F2 entraîne une réduction de la largeur W du filtre numérique. Au début de la période de temps illustrée 87, la largeur du filtre numérique 15 est diminuée en raison d'un contenu spectral plus étroit. Au cours de la période de temps suivante 88, le spectre issu de la 20 transformée de Fourier présente un profil plus étroit, et la largeur W du filtre numérique peut être réduite en conséquence. 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 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 numérique sélectif centré ayant une largeur (W) déduite de la fréquence de croisement FC, à savoir une largeur de 2 x (FO - FC) + FH, FH étant une valeur de correction comprise entre 0 kHz et 10 kHz, /d/ répéter continuellement les étapes /b1/ /b2/ et /c/, jusqu'à un changement de canal fréquentiel, 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 20 d'un second canal adjacent, ou d'interférence en dehors du contenu audio. 2. Procédé selon la revendication 1, 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. 25 3. Procédé selon l'une des revendications 1 ou 2, dans lequel la numérisation et la transformée de Fourier sont répétées au moins toutes les 500 ps. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel le filtre d'entrée (11) présente une largeur d'environ 300 kHz. 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel la 30 largeur (W) du filtre numérique 2 x (FO - FC) + FH prend des valeurs comprises entre 40 kHz et 300 kHz, et est réajusté à une cadence identique à la fréquence de récurrence des étapes /b1/,/b2/ et /c/, de préférence supérieure à 1 kHz. 6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel le minimum local le plus proche est 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. 7. Procédé selon la revendication 5 et la revendication 6, 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. 8. Système audio capable de recevoir et de démoduler des ondes électromagnétiques diffusant des programmes radiophoniques sur la modulation de 10 fréquence, caractérisé en ce qu'il est configuré pour mettre en oeuvre le procédé selon l'une des revendications précédentes.