FR 3030946 A1 20160624 FR 1402939 A 20141219 Procédé de correction des défauts d'intermodulation d'un circuit de conversion DOMAINE TECHNIQUE DE L'INVENTION La présente invention concerne un procédé de correction des défauts d'intermodulation introduits par un circuit de conversion et/ou une chaîne de réception se terminant par un tel circuit. L'invention se rapporte également à un produit programme d'ordinateur et un appareil électronique propres à mettre en oeuvre le procédé de correction. ARR1ERE-PLAN TECHNOLOGIQUE La mise en place d'un traitement par une chaîne de réception unique « partagée » des informations en provenance des nombreux capteurs embarqués sur un avion permet d'augmenter les performances des systèmes tout en réduisant l'encombrement, la consommation d'énergie et le coût de l'électronique embarquée. Pour que ceci soit possible, il convient de numériser le signal très en amont de la chaîne de réception. Mais, lorsqu'un circuit de conversion propre à convertir un signal analogique d'entrée en un signal numérique de sortie est situé en amont dans la chaîne de réception, la conversion analogique-numérique est une fonction qui dimensionne la performance de la chaîne de réception. Pour la mise en oeuvre d'une telle fonction, il convient de prendre en compte les contraintes souvent antinomiques des besoins en dynamique et en largeur de bande des différents capteurs de l'avion. Pour cela, il est souhaitable d'obtenir un circuit de conversion travaillant à la fois sur porteuse élevée et avec une fréquence d'échantillonnage élevée pour traiter les besoins en bande large et présentant des caractéristiques dynamiques excellentes pour satisfaire les besoins en bande étroite. Il est à noter que les deux contraintes précédentes sont antinomiques et que, pour une technologie de réalisation donnée, il n'est possible que d'aboutir à un compromis entre vitesse/bande et dynamique. De plus, du fait qu'il est envisagé de fonctionner avec une fréquence d'échantillonnage élevée, il apparaît que la maîtrise du plan de fréquences est un élément important pour limiter la non-linéarité du circuit de conversion. Pour cela, il est constatable que lorsque que la fréquence de fonctionnement fo est égale au rapport de la fréquence d'échantillonnage fEcH par 5 (c'est-à-dire que fo = fEcH/5), les harmoniques d'ordre 2 et 3 du signal se replient autour de deux fréquences qui sont 2*fEcH/5 et 31EcH/5. En outre, l'harmonique d'ordre 5 se replie autour du continu. Similairement, lorsque la fréquence de fonctionnement fo est égale au rapport de la fréquence d'échantillonnage fEcH par 8 (c'est-à-dire que fo = fEcH/8) ou lorsque la fréquence de fonctionnement fo est égale au rapport du produit de 3 avec la fréquence d'échantillonnage fEcFi par 8 (c'est-à-dire que fo = 3*fEcH/8), l'harmonique d'ordre 2 se replie autour de deux fréquences qui sont fEcH/4 et 3*fEcH/4. De plus, il est constaté que les harmoniques d'ordre 3 et 5 du signal se replient autour d'une fréquence qui est soit 3*fEcFF/8 soit fEcH/8. En outre, l'harmonique d'ordre 4 se replie autour de la fréquence qui est fEcH/2. De plus, il est constaté que l'harmonique d'ordre 6 du signal se replie autour d'une fréquence qui est soit 31EcH/4 soit fEcH/4. Placer le signal d'entrée autour des fréquences de fonctionnement précitées permet de ne pas être gêné par les parasites induits par les mécanismes de distorsion les plus difficiles à combattre. Cependant, il reste des phénomènes d'intermodulation créant des parasites dans la bande qui ne sont pas filtrés par principe. Le problème de l'intermodulation se pose lorsqu'un amplificateur amplifie des signaux de forme complexe. On mesure les caractéristiques d'intermodulation d'un amplificateur en appliquant simultanément à l'entrée de celui-ci deux signaux sinusoïdaux, de niveaux identiques, et de fréquences f1 et f2. Si l'amplificateur est parfaitement linéaire, on retrouve sur sa sortie les deux signaux superposés, de fréquences f1 et f2, et seulement ceux-ci. S'il n'est pas parfaitement linéaire, on retrouvera en sortie, d'autres fréquences en plus de f1 et f2. On retrouvera par exemple les fréquences f1+f2, la fréquence f1 - f2, la fréquence 2f1 - f2, et d'une façon générale, on pourra retrouver des signaux présentant une fréquence à n.f1+m.f2 où n et m sont des entiers relatifs. Si n=0, le produit issu de m est une harmonique du signal utile. En particulier, l'intermodulation d'ordre 3 et l'intermodulation d'ordre 5 sont particulièrement gênantes, l'intermodulation d'ordre 3 présentant en général une amplitude plus importante que l'intermodulation d'ordre 5. Compte tenu du compromis vitesse/dynamique à obtenir, il est difficile de disposer des convertisseurs présentant des performances dynamiques compatibles avec les impératifs associés aux modes radar les plus résolus. Il est donc souhaitable de corriger par post-traitement les défauts d'ordre 3 ou 5 induits non seulement par le convertisseur, mais aussi par toute la chaîne de réception qui le précède, en retranchant au signal une correction calculée comme une fonction de lui-même. Usuellement, une telle fonction de correction est issue d'un modèle de distorsion prenant en compte la complexité des défauts de linéarité. Dans le cas de nombreux composants, la fonction de correction choisie est une fonction monôme de degré 3, à savoir que f(x)=. axx3en notant f la fonction de correction et a un paramètre de proportionnalité. Toutefois, un tel type de modèle est inopérant pour un circuit de conversion analogique-numérique du fait que les mécanismes de distorsion mis en jeu sont beaucoup plus complexes que ceux des composants habituels. En effet, le circuit de conversion introduit une erreur à un instant donné dépendant non seulement de la valeur du signal à l'instant donné mais aussi de ses variations et également de sa fréquence. Pour les circuits de conversion, il est connu du document US 8 825 415 B2 une architecture de correction complexe comportant plusieurs étapes successives de calcul de correction (estimation) et de compensation, chaque estimation s'effectuant sur le signal issu du circuit de conversion. De façon schématique, les calculs afférant à l'une quelconque des estimations s'effectuent comme le calcul parallèle de P corrections (P étant un nombre entier positif) distinctes qui sont ensuite sommées entre chaque estimation et retranchées au signal. Chaque calcul des P corrections comporte successivement un filtre (sélection d'une partie du spectre), une interpolation, une distorsion de type e (n variant d'une correction à l'autre), un filtrage et une décimation. L'objet de la correction a donc pour objet l'intégralité du spectre et porte à la fois, du fait du type de distorsion mise en jeu, de la distorsion harmonique et de l'intermodulation. Cependant, une telle correction n'est pas adaptée pour des signaux en bande étroite sur-échantillonnés. Notamment, une telle correction n'est pas adaptée à un signal dans lequel la distorsion harmonique a été filtrée (séparée de l'intermodulation du fait du plan de fréquence et filtrée par un filtre de décimation). Dans certains cas, l'application d'un tel mécanisme de compensation à un tel signal conduit à générer de la distorsion harmonique non désirée. De plus, mettre en oeuvre une telle compensation implique une charge de calcul très importante du fait des multiples filtrages et interpolations utilisés. Il existe donc un besoin pour un procédé de correction des défauts d'intermodulation introduits par un circuit de conversion qui permette, avec une mise en oeuvre aisée, de diminuer les défauts d'intermodulation sans générer de la distorsion harmonique non désirée. RESUME DE L'INVENTION Pour cela, il est proposé un procédé de correction des défauts d'intermodulation introduits par un circuit de conversion adapté à convertir un signal analogique d'entrée en un signal numérique de sortie. Le procédé comporte les étapes de fourniture du signal numérique de sortie, de génération d'une fonction de correction, la fonction de correction étant une fonction complexe du signal numérique de sortie, et de correction du signal numérique de sortie à l'aide de la fonction de correction, pour obtenir un signal numérique de sortie corrigé. Suivant des modes de réalisation particuliers, le procédé comprend une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : - le signal numérique de sortie corrigé est égal à la différence entre le signal numérique de sortie et le produit du signal numérique de sortie par la fonction de correction. - l'étape de génération comporte les étapes de calcul d'une première fonction polynômiale du carrée du module du signal numérique de sortie, pour obtenir un premier terme final, de calcul d'une deuxième fonction polynômiale du carrée du module du signal numérique de sortie selon une représentation complexe, pour obtenir un deuxième terme intermédiaire, d'application d'un filtre sur le signal numérique de sortie pour obtenir un signal filtré, de calcul du produit du deuxième terme intermédiaire par le produit du complexe conjugué du signal numérique de sortie par le signal filtré, pour obtenir un deuxième terme final, et de calcul de la somme du premier terme final et du deuxième terme final, pour obtenir la fonction de correction. - le filtre appliqué à l'étape d'application d'un filtre restitue la dérivée du signal numérique de sortie. - l'étape de génération comporte l'application d'au moins une opération au signal numérique de sortie, pour obtenir un signal de calcul, l'étape de correction comportant l'application de la fonction de correction au signal de calcul. - au moins une opération est choisie dans le groupe constitué d'un filtrage de Hilbert du signal numérique de sortie, d'une démodulation numérique du signal numérique de sortie, d'un filtre à phase linéaire et d'une décimation. - l'étape de génération comporte les étapes d'application d'un signal de calibration au circuit de conversion, par exemple un signal à deux fréquences distinctes, de mesure du signal obtenue en sortie du circuit de conversion, pour obtenir un signal mesuré et d'analyse des composantes spectrales du signal mesuré. Il est également proposé un produit programme d'ordinateur comportant un support lisible d'informations, sur lequel est mémorisé un programme d'ordinateur comprenant des instructions de programme, le programme d'ordinateur étant chargeable sur une unité de traitement de données et adapté pour entraîner la mise en oeuvre d'un procédé tel que précédemment décrit lorsque le programme d'ordinateur est mis en oeuvre sur l'unité de traitement des données. Il est aussi proposée un appareil électronique comportant un récepteur d'un signal analogique d'entrée, un circuit de conversion adapté à convertir le signal analogique d'entrée en un signal numérique de sortie, et un correcteur adapté à générer d'une fonction de correction, la fonction de correction étant une fonction complexe du signal numérique de sortie, et à corriger le signal numérique de sortie à l'aide de la fonction de correction, pour obtenir un signal numérique de sortie corrigé. Suivant un mode de réalisation particulier, le correcteur est un circuit logique programmable. BREVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui suit de modes de réalisation de l'invention, donnée à titre d'exemple uniquement et en référence aux dessins qui sont : figure 1, une représentation schématique d'appareil électronique, figure 2, une vue schématique d'un exemple de système permettant la mise en oeuvre d'un procédé de correction, et figures 3 à 11, des graphes illustrant des simulations de résultats obtenus après mise en oeuvre d'exemples de procédés de correction. DESCRIPTION DETAILLEE DE MODES DE REALISATION DE L'INVENTION Description d'un exemple d'appareil électronique Un appareil électronique 10 est illustré schématiquement à la figure 1. L'appareil électronique 10 comporte un récepteur 12, un circuit de conversion 14, un correcteur 16 et un émetteur 18. Le récepteur 12 est propre à recevoir un signal analogique d'entrée. Le circuit de conversion 14 est adapté à convertir le signal analogique d'entrée en un signal numérique de sortie. Le circuit de conversion 14 comporte une entrée 14E et une sortie 14S. L'entrée 14E du circuit de conversion 14 est reliée au récepteur 12 et reçoit le signal analogique d'entrée. Le signal numérique de sortie est délivré sur la sortie 14S du circuit de conversion 14. Un tel circuit de conversion 14 comporte un convertisseur analogique-numérique et désigné sous l'acronyme CAN. En variante, le circuit de conversion 14 comprend aussi une chaîne de traitement analogique, par exemple une chaîne dite de réception. Parmi les composants de la chaîne de traitement analogique, se trouvent, par exemple, un ou plusieurs amplificateurs. Le correcteur 16 comporte une entrée 16E et une sortie 16S. L'entrée 16E du correcteur 16 est reliée à la sortie 14S du circuit de conversion 14. La sortie 16S du correcteur 16 est reliée à l'émetteur 18. Le correcteur 16 est adapté à générer une fonction de correction. La fonction de correction est une fonction complexe du signal numérique de sortie. Un signal numérique corrigé est obtenu à partir du signal numérique de sortie et de la fonction de correction. A titre d'exemple, le correcteur 16 est un circuit logique programmable. Un circuit logique programmable, ou réseau logique programmable, est un circuit intégré logique qui peut être reprogrammé après sa fabrication. Notons qu'il serait impropre de parler de programmation au sens logiciel (contrairement à un microprocesseur, il n'exécute aucune ligne de code). Ici, mieux vaudrait parler de « reconfiguration » plutôt que de reprogrammation (on modifie des connexions ou le comportement du composant, on connecte des portes logiques entre elles, etc.). Le verbe programmer est toutefois le plus fréquent, mais au sens de personnaliser. Il s'agit bel et bien de réseaux logiques reconfigurables et modifiables. De préférence, le correcteur 16 est un FPGA, acronyme anglais de l'expression « field-programmable gate array », réseau de portes programmables in situ. Description d'un premier exemple de mise en oeuvre d'un procédé de correction des défauts d'intermodulation introduits par le circuit de conversion 14 Le fonctionnement de l'appareil électronique 10 est maintenant décrit en référence à un premier exemple de mise en oeuvre d'un procédé de correction des défauts d'intermodulation introduits par le circuit de conversion 14. Le procédé de correction comporte une étape a) de fourniture du signal numérique de sortie. L'étape a) de fourniture du signal numérique de sortie est notamment réalisable par fourniture d'un signal analogique d'entrée sur le récepteur 12 et conversion du signal reçu par le circuit de conversion 14. Le procédé de correction comprend également une étape b) de génération d'une fonction de correction. L'étape b) de génération est mise en oeuvre par le correcteur 16. De manière générale, la fonction de correction est une fonction complexe du signal numérique de sortie. Selon le premier exemple de mise en oeuvre du procédé de correction, l'étape b) de génération comporte une étape bli) de calcul d'une première fonction polynômiale du carrée du module du signal numérique de sortie, pour obtenir un premier terme final. L'étape b) de génération comprend aussi une étape b21) de calcul d'une deuxième fonction polynômiale du carrée du module du signal numérique de sortie selon une représentation complexe, pour obtenir un deuxième terme intermédiaire. L'étape b) de génération comporte également une étape b31) de calcul d'une deuxième fonction polynômiale du carrée du module du signal numérique de sortie selon une représentation complexe, pour obtenir un deuxième terme intermédiaire. L'étape b) de génération comprend une étape b41) d'application d'un filtre sur le signal numérique de sortie pour obtenir un signal filtré. Par exemple, le filtre appliqué à l'étape b41) d'application restitue la dérivée du signal numérique de sortie c'est-à-dire restitue préférentiellement l'échantillonnage de la dérivée du signal continu avant échantillonnage. Dans le cas général, par l'expression « restitue l'échantillonnage de la dérivée», il est entendu que le filtre, supposé être réalisable, n'en fournit en fait qu'une approximation et que le filtre ne constitue donc qu'une approximation plus ou moins précise de la fonction idéale interpolation-dérivation-rééchantillonnage du signal numérique de sortie aux instants d'échantillonnage. Dans le cas spécial d'un signal numérique de sortie répétitif, de période multiple de la période d'échantillonnage, le filtre peut alors constituer exactement cette fonction. Selon le premier exemple, l'étape b) de génération comporte aussi une étape b51) de calcul du produit du deuxième terme intermédiaire par le produit du complexe conjugué du signal numérique de sortie par le signal filtré, pour obtenir un deuxième terme final. L'étape b) de génération comporte également une étape b61) de calcul de la somme du premier terme final et du deuxième terme final, pour obtenir la fonction de correction. Le procédé de correction comporte aussi une étape c) de correction du signal numérique de sortie pour obtenir un signal numérique de sortie corrigé. Le signal numérique de sortie corrigé est obtenu à l'aide de la fonction de correction. Selon un mode de réalisation particulier, le signal numérique de sortie corrigé est égale à la différence entre le signal numérique de sortie et le produit du signal numérique de sortie par la fonction de correction. Un tel procédé permet donc d'obtenir un signal numérique de sortie dans lequel les défauts d'intermodulation introduits par le circuit de conversion 14 sont réduits. Cela apparaît à la lecture de la section « démonstration mathématique » qui suit. Autrement formulé, le procédé de correction des défauts d'intermodulation introduits par le circuit de conversion 14 qui permet, avec une mise en oeuvre aisée, de diminuer les défauts d'intermodulation sans générer de la distorsion harmonique non désirée. Un tel procédé permet d'obtenir un circuit de conversion 14 travaillant à la fois sur porteuse élevée et avec une fréquence d'échantillonnage élevée pour traiter les besoins des bandes les plus larges et présentant des caractéristiques dynamiques excellentes pour satisfaire les besoins en bande étroite à haute résolution. Cela ouvre la voie à la réalisation pour l'avionique de chaîne de réception unique « partagée » des informations en provenance des nombreux capteurs embarqués sur un avion. En effet, comme expliqué précédemment, lorsque le circuit de conversion 14 est situé en amont de la chaîne de réception, la fonction de conversion réalisée par le circuit de conversion 14 dimensionne la performance de la chaîne. H convient donc que le circuit de conversion 14 réponde aux besoins différents des radars et des organes de communication de l'avion en termes de dynamique et de largeur de bande. La réalisation de chaîne de réception unique « partagée » des informations permet d'envisager une augmentation des performances du système tout en réduisant l'encombrement, la consommation d'énergie et le coût de l'électronique embarquée. L'invention n'est pas restreinte au seul cas de l'électronique embarquée ou destinée aux seules applications avioniques. L'invention est valable pour toute chaîne de réception et permet d'éliminer les défauts d'intermodulation d'ordres impairs. Il s'agit typiquement des défauts conduisant, pour un ordre d'intermodulation 2N+1 donné, N étant un entier, à la génération de composantes aux fréquences I m f1 - (m+1) f2 I et I (m+1) f1 - m f2 I et avec m positif et inférieur ou égal à l'entier N. L'invention n'a d'intérêt réel que si, par ailleurs, les défauts harmoniques ne se seront pas repliés sur le signal numérique de sortie.35 Description d'un deuxième exemple de mise en oeuvre d'un procédé de correction des défauts d'intermodulation introduits par le circuit de conversion 14 Selon un deuxième exemple de mise en oeuvre du procédé de contrôle, le procédé de correction comporte également une étape a) de fourniture du signal numérique de sortie. L'étape a) de fourniture du signal numérique de sortie est identique au premier exemple de mise en oeuvre. Les remarques précédentes s'appliquent donc également à l'étape a) selon le deuxième exemple. Le procédé de correction comprend également une étape b) de génération de la fonction de correction. L'étape b) de génération comporte une étape b12) de filtrage de Hilbert du signal numérique de sortie, pour obtenir un premier signal. En variante, à l'étape b12), il est effectué une démodulation numérique du signal numérique de sortie pour obtenir le premier signal. L'étape b) de génération comporte aussi une étape b32) d'application d'un premier filtre à phase linéaire sur le premier signal pour obtenir un premier signal filtré. L'étape b) de génération comporte aussi une étape b42) de décimation du premier signal filtré pour obtenir un deuxième signal. En traitement du signal, la décimation est synonyme de sous-échantillonnage : on ne garde qu'un certain nombre d'échantillons par rapport au signal original (un sur 2, un sur 10...). L'étape b) de génération comporte aussi une étape b52) d'application d'un deuxième filtre à phase linéaire sur le deuxième signal pour obtenir un deuxième signal filtré. L'étape b) de génération comporte aussi une étape b62) de décimation du deuxième signal filtré pour obtenir un troisième signal. L'étape b) de génération comporte aussi une étape b72) d'application d'un troisième filtre à phase linéaire sur le troisième signal pour obtenir un troisième signal filtré. L'étape b) de génération comporte aussi une étape b82) de décimation du troisième signal filtré pour obtenir un signal de calcul. En variante, un seul filtre à phase linéaire est appliqué à l'étape b) de génération. Selon une variante, un nombre de filtre à phase linéaire strictement supérieure à 3 est mis en oeuvre. Selon une autre variante, aucune étape de décimation n'est utilisée lors de la mise en oeuvre de l'étape b) de génération. Selon un mode de réalisation particulier, seules certaines étapes d'application d'un filtre à phase linéaire sont suivies d'une étape de décimation. Le procédé de correction comporte aussi une étape c) de correction du signal numérique de sortie. Le signal numérique de sortie corrigé est obtenu par application de la fonction de correction au signal de calcul. Un tel procédé permet donc d'obtenir un signal numérique de sortie dans lequel les défauts d'intermodulation introduits par le circuit de conversion 14 sont réduits. Cela apparaît à la lecture de la section « démonstration mathématique » qui suit. Autrement formulé, le procédé de correction des défauts d'intermodulation introduits par le circuit de conversion 14 qui permet, avec une mise en oeuvre aisée, de diminuer les défauts d'intermodulation sans générer de la distorsion harmonique non désirée. Description d'un troisième exemple de mise en oeuvre d'un procédé de correction des défauts d'intermodulation introduits par le circuit de conversion 14 Selon un troisième exemple de mise en oeuvre du procédé de contrôle, le procédé de correction comporte également une étape a) de fourniture du signal numérique de sortie. L'étape a) de fourniture du signal numérique de sortie est identique au premier exemple de mise en oeuvre. Les remarques précédentes s'appliquent donc également à l'étape a) selon le troisième exemple. Le procédé de correction comprend également une étape b) de génération de la fonction de correction. L'étape b) de génération comporte une étape b13) d'application d'un signal de calibration au circuit de conversion 14. Le signal de calibration est, par exemple, un signal à deux fréquences distinctes. L'étape b) de génération comporte une étape b23) de mesure du signal obtenue en sortie du circuit de conversion 14, pour obtenir un signal mesuré. L'étape b) de génération comporte une étape b33) d'analyse des composantes spectrales du signal mesuré. Une telle étape b33) d'analyse est généralement mise en oeuvre par utilisation d'une transformation de Fourier discrète (TFD) ou d'une transformée de Fourier rapide (FFT). La transformation de Fourier discrète est un outil mathématique de traitement du signal numérique, qui est l'équivalent discret de la transformation de Fourier continue qui est utilisée pour le traitement du signal analogique. La transformation de Fourier rapide (acronyme anglais FFT pour Fast Fourier Transform) est un algorithme de calcul de la transformation de Fourier discrète (TFD). L'étape b) de génération comporte une étape b34) d'utilisation de la valeur des composantes spectrales du signal mesuré pour en déduire la fonction de correction. Le procédé de correction comporte aussi une étape c) d'obtention du signal numérique de sortie corrigé. L'étape c) d'obtention du signal numérique de sortie corrigé est identique au premier exemple de mise en oeuvre. Les remarques précédentes s'appliquent donc également à l'étape c) selon le troisième exemple. Un tel procédé permet donc d'obtenir un signal numérique de sortie dans lequel les défauts d'intermodulation introduits par le circuit de conversion 14 sont réduits. Cela apparaît à la lecture de la section « démonstration mathématique » qui suit. Autrement formulé, le procédé de correction des défauts d'intermodulation introduits par le circuit de conversion 14 qui permet, avec une mise en oeuvre aisée, de diminuer les défauts d'intermodulation sans générer de la distorsion harmonique non désirée. Autres variantes D'autres variantes correspondant à des combinaisons des modes de réalisations précédemment décrits sont également envisageables. Notamment, il est possible de mettre en oeuvre le procédé avec un système 110 et un produit programme d'ordinateur 112. Le système 110 et le produit programme d'ordinateur 112 sont représentés à la figure 2. L'interaction du produit programme d'ordinateur 112 avec le système 110 permet de mettre en oeuvre un procédé d'identification d'une relation entre des éléments physiques. Le système 110 est un ordinateur. Plus généralement, le système 110 est un calculateur électronique propre à manipuler et/ou transformer des données représentées comme des quantités électroniques ou physiques dans des registres du système 110 et/ou des mémoires en d'autres données similaires correspondant à des données physiques dans des mémoires, des registres ou d'autres types de dispositifs d'affichage, de transmission ou de mémorisation. Le système 110 comporte un processeur 114 comprenant une unité de traitement de données 116, des mémoires 118 et un lecteur 120 de support d'informations. Le système 110 comprend également un clavier 122 et une unité d'affichage 124. Le produit programme d'ordinateur 112 comporte un support lisible d'informations 120. Un support lisible d'informations 120 est un support lisible par le système 110, usuellement par l'unité de traitement de données 114. Le support lisible d'informations 120 est un médium adapté à mémoriser des instructions électroniques et capables d'être couplé à un bus d'un système informatique. A titre d'exemple, le support lisible d'informations 120 est une disquette ou disque souple (de la dénomination anglaise de « floppy disk »), un disque optique, un CD-ROM, un disque magnéto-optique, une mémoire ROM, une mémoire RAM, une mémoire EPROM, une mémoire EEPROM, une carte magnétique ou une carte optique. Sur le support lisible d'informations 20 est mémorisé un programme d'ordinateur comprenant des instructions de programme. Le programme d'ordinateur est chargeable sur l'unité de traitement de données 114 et est adapté pour entraîner la mise en oeuvre d'un procédé d'identification d'une relation entre des éléments physiques lorsque le programme d'ordinateur est mis en oeuvre sur l'unité de traitement des données 114. Démonstration mathématique Comme précisé précédemment, pour traduire les défauts induits par un circuit de conversion 14, il convient de prendre une fonction plus compliquée qu'une simple fonction polynomiale du signal. En effet, la fonction représentation de la distorsion traduit le fait que l'erreur à un instant donné dépend non seulement de la valeur instantanée du signal mais également de ses variations. En notant s = Eel: le signal complexe considéré (c'est-à-dire le signal numérique de sortie) où s est le signal, p est le module du signal, j le nombre complexe et cp le déphasage, alors la dérivée logarithmique du signal s'écrit : 1 as 1 ap . -x-.-x-+J2n-xf(t) s at p at où f, est la fréquence instantanée égale à: at En négligeant les variations logarithmiques de l'amplitude du signal par rapport aux variations induites par la fréquence instantanée, ce qui est tout à fait justifié dans le cas d'un fonctionnement à bande étroite sur porteuse, alors, si en notant x et y les parties réelle et imaginaire du signal, les relations approchées suivantes sont obtenues : ax -27rx (t)x y et -ay 2yrx (t)x x at at A l'étude des relations précédentes, il apparaît que les variations de la partie réelle x du signal sont proportionnelles à la partie imaginaire y du signal. Il est donc considéré une traduction mathématique de la distorsion de la partie réelle du signal sous la forme suivante beaucoup plus appropriée : f(x,y).= cr0xx3 + oeixx2y + a2x y2x + ce3x Où: - ao, ai, a2 et a3 sont des coefficients. Comme expliqué précédemment, il est proposé d'exploiter la forme précédente particulier de modèles de défaut. D'un point de vue physique, il est possible de montrer que le modèle de défauts donné par la relation ci-dessus se partage en deux parties distinctes et indépendantes : l'intermodulation proche du signal et la distorsion harmonique et l'intermodulation lointaine. Ainsi qu'expliqué en introduction, il est spécifiquement souhaité se concentrer sur la correction des défauts d'intermodulation puisque les autres défauts sont éliminés grâce au plan de fréquence. Pour cela, il est proposé une reformulation de la relation précédente en fonction d'un seul paramètre noté k(s) qu'il est possible d'acquérir via une phase de calibrage. Pour un signal s n'étant pas le signal réel issu du circuit de conversion 14 mais un signal complexe obtenu par filtrage de Hilbert, il est considéré le signal suivant : s -> [1- k(s)]x s Où: - k(s) est le paramètre et est donné par k(s)= Is12 avec a un nombre. Une telle relation permet de corriger les défauts d'intermodulation d'ordre 3. La fonction k(s)*s correspond dans ce cas à la fonction de correction. Une telle relation peut s'étendre pour les défauts d'intermodulation d'ordre 5 en modifiant la valeur du paramètre k(s). Dans un tel cas, le paramètre k(s) est donné par la relation suivante : kiCO= (ait+ 04.1313).1412 Par généralisation, il s'en suit que pour la correction des défauts d'intermodulation jusqu'à un ordre 2*N+3, N étant un entier positif, le paramètre k(s) est donné par la relation suivante : k(S)= /412 )X IS12 Où: - P est un polynôme d'ordre N Autour d'une fréquence, les coefficients du polynôme peuvent être obtenus à partir d'une phase de calibrage par mesure du niveau des deux tons d'un signal 2-tons et du niveau des parasites d'intermodulation jusqu'à l'ordre 2*N+3. Pour obtenir une correction dans une bande plus large, il convient de prendre en compte le fait que les coefficients du polynôme P dépendent de la fréquence. Il est possible d'utiliser une dépendance linéaire des coefficients en fonction de la fréquence instantanée. Comme la fréquence instantanée est une fonction du signal complexe, il est abouti à une complexification de la fonction de correction qui comporte alors deux termes faisant chacun intervenir un polynôme : k(,)x s. [4,12 )+. Q,s12)x flx 2 /1s12)x s2 + Xs X S = Où: - Q est un polynôme, - S désigne la fonction mathématique partie imaginaire, et - 1-10 désigne la dérivée du filtre interpolateur et a pour réponse aux points d'échantillonnage pour n différent de 0 hD(n)= qui peut être approchée sur une durée finie égale à 2*N par i; D (n)= (-1)5r Nx tan - Pour expliciter le procédé de correction, il est possible de décomposer la fonction f(x,y) représentant la distorsion du signal comme la somme de deux termes. Le premier terme s'exprime mathématiquement comme : a0xx3 + ce,xy3 .9i(a0+ jxajx(x3 +(jx.di =9i(ao+ jxajx(x+ jxy)x x2 -xx(jxy)+Uxyd Où: - 9i désigne la fonction mathématique partie réelle. Le deuxième terme s'exprime mathématiquement comme : Ce, X x2y + Ct2 X y2X = (0(2 + jXCOX ((jXy)2x + x2 X jX)))1 =9i[-(a2+ jxai)x(x+ jxy)xxx(jxy)] En utilisant les deux expressions précédentes du premier terme et du deuxième terme, il en résulte que la fonction f(x,y) représentant la distorsion du signal vaut : f(x,y) =9Î (x+ jxy)xl(ao + ixotjx(x2-xx(jxy)+(jxy)2)-(a2+ jxc)xxx(jxy)}j (ao + jxa3)x9ikx+GXy))21 (x+ jX.y)X (CY2+/X00+(ao-FiXCe3)x jx3kx±(ixy))21 2 L'expression précédente se simplifie en posant z = x+ jxy l'extension analytique de x. Il est alors obtenu les expression mathématiques suivantes : f(x,y) = =9"i[zx = zX zX ixce3)xz2 + z*2 (a2 + jxa,)+(a0+ jxa3)xz2 *211 2 2 2 (ao+ jxa3)-(a2+ jxal) xz2 +3(ao+ jxa3)+(a2+ jxal)xz*2 4 4 (ceo-Fixa3)-(a2-1-ixoei)Xz2+(3(ao-Fixce3)±(ce2±ixcer X z 2} 4 4 Pour simplifier encore plus l'expression précédente, il est posé : 7(4= zx (ao+ jxa3)-(a2+ jxai) xz 2 + (3(a0+ jxa3)+(a2+ jixa)\* x z12 4 4 ) L'expression de la fonction f(x,y) représentant la distorsion du signal s'écrit alors : 7(z)+ 7* (4 f(x,y)= 2 Or, la correction réelle est la suivante : z+ z* z+z* 7(z)+:7*(z) x= 2 2 2 De ce fait, l'expression de la correction du signal complexe s'écrit : Z Z 7(Z) Z Z x [(a0+ jxo(3)-(a2+ jxcej2 lx - +[3(a0 + jxa3)+(a2+ jxai)] x 2 2 2 2 2 Il est alors possible d'effectuer la correction suivante en complexe sur le signal s = z/2 issu du filtrage de Hilbert: s --->[1-k(s)]xs où k(s)=[(ao+ jxa,)-(a2.+ jxai)]xs2 +[3(ao + j x a3)+ (a2 + x ce,)]* xisr Dans le spectre, les deux termes formant la distorsion s x k(s) ont une signification différente : le premier terme représente le cumul de la distorsion harmonique et de l'intermodulation lointaine alors que la deuxième terme représente l'intermodulation proche. Mathématiquement, cela s'écrit ainsi : 2 21 sxk(s)=-sxRce0+ jxcx3)-(a2+ jxcei)]xs2+sx[3(a0+ jx(x3)+(a2+jxa3exis2 distorsion harmonique intermodulation intermodulation lointaine proche porteuse L'opération ainsi décrite s'effectue après passage en complexe. La séparation des termes « lointains » et « proches » permet de réaliser la seule correction des défauts d'intermodulation « proche ». Une telle correction peut s'effectuer après retour du spectre en bande de base puisque la phase introduite transposition n'affecte pas le terme Ie. Une telle opération peut aussi s'effectuer après filtrage de la distorsion harmonique et décimation sans craindre de générer des produits harmoniques non désirés qui dont le repliement est à éviter. Cela n'est bien sûr possible qu'à la condition que la bande du signal soit suffisante faible pour que la distorsion harmonique puisse être filtrée par le plan de fréquence ; pour que cela soit possible après filtrage, il convient également nécessaire que le plan de fréquence ne déforme pas la phase. Dans le cas où le procédé comporte les étapes du troisième exemple, si deux tons à f1 et f2 sont simultanément présents dans la bande, la fonction f(x,y) représentant la distorsion du signal génère à la fois des parasites harmoniques (3 f1 et 3 f2) et des parasites d'intermodulation hors bande (situés à 2 f1 + f2 et 2 f2 + f1) et dans la bande (situés à 2 fl - f2 et 2 f2 - fi). Il est alors possible de poser : x = a, cos çoi + a, cos e2 y= al sin ço, + a2 sin 92 Il est alors possible d'obtenir en calculant les différents termes impliqués dans la fonction f(x,y) représentant la distorsion du signal : 3 3 3 2 X = a, cos3 'o,+ a2cos3 ço' + 3a,2a2cos2 el cos e2 + 3a2 ai cos2 92 C 0 S 91 3 3 X2y = al cos2 91 sin çoi + a2 cos2 92 sin ÇO2 +a,2a2x(cos2 ço, sin 92 + 2 sin çoi cos 91 cos 92 ) + a22a1(cos2 92 sin ço, + 2 sin yo2 cos 92 cos 91) 3 3 - 2 y2x = a, sin2 ço, cos 91 + a2 sm 92 cos 92 3 +a12a2 x (sin2 çoi cos e2 + 2 cos 91 sin 91 sin 92 ) + aa,(sin2 92 cos çoi + 2 cos 92 sM 92 sin 91) 3 - 3 2 2 y3 = a, sin3 91 + a2 sm 92 + 3ai a2 sin2 91 sin 92 + 3a2a1 sin2 92 sin 91 Il est alors obtenu une expression pour l'intermodulation et la distorsion harmonique selon laquelle : 33a12a2 / \ 3a2a x = cos29i 92)+ 2 I COS(402 -9 ) 4 4 a3 a3 3a2a 3a2a + cos 391 + ---a- cos 3ÇO2 H- I 2 cos(29, +92)+ 2 I COS(292 + q)i ) + ... 4 4 4 4 2 \ a2a X2y =ai a2 sin(291 92)+ sin(292 - 9i) 2 1 4 4 a3 +- sin 391 + -2-sin 392 + 3a2 sin(291+92)+ 3a1 sin(292 +91)+... a3 4 4 4 4 ,2, 2 y2x=-1-2 cos(29, e2)+ a2c11 cos(292 -91) 4 4 a3 a3 3a12a2 -- cos 3ço1 ---. cos 392 4 4 4 4 \ y3 = 3a' 2a2 sin(291 92)+ 3a2a 2 1 sin(2ço2 -93 4 4 2 2 a3 a 3 sin 3 - sin 39 3ala2 sin(2çoi ço2 ) 3a2 a sin(292 + e 2 4 4 4 4 Les signaux parasites réels induits par le mécanisme de distorsion cubique considéré sur chacune des six fréquences 2f1 - f2, 2f2 - f1, 3f1, 2f1 + f2, 2f2 + f1, et 3f2 sont alors respectivement les suivants : 2 . 3a2ai cos(2e1+ 92) cos(292 +91)+... 2 a a 1 2 X Pa° + a2)xcos(291 - 4 2 a a x[(3oe0 + a2 )x cos(2(o2 - 2 4 3 X [(a,,- a2)x cos 391 + 4 3a2a 1 2 x Ra° a2)xcos(291 + 4 3a2a 2 1 x K 0 _ a a2)xcos(292+ 4 92)+(a1+3a3)xsin(291 -92)] 9i)+ (''i +36Y3)xsin(492 -91)1 - a3)x sin 391] 92)+ a3)xsin(2ço1 +92)1 91)+ a3)xsin(292 +91)] 3 Ka° - a2)x cos 392 + (a,- a3)x sin 3ço2] 4 Lorsque les deux fréquences f1 et f2 sont des fréquences cohérentes sur un pas de transformée de Fourier rapide (aussi désignée par l'acronyme FFT pour « Fast Fourier Transform ») correspondant à un nombre de NFFT points (NFFT étant un entier positif), il devient possible de déterminer l'ensemble des grandeurs à partir d'une normalisation de la transformée de Fourier rapide (par NFFT) du signal qui permet d'en séparer les différentes contributions spectrales. Pour la suite, il est noté : - 5+1 l'amplitude de la composante obtenue à une première fréquence +fi, - s..1 l'amplitude de la composante obtenue à l'opposé de la première fréquence -f1, - s+2 l'amplitude de la composante obtenue à une deuxième fréquence +f2, - s_2 l'amplitude de la composante obtenue à l'opposé de la deuxième fréquence - f2, - i+12 l'amplitude de la composante obtenue à une troisième fréquence +211 - f2, - i-12 l'amplitude de la composante obtenue à l'opposé de la troisième fréquence -211 + f2, - i+21 l'amplitude de la composante obtenue à une quatrième fréquence +212 - fi, - I-21 l'amplitude de la composante obtenue à l'opposé de la quatrième fréquence -212 + - h+1 l'amplitude de la composante obtenue à une cinquième fréquence +311, - h1 l'amplitude de la composante obtenue à l'opposé de la cinquième fréquence - 311, - h+12 l'amplitude de la composante obtenue à une sixième fréquence +211 + f2, - h_12 l'amplitude de la composante obtenue à l'opposé de la sixième fréquence -211 - f2, - h+21 l'amplitude de la composante obtenue à une septième fréquence +2f2 + fi, - h_21 l'amplitude de la composante obtenue à l'opposé de la septième fréquence -2f2 - fi, - h+2 l'amplitude de la composante obtenue à une huitième fréquence +312, et - h..2 l'amplitude de la composante obtenue à l'opposé de la huitième fréquence - 312. Pour chacune des notations introduites, il est à noter que, comme les signaux sont réels, chacune des amplitudes notées avec un indice négatif, c'est-à-dire les grandeurs négatives que sont Si, -12hh -, i_1 h , - -1, -12, -21 et h-2, sont les conjuguées des amplitudes 10 correspondantes notées avec un indice positif, c'est-à-dire les grandeurs positives que 1, -+12, +21, +1, +12, +21 t h+2. sonts+i i hh h e Or, il est connu que sin (0 = 9î(- j X ei9) , il est alors obtenu pour l'amplitude de la composante obtenue à la première fréquence +fi et l'amplitude de la composante obtenue à la deuxième fréquence +f2 les relations suivantes : { ,s+1 = aieJ9' 12 s+2= a2e°2/2 Similairement, la relation sin ça = 9î(- jxeJ9) permet d'obtenir pour l'amplitude de la composante obtenue à la troisième fréquence +2*f1 - f2 et l'amplitude de la composante obtenue à une quatrième fréquence +2*f2- fi les relations suivantes : { 412 - 8 , _ cia, x[(3a0+ a2)- jx(a1+ 3a3 1 x ei(2e'' -92) = [(3a0 + a2)- j x (a1+3a3)]xss 8 , . a2a, x[(3a0+ a2)- jx(a,+3a,)]xe2-9') =[(3a0+ a2)- jx(a,+ 3a3)1x ,s,si* où la notation « * » en exposant renvoie à l'opération mathématique de conjugaison. Similairement, la relation sinço=9i(- jxel9) permet d'obtenir des relations pour l'amplitude de la composante obtenue à la cinquième fréquence +3*f1, l'amplitude de la composante obtenue à la sixième fréquence +2*fi + f2, l'amplitude de la composante 15 obtenue à la septième fréquence +2f2 + fi et l'amplitude de la composante obtenue à la huitième fréquence +3*f2. Le système suivant donne le détail de chacune des relations mathématiques précédentes : h+1 .-ci3 x[(a0-a2)- jx(al-a,)]x e.'39' =[(oe0- a,)- jx(a,- a, lx si' 8 h+12= h+21 = 3a2a 3a1a2 8 2 x Rao - ce2 )- j(al -a3)1xef(29'+92) =[(a,- a2)- j(a,- a3)1x3s12s2 82 1 x pro _ CY2 )- j(ai - a, lx ef(292 +91) =[(a0- a2)- j(a,-a3)]x3ss, h+2 = -La3 x [(a,- a2)- jx(a,-a,)]x e'3"2 =[(a0-a2)- jx(a,- a,)]xs; 8 Par identification des différents termes, il est obtenu le système d'équations 20 suivant : [3(a0+ jxa3)+(a2+ jxoti)T = /24.1 3 i+12 h+12 2 * Si S2 (a0+ jxot3)-(a2+ jxaj= 3:12'92 i+21 "+21 2 * ,S2S1 h+2 3 S2 A l'analyse du système d'équations obtenu, il apparaît qu'il s'agit d'un système d'équations avec une forte redondance. Dès lors, il est possible de considérer un système d'équations dit « moyen » en utilisant des opérations de moyennages, notamment par emploi de moyenne arithmétique. Selon un premier exemple, il est obtenu le système d'équations dit « moyen » suivant : 1 ( - - \1/2 +12 +21 [3(ao+ jxcr3)+(a2+ jxoel)]*= sis2 s,s2 1 'h +1h+12h+21h+2 (a0+ jxa3)-(a2+ jx(xj= 9s2s2 1 2 Selon un deuxième exemple, il est obtenu le système d'équations dit « moyen » suivant : [3(a0+ jxoe3)+(a2+ jxot,)]* = 1 ; 1 i+12 + 421 ( 2 :s12s; s22s1*1 (a0+ixce3)_(a2+ixtyi)= 1 + h, + h+12 + h+21 ± h+2 4 \. s13 3s1 s23S2S S3 2 1 2 ' Le procédé de correction précédemment décrit s'applique également pour corriger des intermodulations d'ordre différent. Une telle extension s'obtient en utilisant des polynômes d'ordres supérieurs en Ii15 le, c'est-à-dire des termes de correction de la forme k(s)x s =sr x1,si2 )x s . En supposant que s = s, + s2 où s1 et s2 sont deux signaux aux fréquences f1 et f2, il est alors obtenu l'expression mathématique suivante: 12 12 12 12 ( IS1 = + S21 = .5'11 +15'21 )+ \Si*S2 + S;Si Il en résulte que : 20 Irn = [(s112 SI +1,5,211+ (si*s2 ',*23,, 1 )1' = zc811 nk 12 1-1_ 212 )n-k -1- 3131 k=0 où les coefficients Cnk sont les coefficients du binôme. Or, en écrivant le premier signal s1 et le deuxième signal s2 sous forme complexe, c'est-à-dire, si s1 = 1s1 xee' et s2 =Is 2ix el92 , la relation précédente s'écrit: Isl2n 2 211-k Ek UkSi S2S2 S1 2 S2 k Ec:ei(k_2d,,-92) g=0 D'où il résulte que : n 2 S2 2)n-1'1 k lk fr'q i(k-2q)(vi_92) 12n IS1 S 2 k e ISI X S =(IS ile +1s 2jee2 )xEcnk q=0 k=0 Une telle formule s'applique notamment pour l'intermodulation d'ordre 5. Dans un tel cas, l'entier n est égal à 2. Ceci s'écrit mathématiquement comme : c2k + is2 2 )2-k Is1 k 21k É c:e j(k-2q)(91-v2) 1S14 X S = s + s 21e92)x± k=0 q=0 De manière plus explicite, l'expression précédente s'écrit : IS14 X S = 142 X 1S2 2 X (1.51 ei(4°1-292) + S2 ei(3492-291)) 2 ±(2142 31S212 142 xls2jel(2e1-'92)+ (21s2 + 3142 )X 1S212 X 1Silei(292-91) Or, il est rappelé que dans ce cas, nous avons P(1S12 )X1S12 (m lsi2 ) 2 XS =a+a Dans ce qui suit, il est étudié la dépendance par rapport à la fréquence instantanée. Comme expliqué précédemment, la correction à retrancher au signal s'écrit sous la forme mathématique suivante : k(s)x s = /412)x Is12 x s Ce résultat, obtenu sous l'hypothèse que la bande instantanée est faible, peut être conservé dans sa forme dans une bande plus large sous la condition que les coefficients du polynôme P dépendent de la fréquence instantanée. Pour cela, il est considéré une dépendance linéaire : k(s)x s = disilx ix Is2 x s Or la fréquence instantanée est, à un facteur multiplicatif 2-rr près, égale à la partie imaginaire de la dérivée logarithmique du signal : 272- s 2g 1 ,(s'j 1 3(s'xs*) 1S12 Ainsi, l'équation relative au produit k(s)xs s'écrit mathématiquement comme : k(s)x s = H2)+ 1Q012/x jx 1s12 x s = [/412)x + Q's.12)x (s* x s') x s 27z- s1 Dans la relation précédente, supposée être calculée dans le domaine du temps discret (c'est-à-dire sur signal échantillonné), s' représente les valeurs de la dérivée du signal complexe aux instants d'échantillonnage. Il s'agit donc du filtrage du signal s par le filtre dérivateur, qui correspond en fait à l'échantillonnage de la dérivée du filtre interpolateur de réponse temporelle en un sinus cardinal. De ce fait, la relation précédente s'exprime comme : Où: Q,s12) - R 27r n - h D (n) =(-1) du fait que la réponse du filtre interpolateur sin WEcHt) est hr v)= Il s'en déduit que la dérivée du filtre f ECHt W - sinWE''t) interpolateur est h; (t)= H t X COS dont les f ECHt2 un valeurs aux points d'échantillonnage sont pour n différent de 0 données a h1 (n I f Ec = 1 x (-1)n full soit hp (n) =(- 1)n parhp(n) = a(fEcHt) ECH Il est à remarquer que la forme de k(s) (déterminée par le terme entre [J) ne dépend pas de la fréquence instantanée f,. Cela signifie que la forme de k(s) est indépendante de la position dans la chaîne de démodulation. Par contre, la valeur des coefficients des polynômes P et Q dépendent de la fréquence instantanée f, de la porteuse du signal, la démodulation par une fréquence foi_ revenant à transférer le terme Q(1s12) x foL sur P(Is12). Bien entendu, le filtre hl:), qui ne dépend pas du signal, ne dépend pas de la valeur de la fréquence instantanée f,. Il est donc constaté que, malgré la dépendance par rapport à la fréquence instantanée, la forme de la correction reste invariante dans la transposition du spectre (même si les coefficients varient) si bien qu'il reste possible d'effectuer la correction après retour en bande de base et décimation, ce qui facilite grandement les opérations. Par exemple, la connaissance des coefficients des deux polynômes P et Q s'obtient par calibrage. Pour cela, il convient d'effectuer deux mesures de calibrage sur deux fréquences instantanées différentes f1 et f2. Pour les deux mesures, il est obtenu par la méthode précédente deux polynômes P1 et P2 a priori différents : JI =P+Qxfi P2 =P+Qxf2 ce qui permet d'obtenir les valeurs de a et de 13 en utilisant le système : = P2 f2 PI fl f2 QP2_Pl = f225 Pour un couple de fréquences fA et fB donné, il est obtenu : , x = a cos e4 + aB cos (9,9 f avec {9A = DfAt ly = a, sin yoA + aB sin 9 B VB = 2111.13t Il en résulte alors que la fréquence instantanée vaut : = 2 27e1S (a, cos (DA + aB cos çoB)(a,f, cos e,+aBfB cos ÇOB)+ (aA sin ço,4 + aB sin 0.8)(a A f A sin ço'+aBfB sin çoB) (a Acos ç 0 A + aB COS 9B) +(a, sin 9A + aB sin çoj C'est-à-dire : a ,2 f A+ a2B + a ,a B(f + fB)cosk- 9B) Ji= 2 2 a,+aB+2a,aB coskoA - çoB) Comme la fréquence instantanée f, varie avec le temps, il est considéré la valeur moyenne de la fréquence instantanée f, pondérée par la puissance du signal pour favoriser les instants où le signal présente l'amplitude la plus grande. En remarquant que la valeur moyenne de cos((pA - cpB) est nulle étant donné que, par hypothèse, la fréquence fA est différente de la fréquence fB, il est obtenu la formule suivante pour la valeur moyenne de la fréquence instantanée f' valeur moyenne notée : (f)= -AJ A ' B) xa2' + a2B + 2a ,a B cos )d ço f ,B2fB AaB(fA + fB)cos(ço, - A B "f , A B 2 0 a2 + a2 + 2a ,aB cos ço 27,( 2 2 a,+aB +2a,aB cos cl)c/ço 2 2 a,+aB Dans le cas où l'amplitude des deux signaux à la fréquence fA et à la fréquence fB sont les mêmes, on a alors aA = aB, ce qui résulte en = (fA + fB) / 2. Les amplitudes de chacune des composante peuvent être calculés par transformée de Fourier rapide ou, compte tenu du faible nombre de points à calculer, par simple transformation de Fourier discrète. Pour des questions de simplicité, il est généré des fréquences sur un pas de transformée de Fourier rapide le plus large possible. Le nombre de points dépend de l'ordre maximum à corriger. Pour une correction jusqu'à l'ordre 2*N+1, le nombre total nFFT de composantes à mesurer est de 2*(N+1). Il est souhaitable que les fréquences des deux calibrages soient le plus éloignées l'une de l'autre tout en conservant un nombre minimum de points de transformée de Fourier rapide. Le cas extrême consiste à prendre : - nFFT = 2(N+1) d'une part, et - fA = ±fEcH/2 et fB = ±N/(N+1) x fEcH/2 d'autre part En variante, il est également possible de rapprocher les fréquences pour être moins sensible à la réponse du filtre de décimation qui atténue inévitablement les bords de la bande d'échantillonnage. Ainsi, lorsque seuls les termes d'ordre 3 sont à corriger, il est possible de considérer une transformée de Fourier rapide à 4 points dont le calcul est extrêmement simple avec fA = ±fEcHi4 et fB = O. Dans ce cas les produits d'intermodulation sont situés à ±fEcH/2 et -fA. Simulations Des simulations sont représentées aux figures 3 à 10. Les simulations qui suivent ont été réalisées à partir de séries de mesures 2 tons d'amplitudes respectives voisines de -6 dBFS (DéciBell à pleine échelle) chacune. De telles mesures ont été effectuées sur un convertisseur analogique-numérique autour de trois points de fonctionnement fo différents voisins des fréquences 2,1*fEcH, 2.2*fEcH et 2.3 fEcH. Les deux tons séparés de fEcH/256 étaient situés respectivement à fA = fo - fEcEi/512 et fB = fo + fEcF1/512 avec, suivant les cas : - fo = ( 2 + 103/1024 ) fEcH , premier point voisin de 2.1 fEcH, - fo = ( 2 + 205/1024 ) x fEcH , second point voisin de 2.2 fEcH, et - fo = ( 2 + 307/1024 ) x fEcH , troisième point voisin de 2.3 fEcH Pour chaque point de fonctionnement, quatre mesures différentes ont été effectuées correspondant à quatre conditions de phase différentes entre signaux et entre signal et horloge. Une mesure du premier point de fonctionnement et une mesure du troisième point de calibrage sont utilisées comme mesures de calibrage pour un traitement visant à corriger les défauts jusqu'à l'ordre 5. Plusieurs simulations du traitement calibrage/correction ont été effectuées à partir de ces mesures. Pour les cas des figures 3 à 5, il a été effectué un traitement calibrage/correction à pleine vitesse après filtrage de Hilbert sans translation du spectre ni décimation. Pour chaque cas, les courbes représentent le maximum de l'amplitude obtenue sur les quatre mesures. Chacune des figures 3 à 5 illustrent respectivement la variation de l'amplitude en fonction de la fréquence pour les trois points de fonctionnement différents. Pour les cas des figures 6 à 8, il a été effectué un traitement calibrage/correction à pleine vitesse après filtrage de Hilbert puis translation de -205/1024 fEcH mais sans décimation. Pour chaque cas, les courbes représentent le maximum de l'amplitude k2z2 obtenue sur les quatre mesures. Chacune des figures 6 à 8 illustrent respectivement la variation de l'amplitude en fonction de la fréquence pour les trois points de fonctionnement différents. Pour les cas des figures 9 à 11, il a été effectué un traitement calibrage/correction à pleine vitesse après filtrage de Hilbert puis translation de -205/1024 fEcH avec décimation par 4. Pour chaque cas, les courbes représentent le maximum de l'amplitude obtenue sur les quatre mesures. Chacune des figures 9 à 11 illustrent respectivement la variation de l'amplitude en fonction de la fréquence pour les trois points de fonctionnement différents. Pour les simulations, les filtres de Hilbert et de décimation utilisés sont des filtres idéaux. Ainsi, dans le cas d'un signal répétitif, donc mesuré sur une durée finie multiple de sa période (1024 points en l'occurrence), il y a lieu de replier le filtre de dérivation dont la durée est infinie ; si N est la durée totale. Cela s'exprime mathématiquement comme : hp(n)= (-1)n devient '171,0(n)= L , n+kxN C'est-à-dire, en utilisant =-n que l'expression précédente devient: îlio(n)_ 4. -11 (-1)n (-1)n (-1)n + 27Z- *-c-; n k+4...con+lexAr k=1 n+lcxN n N k=1 -1z-2 n N k=171-2e Or: 1 = 1+ 2u (-1)" Li tan u u u2 k2z2 D'où : ^' (-1)n (-1)r ( 1 1 (-1)n Tt- hp(n)= + x n N tan ral Nx tan - Conclusion Il est proposé un procédé de correction des défauts d'intermodulation d'une chaîne de réception comportant un convertisseur analogique-numérique et une chaîne de réception mettant en oeuvre un tel procédé. Le procédé permet de trouver un mécanisme de compensation travaillant sur signaux démodulés, filtrés, décimés et ne compensant que l'intermodulation sans introduire de distorsion harmonique non désirée. Ainsi, le procédé de correction des défauts d'intermodulation introduits par un circuit de conversion permet, avec une mise en oeuvre aisée, de diminuer les défauts d'intermodulation sans générer de la distorsion harmonique non désirée. REVENDICATIONS1.- Procédé de correction des défauts d'intermodulation introduits par un circuit de conversion (14) adapté à convertir un signal analogique d'entrée en un signal numérique de sortie, le procédé comportant les étapes de : - fourniture du signal numérique de sortie, - génération d'une fonction de correction, la fonction de correction étant une fonction complexe du signal numérique de sortie, et - correction du signal numérique de sortie à l'aide de la fonction de correction, ptur obtenir un signal numérique de sortie corrigé. 2.- Procédé selon la revendication 1, dans lequel, le signal numérique de sortie corrigé est égal à la différence le signal numérique de sortie et le produit du signal numérique de sortie par la fonction de correction. 3.- Procédé selon la revendication 1 ou 2, dans lequel, l'étape de génération comporte les étapes de : - calcul d'une première fonction polynômiale du carrée du module du signal numérique de sortie, pour obtenir un premier terme final, - calcul d'une deuxième fonction polynômiale du carrée du module du signal numérique de sortie selon une représentation complexe, pour obtenir un deuxième terme intermédiaire, - application d'un filtre sur le signal numérique de sortie pour obtenir un signal filtré, - calcul du produit du deuxième terme intermédiaire par le produit du complexe conjugué du signal numérique de sortie par le signal filtré, pour obtenir un deuxième terme final, et - calcul de la somme du premier terme final et du deuxième terme final, pour obtenir la fonction de correction. 4.- Procédé selon la revendication 3, dans lequel, le filtre appliqué à l'étape d'application d'un filtre restitue la dérivée du signal numérique de sortie. 5.- Procédé selon l'une quelconque des revendications 1 à 4, dans lequel l'étape de génération comporte l'application d'au moins une opération au signal numérique desortie, pour obtenir un signal de calcul, l'étape de correction comportant l'application de la fonction de correction au signal de calcul. 6.- Procédé selon la revendication 5, dans lequel au moins une opération est choisie dans le groupe constitué de: - un filtrage de Hilbert du signal numérique de sortie, - une démodulation numérique du signal numérique de sortie, - un filtre à phase linéaire, et - une décimation. 7.- Procédé selon l'une quelconque des revendications 1 à 6, dans lequel l'étape de génération comporte les étapes de : - application d'un signal de calibration au circuit de conversion (14), par exemple un signal à deux fréquences distinctes, - mesure du signal obtenue en sortie du circuit de conversion (14), pour obtenir un signal mesuré et - analyse des composantes spectrales du signal mesuré. 8.- Produit programme d'ordinateur comportant un programme d'ordinateur comprenant des instructions de programme, le programme d'ordinateur étant mémorisé sur un support d'informations lisible par un ordinateur et le programme d'ordinateur étant chargeable sur un ordinateur et adapté pour entraîner la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à 7 lorsque le programme d'ordinateur est mis en oeuvre sur l'ordinateur. 9.- Appareil électronique (10) comportant : - un récepteur (12) d'un signal analogique d'entrée, - un circuit de conversion (14) adapté à convertir le signal analogique d'entrée en un signal numérique de sortie, et - un correcteur (16) adapté à générer d'une fonction de correction, la fonction de correction étant une fonction complexe du signal numérique de sortie, et à corriger le signal numérique de sortie à l'aide de la fonction de correction, pour obtenir un signal numérique de sortie corrigé. 10.- Appareil selon la revendication 9, dans lequel le correcteur (16) est un circuit logique programmable.