Come convertire i parametri S in parametri Y e Z per l'analisi del circuito

dove Z₀] è la caratteristica matrice di impedenza, io sono la matrice di identità, e la radice quadrata indica la radice quadrata della matrice (solitamente banale per una diagonale Z₀]]]]). Quando tutti i porti usano la stessa impedenza Z_mp

Z = Z₀ (I + S) (I - S)^]−1]

Questa forma è ampiamente utilizzata e deriva dal presupposto che Z₀]] è scalare. Per evitare confusione, verificare sempre se i tuoi parametri S sono normalizzati a una singola impedenza scalare o a una matrice diagonale con valori variabili. Nei file commerciali touchstone, l'impedenza di riferimento per ogni porta è spesso fornita nell'intestazione; se non, il default è di solito 50Ω.

Procedura dettagliata

1. Costruire la S‐matrix:[ Assicurare che i dati del parametro S siano disposti come una matrice complessa N×N, dove N è il numero di porte. Rappresentare ogni voce in forma rettangolare (reale + immaginaria) per il calcolo numerico; la forma polare può essere utilizzata per il display ma deve essere convertita per l'aritmetica.
2. Form the I - S matrix:[] Sottrarre la S-matrix dalla matrice di identità della stessa dimensione.
3. Computa l'inverso (I - S)⁻¹]: Questo passaggio è fondamentale. La matrice deve essere non-singolare. Controllare il numero di condizione o determinante; se è ill-condizionato, la conversione può fallire o produrre risultati imprecisi.
4. Form I + S:[] Aggiungi la matrice di identità alla S-matrix.
5. Multiply (I + S) per l'impedenza di riferimento: Se Z_0[ è scalare, moltiplicare ogni elemento di (I + S) per quel valore. Per una matrice di impedenza diagonale, eseguire moltiplicazione di ogni riga in senso di elemento secondo l'impedenza corrispondente (o diagonale moltiplicazione Ztrix)
6. Multiply il risultato di (I - S)^]−1]:[] Il prodotto matrice finale produce la matrice Z-parametro. Preferire risolvere il sistema lineare (I - S) * X = (I + S) e poi moltiplicare per Z_0[FFT:

Matematicamente, l'operazione è:
]]Z = Z₀ (I + S) ^]−1[]]

Per le reti a due rapporti, il calcolo manuale è gestibile, ma per i conti delle porte più elevate è consigliato il calcolo numerico.

Renormalizzazione Prima della Conversione

Una comune sovrapposizione si verifica quando i parametri S sono misurati o simulati con diverse impedenze di riferimento in ogni porta. Ad esempio, un dispositivo differenziale potrebbe avere parametri S normalizzati a 100 Ω in porte differenziali e 50 Ω in porte single-ended.

S' = R[⁻¹ (S - Γ)(I - Γ S)⁻¹ R]

dove la figura di questa matrice diagonale dei coefficienti di riflessione calcolati dalle vecchie e nuove impedenze di riferimento, R è una matrice diagonale che normalizza le ampiezza dell'onda. Questa trasformazione viene implementata nella maggior parte degli strumenti di simulazione RF. Quando si eseguono conversioni manuali, sempre confermare che tutte le porte condividono una comune impedenza di riferimento prima di usare le formule semplificate.

Convertire S‐Parameters in Y‐Parameters

Ci sono due percorsi principali per ottenere i parametri Y da parametri S. Il primo è una formula diretta analoga alla conversione Z‐:

Y = Y₀ (I - S) (I + S)^]−1]

dove Y_0[] è la matrice di riconoscimento di riferimento (ingressi diagonali = 1 / Z[_0k]]]]]. Per una scalar Z[_]]]0[[]], questo diventa:

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

Il secondo metodo, spesso più sicuro, sfrutta il rapporto inverso tra Z‐ e Y‐matrices: una volta che la Z‐matrix è calcolata in modo affidabile, l'inversione semplice dà la Y‐matrix.

1. Convertire S in Z utilizzando il metodo sopra.
2. Computo Y = Z⁻¹.

Questo approccio a due fasi può essere più robusto perché la matrice per invertire, (I + S), può essere meglio condizionato di (I - S) in alcuni casi. Tuttavia, entrambi i passaggi di inversione devono essere gestiti con cura. La stabilità può essere migliorata utilizzando pseudo-inversione o inversione regolarizzata se le conversioni sono quasi singolari.

Esempio di lavoro per una rete a due porte

Considerare una rete a due rapporti misurata ad una singola frequenza con parametri S:

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

Assumere un'impedenza di riferimento Z[₀ = 50 Ω in entrambi i porti. La S-matrix in forma rettangolare (arrotondata a tre decimali) è:

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

La matrice di identità che sono [[1,0],[0,1]]. Computo I - S e 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]]

Ora risolvi (I - S) * X = (I + S) per X. Utilizzando un risolutore numerico (vedere il codice qui sotto), otteniamo:

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

Poi Z = Z₀] * X = 50 * X, che produce:

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

Da questa Z‐matrix, la Y‐matrix è ottenuta per inversione:

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

Questo esempio dimostra che anche con modeste magnitudine S‐parameter, i valori di immittanza risultanti possono essere significativi e rivelare le caratteristiche di accoppiamento, perdita e reattività. La parte negativa reale di Z₁₁] indica una potenziale instabilità a questa frequenza, un'osservazione critica per il design degli amplificatori.

Considerazioni avanzate e stabilità numerica

Gli ingegneri che effettuano ripetute conversioni attraverso le spazzate di frequenza devono essere consapevoli delle insidie numeriche. Il numero di condizione di (I ± S) può deteriorarsi vicino alle frequenze in cui il dispositivo mostra una riflessione estrema (|S_]ii[]]| vicino a 1) o quando la rete è scarsamente abbinata. In tali casi, l'inversione diretta può amplificare il rumore.

• Controllo il numero di condizione prima dell'inversione: Molte librerie di algebra lineari forniscono stima del numero di condizione. Se il numero di condizione supera, diciamo, 10⁶, consideri il rifinanziamento dei dati o l'applicazione di regolarizzazione.
• Rinormalizzazione della terminazione del porto:[] Se i parametri S sono presi con diverse impedenze di riferimento, utilizzare formule di renormalizzazione per portare tutti i porti ad un comune Z₀] prima di applicare la conversione standard. L'espressione generale coinvolge matrici di trasformazione dell'onda e viene documentata in [FLT: formalismo.
• Regolarizzazione: Per le matrici ill-condizionate, la regolarizzazione di Tikhonov o il SVD troncato possono produrre parametri Z utilizzabili, anche se questo introduce un trade-off di bias-variance. In pratica, una piccola tecnica di caricamento diagonale (ad esempio, 10⁻⁶]]
• Altre formule dirette:] La conversione diretta del parametro Y (I - S)(I + S)[[⁻¹ può essere meglio condizionata quando il dispositivo ha una bassa perdita di inserimento, perché (I + S) può diventare vicino a un'interfaccia numerica per reti ben mate.
• L'applicazione passiva:[] I parametri misurati violano spesso la passività a causa di errori di rumore o di calibrazione. Prima di conversione, applicare la passività utilizzando algoritmi come la riduzione dell'autovalore o metodi basati su Hamiltonian per prevenire i valori di simulazione non fisici.

Attuazione software: MATLAB, Python e Julia Scripts

Per automatizzare la conversione attraverso centinaia di punti di frequenza, gli ingegneri spesso si affidano alla scrittura. Di seguito sono frammenti di codice concisi e riutilizzabili. Entrambi gli esempi assumono S è una matrice complessa N×N e Z0 è uno scalare o un vettore di N-element di impedenze di riferimento (uno per porta). Queste implementazioni evitano l'inversione di matrice esplicita utilizzando risolutori di sistema lineari, che migliora l'accuratezza numerica.

Funzione 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

Per i parametri Y, si può semplicemente chiamare `Y = inv(s2z(S, Z0) o implementare una funzione `s2y` dedicata utilizzando la formula alternativa: `Y = (1/Z0) * ((I - S) / (I + S)`)` (utilizzando la destra-divisione).

Python con 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

Per i parametri Y, un approccio diretto è:

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

Julia Attuazione

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

Per i dati sintetici senza rumore, queste funzioni producono risultati esatti. Quando si lavora con parametri misurati, è consigliabile lisciare i dati e rimuovere i outlier prima della conversione per evitare spurie punte di impedenza. La libreria open-source scikit‐rf fornisce routine di impedenza per queste conversioni, comprese le revisioni e i controlli di passività.

Interpretazione fisica e utilizzo dei parametri convertiti

Una volta che avete la Y‐ o Z‐matrix, una ricchezza di analisi diventa possibile. I parametri Z consentono una valutazione semplice dell’impedenza di ingresso della rete in condizioni di carico arbitrarie. Ad esempio, l’impedenza di ingresso della porta 1 quando la porta 2 è terminata in Z_L] è:

] _] = Z_]11] [₁₂ Z₂₁]]] _{]] [[FLT]]]]]]] [[FLT]]]]]]]]]]]]] [[FLT]]]]] [[FLT]]]]]]]]]]] [[[[[[F[FLT]]]]]]]]]]]]]]]]]] [[F[F[F[F[F[F[FLT]]]]]]]]]]]]]]]]]]]]]]]]]]]]][[[F[F[F[F[F[F[F[FLT]]]]]]]]]]]]]]]]]]]]]]]]]]]]]}

S = (Z - Z₀]] [Z + Z[_]0]] ⁻¹[]]]]]

Se i parametri S-parametri a giro corrispondono all'originale in precisione numerica (ad esempio, errore max [< 10^-12]), la conversione è corretta. Allo stesso modo, la conversione S → Z → Y e poi Y di nuovo a Z dovrebbe recuperare la Z-matrix originale all'interno della macchina epsilon.

Impedenza di riferimento dipendente dalla frequenza

Nelle misurazioni pratiche dell'analizzatore di rete vettoriale, l'impedenza di riferimento non può essere costante con la frequenza, in particolare quando si verificano effetti di fissaggio o si utilizzano linee di trasmissione non-50-Ω.

Z = (I - S R[_]w[]]]^]−1[ [I + S R_]w]]]] _]]] ]][FLT]]]]][FLT]]]][FLT]]]]] [[FLT]]]]]]]]]]] [[FLT]]]]]]]] [[FLT]]]]]]]]] [[FLT]]] [[FLT]]]]]]]]]]]]]]]]]]]]] [[FLT]]]] [[[FLT]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]] [[[[[[[[[