FR 2766992 A1 19990205 FR 9710052 A 19970801 PROCEDE DE SIMULATION D'AMPLIFICATEUR NON LINEAIRE A MEMOIRE D'ENVELOPPE La présente invention concerne le domaine de la simulation de réponse en signal d'amplificateurs non linéaires et a pour objet un système de simulation du signal de réponse d'un amplificateur non linéaire présentant un effet de mémoire. Une application d'un tel système est la simulation de l'amplification hyperfréquence à haut rendement, notamment en classe AB, B ou C d'amplification et plus particulièrement la simulation de réponse d'amplificateurs de puissance à l'état solide, de type SSPA, et de dispositifs amplificateurs à base de tubes à ondes progressives, de type TOP ou TWT, utilisés dans les liaisons de transmission par faisceau hertzien terrestre ou satellite. A haute fréquence et à haut rendement, de tels dispositifs amplificateurs présentent une caractéristique de réponse non-linéaire. La figure 1 illustre un exemple de courbe de réponse entrée-sortie d'amplificateur non linéaire ANL. La courbe donnant le niveau de signal de sortie g de l'amplificateur ANL en fonction du niveau de signal d'entrée x est typiquement infléchie aux fortes amplitudes A du signal d'entrée x, à cause de phénomènes de saturation. Lorsque l'amplificateur est utilisé dans des conditions telles que le gain n'est pas constant en fonction du niveau de signal d'entrée x, on parle d'amplificateur fonctionnant en régime non linéaire ou plus simplement d'amplification non linéaire. Parmi les dispositifs non-linéaires, on distingue des dispositifs sans mémoire, quasi sans mémoire et à mémoire. Les amplificateurs sans mémoire présentent une forte non-linéarité en amplitude et une moindre distorsion de phase. La caractéristique de réponse entree/sortie d'un amplificateur sans mémoire ANL peut alors se réduire, comme illustrée figure 1, à une seule courbe g(x). On peut modéliser et simuler la réponse en fonction du temps t d'un amplificateur sans mémoire ANL à un signal d'entrée x sinusoïdal, de fréquence f0, modulé en amplitude A et en phase ó, signal ayant donc la forme suivante: x(t) = A(t).cos(2n.f0.t+ó(t)) (1) A(t) représente l'enveloppe du signal d'entrée, qui est définie par les limites en amplitude dans lequel le signal sinusoïdal x évolue, l'enveloppe variant en fonction du temps. La figure 3 illustre un chronogramme d'un signal x(t) ayant une enveloppe A constante. Le signal de sortie g de l'amplificateur sans mémoire est alors de la forme suivante: g(t) = C(A(t)).cos(2n.f0.t+P(t)) (2) Il est intéressant de donner une enveloppe complexe à ces signaux x et g, en abandonnant toute référence au temps t puisque les amplificateurs sans mémoire ont une réponse instantanée. Le signal d'entrée x de l'équation (1) a une enveloppe complexe X de la forme suivante: X = A.exp(j.p) (3) Toute l'information utile de modulation d'amplitude A(t) et de phase 0(t) est inscrite dans cette enveloppe complexe X. Le signal de sortie g de l'équation (2) a de même une enveloppe complexe G de la forme suivante: G = C(A).exp(j.ó) (4) Pour un amplificateur ANL non linéaire sans mémoire, on démontre que C(A) est la transformée de Chebyshev de la courbe de réponse entrée/sortie g(x). La réponse d'un amplificateur ANL non-linéaire sans mémoire à un signal x modulé peut donc être modélisée et simulée simplement par une seule courbe C(A), dont un exemple est illustré en trait plein sur la figure 6A. Une telle courbe donnant l'amplitude C de sortie en fonction de l'amplitude A d'entrée est appelée courbe de non- linéarité en amplitude et notée en abrégé courbe AM/AM. La figure 2 illustre une caractéristique de réponse d'un amplificateur non-linéaire à mémoire ANLAM, sur laquelle apparaît un phénomène d'hystérésis provoqué par effet de mémorisation. On voit que les courbes de montée m et m' d'hystérésis ne se superposent pas lorsque les amplitudes A et A' respectives du signal d'entrée x sont différentes. La variation du temps de mémorisation liée à la variation d'amplitude empêche les courbes m et m' de se superposer. Lorsque le temps de mémorisation de l'amplificateur ANLAM est négligeable par rapport à la période de variation d'amplitude A(t), on peut encore considérer que l'amplitude A est stable et on parle d'amplificateur quasi sans mémoire. La figure 4 représente un chronogramme de signal de sortie y d'un amplificateur quasi sans mémoire auquel est appliqué le signal d'entrée x illustré sur le chronogramme de la figure 3 dont la formule est rappelée ci-après: x(t) = A(t).cos(27.f0.t+ó(t)) (1) Pour un amplificateur non-linéaire quasi sans mémoire auquel on applique le signal x de l'équation (1), le signal de sortie y prend la forme suivante: y(t) = C(A).cos(2.fo.t-D(A)+O(t)) (5) o C(A) est l'amplitude du signal de sortie y, et $(A) est le déphasage du signal de sortie y, qui dépendent de l'amplitude A(t) du signal d'entrée x. Ainsi, à un instant t donné, l'amplitude C(A) et le déphasage $(A) du signal de sortie y dépendent uniquement de l'amplitude A du signal d'entrée x à cet instant t. On peut donc négliger les variations d'amplitude A(t) en fonction du temps, et considérer quasiment que l'amplitude A est constante comme visible figure 3. Il est également intéressant d'écrire sous forme complexe les enveloppes des signaux x et y exprimés précédemment, enveloppes X et Y qui prennent les formes respectives suivantes: X = A.exp(j.ó) (3) (enveloppe A considérée comme constante dans le temps) Y = C(A).exp(j.0-D(A)) (6) La réponse d'un amplificateur non-linéaire quasi sans mémoire peut donc être modélisée et simulée simplement en connaissant deux courbes caractéristiques suivantes: - une courbe C(A) donnant l'amplitude C du signal de sortie y en fonction de l'amplitude A du signal d'entrée, illustrée par exemple par la courbe AM/AM en trait plein de la figure 6A (courbe de conversion amplitude/amplitude). - une courbe 4(A) donnant le déphasage b du signal de sortie y par rapport au signal d'entrée x, en fonction de l'amplitude A du signal x, dite courbe de conversion amplitude/phase, en abrégé AM/PM, dont un exemple est illustré en trait plein sur la figure 6B. On démontre que la courbe C(A) est, de façon similaire, la norme de la transformée complexe de Chebychev de la caractéristique de réponse y(x), la courbe D(A) étant l'argument de la transformée complexe. L'enveloppe complexe Y du signal de sortie peut également être écrite sous forme de deux parties, réelles et imaginaires, correspondant à une composante P en phase et à une composante Q en quadrature de phase, composantes P et Q de formes suivantes: P(A) = C(A). cos(b(A)) (7') O(A) =+ (A).sin(sb(A)) (7"1) La figure 7 illustre un exemple de courbes P(A) et Q(A) équivalentes aux courbes C(A) et et 6B. Les modèles connus de simulation de réponse d'amplificateur quasi sans mémoire préfèrent généralement utiliser des caractéristiques sous forme de paire de courbes P(A) et Q(A) plutôt que sous forme de paire de courbes C(A) et @(A), bien que ces paires de courbes soient strictement équivalentes. Un principe connu de simulation d'amplificateur non- linéaire consiste à précaractériser sur un banc de test l'amplificateur réalisé, la précaractérisation étant effectuée avec un signal de fréquence et d'amplitude déterminées, afin de simuler ensuite la réponse à un signal de fréquence et d'amplitude quelconques. Comme schématisé figure 5, un signal monofréquence f0 prenant différentes amplitudes A', A", A"'. est appliqué à l'amplificateur testé pour obtenir ses caractéristiques, illustrées par exemple en figure 6 ou 7. Cependant, pour des amplificateurs présentant une certaine quantité de mémoire, on constate que les caractéristiques varient de façon importante selon la fréquence f0-, f0 ou f0+ du signal à amplifier. Un système connu de simulation de tels amplificateurs a été exposé par H.B. Poza dans un article intitulé "A Wideband data link computer simulation model", publié dans les actes de conférence "NAECON'75 record", page 71. Il propose d'établir plusieurs paires de courbes AM/AM et AM/PM pour plusieurs fréquences f0-, f0 ou f0+ de fonctionnement de l'amplificateur. Sur les figures 6A et 6B, on a ainsi représenté un exemple de trois paires de courbes AM/AM et AM/PM, obtenues + respectivement à trois fréquences f0-, f0 et f0 situées dans la bande utile BU d'un amplificateur de liaison hertzienne. Le simulateur de H.B. Poza stocke seulement une paire de courbes AM/AM et AM/PM, par exemple la paire de courbes obtenues à la fréquence f0, et reconstitue les autres paires de courbes (non stockées) correspondant aux -+ autres fréquences f0, f+0 ou à des fréquences intermédiaires. Une courbe non stockée est déduite simplement en translatant la courbe stockée d'un vecteur adéquat. Le simulateur calcule les composantes selon l'axe (A) et selon l'axe (C) ou (ó) du vecteur de translation pour rapprocher au mieux la courbe de fréquence f0, représentée en trait plein à la figure 6A ou 6B, de la courbe de fréquence f0o- ou f0+, représentée en pointillés. Un tel système de simulation est trop approximatif pour simuler les distorsions apparaissant sur un amplificateur quasi sans mémoire à différentes fréquences. Un autre inconvénient d'un tel système est qu'il ne permet pas de simuler la réponse d'un amplificateur à mémoire. La figure 9 illustre un autre système connu de simulation, selon le modèle de A.A.M. Saleh, décrit dans un article intitulé "Frequency-independent and frequency- dependent models of TWTA Amplifiers" paru en novembre 1981 dans la revue "IEEE transactions on communication", volume com-29, n 11, page 1715. Le calcul de la réponse à un signal d'entrée x d'un amplificateur non-linéaire quasi sans mémoire est décomposé en deux étapes P(A,f) et Q(A,f) de calcul de deux composantes yp et yq respectives du signal de sortie y. La composante yp est en phase avec le signal d'entrée x, tandis que la composante yq est en quadrature de phase. Le système de simulation de Saleh utilise des résultats de caractérisation de l'amplificateur à plusieurs fréquences f, plusieurs paires de courbes P(A) et Q(A) étant établies à plusieurs fréquences f pour calculer des fonctions de transfert P(A,f) et Q(A, f) correspondant à chaque branche de calcul des composantes yp et yq du signal de sortie y. Chacune des fonctions de transfert P(A,f) ou Q(A,f) correspond au calcul d'une réponse non linéaire sans effet de mémoire, aucun déphasage n'étant introduit dans chacune de ces branches de calcul. L'article de A.A.M. Saleh signale que le modèle s'applique à l'amplification de signaux mono-fréquence et, par conjecture, suppose qu'il serait applicable à la simulation de signaux de forme quelconque. Un autre inconvénient de ce modèle est qu'il ne s'applique pas aux amplificateurs à mémoire, le système ne rendant pas compte des distorsions apparaissant à différentes fréquences lorsque l'amplitude du signal varie rapidement par rapport aux constantes de temps de mémorisation. Pour les amplificateurs non linéaires à mémoire ANLAM, les effets de mémorisation sont plus importants, les temps de mémorisation n'étant pas négligeables par rapport aux temps de variation de l'amplitude A(t). Dans de tels cas, une méthode plus complexe telle que la résolution d'équations différentielles non linéaires ou le développement des caractéristiques en suite de fonctions est nécessaire pour simuler la réponse de l'amplificateur à mémoire ANLAM. Un modèle connu plus perfectionné est proposé par M.T. Abuelma'atti dans un article intitulé "Frequency- dependent nonlinear quadrature model for TWT amplifiers", paru en août 1984, dans la revue IEEE transaction on communication, volume com-32, N 8, page 982. La modélisation utilise un développement des courbes caractéristiques d'un amplificateur non linéaire en suite de fonctions de Bessel, ce qui permet de simuler la réponse d'amplificateurs à mémoire. La figure 10 illustre le système de simulation décrit par M.T. Abuelma'atti qui comporte comme précédemment deux branches de calcul des composantes en phase yp et en quadrature yq du signal de sortie y, chaque branche calculant la contribution d'une suite de N fonctions de Bessel J1 dépendant de l'amplitude A du signal d'entrée x, en pondérant chaque fonction J1 d'un coefficient a et d'un facteur G(f) de correction en fréquence. Les coefficients an et les facteurs Gn(f) sont calculés après avoir établi plusieurs paires de courbes P(A) et Q(A) caractéristiques de l'amplificateur non- linéaire, les courbes étant établies à plusieurs fréquences f de test de l'amplificateur. Avec la figure 8, on peut voir que deux sommes pondérées de N fonctions de Bessel du premier type d'ordre 1, notées Jl(nzA/D) permet d'interpoler très précisément des courbes P(A) et Q(A), comme celles de la figure 7, en ajustant les coefficients anp et anq de pondération de ces fonctions J1. On peut noter que les fonctions de Bessel sont les transformées de Chebychev de fonctions sinusoïdales et que le développement en suite de fonctions de Bessel des courbes P(A) et Q(A) correspond à un développement de Fourier en suite de fonctions sinusoïdales des courbes y(x) de non linéarité de l'amplificateur, ce qui convient élégamment à l'allure sinusoïdale des courbes d'hystérésis y(x), tel qu'illustrées figure 2. Théoriquement, un tel système devrait permettre la simulation de l'amplification de signaux multiporteuses, c'est-à-dire de signaux comportant plusieurs composantes sinusoïdales de fréquence distinctes. Toutefois, l'amplification non-linéaire d'un signal d'entrée multiporteuses est compliqué par l'apparition de distorsions appelées phénomènes d'intermodulation. Lorsqu'un amplificateur non-linéaire reçoit en entrée plusieurs fréquences porteuses, on obtient en sortie, en plus des porteuses amplifiées, des harmoniques non désirées appelées produits d'intermodulation, chaque harmonique ayant une fréquence distincte des fréquences des porteuses. Les figures 11-14 illustrent l'apparition du phénomène d'intermodulation lors de l'amplification non- linéaire d'un signal à deux porteuses. La figure 11 représente un diagramme fréquentiel de signal d'entrée x biporteuse, sur lequel on voit par exemple, deux composantes porteuses de fréquences respectives f-1 et fl ayant une égale amplitude A d'entrée. Les deux fréquences porteuses f-1 et fl, situées dans la bande utile BU de fréquence de l'amplificateur, sont séparées par un écart fréquentiel df. La figure 14 représente un diagramme fréquentiel de signal de sortie y qui correspond à l'amplification ANLAN avec effet de mémoire du signal d'entrée x biporteuse de la figure 11. On voit que le signal de sortie y comporte une série d'harmoniques de diverses fréquences f et d'amplitudes C différentes. La fréquence de chacune des composantes d'intermodulation est une combinaison entière des fréquences porteuses en entrée. Le détail des harmoniques de sortie comprises dans la bande utile BU, illustrées figure 14, est le suivant: - les porteuses, de fréquence f-1 et fl et d'amplitude de sortie respectives C_1 et C1. - les composantes d'intermodulation d'ordre 3, de fréquences f-3, f3, et d'amplitudes de sortie respectives C-3, C3- La figure 15 reproduit des résultats d'après le modèle d'Abuelma'atti qui a permis une estimation des amplitudes C_1, C1, et C_3, C3 de composantes porteuses et de composantes d'intermodulation d'ordre 3. Toutefois, l'estimation de l'amplitude des composantes d'intermodulation par le modèle d'Abuelma'atti correspond mal à la réalité des mesures de distorsions d'intermodulation sur les amplificateurs non- linéaires en fonctionnement multiporteuses à haut rendement. Il se produit en effet un phénomène général, appelé mémoire d'enveloppe lorsqu'on amplifie un signal multiporteuses. Pour ce type de signal, l'enveloppe, définie par les limites positive et négative en amplitude du signal, n'est pas constante. Dans ce cas, on ne peut plus considérer que l'enveloppe est constante comme dans les modèles quasi- sans mémoire. En fait, les constantes de temps de mémorisation ne sont plus négligeables vis-à-vis du temps de variation de l'enveloppe. La figure 12 qui illustre l'évolution temporelle du signal d'entrée x biporteuses de la figure 11, montre par exemple que l'enveloppe X(t) et -X(t) du signal biporteuse varie très rapidement et dans des proportions importantes, alors que l'amplitude A de chaque porteuse fl et f-1 porteuse est supposée constante. La figure 13 qui illustre l'évolution temporelle du signal de sortie y correspondant au signal d'entrée x biporteuse précédent, montre ainsi que l'enveloppe Y(t) et -Y(t) du signal de sortie y, apparaît alors très déformée par les distorsions d'intermodulation. On désigne sous le nom d'effet de mémoire d'enveloppe de telles déformations d'amplitude du signal par un amplificateur non-linéaire à mémoire. Les modèles connus ne rendent pas compte, en particulier, de deux effets illustrés aux figures 16 et 17. Le premier effet observé est que le ratio C1/C3 illustré figure 16, comparant l'amplitude Cl de porteuses par rapport à l'amplitude C3 des composantes d'intermodulation, varie considérablement selon l'écart fréquentiel df des porteuses et selon la valeur précise de l'amplitude A des porteuses du signal d'entrée. En appliquant le modèle d'Abuelma'atti, le ratio C1/C3 ne dépend pas de l'écart fréquentiel et varie continûment selon l'amplitude A d'entrée. Le second effet observé est que la présence d'une seconde porteuse influence l'amplitude de sortie de la première porteuse. La figure 17 montre par exemple des pics d'affaiblissement ou de résonnance d'amplitude Cl ou C-1 de sortie de chacune des porteuses de fréquence fi ou f_ 1, selon l'écart fréquentiel df des deux porteuses (df valant f1-f-1) Un inconvénient des modèles connus est qu'ils ne rendent pas compte de tels effets. De, façon générale, l'inconvénient des modèles connus est qu'ils ne simulent pas d'effet de mémoire d'enveloppe. Un but de l'invention est de prévoir une simulation précise de l'effet de mémoire d'enveloppe sur des amplificateurs non-linéaires. Un but particulier de l'invention est de prévoir un procédé et un dispositif permettant de simuler de façon précise la réponse d'un amplificateur non-linéaire à des signaux multiporteuses. Un autre but de l'invention est de prévoir un procédé et un dispositif simples de simulation d'amplificateurs non-linéaires à mémoire. Selon l'invention ces buts sont atteints en prévoyant une mesure dynamique des caractéristiques de l'amplificateur dans les conditions ou apparaissent les distorsions. Il est prévu, selon l'invention, de mesurer les caractéristiques de distorsion de modulation d'amplitude, en prévoyant, à chaque mesure, d'appliquer un signal modulé en amplitude, ce qui permet de faire varier l'enveloppe du signal. Les caractéristiques de distorsions permettent de corriger les réponses simulées à partir de caractéristiques établies à amplitude constante. Il est prévu notamment une mesure de caractéristiques en appliquant un signal multiporteuse, signal qui est donc modulé en amplitude. L'invention est réalisée en mettant en oeuvre un procédé de simulation de réponse en signal d'un amplificateur non-linéaire présentant un effet de mémoire, le procédé comportant des étapes parmi les étapes consistant à: mesurer des caractéristiques de conversion amplitude/amplitude et de conversion amplitude/déphasage de l'amplificateur, chaque mesure étant effectué à amplitude constante en entrée, et - développer les caractéristiques en suite de fonction de transfert directes, le procédé étant caractérisé par des étapes consistant à: - mesurer des caractéristiques de distorsion de modulation d'amplitude, chaque mesure étant effectuée en modulant l'amplitude d'entrée, - calculer des fonctions de transfert de modulation reproduisant en sortie les amplitudes de distorsion selon les amplitudes de modulation d'entrée et corrigeant les fonctions de transfert directes lorsque l'amplitude d'entrée est modulée, afin de simuler l'effet de mémoire d'enveloppe. De préférence, le procédé selon l'invention comporte des étapes intercalaires consistant à: mesurer les caractéristiques à différentes fréquences, et - calculer des correcteurs en fréquence pour les fonctions de transfert directes. Il est prévu, selon l'invention, qu'une mesure de caractéristiques de distorsion de modulation d'amplitude consiste à mesurer des distorsions d'intermodulation en appliquant en entrée un signal comportant au moins deux porteuses de fréquences différentes. Une première réalisation du procédé selon l'invention prévoit de: mesurer des caractéristiques de réjection de composantes d'intermodulation en sortie, chaque mesure étant effectuée avec un groupe de porteuses en entrée d'amplitudes déterminées. Une deuxième réalisation du procédé selon l'invention prévoit de: mesurer des caractéristiques d'interaction de porteuses en sortie, chaque mesure étant effectuée avec un groupe de porteuses en entrée d'amplitudes déterminées. Une troisième réalisation du procédé selon l'invention prévoit de.: mesurer des caractéristiques de bruit de modulation en sortie, chaque mesure étant effectuée avec des amplitudes déterminées de modulation en entrée. De préférence, il est prévu selon l'invention de: - mesurer les caractéristiques de distorsion avec plusieurs écarts fréquentiels de modulation en entrée. De préférence, il est prévu selon l'invention de: - mesurer les caractéristiques de distorsions d'intermodulation avec plusieurs valeurs d'amplitudes de modulation en entrée. Selon l'invention, il est prévu en outre un dispositif de simulation de réponse en signal d'un amplificateur non-linéaire présentant un effet de mémoire, le dispositif comportant: - un module de calcul de la réponse d'un amplificateur non-linéaire à un signal monoporteuse, le module comportant au moins un filtre ayant une fonction de transfert établie à partir de caractéristiques réelles de l'amplificateur, notamment des caractéristiques de non-linéarité d'amplitude et de conversion amplitude/déphasage, la fonction de transfert étant corrigée en fréquence par un correcteur calculé à partir de caractéristiques mesurées à différentes fréquences de signal monoporteuse, caractérisé en ce que le dispositif comporte des moyens de génération d'harmoniques en réponse à un signal multiporteuses. Les moyens de génération d'harmoniques du dispositif comportent de préférence selon l'invention un module de décomposition de Fourier d'un signal en composantes fréquentielles et/ou un module combinatoire calculant les combinaisons de fréquence de composantes d'intermodulation et/ou un module normatif calculant une valeur efficace d'amplitude de signal. Un premier mode de réalisation de dispositif selon l'invention prévoit un filtre de modulation de la réponse du module de calcul, le filtre de modulation ayant une fonction de transfert établie à partir de caractéristiques de réjection de composantes d'intermodulation lorsque l'amplification d'un signal multiporteuses est simulée par le dispositif. Un deuxième mode de réalisation de dispositif selon l'invention prévoit un filtre de correction de la réponse du module de calcul, le filtre ayant une fonction de transfert établie à partir de caractéristiques d'interaction de porteuses, de sorte que l'amplitude d'une porteuse fournie par le module de calcul est corrigée par le filtre de correction lorsque l'amplification d'un signal multiporteuse est simulée par le dispositif. De façon avantageuse, l'invention s'applique à la simulation de réponse en signal d'un amplificateur hyperfréquence fonctionnant à haut rendement. D'autres caractéristiques, buts et avantages de la présente invention apparaîtront à la lecture de la description détaillée qui va suivre, et des dessins annexés, donnés à titre d'exemples non limitatifs et sur lesquels: - les figures 1 et 2, précédemment décrites, illustrent des courbes de réponse d'amplificateurs non- linéaires; - les figures 3 et 4, précédemment décrites, illustrent des diagrammes temporels de réponse en signal d'amplificateur non-linéaire sans mémoire ou quasi sans mémoire; - la figure 5, précédemment décrite, schématise un diagramme fréquentiel de mesures de caractéristiques selon un principe connu; - les figures 6A et 6B, précédemment décrites, illustrent des caractéristiques d'amplificateur non- linéaire sous forme de courbes de conversion amplitude/amplitude et de courbes de conversion amplitude/déphasage, obtenues selon un procédé connu; - la figure 7, précédemment décrite, illustre des caractéristiques d'amplificateur non- linéaire, sous forme de courbes de non-linéarité en phase et en quadrature, obtenues selon un procédé équivalent connu; - la figure 8, précédemment décrite, illustre une suite de fonction de Bessel de décomposition des caractéristiques mises en oeuvre dans un procédé de simulation connu; - les figures 9 et 10, précédemment décrites, représentent deux dispositifs de simulation connus; - les figures 11 et 14, précédemment décrites, illustrent des diagrammes fréquentiels d'amplitude de signaux multiporteuses, respectivement en entrée et en sortie d'un amplificateur non-linéaire à mémoire; - les figures 12 et 13, précédemment décrites, illustrent des diagrammes temporels d'amplitude de signaux multiporteuses, respectivement en entrée et en sortie d'un amplificateur non-linéaire à mémoire; - la figure 15, précédemment décrite, reproduit des résultats de simulation d'amplitudes de composantes fréquentielles de sortie d'un amplificateur non-linéaire selon l'amplitude des porteuses d'un signal d'entrée biporteuse, simulation effectuée selon un procédé connu; - la figure 16 illustre des mesures de réjection de composante d'intermodulation en sortie par rapport aux composantes porteuses en entrée d'un amplificateur non- linéaire à mémoire exprimé sous forme de ratio C1/C3 d'amplitude de sortie de porteuse/amplitude de composante d'intermodulation, en fonction de l'amplitude A en entrée des porteuses, pour divers écarts fréquentiels df des porteuses, mesures mises en oeuvre selon l'invention; - la figure 17 illustre des mesures d'interaction de porteuses en sortie d'un amplificateur non-linéaire à mémoire, sous forme d'amplitude de sortie y d'une porteuse en fonction de l'écart fréquentiel des porteuses, mesures mises en oeuvre selon l'invention; - la figure 18 représente un diagramme fréquentiel de mesures de caractéristiques mis en oeuvre selon l'invention; - Les figures 19A et 19B représentent des exemples de densité spectrale de puissance (DSP) du signal utilisé en entrée de l'amplificateur pour la mesure du NPR; - la figure 20 représente une réalisation de dispositif de simulation d'amplification non-linéaire à mémoire comportant des moyens de génération d'harmoniques selon l'invention; - la figure 21 représente un premier mode de réalisation dedispositif de simulation d'amplification non-linéaire à mémoire selon l'invention; - la figure 22 représente une variante du premier mode de réalisation de dispositif de simulation d'amplification non-linéaire à mémoire selon l'invention; - la figure 23 représente un deuxième mode de réalisation de dispositif de simulation d'amplification non-linéaire à mémoire selon l'invention; - la figure 24 représente un banc de mesure de caractéristiques de distorsions d'intermodulation utilisé par le procédé ou le dispositif selon l'invention; - la figure 25 représente un diagramme de fichiers de caractéristiques mis en oeuvre par un processus de simulation de réponse d'un amplificateur non-linéaire selon l'invention; et, - la figure 26 représente un schéma d'étapes de processus de simulation de réponse d'un amplificateur non-linéaire selon l'invention. Les procédé et dispositif de simulation selon l'invention sont développés à partir des modélisations connues. L'invention se propose de simuler l'amplification de signaux de forme quelconque et notamment de signaux comportant plusieurs composantes fréquentielles. Selon une première caractéristique de dispositif selon l'invention, le dispositif de simulation comporte des moyens de génération d'harmoniques en réponse à un signal multiporteuses. Ainsi les diverses structures, illustrées aux figures 21 à 23, de dispositif selon l'invention, comportent systématiquement un bloc 11 ou 21 ou 31 de calcul de signal de sortie y selon un modèle connu et un bloc parallèle 12 ou 22, 23 ou 32 de génération d'harmonique qui calcule les fréquences des composantes d'intermodulation avant le calcul de leur amplitude précise pour corriger le signal de sortie y. Les blocs de calcul de signal de sortie 11, 21 et 31 sont prévus pour mettre en oeuvre un modèle de simulation connu, par exemple et respectivement le modèle non linéaire sans mémoire MNLSM de H.B. Poza, le modèle non linéaire sans mémoire MLLSM quadratique de A.A.M. Saleh et le modèle à mémoire MM de M.T. Abuelma'atti. D'autres modèles peuvent être mis en oeuvre, et la description de la mise en oeuvre de l'invention sera développée essentiellement à partir du modèle perfectionné de M.T. Abuelma'atti. On notera que la structure en parallèle des blocs de calcul de sortie MNLSM ou MM et des blocs Harm de génération d'harmoniques est indicative, la structure des blocs pouvant être en série, et la réalisation effective des blocs de calculs étant de préférence effectuée par des programmes d'ordinateur, les blocs schématisant dans ce cas les fonctions de sous-programmes de calcul. La figure 23 décrit le détail d'une réalisation alternative de moyens Harm de génération d'harmoniques et de moyens MM de calcul des composantes d'intermodulation perfectionnée à partir du modèle d'Abuelma'atti simple illustré figure 10. Les modules Harm et MM permettent d'effectuer le' calcul des amplitudes de composantes d'intermodulation..., Yo1o, Y1-11, Y-120, Y-111, Y0-12, Y-102, Y20-1 Y2Y2-10, Y1l-1, Y02-1,... lorsque le signal d'entrée x comporte une multitude de composantes fréquentielles A, B,..., M. Cette réalisation complexe de moyens Harm et MM est destinée à la simulation de l'amplification de signaux de forme quelconque, en particulier de signaux multiporteuses. On note sur la figure 20 que les moyens Harm de génération d'harmoniques comportent un filtre FFT de décomposition en série de Fourier pour décomposer le signal d'entrée x multifréquence en composantes fréquentielles d'amplitude A, B,..., M. Ces composantes sont appliquées en entrée de filtres élémentaires groupés par tableaux correspondant au calcul d'une composante d'intermodulation, par exemple Yo10. Chaque filtre élémentaire a une fonction de transfert issue d'une fonction de Bessel J généralisée sur l'exemple des filtres de la figure 10. Le détail des fonctions de transfert des tableaux de filtre et des fonctions de Bessel généralisé sera détaillé suite au développement ci-après. On retiendra que l'architecture d'un dispositif selon l'invention comporte des moyens de génération d'harmoniques avec de préférence un module de décomposition en série de Fourier et un modèle de simulation connu. Un développement des formules proposées par le modèle connu le plus perfectionné, le modèle d'Abuelma'atti aboutit à des erreurs de simulation, comme démontré ci-après, erreurs qui seront corrigées selon l'invention. En effet, suivant les équations précédentes 7' et 7", Abuelma'atti décompose isolément les caractéristiques de non-linéarité en phase P et en quadrature Q, en suite de N fonctions de Bessel, selon les formules suivantes: N P =n anp Jl (n.n.A/D) (8') Q =nJ %nq J1 (n.n.A/D) (8") o 2D est la dynamique en entrée de l'amplificateur et anp, anq sont des coefficients calculés par une méthode de minimisation de l'erreur quadratique présentée dans l'article cité. Cette décomposition et le dispositif de simulation d'Abuelma'atti correspondant sont illustrés figure 10. Supposons maintenant que le signal d'entrée x comporte deux porteuses de fréquence f-1 et fl et d'amplitudes égales notée A, c'est-à-dire que le signal x a la forme suivante: x(t) = A.cos(27T.f_1. t) t A.cos(27.f1.t) (9) Le signal d'entrée x a alors une enveloppe X fortement modulée en amplitude, comme le montre la figure 12 et la formule suivante: X(t) = 2A cos(n.df.t) (10) o df = f1-fl représente l'écart fréquentiel des deux porteuses f-1 et fl. D'après les modèles quadratiques quasi sans mémoire, l'enveloppe Y du signal de sortie correspondant au signal d'entrée x des équations 9 et 10, prendrait la forme quadratique suivante: Yp(t) = C[2A.cos(7r.df. t)].cos($[2A.cos(7.df.t)]) (11') YQ(t) = C[2A.cos(K.df.t)].sin(b[2A. cos(ff.df.t)]) (11") o Yp serait la composante en phase de l'enveloppe Y du signal de sortie y correspondant à l'entrée x précédente, et Yq serait la composante en quadrature de l'enveloppe Y. Or les formules 11 et 11' ne peuvent pas être appliquées directement pour la simulation d'amplificateurs ANLAM non- linéaires à mémoire car l'enveloppe d'entrée X est sinusoïdale de fréquence moitié de l'écart df et varie trop rapidement dans le temps par rapport aux constantes de temps de mémorisation. D'après le modèle plus perfectionné d'Abuelma'atti, l'enveloppe Y du signal de sortie correspondant au signal d'entrée x biporteuse des équations 9 et 10 comporte une série de composantes de fréquence notée fi et d'amplitude notée Ci, i étant un entier positif ou négatif. Dans le détail, comme illustré figure 14, les composantes fréquentielles du signal de sortie y ont les fréquences et les amplitudes suivantes: une composante porteuse, de fréquence fi d'amplitude de sortie C1; - une composante porteuse, de fréquence f1 d'amplitude de sortie C_1; - des composantes d'intermodulation d'ordre 2, de fréquences valant respectivement 2.f_1, f-l+fl, 2.f1 ainsi que f_1-fl, d'amplitudes non référencées; - une composante d'intermodulation d'ordre 3, de fréquence notée f3 valant 2.f1-f_1 et d'amplitude C3, - une composante d'intermodulation d'ordre 3, de fréquence notée f-3 valant 2.f_1-f1 et d'amplitude C_3, - une composante d'intermodulation d'ordre 5, de fréquence notée f5 valant 3.f1-2.fl1 et d'amplitude C5, - une composante d'intermodulation d'ordre 5, notée f-5 valant 3.f_1-2.f1 et d'amplitude C_5, et - d'autres composantes d'ordre impair d'importance moindre et d'autres composantes d'ordre pair généralement rejetées hors de la bande utile BU. Seules les composantes d'intermodulation d'ordre 3 et éventuellement d'ordre 5, sont situées dans la bande utile BU et ont des amplitudes C_3, C3, C_5, C5 suffisantes pour être étudiés dans la réponse de l'amplificateur. On notera que l'ordre des composantes d'intermodulation est la somme des valeurs absolues des coefficients entiers de combinaison des fréquences porteuses. Le signal de sortie y correspondant est la somme des amplitudes des composantes fréquentielles selon la formule suivante: y(t) = Z 2.Ci(A).cos(24.fi.t) (12) i o i prend toutes les valeurs entières impaires négatives et positives, et Ci est une amplitude complexe qui se décompose en partie réelle Cip (en phase) et en partie imaginaire Ciq (en quadrature), selon l'écriture suivante: Ci(A) = Cip(A) + Ciq(A) (13) Abuelma'atti a montré, dans l'article précité, que le calcul des amplitudes Ci des composantes d'intermodulation issues de l'amplification d'un signal biporteuse pouvait être effectué à partir d'une sommation de fonctions de Bessel Jn du premier type mais d'ordres n divers, selon les formules suivantes: N C1p(A) =nZi1 np J0(n.T.A/D).Jl(n.-.A/D). Gnp(fl) (14') N Cjq(A) =nz1 anq- J0(n.'.A/D).Jl(n.n.A/D).Gnq(fl) (14") N C-lp(A) =nz1 enp.-Jl(n..A/D).J0(n.1.A/D).Gnp(f-l) (15') N Clq(A) =nZ1 anq.-Jl(n.n.A/D).JO(n.n.A/D).Gnq(fl) (15") N C3p(A) =nz1 anp.Jl(n.z.A/D). J2(n.n.A/D).Gnp(f3) (16') N C3q(A) =nZ1anq-Jl(n. n.A/D). J2(n.n.A/D).Gnq(f3) (16") N C_3p(A) =nZl anp.J2(n..A/D).Jl(n.n. A/D).Gnp(f-3) (17') N C_3q(A) =nz1 anq.J2(n..A/D).Jl(n.IT.A/D).Gnq(f_3) (17") o anp, anq sont les coefficients obtenus aux équations 8' et 8" en développant les caractéristiques P et Q de non-linéarité en fonctions de Bessel J1 d'ordre 1 uniquement. et J0, Jl, J2 sont des fonctions de Bessel du premier type d'ordre 0, 1, 2, respectivement. L'ordre 0, 1, 2 ou N d'une fonction de Bessel J0, J1, J2 ou JN est choisi de façon à correspondre à la valeur absolue du coefficient entier de combinaison d'une fréquence porteuse f-1 ou f1 dans la fréquence fi de la composante d'intermodulation d'amplitude Ci. L'ordre |i| de chaque composante d'intermodulation Ci correspond donc à la somme des ordres 0, 1, 2 ou N des fonctions de Bessel J0, J1, J2 ou JN impliquées dans chaque formule. Dans le cas général o le signal d'entrée x comporte m composantes fréquentielles (porteuses), d'amplitudes quelconques, notées A, B,... , M, on obtient une multitude de produits d'intermodulation ayant chacun une fréquence donnée par la formule suivante: fi = a-fA + - fB + --.. + 4-fM o a, 13,...,, sont des entiers relatifs, et fA, fB,.- -, fM sont les fréquences des porteuses A, B,..., M respectivement présentes dans le signal d'entrée x. L'amplitude Ci de chaque composante d'intermodulation de fréquence fi suit les formules générales suivantes: N N Ciq(AB,...,M)=n-1 anq.Jlel(nA/D).Jl6l1(nzB/D)-Jlgl(nM/D) o JinI est la fonction de Bessel du premier type et d'ordre valeur absolue de n (généralisation). La figure 20 illustre ainsi une réalisation sophistiquée de moyens de calcul MM de composantes d'intermodulation développée à partir de ces formules de modèle à mémoire. On a représenté schématiquement sous forme de tableau de filtres élémentaires connectés en série horizontalement pour calculer le produit des fonctions de Bessel JIal généralisées, modulées par l'amplitude des composantes A, B,..., M du signal multiporteuse. Chaque résultat horizontal est ensuite additionné avec les autres résultats horizontaux du tableau et avec les résultats d'un tableau similaire pour les composantes en quadratures avant de fournir en sortie l'amplitude d'une harmonique Ye,,..., l Sur cet exemple de réalisation sophistiqué, on voit que à partir d'une porteuse centrale B et de deux porteuses latérales A et M on obtient, outre la porteuse amplifiée Y010 de fréquence fB, au moins neuf produits d'intermodulation Y20-1, Y2-10, Yll-1, Y02-1, Yl-1l, Y- , Y-111, Y0-12 et Y-102 aux environs immédiats de la fréquence centrale fB. On voit donc que le nombre de produits d'intermodulation croît exponentiellement avec le nombre m de porteuses, que la réalisation d'une simulation pour un signal à plus de deux porteuses ne peut être détaillée dans la présente, l'homme de métier pouvant la concevoir en généralisant les calculs décrits pour deux porteuses. Un facteur de qualité déterminant dans là conception et le choix d'un système d'amplification hyperfréquence est le rapport de l'amplitude C3 de la plus forte composante d'intermodulation en sortie et de l'amplitude C1 en sortie d'une composante porteuse. Dans la présente ce facteur, appelé facteur de réjection d'intermodulation, est exprimé par la formule suivante: P/I = C1/C3 (18) Plus le ratio P/I est élevé, plus la qualité du système est élevée, la porteuse Cl se distinguant mieux du bruit d'intermodulation représenté par C3. Dans le domaine hyperfréquence, le facteur de réjection est généralement exprimé en puissance, ce qui correspond au ratio suivant: C12/C32 noté ici P/I2. Sous forme quadratique le facteur de réjection d'intermodulation s'écrit sous la forme suivante: p/I2= (C1 2 + Clq2) (19) (u3p' + u3qz) En utilisant les résultats précédents des équations 14', 14" et 16', 16" établies par Abuelma'atti, on s'aperçoit que le facteur de réjection d'intermodulation simulé pour un amplificateur hyperfréquence varie très peu selon la fréquence. En effet les correcteurs en fréquence Gn(fi) ont un faible effet de correction, la courbe des correcteurs Gn étant calculée à partir de mesures avec une seule porteuse, donc sans modulation de l'enveloppe d'entrée. En particulier, les facteurs de réjection d'intermodulation P/I ainsi calculés ne sont pas fonction de l'écart fréquentiel df des porteuses, contrairement aux effets de résonnance constatés et illustrés par exemple sur la figure 16. Selon un premier mode de réalisation de l'invention, il est prévu précisément de mesurer des caractéristiques de réjection d'intermodulation et d'utiliser ces mesures pour corriger le signal y de sortie calculé à partir d'un modèle connu. Dans le développement qui va suivre, le signal y(t) calculé selon le modèle d'Abuelma'atti est comparé à un signal z(t) comportant des composantes d'intermodulation telles que réellement mesurées. Une telle comparaison fournit la valeur d'un signal u(t) de modulation qui sera simulée par un filtre de modulation selon l'invention. La figure 21 schématise un premier mode de réalisation de dispositif selon l'invention, qui comporte un tel filtre 14 de modulation. Le filtre 14 a une fonction de transfert E(f, H) qui calcule le signal u en fonction du signal d'entrée. Le signal y du modèle connu 11 est alors corrigé par le signal modulateur u pour fournir un signal z qui reproduit les distorsions de modulation effectivement mesurées. Le développement suivant permet d'établir un exemple de relations permettant de construire la fonction de transfert E(f, H) du filtre IMF de calcul des intermodulations. On constate d'abord que le signal modulateur u est utile seulement si le signal d'entrée a une enveloppe X modulée en amplitude. En effet, lorsque le signal d'entrée x n'est pas modulé, le modèle connu 11 a un signal de sortie y qui simule convenablement la réponse de l'amplificateur. L'invention prévoit précisément d'obtenir l'enveloppe du signal en prenant le carré de la norme du signal x. En considérant l'équation 10 donnant l'enveloppe X du signal biporteuse d'amplitude A et d'écart fréquentiel df, on obtient ainsi le carré X2 de l'enveloppe suivant: X2(t) = 4.A2.cos2(n.df.t) (20) Le carré X2 donné par l'équation 20 se décompose en une amplitude A2 constante et un facteur de modulation donnés par l'équation suivante: X2(t) = 2.A2. + 2.A2.cos(2n.df.t) (20') La figure 21 montre ainsi que le filtre 14 reçoit en entrée la valeur efficace X2(t) de l'enveloppe du signal d'entrée x. Selon l'invention, on considère que le signal modulateur u issu du filtre 14 doit permettre d'écrire une formule de correction suivante: Z(t) = Y(t).[l+U(t)] (21) o Z est l'enveloppe du signal z de sortie avec les intermodulations réellement mesurées, et Y est l'enveloppe du signal y calculé par un modèle connu, et U est l'enveloppe du signal u. De la sorte, en l'absence de modulation en entrée, le signal z simulé selon l'invention a une enveloppe Z qui correspond à l'enveloppe Y du signal de sortie y simulé par le modèle connu MNLSM ou MM, et le signal modulateur u est avantageusement nul. Par hypothèse, on exprime l'enveloppe U du signal modulateur, de façon analogue à l'équation 20', comme le résultat de la transformation de la valeur efficace X2 par la fonction de transfert E(f,H) qui dépend de la fréquence f et de l'amplitude H des harmoniques h en entrée, c'est- à-dire sous la forme suivante: U(t) = 2A2.E(0,A) + 2A2.E(df,A).cos(2n.df. t) (22) o E(df,A) est la réponse de la fonction de transfert à un signal de fréquence df et d'amplitude A, et E(O,A) est la réponse de la fonction de transfert à écart fréquentiel nul. Conventionnellement, E(0,A) a une valeur nulle, puisqu'il est inutile d'avoir un signal modulateur u quand il n'y a pas d'écart de fréquence de porteuses. On note K l'inconnue suivante: K = A2.E(df,A) Il est préférable d'écrire l'enveloppe U du signal modulateur donné par l'équation 22, sous la forme complexe et exponentielle suivante: l+U(t) = K.exp(-j.2n.df.t) + 1 + K.exp(j.2n.df.t) (23) De façon équivalente à l'équation (12), on préfère ici écrire l'enveloppe Y du signal de sortie y du modèle connu sous la forme complexe et exponentielle suivante: Y(t) = S Ci(A).exp(j.i.r.df.t) (24) i o i est un entier impair prenant toutes les valeurs négatives et positives. Y(t) =...+C_5.exp(-j.5.df.t)+C_3.exp(-j.37.df.t) +C_1.exp(-j.i.df.t)+C1. exp(j.i.df.t) (24') +C3.exp(j.3v.df.t)+C5.exp(j.5f.df.t)+... On déduit donc des équations (21), (23) et (24'), le développement suivant: Z(t) =...+(K.C_5+C_3+K.C_1).exp(-j.37.df.t) +(K.C_3+C_1+K.C1). exp(-j.i.df.t) +(K.C_1+Ci+K.C3).exp(j.ff.df.t) (25) +(K.Cl+C3+K.C5).exp(j.3i.df.t)+... Dans le premier mode de réalisation on cherche à moduler uniquement l'amplitude des composantes d'intermodulation d'ordre 3 en modifiant au minimum l'amplitude des composantes porteuses (ordre 1) calculée par le modèle connu. Ceci s'obtient en minimisant la valeur de K, de sorte que les termes K.C_3, K.C1, K.C-1 et K.C3 soient négligeables. Dans ce cas, les deux lignes centrales de l'équation (25) indiquent que les composantes porteuses de fréquence f-1 et f, ont bien une amplitude sensiblement égale à C_1 et C1, respectivement, dans le signal de sortie z modulé. Dans un premier temps, on peut considérer que les amplitudes Ci de composantes d'intermodulation ne dépendent pas de la fréquence fi en supprimant les corrections en fréquence Gnp,q(fi) prévues dans les formules de calcul 12' à 15". Le calcul est alors simplifié, car les composantes symétriques en fréquence f-i et fi autour des porteuses f-1 et f1 ont des amplitudes égales, selon l'identité suivante: C-i = Ci Les deux dernières lignes de l'équation (23) représentent les composantes de fréquence fl et f3 respectivement, ce qui permet de calculer le facteur de réjection P/I qui a finalement la valeur suivante: P/I = C1+K(C1+C3) (26) Ic3+K (c1+5)| Numériquement, après avoir déterminé les amplitudes C1(A), C3(A) et C5(A) des composantes selon les méthodes connues, et après avoir mesuré la valeur du facteur de réjection P/I pour cette valeur d'amplitude A des porteuses, on résoud l'équation (26) pour déterminer la valeur de l'inconnue K. Pour la résolution, on écrit de préférence l'inconnue K sous la forme complexe suivante K = IKI. exp(j.8) On fait varier l'argument e pour calculer les valeurs de la norme IKI en fonction des valeurs de *l'argument e. La minimisation recherchée précédemment de la modification d'amplitude C_1 et C1 des porteuses f-1 et fi est obtenue en retenant l'argument 8 permettant d'obtenir la valeur minimale de la norme IKI. La figure 18 permet de schématiser la systématique de mesures effectuées selon l'invention. On voit que les caractéristiques sont déterminées en effectuant une série de mesures avec un signal multiporteuses f-1 et fi, en faisant varier à la fois l'amplitude A, A', A" des porteuses mais aussi l'écart fréquentiel df' = f1-fl des porteuses. En outre, on peut faire varier la fréquence de base f de la porteuse fl ou f-1 servant de référence à l'écart fréquentiel df. Pour chaque écart fréquentiel df, comme illustré sur la figure 16, on peut mesurer une courbe P/I caractéristique en fonction de l'amplitude A ou de la puissance de porteuse en entrée. En appliquant le calcul précédent, on obtient alors pour chaque écart fréquentiel df une courbe des valeurs d'inconnue K en fonction de l'amplitude A des porteuses. En faisant varier les écarts fréquentiels df, df', df", df"', df"" des porteuses, on obtient ainsi un tableau de courbes de valeurs d'inconnues K. Ce tableau de valeurs K établit la trame de la fonction de transfert E(df,H) du filtre 14. La continuité de la fonction de transfert E(df,H) selon l'amplitude H des harmoniques h en entrée ou selon la puissance X2 de l'enveloppe du signal d'entrée x est obtenue en interpolant les valeurs discrètes du tableau des valeurs d'inconnues K. Les calculs et interpolations sont effectuées également pour toutes les valeurs d'écart fréquentiel df, les valeurs d'inconnue K pouvant également être interpolées entre les valeurs choisies df', df", df"', df"" d'écart fréquentiel df. Un programme informatique exécute de préférence de tels calculs et stocke les valeurs d'inconnue K sous forme de matrice. Dans un second temps, on considère que les amplitudes Ci des composantes d'intermodulation dépendent de la fréquence fi selon les formules complètes des équations 14' à 17". I1 n'y a plus alors d'identité entre Ci et Ci. Les dissymétries en fréquence de la réponse de l'amplificateur peuvent alors être simulées par cette variante du premier mode de réalisation. Dans ce cas, l'invention prévoit que la fonction de transfert E'(df,H) du filtre 14 peut être dissymétrique. Les hypothèses de l'équation 22 et du calcul précédent sont alors modifiés selon la formule complexe suivante: U'(t) = A2.E'(-df,A).exp(-j. 2n.df.t) +2A2.E' (O,A)+ A2.E'(df,A).exp(j.2n.df.t) (27) o E'(O,A) a une valeur nulle car le filtre ne transfère pas de composante continue de signal, et E' (- df,A) est une réponse de la fonction de transfert E' à un écart fréquentiel -df, réponse distincte de la réponse E' (df,A) à un écart fréquentiel dl, en raison de la dissymétrie du filtre E'. Notons les termes de l'équation 22', sous la forme simplifiée suivante: K+ = A2.E'(df,A) et K = A2.E'(-df,A) L'équation 23 est alors modifiée selon la formule suivante: l+U'(t) = K.exp(-j.2K.df. t)+1 + K+.exp(j.2f.df.t)(28) Notons Ci' les composantes modulées du signal de sortie z' modulé obtenu selon cette variante du premier mode de réalisation. Le signal de sortie z' est encore obtenu en modulant le signal y d'un modèle connu par un signal modulateur u' calculé par la fonction de transfert E' du filtre 34 de la figure 23. Ceci permet de distinguer les composantes C'i modulées selon l'invention, des composantes Ci calculées par le modèle sans mémoire MNLSM connu. Selon les équations 21, 24' et 28 précédentes, les composantes modulées ont des amplitudes C_3', C_1', Cl', C3' déterminées par les formules suivantes: C1' = K.C_1 + C1 + K.C3 (29) C3' = K+.C1 + C3 + K.C5 (30) C-1' = K.C_3 + C_1 + K.C1 (31) C_31 = K+.C_5 + C_3 + K.C-1 (32) On obtient alors deux résultats plus complexes de facteurs de réjection d'intermodulation, notés P/I+ et P/I-, selon que le facteur de réjection est mesuré sur les fréquences fi et f3 supérieures ou sur les fréquence f!1 et f-3 inférieures de porteuse et d'intermodulation. Les deux facteurs de réjection ont les formes suivantes: P/I+ = ICî'/C3'| = K.C-1 + C1 + C + K C3 (33) K+.C1 + C3 + K C5 P/I- = C_1'/C_3'l = IK+.C_3 + C-1 + K-.Clll (34) K+.C_5 + C_3 + K.C-l On résoud ensuite, de façon connue, le système des deux équations 33 et 34 à deux inconnues K et K+. Ainsi quand l'amplificateur a une réponse dissymétrique autour des fréquences porteuses, il est prévu selon l'invention de mesurer des caractéristiques de réjection d'intermodulation, en prévoyant de mesurer les deux facteurs de réjection P/I+ supérieur et P/I inférieur pour chaque amplitude A et écart fréquentiel df de porteuses, afin de simuler les distorsions de modulation en sortie. Un second mode de réalisation de l'invention permet de simuler le second effet de distorsion d'intermodulation. Le second effet observé est la variation du gain des porteuses en fonction de la fréquence d'enveloppe df/2, illustré par exemple sur la figure 17. En mesurant le gain de l'amplificateur pour une porteuse fl d'amplitude C1 de sortie, en présence d'une autre porteuse f-l, on se rend compte que le gain peut avoir des variations importantes, en fonction de l'écart fréquentiel df des deux porteuses f-, et fl. L'invention prévoit d'ajuster les amplitudes Cl et C_1 de sortie des porteuses calculée par un modèle connu en tenant compte de mesures de caractéristiques Cl/C-1 d'interaction de porteuses. De préférence, l'ajustement du second mode de réalisation se combine à l'ajustement des principales composantes d'intermodulation C3 et C-3 prévu dans le premier mode de réalisation de l'invention, précédemment exposé. On peut prévoir en outre un ajustement des composantes C5,C-5 d'intermodulation secondaires. Dans le second mode de réalisation de dispositif selon l'invention, l'ajustement est effectué en corrigeant le signal de sortie y, ou de préférence le signal modulé z ou z', par un signal de correction v issu d'un filtre 37 ayant une fonction de transfert F de correction. La fonction de transfert F est calculée, de façon analogue à la fonction E, en supposant une correction des composantes modulées Ci' formulées précédemment aux équations 29 à 32 selon les formules suivantes: C3 = K.C1 + C3 + K.C5 (35) C_3 i= K.C_5 + C_3 + K.C_1 (36) Ci = K.C_1 + C1 + K-.C3 +L+ (37) C_1t = K+.C_3 + C_1 + K.C1 + L- (38) dans lesquelles L+ = A.F(df,A) L = A.F(df,A) F étant la fonction de transfert de correction exprimée en fonction de la fréquence de modulation de l'enveloppe, c'est-à-dire de l'écart fréquentiel df des porteuses d'amplitude A en entrée. Les composantes C3", C_3", C1i " et C_1I" sont appelées composantes corrigées. Les équations 35 à 38 forment un système de quatre équations à quatre inconnues K+, K, L+, L qui est résolu de façon connue. La résolution des deux dernières équations 37 et 38 fournit les valeurs des inconnues de correction L+ et L -. La mesure précise des amplitudes réelles des composantes f3, f- 3 d'intermodulation d'ordre 3 et des porteuses fl, f-1 permet ainsi, après calcul des composantes C_5, C_3, C_1, C1, C3, C5 d'ordre 1, 3 et 5 selon un modèle connu, de calculer les fonctions de transfert de modulation E(df,A) et de correction F(df,A) en fréquence. De telles mesures complètes d'amplitudes de composantes en sortie d'un amplificateur comportant à la fois des mesures de réjection de composante d'intermodulation d'ordre 3 et 5 en sortie et des mesures d'interaction de porteuses permettent ainsi de simuler complètement la réponse en signal avec les distorsions dues à l'effet de mémoire d'enveloppe d'un amplificateur non linéaire à mémoire. On notera cependant que des mesures partielles comme les mesures de réfection Cl/C3 d'intermodulation d'ordre 3 ou les mesures d'interaction de porteuses C1/C-1 suffisent à caractériser certains effets de distorsion d'intermodulation. Un deuxième mode de réalisation permet la validation du modèle par la mesure du rapport signal à bruit en présence d'un nombre élevé de porteuses dont la densité spectrale de puissance DSP est modélisée par un bruit blanc. Les figures 19A et 19B illustrent ainsi un diagramme fréquentiel de bruit, respectivement appliqué en entrée et mesuré en sortie d'un amplificateur non-linéaire à mémoire pour le caractériser. Selon l'invention, comme illustré sur la figure 19A, on applique un bruit blanc ayant une densité spectrale de puissance DSP de niveau NP constant dans une bande de fréquence BW sauf dans une fenêtre BH de fréquence dans laquelle le bruit est absent. Ce signal d'entrée bien qu'ayant une densité spectrale constante est d'amplitude variable. Il permet donc de mesurer les distorsions de modulation apparaissant en sortie de l'amplificateur. Les distorsions apparaissent en sortie de l'amplificateur sous forme de report de bruit dans la fenêtre BH précédemment définie. La caractérisation de l'amplificateur consiste à mesurer sa capacité à limiter le report de bruit dans la fenêtre BH en mesurant la différence de niveau de bruit NPR en sortie, dans la bande de fréquence BW et dans la fenêtre BH. Ces mesures de caractéristiques sont effectuées à différent niveau NP de bruit et sont introduites dans le simulateur afin qu'il reproduise ses caractéristiques de bruit selon les conditions de fonctionnement. Un troisième mode de contrôle de la validité du modèle consiste à mesurer le taux d'erreur binaire lors d'une amplification de porteuses modulées. La figure 24 illustre un banc de test utilisé pour établir selon l'invention des mesures de caractéristiques d'amplificateur DUT nonlinéaire à mémoire. Ce banc de test dont le détail de fonctionnement est à la portée de l'homme de l'art comporte notamment un bloc 10 de génération de signal modulé en amplitude avec par exemple deux générateurs SG1 et SG2 de fréquences pouvant être distinctes pour appliquer un signal multiporteuse au dispositif DUT amplificateur à caractériser. Il comporte en outre un bloc 20 de mesure de bruit avec un générateur NG de bruit blanc à fenêtre et éventuellement un bloc 30 de mesure de taux d'erreur binaire avec un générateur MOD de signaux modulés par un message binaire aléatoire et un comparateur Comp du signal d'origine et du signal de retour de l'amplificateur. La figure 25 illustre enfin un diagramme des fichiers de caractéristiques mis en oeuvre par un programme informatique pour simuler l'amplification non- linéaire selon le procédé de l'invention. Le programme Pro fait ainsi appel à un fichier 203 stockant des résultats de mesures de caractéristiques AP statiques, c'est-à-dire de courbes C,0 ou de courbes P,Q établies à amplitude constante (de façon connue). Le programme Pro fait encore appel à une série de fichiers 205, 206, 207 stockant des résultats de mesures de caractéristiques P/I dynamiques comme des courbes Cl/C3 de réjection d'intermodulation ou des courbes C1/C- 1 d'interaction de porteuses. Chaque fichier 205 ou 206 ou 207 de la série est établi par exemple pour un écart fréquentiel df de porteuses. Le programme est lancé après avoir introduit des données générales Dat sur l'amplificateur à simuler. En sortie, le programme Pro fournit un fichier 204 des coefficients a de suite de fonctions de transfert directes, notamment de fonctions de Bessel. Il fournit surtout selon l'invention des fichiers 208, 209 contenant les données des fonctions de transfert de modulation E et de correction F précitées. La figure 26 représente un exemple de schéma d'étape de processus mis en oeuvre par un tel programme Pro. La réalisation de dispositif selon l'invention est effectuée de préférence sous forme de programme informatique dont les fonctions correspondent aux blocs schématisés à titre d'exemple sur les figures 20 à 24. D'autres modes de réalisation, caractéristiques et avantages, peuvent être développés par l'homme de l'art sans sortir du cadre de l'invention, défini notamment par les revendications ci-après. REVENDICATIONS 1. Procédé de simulation de réponse en signal d'un amplificateur nonlinéaire présentant un effet de mémoire, le procédé comportant des étapes consistant à: - mesurer des caractéristiques (C, à, P, Q) de conversion amplitude/amplitude (C[A]) et de conversion amplitude/déphasage (4[A]) de l'amplificateur, chaque mesure étant effectué à amplitude (A) constante en entrée, et - développer les caractéristiques (P, Q) en suite de fonctions (J1) de transfert directes, le procédé étant caractérisé par des étapes consistant à: - mesurer des caractéristiques de distorsion de modulation d'amplitude (Y), chaque mesure étant effectuée en modulant l'amplitude (X) d'entrée, - calculer des fonctions (H) de transfert de modulation reproduisant en sortie (a) les amplitudes de distorsion selon les amplitudes de modulation d'entrée (h) et corrigeant les fonctions de transfert directes lorsque l'amplitude d'entrée (x) est modulée, afin de simuler l'effet de mémoire d'enveloppe. 2. Procédé selon la revendication 1, caractérisé par des étapes intercalaires consistant à: - mesurer les caractéristiques à différentes fréquences (f0-, fo, fo+), et - calculer des correcteurs (G[f]) en fréquence pour les fonctions de transfert directes. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce qu'une mesure de caractéristiques de distorsion de modulation d'amplitude consiste à mesurer des distorsions d'intermodulation en appliquant en entrée un signal comportant au moins deux porteuses de fréquences (fl, f-1) différentes. 4. Procédé selon l'une des revendications 1 à 3, caractérisé par une étape consistant à: - mesurer des caractéristiques de réjection (C1/C3) de composantes d'intermodulation en sortie, chaque mesure étant effectuée avec un groupe de porteuses (fl, f-1) en entrée d'amplitudes (A) déterminées. 5. Procédé selon l'une des revendications 1 à 4, caractérisé par une étape consistant à: - mesurer des caractéristiques d'interaction (Cl, C- 1) de porteuses en sortie, chaque mesure étant effectuée avec un groupe de porteuses (fl, f-1) en entrée d'amplitudes (A) déterminées. 6. Procédé selon l'une des revendications 1 à 5, caractérisé par une étape consistant à: - mesurer des caractéristiques de bruit (DSP) de modulation en sortie, chaque mesure étant effectuée avec des amplitudes (A) déterminées de modulation en entrée. 7. Procédé selon l'une des revendications 1 à 6, caractérisé par une étape consistant à: - mesurer les caractéristiques de distorsion avec plusieurs écarts fréquentiels (df', df", df"') de modulation en entrée. 8. Procédé selon l'une des revendications 1 à 7 caractérisé par une étape consistant à: - mesurer les caractéristiques de distorsion avec plusieurs valeurs d'amplitudes (A', A", A"') de modulation en entrée. 9. Procédé selon l'une des revendications 1 à 8, caractérisé par une étape consistant à: - mesurer des caractéristiques de report (NPR) de bruit dans une fenêtre (BH) de fréquence, chaque mesure étant effectuée avec un bruit blanc (DSP) à fenêtre (BH) en entrée. 10. Procédé selon l'une des revendications 1 à 9, caractérisé par une étape consistant à: - décomposer l'amplitude d'entrée en série de Fourier (FFT) de composantes fréquentielles. 11. Procédé selon l'une des revendications 1 à 10, caractérisé en ce qu'il est appliqué à la simulation de réponse en signal d'un amplificateur hyperfréquence fonctionnant à haut rendement. 12. Dispositif de simulation de réponse en signal d'un amplificateur nonlinéaire présentant un effet de mémoire, le dispositif comportant: - un module (MNLSM, MM) de calcul de la réponse (y) d'un amplificateur nonlinéaire à un signal monoporteuse, le module comportant: - au moins un filtre direct (P0, Q0, J1) ayant une fonction de transfert établie à partir de caractéristiques (P, Q) de conversion amplitude/amplitude (C[A]) et de conversion amplitude/déphasage (ô[A]), mesurées réellement sur l'amplificateur, chaque fonction de transfert directe étant établie avec un signal monoporteuse (fo) prenant différentes amplitudes (A',A"), caractérisé en ce que le dispositif comporte: - des moyens de génération d'harmoniques (Harm) en réponse à un signal multiporteuse (f-1, fl), - au moins un filtre (IMF, CF, J2, JO) de modulation de la réponse (y, Y010) du module (MNLSM, MM) de calcul, le filtre de modulation ayant une fonction de transfert (E[f,H], F[f,H]) établie à partir de caractéristiques de distorsions de modulation d'amplitude mesurées réellement sur l'amplificateur, chaque fonction de transfert de modulation étant établie avec des signaux multiporteuses (f-l,fl). 13. Dispositif selon la revendication 12 caractérisé en ce que le module comporte en outre: - au moins un correcteur (Gp[f]) en fréquence pondérant la fonction de transfert d'un filtre direct (Po[A], Qo[A]), chaque correcteur ayant une fonction de correction établie à partir de caractéristiques de conversion mesurées sur l'amplificateur à différentes fréquences (fo-, f0, fo0+) de signal monoporteuse. 14. Dispositif selon la revendication 12 ou 13, caractérisé en ce que les moyens de génération d'harmoniques comportent un module de décomposition de Fourier (FFT) d'un signal en composantes (A, B, M) fréquentielles. 15. Dispositif selon l'une des revendications 12 à 14, caractérisé en ce que les moyens de génération d'harmoniques (Harm) comportent un module combinatoire calculant les combinaisons de fréquence (2f_1-fl,2fl-fl) de composantes d'intermodulation (f-3, f3). 16. Dispositif selon l'une des revendications 12 à , caractérisé en ce que les moyens de génération d'harmoniques comportent un module (Squa) normatif calculant une valeur efficace d'amplitude de signal. 17. Dispositif selon l'une des revendications 12 à 16, caractérisé en ce que le filtre (IMF) de modulation a une fonction de transfert (E[f,H]) établie à partir de caractéristiques de réjection (C1/C3) de composantes d'intermodulation, de sorte que des composantes d'intermodulation (f-3, f3) d'amplitude (C_3, C3) calculée par le filtre de modulation (IMF) sont ajoutées à la réponse (y) du module de calcul lorsque l'amplification d'un signal multiporteuse (f-l, fl) est simulée par le dispositif. 18. Dispositif selon l'une des revendications 12 à 17, caractérisé en ce qu'il comporte en outre un filtre de correction (CF) de la réponse du module de calcul (MM), le filtre ayant une fonction de transfert (F) établie à partir de caractéristiques d'interaction (Cl, C-i) de porteuses, de sorte que l'amplitude (C-1) d'une porteuse (f-l) fournie par le module de calcul (MM) est corrigée par le filtre de correction (CF) lorsque l'amplification d'un signal multiporteuse (f-1, fl) est simulée par le dispositif. 19. Dispositif selon l'une des revendications 12 à 18, caractérisé en ce qu'il est appliqué à la simulation de réponse en signal d'un amplificateur hyperfréquence fonctionnant à haut rendement.