Ex1: Ring SlotΒΆ

The Vector Fitting feature is demonstrated using the ring slot example network from the scikit-rf data folder. Additional explanations and background information can be found in the Vector Fitting tutorial.

[1]:
import skrf
import numpy as np
import matplotlib.pyplot as mplt

To create a VectorFitting instance, a Network containing the frequency responses of the N-port is passed. In this example the ring slot is used, which can be loaded directly as a Network from scikit-rf:

[2]:
nw = skrf.data.ring_slot
vf = skrf.VectorFitting(nw)

Now, the vector fit can be performed. The number of poles has to be specified, which depends on the behaviour of the responses. A smooth response would only require very few poles (2-5). In this case, 3 real poles are sufficient:

[3]:
vf.vector_fit(n_poles_real=3, n_poles_cmplx=0)

As printed in the logging output (not shown), the pole relocation process converged quickly after just 5 iteration steps. This can also be checked with the convergence plot:

[4]:
vf.plot_convergence()
../../_images/examples_vectorfitting_vectorfitting_ex1_ringslot_8_0.png

The fitted model parameters are now stored in the class attributes poles, residues, proportional_coeff and constant_coeff for further use. To verify the result, the model response can be compared to the original network response. One option is to analyze the rms error magnitude, which should be smaller than 0.05 for a good fit:

[5]:
vf.get_rms_error()
[5]:
0.003855803397460313

As the model will return a response at any given frequency, it makes sense to also manually check its response outside the frequency range of the original samples by plotting it at lower and higher frequencies:

[6]:
freqs1 = np.linspace(0, 200e9, 201)
fig, ax = mplt.subplots(2, 2)
fig.set_size_inches(12, 8)
vf.plot_s_mag(0, 0, freqs1, ax=ax[0][0]) # plot s11
vf.plot_s_mag(1, 0, freqs1, ax=ax[1][0]) # plot s21
vf.plot_s_mag(0, 1, freqs1, ax=ax[0][1]) # plot s12
vf.plot_s_mag(1, 1, freqs1, ax=ax[1][1]) # plot s22
fig.tight_layout()
mplt.show()
../../_images/examples_vectorfitting_vectorfitting_ex1_ringslot_12_0.png

To use the model in a circuit simulation, an equivalent circuit can be created based on the fitting parameters. This is currently only implemented for SPICE, but the structure of the equivalent circuit can be adopted to any kind of circuit simulator.

vf1.write_spice_subcircuit_s('/home/vinc/Desktop/ring_slot.sp')

For a quick test, the subcircuit is included in a schematic in QUCS-S for AC simulation and S-parameter calculation based on the port voltages and currents (see the equations): 8d9ae8833cc54743abe5170cbd603805

The simulation outputs from ngspice compare well to the plots above: 5e7a3a19d10b4f688f157b8974c4ffd1

e594f36152f84dd8955da1f508ccd722