Comment convertir les paramètres S en paramètres Y et Z pour l'analyse de circuit

Sur le plan mathématique, l'opération est:
Z = Z[₀ (I + S) (I - S)⁻¹

Pour les réseaux à deux ports, le calcul manuel est gérable, mais pour les comptes de port plus élevés, le calcul numérique est recommandé. Nous examinerons un exemple pratique sous peu, y compris l'extraction des paramètres Y par inversion de Z.

Rénormalisation avant conversion

Une surveillance commune se produit lorsque les paramètres S sont mesurés ou simulés avec différentes impédances de référence à chaque port. Par exemple, un dispositif différentiel peut avoir des paramètres S normalisés à 100 γ dans les ports différentiels et 50 γ dans les ports à une seule extrémité. L'application directe de la formule scalaire Z[₀ donnerait des valeurs d'immittance incorrectes. La renormalisation transforme la S‐matrix en un ensemble commun d'impédances de référence avant la conversion. La formule générale de renormalisation est :

S' = R⁻¹ (S - ι)(I - ι S)⁻¹ R

où Ι est une matrice diagonale des coefficients de réflexion calculés à partir des anciennes et nouvelles impédances de référence, et R est une matrice diagonale qui normalise les amplitudes d'onde. Cette transformation est mise en œuvre dans la plupart des outils de simulation RF. Lors de la conversion manuelle, confirmez toujours que tous les ports partagent une impédance de référence commune avant d'utiliser les formules simplifiées. Le Marks , papier sur l'étalonnage multiligne TRL fournit un traitement approfondi de ce sujet, en particulier pour les mesures traçables sur l'échafaudage où les impédances de ligne peuvent varier en fréquence.

Conversion des paramètres S en paramètres Y

Il existe deux voies principales pour obtenir des paramètres Y à partir de paramètres S. La première est une formule directe analogue à la conversion Z :

Y = Y₀ (I - S) (I + S)⁻¹

où Y₀ est la matrice d'admission de référence (entrées diagonales = 1 / Z_0k). Pour un Z₀ scalaire, ceci devient:

Y = (1/Z₀) (I - S) (I + S)⁻¹

La seconde méthode, souvent plus sûre, exploite la relation inverse entre Z‐ et Y‐matrices : une fois la matrice Z‐matrice calculée de façon fiable, une inversion simple donne la matrice Y. C'est-à-dire :

1. Convertir S en Z en utilisant la méthode ci-dessus.
2. Calculer Y = Z⁻¹.

Cette approche en deux étapes peut être plus robuste car la matrice à inverser (I + S) peut être mieux conditionnée que (I - S) dans certains cas. Cependant, les deux étapes d'inversion doivent être manipulées avec soin. La stabilité peut être améliorée par l'utilisation de pseudo-inversion ou d'inversion régularisée si les matrices deviennent presque singulières. Pour les dispositifs bien appariés (=S_ij), les deux chemins produisent des résultats presque identiques; pour les réseaux très réfléchissants, la conversion directe en Y peut être préférable parce que (I + S) tend à rester diagonalement dominante.

Exemple de travail pour un réseau à deux ports

Considérez un réseau à deux ports mesuré à une seule fréquence avec des paramètres S :

S11 = 0.1∠40° S12 = 0.02∠-60°
S21 = 4.5∠120° S22 = 0.3∠-30°

Supposons une impédance de référence Z[_0[ = 50 ↓ aux deux ports. La matrice S en forme rectangulaire (arrondie à trois décimales) est:

S = [[0.0766 + j0.0643, 0.0100 - j0.0173],
 [ -2.2500 + j3.8971, 0.2598 - j0.1500]]

La matrice d'identité I est [[1,0],[0,1]]. Calculer I - S et I + S:

I - S = [[0.9234 - j0.0643, -0.0100 + j0.0173],
 [ 2.2500 - j3.8971, 0.7402 + j0.1500]]

I + S = [[1.0766 + j0.0643, 0.0100 - j0.0173],
 [-2.2500 + j3.8971, 1.2598 - j0.1500]]

Maintenant résolvez (I - S) * X = (I + S) pour X. En utilisant un solveur numérique (voir le code ci-dessous), nous obtenons:

X = [[-0.2475 - j0.5493, 0.0011 + j0.0064],
 [ 1.9736 - j0.3036, 0.9162 - j0.1786]]

Ensuite, Z = Z₀ * X = 50 * X, donnant:

Z11 ≈ -12.37 - j27.47 Ω
Z12 ≈ 0.057 + j 0.32 Ω
Z21 ≈ 98.68 - j15.18 Ω
Z22 ≈ 45.81 - j 8.93 Ω

À partir de cette Z‐matrix, la Y‐matrix est obtenue par inversion:

Y = Z−1 ≈ [[-0.0122 + j0.0272, -0.0001 - j0.0007],
 [ 0.0216 + j0.0033, 0.0068 + j0.0048]] siemens

Cet exemple démontre que même avec des grandeurs de S-paramètres modestes, les valeurs d'immittance qui en résultent peuvent être significatives et révéler des caractéristiques de couplage, de perte et de réaction. La partie réelle négative de Z[₁₁ indique une instabilité potentielle à cette fréquence, une observation critique pour la conception de l'amplificateur.

Considérations avancées et stabilité numérique

Les ingénieurs effectuant des conversions répétées à travers les balayages de fréquence doivent être conscients des écueils numériques. Le nombre de conditions (I ± S) peut se détériorer près des fréquences où l'appareil présente une réflexion extrême (=S_ii== ou lorsque le réseau est mal adapté.

• Vérifiez le numéro de condition avant l'inversion:[ De nombreuses bibliothèques d'algèbre linéaire fournissent une estimation du numéro de condition. Si le numéro de condition dépasse, par exemple, 10⁶, envisagez de raffiner les données ou d'appliquer une régularisation. Un numéro de condition de l'ordre de 10⁸ ou plus indique généralement que les résultats de conversion peuvent être peu fiables.
• Renormalisation de la terminaison des ports:[ Si les paramètres S sont pris avec différentes impédances de référence, utilisez des formules de renormalisation pour amener tous les ports à un Z₀ commun avant d'appliquer la conversion standard. L'expression générale implique des matrices de transformation des vagues et est documentée dans Kurokawa=s power wave formalisme.
• Regularisation: Pour les matrices mal conditionnées, la régularisation de Tikhonov ou la SVD tronquée peut donner des paramètres Z utilisables, bien que cela introduit un compromis de biais-variance. Dans la pratique, une petite charge diagonale (p. ex. 10^]−6 * max=S=) peut stabiliser l'inversion sans modifier significativement le résultat. Cette technique est particulièrement utile pour traiter des données mesurées qui contiennent des pics sonores.
• Autres formules directes:[ La conversion directe par paramètre Y (I - S)(I + S)[⁻¹ peut être mieux conditionnée lorsque l'appareil a une faible perte d'insertion, car (I + S) peut devenir quasi-diagonale pour des réseaux bien appariés. Inversement, pour les appareils à gain élevé, (I - S) peut être plus facile à inverser. Une approche pragmatique consiste à calculer Z et Y par les deux voies et à comparer les résultats; des différences significatives indiquent des problèmes numériques.
• Application de la passivité:[ Les paramètres S mesurés violent souvent la passivité due à des erreurs de bruit ou d'étalonnage. Avant la conversion, appliquer la passivité à l'aide d'algorithmes tels que la réduction de la valeur propre ou des méthodes basées sur Hamilton pour empêcher des valeurs d'immittance non physiques.

Mise en œuvre des logiciels: MATLAB, Python et Julia Scripts

Pour automatiser la conversion sur des centaines de points de fréquence, les ingénieurs comptent souvent sur le scripting. Ci-dessous sont des extraits de code concis et réutilisables. Les deux exemples supposent que S est une matrice complexe N×N et que Z0 est soit un vecteur scalaire ou N-element d'impédances de référence (un par port).

Fonction MATLAB

function Z = s2z(S, Z0)
 % S2Z Convert S-parameters to Z-parameters
 % S : NxN scattering matrix
 % Z0 : Nx1 vector or scalar reference impedance
 N = size(S,1);
 if isscalar(Z0)
 Z0 = Z0 * eye(N);
 else
 Z0 = diag(Z0);
 end
 I = eye(N);
 Z = sqrt(Z0) * ((I - S) \ (I + S)) * sqrt(Z0);
end

Pour les paramètres Y, on peut simplement appeler `Y = inv(s2z(S, Z0))` ou mettre en œuvre une fonction `s2y` dédiée en utilisant la formule alternative: `Y = (1/Z0) * ((I - S) / (I + S))` (en utilisant l'opérateur de division de droite `/`). Notez que la division de droite est équivalente à se multiplier par l'inverse de (I + S) à partir de la droite.

Python avec NumPy

import numpy as np

def s_to_z(S, Z0):
 N = S.shape[0]
 if np.isscalar(Z0):
 Z0 = Z0 * np.eye(N)
 else:
 Z0 = np.diag(Z0)
 I = np.eye(N)
 Z0_sqrt = np.sqrt(Z0) # element-wise; Z0 diagonal
 Z = Z0_sqrt @ np.linalg.solve(I - S, I + S) @ Z0_sqrt
 return Z

Pour les paramètres Y, une approche directe est :

def s_to_y(S, Z0):
 N = S.shape[0]
 if np.isscalar(Z0):
 Y0 = (1/Z0) * np.eye(N)
 else:
 Y0 = np.diag(1 / Z0)
 I = np.eye(N)
 Y = Y0 @ np.linalg.solve(I + S, I - S)
 return Y

Mise en œuvre de Julia

function s2z(S, Z0)
 N = size(S,1)
 I = Matrix{ComplexF64}(I, N, N)
 Z0_mat = isa(Z0, Number) ? Z0 * I : diagm(Z0)
 Z0_sqrt = sqrt.(Z0_mat)
 return Z0_sqrt * ((I - S) \ (I + S)) * Z0_sqrt
end

Pour les données synthétiques sans bruit, ces fonctions donnent des résultats exacts. Lorsqu'on travaille avec des paramètres S mesurés, il est conseillé de lisser les données et de supprimer les valeurs aberrantes avant la conversion pour éviter les pics d'impédance fallacieux. La bibliothèque open-source scikit‐rf fournit des routines testées par la communauté pour ces conversions, y compris des contrôles de renormalisation et de passivité. scikit‐rf gère également les réseaux multiports et peut lire automatiquement les fichiers Touchstone avec des impédances de référence variables.

Interprétation physique et utilisation des paramètres convertis

Une fois que vous avez le Y ou Z-matrix, une abondance d'analyse devient possible. Les paramètres Z permettent une évaluation simple de l'impédance d'entrée du réseau dans des conditions de charge arbitraires. Par exemple, l'impédance d'entrée du port 1 lorsque le port 2 est terminé dans Z_L est:

Z_en = Z[₁₁ - (Z[₁₂ Z₂₁) / [Z₂₂ + Z_L]]]

Cette expression est fondamentale dans la conception de l'amplificateur pour calculer les cercles de stabilité et les réseaux correspondants. Les paramètres Y alimentent de la même façon la conception des oscillateurs et des amplificateurs de résistance négative, où la détermination des régions d'instabilité potentielle repose sur la matrice d'admission.

Les paramètres convertis permettent également l'extraction de circuits équivalents à éléments remorqués. En adaptant les éléments Z ou Y dépendants de la fréquence aux fonctions rationnelles, on peut synthétiser des modèles compacts adaptés à la simulation du domaine temporel dans SPICE. Ce processus est central pour la modélisation de l'intégrité des signaux et de la compatibilité électromagnétique. Par exemple, une matrice Z qui culmine à une fréquence résonante peut être modélisée comme un circuit-citerne RLC dont les valeurs de composants sont dérivées directement des données d'immittance.

Pièges communs et stratégies de vérification

Malgré la simplicité apparente des formules, des erreurs peuvent survenir de plusieurs sources :

• Impédances de référence différenciées: Si le fichier S-paramètre n'indique pas l'impédance de référence, ou si elle varie selon la fréquence, vous devez renormaliser avant la conversion. De nombreux fichiers de touches modernes incluent l'impédance de référence pour chaque point de données; vérifiez attentivement la ligne d'en-tête.
• Passivité et causalité:[ Le processus de conversion peut créer des matrices non passives Z ou Y si les paramètres S originaux violent la passivité (p. ex., en raison du bruit de mesure). Appliquer la passivité par des algorithmes comme le quadrage matriciel ou des méthodes basées sur Hamilton si nécessaire. Un Z‐matrix passif doit avoir une partie réelle définie positive. De même, les paramètres d'immittance résultants devraient satisfaire les relations Kramers‐Kronig si le réseau est causal; l'échec indique des problèmes de qualité des données.
• Unités de confusion : Les paramètres S sont inunité, mais les paramètres Z résultants sont en ohms et en Y-paramètres dans les siemens. Confirmez toujours que votre impédance de référence est constante dans les ohms. Si vos paramètres S sont en magnitude dB, convertissez-les en magnitude linéaire avant la conversion.
• Ordre des ports: Pour les réseaux multiports, maintenir une numérotation cohérente des ports. Un indice échangé donnera un couplage interparamètre incorrect. Il est de bonne pratique d'étiqueter les ports et de les comparer aux mesures connues en ligne ou en circuit ouvert. Certains appareils de test de ré-embudding routines réorder les ports, donc toujours vérifier la cartographie.
• Echantillonnage fréquent : La conversion à chaque point de fréquence est indépendante. Évitez d'interpoler les paramètres S avant la conversion; au lieu de cela, interpolez les paramètres Z ou Y résultants si une grille de fréquence uniforme est nécessaire. L'interpolation dans le domaine complexe de S peut introduire des artefacts non causaux qui corrompent les résultats de l'immittance.

Pour vérifier une conversion, retransformer les paramètres Z en paramètres S calculés en utilisant la formule inverse :

S = (Z - Z₀) (Z + Z₀)⁻¹

Si les paramètres S aller-retour correspondent à l'original en précision numérique (par exemple, erreur max < 10^-12]), la conversion est correcte. De même, convertir S → Z → Y et ensuite Y de retour en Z devrait récupérer la Z-matrix originale dans l'epsilon machine. Ce test aller-retour est une pratique exemplaire qui devrait être automatisée dans tout script de conversion.

Impédances de référence de la fréquence

Dans les mesures pratiques des analyseurs de réseau vectoriels, l'impédance de référence peut ne pas être constante avec la fréquence, en particulier lorsque les effets de montage de dé-embudding ou l'utilisation de lignes de transmission non-50---. Formules de conversion généralisées qui permettent de varier la fréquence Z₀ utilisent les relations pleines d'onde:

Z = (I - S R_w[)⁻¹ (I + S R_w) Z[₀

Cette approche, détaillée dans Marques du papier sur l'étalonnage multiligne TRL, est essentielle pour maintenir la traçabilité aux normes de référence non-50‐=. Par exemple, dans les systèmes de mesure cryogénique, l'impédance de la ligne peut varier considérablement en fonction de la température, nécessitant une renormalisation fréquence par fréquence. Heureusement, pour la plupart des applications courantes où l'impédance de référence est 50 = Résistive et constante, les matrices simplifiées tiennent parfaitement. Lorsque l'on travaille avec des mesures de niveau wafer où l'impédance de la ligne varie, la mise en œuvre de la conversion généralisée complète est obligatoire pour éviter des erreurs systématiques dans les données d'immittance extraites.

Connexion à des flux de travail plus larges

La conversion entre types de paramètres est rarement une fin en soi. En général, elle est un précurseur de tâches telles que:

• Analyse de stabilité: Calcul du facteur de stabilité Rollett K et Δ à partir des paramètres S, puis conversion en Z ou Y pour dessiner des cercles de stabilité sur le graphique Smith. Les paramètres Z convertis fournissent le gain de boucle ouverte nécessaire à l'analyse Nyquist. De nombreux manuels fournissent des formules pour les cercles de stabilité directement en termes de paramètres Z ou Y.
• Résistance des filtres: En utilisant la matrice Y extraite pour déterminer les coefficients de couplage et les fréquences de résonance pour les filtres planaires et coaxiaux. Les paramètres Y hors-diagonal représentent directement l'admission mutuelle entre les résonateurs, qui peut être utilisée pour concevoir des filtres couplés avec des zéros de transmission prescrits.
• Dé-embuding:[ Suppression des effets des appareils d'essai en convertissant les paramètres S en paramètres Z, en soustrayant le montage Z-matrix (obtenu d'une mesure ou simulation séparée), et en se convertissant à S. C'est une pratique courante dans la caractérisation du transistor sur le wafer, où le parasite du pad affecte de façon significative les données mesurées.
• Extraction de paramètres de modèle: Monter les données du paramètre Z sur la fréquence à des modèles de circuits équivalents basés sur la physique, permettant d'optimiser les configurations des transistors et les composants passifs. Les parties réelles et imaginaires de Z_ij limitent la résistance et la réactivité de chaque branche.

Des ressources en ligne telles que l'encyclopédie Mibrowaves101 S-parameter et la documentation scikit-rf fournissent d'autres exemples pratiques et des outils prêts à l'emploi pour ces workflows. Les ingénieurs sont encouragés à explorer ces références pour approfondir leur compréhension du processus de conversion et de ses applications.

Résumé et pratiques exemplaires

La conversion des paramètres S en paramètres Y et Z est une pierre angulaire de la conception des circuits RF et micro-ondes. Les étapes mathématiques sont bien définies, mais une application réussie exige une attention à l'impédance de référence, au conditionnement matriciel et à la sélection des méthodes numériques.

• Utiliser la formule appropriée:[ Z = Z[₀ (I + S) (I - S)[⁻¹ pour les scalaires Z[₀; pour les références non uniformes, utiliser la transformation de la variable d'onde complète ou rénormaliser en premier.
• Préférence des résolveurs linéaires de système par rapport aux inverses explicites: Les méthodes comme `mldidivid` dans MATLAB ou `numpy.linalg.solve` dans Python sont plus stables et plus rapides.
• Vérifier les conversions par des tests aller-retour:[ Back‐transformer en paramètres S pour assurer la cohérence. Une erreur maximale inférieure à 10⁻¹⁰ indique un résultat digne de confiance.
• Numéros de condition de moniteur:[ Détecter les fréquences problématiques tôt et appliquer lissage ou régularisation. Le numéro de condition de (I - S) doit être calculé à chaque point de fréquence.
• Inforcez la passivité au besoin:[ Utilisez des algorithmes établis pour restaurer la passivité dans le domaine de l'immittance avant d'utiliser les données pour l'analyse de synthèse ou de stabilité.
• Intégrer la conversion en un flux de conception plus large:[ Utiliser des bibliothèques établies comme scikit‐rf pour la vitesse et la fiabilité des véhicules communautaires. Ces bibliothèques gèrent également les conversions de renormalisation et de multiports qui sont susceptibles d'erreur pour être mises en œuvre à partir de zéro.

En internalisant ces principes, les ingénieurs peuvent combler avec confiance l'écart entre les données mesurées de diffusion et l'analyse basée sur l'immittance qui sous-tend la conception pratique de haute fréquence. La capacité de se déplacer couramment entre les ensembles de paramètres réduit non seulement la dépendance à un seul type de mesure, mais approfondit également la compréhension du comportement du réseau dans différentes conditions de terminaison.