def setUp(self):
self.n_ports = 2
self.wg = WG
wg = self.wg
self.X = wg.random(n_ports=2, name='X')
self.Y = wg.random(n_ports=2, name='Y')
self.If = wg.random(n_ports=1, name='If')
self.Ir = wg.random(n_ports=1, name='Ir')
self.gamma_f = wg.random(n_ports=1, name='gamma_f')
self.gamma_r = wg.random(n_ports=1, name='gamma_r')
actuals = [
wg.thru(),
rf.two_port_reflect(wg.load(-.98 - .1j), wg.load(-.98 - .1j)),
rf.two_port_reflect(wg.load(.99 + 0.05j), wg.load(.99 + 0.05j)),
wg.line(100, 'um'),
wg.line(200, 'um'),
wg.line(900, 'um'),
]
self.actuals = actuals
measured = [self.measure(k) for k in actuals]
self.cal = NISTMultilineTRL(measured=measured,
isolation=measured[1],
Grefls=[-1, 1],
l=[0, 100e-6, 200e-6, 900e-6],
er_est=1,
switch_terms=(self.gamma_f, self.gamma_r),
gamma_root_choice='real')
def setUp(self):
global NPTS
self.n_ports = 2
self.wg = WG
wg = self.wg
r = npy.random.uniform(10, 100, NPTS)
l = 1e-9 * npy.random.uniform(100, 200, NPTS)
g = npy.zeros(NPTS)
c = 1e-12 * npy.random.uniform(100, 200, NPTS)
rlgc = rf.media.DistributedCircuit(frequency=wg.frequency,
z0=None,
R=r,
L=l,
G=g,
C=c)
self.rlgc = rlgc
self.X = wg.random(n_ports=2, name='X')
self.Y = wg.random(n_ports=2, name='Y')
self.gamma_f = wg.random(n_ports=1, name='gamma_f')
self.gamma_r = wg.random(n_ports=1, name='gamma_r')
actuals = [
rlgc.thru(),
rlgc.short(nports=2),
rlgc.line(10, 'um'),
rlgc.line(100, 'um'),
rlgc.line(500, 'um'),
]
self.actuals = actuals
measured = [self.measure(k) for k in actuals]
self.measured = measured
self.cal = NISTMultilineTRL(measured=measured,
Grefls=[-1],
l=[0, 10e-6, 100e-6, 500e-6],
switch_terms=(self.gamma_f, self.gamma_r),
ref_plane=50e-6,
c0=c,
z0_ref=50,
gamma_root_choice='real')
def setUp(self):
self.n_ports = 2
self.wg = WG
wg = self.wg
self.X = wg.random(n_ports=2, name='X')
self.Y = wg.random(n_ports=2, name='Y')
self.gamma_f = wg.random(n_ports=1, name='gamma_f')
self.gamma_r = wg.random(n_ports=1, name='gamma_r')
# make error networks have s21,s12 >> s11,s22 so that TRL
# can guess at line length
self.X.s[:, 0, 0] *= 1e-1
self.Y.s[:, 0, 0] *= 1e-1
self.X.s[:, 1, 1] *= 1e-1
self.Y.s[:, 1, 1] *= 1e-1
actuals = [
wg.thru(),
rf.two_port_reflect(wg.load(-.9 - .1j), wg.load(-.9 - .1j)),
rf.two_port_reflect(wg.load(.92 + 0.05j), wg.load(.92 + 0.05j)),
wg.line(100, 'um'),
wg.line(200, 'um'),
wg.line(500, 'um'),
]
self.actuals = actuals
measured = [self.measure(k) for k in actuals]
self.cal = NISTMultilineTRL(
measured=measured,
Grefls=[-1, 1],
l=[0, 100e-6, 200e-6, 500e-6],
er_est=1,
switch_terms=(self.gamma_f, self.gamma_r),
gamma_root_choice='imag' #Losslessness of the line causes problems
#without this option
)
D_DUT_Measure = f'{Home_1}/{N_DUT}.s2p'
reflect = rf.Network(D_reflect)
T = rf.Network(D_thru)
R = rf.two_port_reflect(reflect, reflect)
L1 = rf.Network(D_Line_1)
L2 = rf.Network(D_Line_2)
measured = [T, R, L2]
cal = NISTMultilineTRL(
measured = measured,
Grefls=[-1],
l = [0, 7.25e-3],
gamma_root_choice = 'auto',
)
#Ideal DUT characteristic for compare
DUT_IDEAL = rf.Network(D_DUT_Ideal)
dut_raw = rf.Network(D_DUT_Measure)
dut_corrected = cal.apply_cal(dut_raw)
plt.figure()
dut_raw.plot_s_smith(m=1 - 1, n=1 - 1, label ='Measure DUT', linewidth='3')
dut_corrected.plot_s_smith(m=1 - 1, n=1 - 1, label ='Calculated DUT', linewidth='3')
DUT_IDEAL.plot_s_smith(m=1 - 1, n=1 - 1, label ='Ideal DUT', linewidth='3', linestyle='--')
plt.ylabel('S11')