FR 3025320 A1 20160304 FR 1457980 A 20140826 Procédé de détermination de paramètres linéiques d'une ligne de transmission La présente invention concerne le domaine des systèmes et procédés 5 de diagnostic pour câble et en particulier le domaine des méthodes de réflectométrie permettant de détecter et/ou localiser les défauts impactant un câble. L'invention porte plus précisément sur un procédé permettant de déterminer les paramètres linéiques d'une ligne de transmission tels que la 10 résistance linéique, l'inductance linéique, la capacitance linéique et la conductance linéique. De tels paramètres présentent une importance pour caractériser un câble et peuvent être utiles dans le cadre du diagnostic filaire basé sur une 15 méthode de réflectométrie. En particulier les paramètres linéiques déterminés grâce à la présente invention peuvent être utilisés comme données d'entrée du procédé de réflectométrie décrit dans la demande de brevet française du Demandeur déposée sous le numéro FR 1355377. 20 Mais ils peuvent également être utilisés pour d'autres applications qui nécessitent la connaissance des ces paramètres, par exemple ils peuvent être employés pour simuler un réseau de câbles complexe pour permettre le diagnostic de défauts à partir de cet environnement simulé. L'invention s'applique à tout type de câble électrique, en particulier 25 des câbles de transmission d'énergie ou des câbles de communication, dans des installations fixes ou mobiles. Les câbles concernés peuvent être coaxiaux, bifilaires, en lignes parallèles, en paires torsadées, en toron de câble ou autre. L'invention peut également s'appliquer à des câbles mécaniques, par exemple des câbles de soutien d'infrastructures telles un 30 ascenseur ou un pont. 3025320 2 Les paramètres linéiques d'une ligne de transmission ne sont pas toujours connus ni fournis dans une documentation du fabricant de câble. En outre, les données fournies par les constructeurs de câble ne permettent pas toujours de déterminer avec une précision suffisante ces paramètres. On connait par ailleurs les méthodes décrites dans les publications [1],[2] et [3]. Ces méthodes présentent les limitations suivantes. Elles nécessitent le plus souvent au moins deux points de mesure sur le câble à tester, par exemple un point de mesure en entrée et un point de mesure en sortie ou deux points de mesure en entrée en modifiant, entre les deux mesures, la charge en entrée. En outre ces méthodes présentent l'inconvénient d'être sensibles aux défauts électriques qui viennent perturber les mesures des paramètres linéiques obtenus. A partir d'une mesure en réflectométrie temporelle ou fréquentielle, la présente invention permet d'extraire, par un ensemble de calcul, les paramètres linéiques intrinsèques et génériques au câble et représentatifs de la propagation d'un signal électromagnétique au sein dudit câble, sur une plage de fréquence recherchée qui peut être plus large que la bande de fréquence du signal sonde. Les paramètres linéiques d'un câble sont la résistance linéique R, l'inductance linéique L, la capacitance linéique C et la conductance linéique G. Dans la suite de la description ces paramètres seront désignés par l'appellation paramètres RLCG. L'invention a pour objet un procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission, ledit procédé comprenant les étapes suivantes : - Déterminer au moins une mesure du facteur de propagation complexe y en fonction de la fréquence à partir d'au moins une mesure réalisée sur la ligne de transmission, 3025320 3 - Déterminer au moins une mesure de l'atténuation linéique a de la ligne de transmission égale à la partie réelle de la mesure du facteur de propagation complexe y et/ou au moins une mesure du facteur de phase 13 de la ligne de transmission égale à la partie 5 imaginaire de la mesure du facteur de propagation complexe y, - Filtrer la mesure de l'atténuation linéique a et/ou la mesure du facteur de phase 13 à partir d'un modèle de régression polynomial en fréquence dépendant des caractéristiques physiques de la ligne de transmission. 10 Selon un aspect particulier de l'invention, le modèle de régression pour l'atténuation linéique est du type a(f) = aoj + a1 f , où ao et al sont deux nombres réels. Selon un aspect particulier de l'invention, le modèle de régression pour le facteur de phase est du type /3(f) = bof où bo est un nombre réel. 15 Selon un aspect particulier, le procédé selon l'invention comporte en outre la détermination d'une estimée d'au moins un paramètre linéique de la ligne de transmission parmi la résistance linéique, l'inductance linéique, la capacitance linéique ou la conductance linéique, à partir d'au moins une mesure parmi la mesure filtrée de l'atténuation linéique, la mesure filtrée du 20 facteur de phase et la mesure du coefficient de réflexion pe à une première extrémité de la ligne de transmission, la dite estimée étant déterminée comme une fonction polynomiale de la fréquence. Selon un aspect particulier, le procédé selon l'invention comporte en outre les étapes suivantes : 25 - Déterminer une nouvelle estimée Xi d'un paramètre caractéristique de la propagation dans la ligne de transmission à partir de ladite estimée d'au moins un paramètre linéique de la ligne de transmission, - Mesurer le degré de similitude entre la nouvelle estimée Xi; et une 30 mesure X du même paramètre caractéristique. 3025320 4 Selon un aspect particulier de l'invention, le paramètre caractéristique de la propagation dans la ligne de transmission est le facteur de propagation complexe y ou une fonction de ce paramètre. Selon un aspect particulier de l'invention, la mesure du facteur de 5 propagation complexe y est déterminée à l'aide des étapes suivantes : - Déterminer une mesure H de la fonction de transfert fréquentielle de la ligne de transmission, - Déterminer une mesure du coefficient de réflexion pe à une première extrémité de la ligne de transmission, 10 - Déterminer une mesure du facteur de propagation complexe y à partir de la mesure H de la fonction de transfert fréquentielle, de la mesure du coefficient de réflexion pe à une première extrémité de la ligne de transmission et d'une estimée de la longueur de la ligne de transmission, 15 Selon un aspect particulier de l'invention, la mesure H de la fonction de transfert fréquentielle de la ligne de transmission est déterminée à l'aide d'une méthode de réflectométrie fréquentielle à partir de l'injection, à une première extrémité de la ligne, d'un signal de réflectométrie. Selon un aspect particulier de l'invention, la mesure du coefficient de 20 réflexion pe en entrée de la ligne de transmission est déterminée par une moyenne, sur une bande fréquentielle donnée, de la mesure H de la fonction de transfert fréquentielle de la ligne de transmission. Selon un aspect particulier de l'invention, la mesure H de la fonction de transfert fréquentielle de la ligne de transmission est déterminée à l'aide 25 d'une méthode de réflectométrie temporelle à partir de l'injection, à une première extrémité de la ligne, d'un signal de réflectométrie. Selon un aspect particulier de l'invention, la mesure du coefficient de réflexion pe en entrée de la ligne de transmission est déterminée par la mesure de l'amplitude du pic du signal de réflectométrie réfléchi au point 30 d'injection. 3025320 5 Selon un aspect particulier, le procédé selon l'invention comporte en outre les étapes suivantes : - Déterminer une mesure de l'impédance caractéristique Zc de la ligne de transmission ou une mesure du coefficient de réflexion à 5 une seconde extrémité de la ligne de transmission à partir de la mesure du coefficient de réflexion pe à une première extrémité de la ligne de transmission. - Déterminer une mesure du facteur de propagation complexe y à partir de la mesure H de la fonction de transfert fréquentielle et de 10 la mesure de l'impédance caractéristique Zc de la ligne de transmission ou de la mesure du coefficient de réflexion à une seconde extrémité de la ligne de transmission. Selon un aspect particulier de l'invention, le paramètre caractéristique de la propagation dans la ligne de transmission est la fonction de transfert 15 fréquentielle H ou l'impédance caractéristique Zc ou une fonction de l'un de ses paramètres. Selon un aspect particulier de l'invention, le degré de similitude est mesuré par le calcul d'un coefficient de détermination entre la nouvelle estimée d'un paramètre caractéristique Xi de la propagation dans la ligne de 20 transmission et la mesure X du même paramètre caractéristique. Selon un aspect particulier, le procédé selon l'invention comporte en outre les étapes suivantes : - Déterminer une nouvelle estimée y; du facteur de propagation complexe de la ligne de transmission à partir des estimées de la 25 résistance linéique, de l'inductance linéique, de la capacitance linéique et de la conductance linéique de la ligne de transmission, - exécuter les étapes du procédé permettant de déterminer une estimée d'au moins un paramètre linéique de la ligne de transmission à partir de la nouvelle estimée y; du facteur de 30 propagation complexe. 3025320 6 Selon un aspect particulier de l'invention, les étapes du procédé sont itérées un nombre donné de fois et les paramètres linéiques retenus sont ceux déterminés à l'itération correspondant au degré de similitude le plus élevé. 5 Selon un aspect particulier de l'invention, les étapes du procédé sont itérées tant que le degré de similitude n'a pas atteint une valeur donnée. Selon un aspect particulier de l'invention, la ligne de transmission est un câble coaxial ou un câble bifilaire ou un câble torsadé. Selon un aspect particulier, le procédé selon l'invention comporte en 10 outre une étape de décimation en fréquence appliquée à une pluralité de mesures de l'atténuation linéique a et/ou de mesures du facteur de phase 13 avant l'étape de filtrage. L'invention a aussi pour objet un programme d'ordinateur comportant des instructions pour l'exécution du procédé de détermination d'au moins un 15 paramètre linéique d'une ligne de transmission selon l'invention, lorsque le programme est exécuté par un processeur et un dispositif comprenant des premiers moyens de mesure connectés à une ligne de transmission et des seconds moyens de calcul configurés pour exécuter le procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission 20 selon l'invention. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent : 25 - La figure 1, un schéma de principe d'un câble à tester dont on souhaite déterminer les paramètres linéiques par application du procédé selon l'invention, - La figure 2, un organigramme décrivant les étapes du procédé selon l'invention, 30 3025320 7 Les paramètres linéiques d'une ligne de transmission R, L, C et G apparaissent dans les équations bien connues dites des télégraphistes qui permettent de décrire l'évolution de la tension et du courant sur une ligne électrique en fonction de la distance et du temps. Une portion de ligne 5 électrique peut être représentée par un quadripôle comprenant une résistance linéique R, une inductance linéique L, une capacitance linéique C et une conductance linéique G. L'usage de ce modèle à constantes réparties R, L, C et G possède de nombreux avantages pour modéliser la propagation d'une onde sur un câble. 10 Ce modèle offre en effet une description fréquentielle indépendante de la longueur lié à quatre paramètres définis par la géométrie du câble ainsi que sa composition matérielle. La résistance R et l'inductance L sont directement liées aux conducteurs métalliques du câble et à leurs dimensions tandis que la capacitance C et la conductance G décrivent la nature du diélectrique 15 utilisé pour séparer ces conducteurs. Les paramètres R,L,C,G se calculent à partir des équations de Maxwell. Cette modélisation permet de caractériser un câble en déduisant par exemple la vitesse de propagation ainsi que l'atténuation ou encore 20 l'impédance caractéristique du câble qui régissent la propagation de n'importe quel signal dans ce câble. Cependant, s'il existe un lien direct entre ces trois grandeurs précitées et les quatre paramètres RLCG, il n'est pas forcément évident de pouvoir les déduire d'une mesure issue de méthodes de réflectométrie. La faute en 25 revient à l'estimation peu précise des trois grandeurs que sont la vitesse de phase, l'atténuation et principalement l'impédance caractéristique, dont la partie imaginaire est particulièrement compliquée à estimer correctement. L'erreur générée par les méthodes de calcul usuelles basées sur l'isolation du pic de désadaptation (en réflectométrie temporelle) et l'étude de 30 l'amplitude de sa transformée de Fourier, est suffisante pour provoquer un écart important dans le calcul des paramètres RLCG, particulièrement pour 3025320 8 la conductance G qui peut accuser un écart de plus de 100% par rapport à la théorie. De plus, dans le domaine de la réflectométrie fréquentielle, les méthodes existantes sont très contraignantes car elles nécessitent d'avoir un 5 accès à chaque extrémité du câble. Elles souffrent d'ailleurs de la même imprécision que les méthodes de réflectométrie temporelle. L'invention est basée sur le principe que le câble étudié peut être décrit par l'ensemble des quatre paramètres RLCG, ce qui implique qu'il est 10 quasi homogène sur toute sa longueur. De plus, ces paramètres RLCG respectent certaines propriétés mathématiques comme la continuité en fréquence, et des propriétés physiques comme le fait qu'ils sont forcément positifs. L'invention présente de bons résultats également dans le cas où le câble est impacté par des défauts non francs qui induisent une variation 15 locale d'impédance de l'ordre de quelques % par rapport à l'impédance caractéristique du câble et qui sont de faibles dimensions. On définit les paramètres suivants utilisés dans la suite de la description et en référence à la figure 1 qui donne un schéma de principe 20 d'un câble sous test CA relié à son entrée E à un appareil de mesure AM adapté à effectuer une mesure de réflectométrie en vue d'obtenir une mesure de la fonction de transfert du câble. L'extrémité opposée S, de sortie, du câble est reliée à une charge en bout de ligne CL. 25 Z, est l'impédance caractéristique du câble CA, y est le facteur de propagation complexe du câble CA, a est l'atténuation linéique du câble CA qui est égale à la partie réelle du facteur de propagation complexe, le est le facteur de phase du câble CA qui est égal à la partie imaginaire du 30 facteur de propagation complexe, 3025320 9 ZG est l'impédance en entrée E du câble CA, c'est-à-dire en sortie de l'appareil de mesure AM connecté en entrée E du câble pour réaliser une mesure de réflectométrie, ZF est l'impédance de la charge CL en bout de ligne, c'est à dire à l'extrémité 5 S opposée du câble, zc -zG Pe = - est le coefficient de réflexion dû à la désadaptation entre le câble zc-EzG CA et l'appareil de mesure AM, autrement dit le coefficient de réflexion en entrée E du câble, zf -zc pf = est le coefficient de réflexion dû à la charge CL en bout de zf+zc, 10 ligne, c'est-à-dire le coefficient de réflexion en sortie S du câble, I est la longueur du câble et H est la fonction de transfert fréquentielle du câble. Les équations des télégraphistes, issues de la modélisation d'une 15 ligne de transmission par constante répartie, permettent de relier les paramètres RLCG d'une ligne à son impédance caractéristique Zc et au facteur de propagation y dans le sens dit direct. f est la fréquence du signal se propageant le long de la ligne. 20 Zc = -\/(R + i2n-fL)/(G + i2n-fC) (i) y = \ I (R + i2n-fL) * (G + i2n-fC) (ii) Les équations (i) et (ii) montrent qu'il est possible de déterminer les paramètres RLCG à partir de la connaissance de l'impédance caractéristique 25 et du facteur de propagation : {R =9U{YZc} G =9U{Y/Zc} L = Z'm{yZc}/2n-f C = Z'm{y/Zc}/2n-f 3025320 10 Re() désigne la fonction partie réelle et lm() désigne la fonction partie imaginaire. La figure 2 décrit, sur un organigramme, les étapes de mise en oeuvre 5 du procédé selon l'invention. Selon une première étape 101, on effectue une mesure de la fonction de transfert fréquentielle H du câble à analyser. Cette étape est réalisée à l'aide d'un appareil de mesure AM adapté à effectuer une mesure de réflectométrie. Cette mesure consiste en l'injection d'un signal de 10 réflectométrie à l'entrée E du câble connectée à l'appareil de mesure AM puis une mesure du signal réfléchi. De nombreuses méthodes de réflectométrie existent dans l'état de l'art, elles ne sont pas décrites ici car elles ne font pas l'objet de l'invention. On peut cependant distinguer d'une part les méthodes de réflectométrie temporelles qui analysent le signal 15 réfléchi dans le domaine temporel et les méthodes de réflectométrie fréquentielles qui analysent le signal réfléchi dans le domaine fréquentiel. Si une méthode fréquentielle est utilisée, l'appareil de mesure AM délivre directement un réflectogramme fréquentiel qui est égal à la fonction de transfert en fréquence du câble sur toute la gamme de fréquence 20 analysée. Si une méthode temporelle est utilisée, l'appareil de mesure AM délivre un réflectogramme temporel auquel il faut appliquer une transformée en fréquence, par exemple une transformée de Fourier, puis diviser le résultat par la transformée de Fourier du signal de réflectométrie injecté dans 25 le câble. Sans sortir du cadre de l'invention, toute autre méthode de réflectométrie ou toute méthode d'analyse permettant de déterminer, à partir de mesures faites sur un câble, la fonction de transfert en fréquence de ce câble entrent dans la portée de l'invention. 30 3025320 11 Dans une deuxième étape 102 du procédé selon l'invention, on détermine ensuite une mesure du coefficient de réflexion pe dû à la désadaptation entre le câble et l'appareil de mesure AM, autrement dit le coefficient de réflexion en entrée E du câble. 5 Plus précisément, on mesure la partie réelle du coefficient de réflexion Pe. Cette mesure peut être réalisée à partir du réflectogramme mesuré par l'appareil de mesure AM. Si une méthode de réflectométrie temporelle est utilisée, la partie réelle du coefficient de réflexion pe est obtenue en mesurant, au point d'injection E du signal de réflectométrie, l'amplitude du pic 10 de désadaptation dans le réflectogramme correspondant au signal réfléchi. Si une méthode de réflectométrie fréquentielle est utilisée, alors cette mesure est obtenue en effectuant la moyenne de la partie réelle de la fonction de transfert fréquentielle H sur la bande de fréquence considérée. 15 Il est connu que la fonction de transfert fréquentielle est liée aux coefficients de réflexion en entrée et en sortie du câble ainsi qu'au facteur de propagation par la relation suivante, où I est la longueur du câble : (1-ePfe-2Y1 (iv) 1-Me) = Pe + 1+pepfe-2" 20 Selon une troisième étape 201, il est possible de déterminer, à partir de la mesure de la partie réelle du coefficient de réflexion pe , une mesure de la partie réelle de l'impédance caractéristique Zc et/ou une mesure de la partie réelle du coefficient de réflexion pf en bout de ligne S. Ces mesures 25 peuvent être obtenues à l'aide des relations introduites ci-dessus reliant ces différents paramètres entre eux et à partir de la connaissance de l'impédance Z0 en sortie de l'appareil de mesure AM et de l'impédance ZF de la charge CL en bout de ligne. Cette troisième étape 201 est optionnelle dans la mesure où le calcul 30 d'une estimée du facteur de propagation peut être fait à partir de la seule 3025320 12 connaissance de l'une des trois mesures pe, Zc ou pf comme cela se déduit de la relation (vi) et sachant que ces trois grandeurs sont reliées entre elles par des relations linéaires. 5 A partir de la relation (iv) on peut déterminer la grandeur H1 en retranchant à la mesure H de la fonction de transfert, la mesure de la partie réelle du coefficient de réflexion pe. (1- ep f e2Y1 H1 = H - pe = r + où r est le reste provenant de l'estimation 1+ p ep f e-2YI incomplète de pe. (v) 10 A partir des relations (iv) et (v) on peut déterminer le facteur de propagation grâce à la relation (vi) : Y = In (vi) 1 (1-epf-pepfH1 21 H1 15 A partir de la relation (vi), on voit qu'il est possible d'estimer 103 une mesure du facteur de propagation à partir de H1, Pe, et pf et de la connaissance de la longueur du câble I. Alternativement, pf peut être remplacé par son expression en fonction de pe. Il est aussi possible de faire intervenir 20 l'impédance caractéristique Zc dans l'expression (vi). Dans une étape supplémentaire 104, on détermine ensuite les valeurs de l'atténuation linéique a et du facteur de phase 13 qui sont respectivement égales à la partie réelle et à la partie imaginaire du facteur de propagation y. 25 On précise que toutes les grandeurs calculées à partir de la fonction de transfert H dépendent de la fréquence et sont ainsi mesurées pour plusieurs valeurs de fréquences dans une bande fréquentielle de mesure prédéfinie. 3025320 13 Sans sortir du cadre de l'invention, d'autres méthodes sont envisageables pour déterminer, à partir d'au moins une mesure réalisée sur un câble à tester, une estimée du facteur de propagation y. Par exemple, les étapes précitées 101,102,201,103 peuvent être remplacées par la méthode 5 décrite dans la référence [1] qui se base sur la mesure d'un réflectogramme temporel. La méthode décrite dans la référence [6] peut également être utilisée. Du fait des imperfections liées à une mesure non parfaite du 10 coefficient de réflexion pe ,notamment du fait qu'on ne mesure que sa partie réelle, les estimées des paramètres a et [3 sont entachées d'erreurs qui apparaissent sous la forme de variations ou d'oscillations périodiques autour de la valeur réelle. Pour lutter contre ces variations, il est nécessaire d'opérer un filtrage 15 105 des paramètres a et [3. Pour réaliser cette étape de filtrage 105, plusieurs solutions sont possibles. Une solution consiste à effectuer un filtrage par régression linéaire ou 20 non linéaire. Pour cela, on détermine un modèle de variation, en fonction de la fréquence, des paramètres a et [3. Ce modèle dépend du type de câble visé et en particulier de ses caractéristiques physiques. Pour un câble du type câble coaxial ou câble bifilaire ou encore câble torsadé, un modèle possible est donné par les relations (vii) et (viii) : 25 a(Î) = aoj + ctif /3(f) = bof Les nombres réels al:), a1 et bo sont les paramètres de la régression. Un 30 modèle linéaire est choisi pour le facteur de phase ce qui permet de représenter l'aspect faiblement dispersif des câbles précités. Un modèle non 3025320 14 linéaire utilisant un polynôme en racine carré de la fréquence est choisi pour l'atténuation linéique afin de tenir compte de deux types de pertes principales : les pertes métalliques et les pertes diélectriques. De façon générale, le modèle d'évolution des grandeurs a et p est établi à 5 partir de connaissances des caractéristiques physiques du câble à analyser. Sans sortir du cadre de l'invention, les modèles d'évolution des paramètres a et p peuvent prendre la forme de toute fonction polynomiale de la fréquence. 10 A partir des mesures obtenues à l'issu de l'étape 104 et des modèles donnés par les relations (vii) et (viii), on effectue une régression selon le critère des moindres carrés de sorte à approximer au mieux les paramètres de la régression à partir des mesures effectuées à l'étape 104. 15 D'autres méthodes de filtrage sont envisageables, par exemple un filtrage par moyenne glissante ou un filtrage par transformée de Fourier ou toute autre méthode permettant de filtrer les variations des mesures des paramètres a et p autour d'une valeur cible. 20 Pour améliorer encore la précision des résultats, selon une autre variante de réalisation de l'invention, une étape supplémentaire de décimation fréquentielle est réalisée avant l'opération de filtrage effectuée sur les grandeurs a et p. Autrement dit, si les mesures des paramètres a et p sont réalisées pour un nombre N de valeurs de fréquences dans une bande de 25 fréquence donnée [Fmin Fmaxl, on sélectionne un sous ensemble de ces valeurs avant d'appliquer l'opération de filtrage. Par exemple, une solution consiste à retirer au moins un point de mesure dans les basses fréquences de la bande considérée ou dans les hautes fréquences ou alternativement dans les basses fréquences puis dans les 30 hautes fréquences. 3025320 15 Comme cela est explicité plus en avant de la description, le procédé selon l'invention peut être exécuté de façon itérative. Dans ce cas, à chaque itération, au moins un point de mesure différent peut être retiré de l'ensemble des mesures disponibles pour a et 13. 5 Ce processus permet d'éliminer les valeurs de a et 13 qui présentent des distorsions trop importantes et d'assurer une convergence vers des valeurs stables pour les paramètres de la régression linéaire. Le procédé selon l'invention peut être interrompu après le calcul des 10 paramètres a et 13 qui peuvent être exploités en tant que tels pour l'information qu'ils contiennent. Si l'on souhaite déterminer une estimée des paramètres linéiques RLCG, une étape supplémentaire 106 doit être exécutée. 15 Selon cette étape supplémentaire 106, on détermine une estimée des paramètres linéiques RLCG à partir des mesures filtrées des paramètres a et P. Ces estimées peuvent être obtenues à partir des équations (iii) et des paramètres de la régression linéaire donnés par les équations (vii) et (viii). 20 , (1-F pe) z. G 2n(1- pe) bo 4a0 In r= bo = q, (ix) 4144 / 2 +bP 2n2 \ ai bî, I(4a.î+b(î)127r {L = 1 R = rj C= c G = gf (x) 3025320 16 Les relations (x) donnent également un modèle d'évolution des paramètres RLCG en fonction de la fréquence qui dépend des caractéristiques physiques du câble. C'est pourquoi les expressions (x) sont 5 valables pour des câbles du type câble coaxial ou câble bifilaire ou câble torsadé. Les relations (x) peuvent prendre la forme de fonctions polynomiales de la fréquence d'ordre supérieur à celles décrites ci-dessus dans le cas de câble différents, par exemple dans le cas de câbles en toron. 10 Selon une variante particulière de réalisation de l'invention, le procédé peut être itéré plusieurs fois afin d'améliorer le filtrage 105 des erreurs sur les paramètres a et [3. Pour cela, dans une étape supplémentaire 202, on effectue une 15 reconstruction Hi de la fonction de transfert fréquentielle du câble à partir des paramètres RLCG estimés à l'étape 106 grâce aux relations (x). Cette reconstruction peut être obtenue à l'aide d'une méthode basée sur une matrice chaine telle que décrite par exemple dans l'un des documents [4] ou [5]. 20 Un test de similitude 203 peut être effectué pour mesurer le degré de similitude entre la fonction de transfert reconstruite Hi à l'issue de l'étape 202 et la fonction de transfert H mesurée initialement à l'étape 101. Ce test peut par exemple être réalisé en calculant le coefficient de 25 détermination R entre les deux fonctions de transfert Hi et H ou à partir d'un coefficient de vraisemblance V défini par la relation V2 = -log(1-R2). Dans une variante de réalisation de l'invention, les étapes 202,203 peuvent être réalisées en remplaçant la fonction de transfert fréquentielle H 30 par tout autre paramètre caractéristique X de la propagation dans la ligne de transmission. En particulier, la fonction de transfert H peut être remplacée 3025320 17 par le facteur de propagation y ou par l'impédance caractéristique Zc ou encore par une fonction de l'un de ces paramètres. Dans ce cas le degré de similitude est déterminé en comparant l'estimée et la mesure du paramètre X de la même façon que celle introduite ci-dessus pour la fonction de transfert 5 H. A chaque itération du procédé, on recalcule, à partir des paramètres R, L,C,G estimés à l'étape 106, l'impédance caractéristique Zc à l'aide de la relation (i) puis le coefficient de réflexion en entrée et/ou le coefficient de 10 réflexion en sortie à partir des relations qui les lient à l'impédance caractéristique. Les étapes 103,104,105,106 sont ensuite exécutées à partir de ces nouvelles valeurs. Alternativement, il est également possible de déterminer directement une nouvelle valeur du facteur de propagation y à partir de la relation (ii) puis 15 d'exécuter ensuite les étapes 104,105,106 pour déterminer de nouveaux paramètres linéiques. A chaque itération du procédé, on mémorise les paramètres linéiques obtenus associés au degré de similitude calculé à l'étape 203. Au bout d'un 20 nombre donné d'itérations, on conserve in fine les paramètres correspondants au degré de similitude le plus élevé, c'est-à-dire pour lequel la fonction de transfert reconstruite à l'étape 202 est la plus proche de celle mesurée à l'étape 101. Le fonctionnement itératif du procédé permet d'affiner les paramètres 25 R, L, C, G estimés notamment du fait que le filtrage des paramètres a, 13 est amélioré à chaque itération. Alternativement, il est également possible de paramétrer le nombre d'itérations du procédé en fonction d'un critère de convergence qui peut 30 consister à comparer le degré de similitude à un seuil de convergence et à arrêter les itérations lorsque la convergence est obtenue, autrement dit 3025320 18 lorsque le degré de similitude est supérieur au seuil de convergence. Cette variante de fonctionnement permet d'optimiser le nombre d'itérations du procédé. Comme explicité ci-dessus, à chaque itération, il est aussi possible 5 d'effectuer une décimation en fréquence des mesures des paramètres a, 13 avant filtrage. Les applications du procédé selon l'invention sont nombreuses. On peut citer sans être exhaustif ni limitatif les applications suivantes : le suivi du 10 vieillissement d'un câble à partir du suivi de l'évolution de ses paramètres linéiques, la définition d'une base de données contenant ces paramètres ou l'amélioration de méthodes de réflectométrie connues. 3025320 19 Références [1] "Analysis of the parameters of a lossy coaxial cable for cable fault location", Shi, Troltzsch, Kanoun 5 [2] "Inverse Scattering Experiments for Electric Cable Soft Fault Diagnosis and Connector Location", Loete, Zhang, Sorine [3] "S-Parameter-Based IC Interconnect Transmission Line Characterization", Eisenstadt [4] "A Mathematical and Lumped-Element Model for Multiple Cascaded Lossy 10 transmission Lines with Arbitrary impedances and discontinuities", Scott Huss [5] "System simulation of network analysis for a lossy cable system", Qinghai Shi, Olfa Kanoun [6] "Extraction of high frequency power cable characteristics from S15 parameter measurements", R. Papazyan et al. REVENDICATIONS1. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission, ledit procédé comprenant les étapes suivantes : - Déterminer (103) au moins une mesure du facteur de propagation complexe y en fonction de la fréquence à partir d'au moins une mesure réalisée sur la ligne de transmission, - Déterminer (104) au moins une mesure de l'atténuation linéique a de la ligne de transmission égale à la partie réelle de la mesure du facteur de propagation complexe y et/ou au moins une mesure du facteur de phase 13 de la ligne de transmission égale à la partie imaginaire de la mesure du facteur de propagation complexe y, - Filtrer (105) la mesure de l'atténuation linéique a et/ou la mesure du facteur de phase 13 à partir d'un modèle de régression polynomial en fréquence dépendant des caractéristiques physiques de la ligne de transmission. 2. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon la revendication 1 dans lequel le modèle de régression pour l'atténuation linéique est du type a(f) = aoj + aif , où ao et a1 sont deux nombres réels. 3. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une des revendications précédentes dans lequel le modèle de régression pour le facteur de phase est du type /3(f) = bof où bo est un nombre réel. 4. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une des revendications précédentes comprenant en outre la détermination (106) d'une estimée d'au moins un paramètre linéique de la ligne de transmission parmi la résistance linéique, l'inductance linéique, la capacitance linéique ou la conductance linéique, à partir d'au moins une mesure parmi la mesure filtrée de l'atténuation linéique, la mesure filtrée du facteur de phase et la mesure du coefficient 5 de réflexion pe à une première extrémité de la ligne de transmission, la dite estimée étant déterminée comme une fonction polynomiale de la fréquence. 5. Procédé de détermination d'au moins un paramètre linéique d'une ligne 10 de transmission selon l'une des revendications précédentes comprenant en outre les étapes suivantes : - Déterminer (202) une nouvelle estimée Xi d'un paramètre caractéristique de la propagation dans la ligne de transmission à partir de ladite estimée d'au moins un paramètre linéique de la 15 ligne de transmission, - Mesurer (203) le degré de similitude entre la nouvelle estimée Xi et une mesure X du même paramètre caractéristique. 6. Procédé de détermination d'au moins un paramètre linéique d'une ligne 20 de transmission selon la revendication 5 dans lequel le paramètre caractéristique de la propagation dans la ligne de transmission est le facteur de propagation complexe y ou une fonction de ce paramètre. 7. Procédé de détermination d'au moins un paramètre linéique d'une ligne 25 de transmission selon l'une des revendications précédentes dans lequel la mesure du facteur de propagation complexe y est déterminée à l'aide des étapes suivantes : - Déterminer (101) une mesure H de la fonction de transfert fréquentielle de la ligne de transmission, 30 - Déterminer (102) une mesure du coefficient de réflexion pe à une première extrémité de la ligne de transmission, - Déterminer (103) une mesure du facteur de propagation complexe y à partir de la mesure H de la fonction de transfert fréquentielle, de la mesure du coefficient de réflexion pe à une première extrémité de la ligne de transmission et d'une estimée de la 5 longueur de la ligne de transmission, 8. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon la revendication 7 dans lequel la mesure H de la fonction de transfert fréquentielle de la ligne de transmission est 10 déterminée (101) à l'aide d'une méthode de réflectométrie fréquentielle à partir de l'injection, à une première extrémité de la ligne, d'un signal de réflectométrie. 9. Procédé de détermination d'au moins un paramètre linéique d'une ligne 15 de transmission selon la revendication 8 dans lequel la mesure du coefficient de réflexion pe en entrée de la ligne de transmission est déterminée (102) par une moyenne, sur une bande fréquentielle donnée, de la mesure H de la fonction de transfert fréquentielle de la ligne de transmission. 20 10. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon la revendication 7 dans lequel la mesure H de la fonction de transfert fréquentielle de la ligne de transmission est déterminée (101) à l'aide d'une méthode de réflectométrie temporelle à 25 partir de l'injection, à une première extrémité de la ligne, d'un signal de réflectométrie. 11. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon la revendication 10 dans lequel la mesure du 30 coefficient de réflexion pe en entrée de la ligne de transmission est déterminée (102) par la mesure de l'amplitude du pic du signal de réflectométrie réfléchi au point d'injection. 12. Procédé de détermination d'au moins un paramètre linéique d'une ligne 5 de transmission selon l'une des revendications 7 à 11 comprenant en outre les étapes suivantes : - Déterminer (201) une mesure de l'impédance caractéristique Zc de la ligne de transmission ou une mesure du coefficient de réflexion à une seconde extrémité de la ligne de transmission à partir de la 10 mesure du coefficient de réflexion pe à une première extrémité de la ligne de transmission. - Déterminer (103) une mesure du facteur de propagation complexe y à partir de la mesure H de la fonction de transfert fréquentielle et de la mesure de l'impédance caractéristique Zc de la ligne de 15 transmission ou de la mesure du coefficient de réflexion à une seconde extrémité de la ligne de transmission. 13. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une des revendications 5 à 12 dans lequel le 20 paramètre caractéristique de la propagation dans la ligne de transmission est la fonction de transfert fréquentielle H ou l'impédance caractéristique Zc ou une fonction de l'un de ses paramètres. 14. Procédé de détermination d'au moins un paramètre linéique d'une ligne 25 de transmission selon l'une des revendications 5 à 13 dans lequel le degré de similitude est mesuré (202) par le calcul d'un coefficient de détermination entre la nouvelle estimée d'un paramètre caractéristique Xi de la propagation dans la ligne de transmission et la mesure X du même paramètre caractéristique. 15. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une des revendications 4 à 14 comprenant en outre les étapes suivantes : - Déterminer (103) une nouvelle estimée y; du facteur de 5 propagation complexe de la ligne de transmission à partir des estimées de la résistance linéique, de l'inductance linéique, de la capacitance linéique et de la conductance linéique de la ligne de transmission, - exécuter les étapes (104,105,106) du procédé permettant de 10 déterminer une estimée d'au moins un paramètre linéique de la ligne de transmission à partir de la nouvelle estimée y; du facteur de propagation complexe. 16. Procédé de détermination d'au moins un paramètre linéique d'une ligne 15 de transmission selon l'une des revendications 5 à 15 dans lequel les étapes du procédé sont itérées un nombre donné de fois et les paramètres linéiques retenus sont ceux déterminés à l'itération correspondant au degré de similitude le plus élevé. 20 17. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une des revendications 5 à 15 dans lequel les étapes du procédé sont itérées tant que le degré de similitude n'a pas atteint une valeur donnée. 25 18. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une des revendications précédentes dans lequel la ligne de transmission est un câble coaxial ou un câble bifilaire ou un câble torsadé. 30 19. Procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une des revendications précédentes comprenant en outre une étape de décimation en fréquence appliquée à une pluralité 3025320 25 de mesures de l'atténuation linéique a et/ou de mesures du facteur de phase 13 avant l'étape de filtrage (105). 20. Programme d'ordinateur comportant des instructions pour l'exécution du 5 procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une quelconque des revendications 1 à 19, lorsque le programme est exécuté par un processeur. 21. Dispositif comprenant des premiers moyens de mesure connectés à une 10 ligne de transmission et des seconds moyens de calcul configurés pour exécuter le procédé de détermination d'au moins un paramètre linéique d'une ligne de transmission selon l'une quelconque des revendications 1 à 19. 15