from lmfit import minimize, Parameters, report_fit # , Parameter import matplotlib.pyplot as plt import PyGnuplot as gp plt.ion() # load hdf5 data # filename = 'S1_203.hdf5' filename = 'S1_945_S11_4p1_BPF7.hdf5' # filename = 'S1_945_S11_4p8_BPF7.hdf5' measdata = get_hdf5data(filename) flux0 = 2.07e-15 # Tm^2; Flux quanta: flux0 = h / (2*charging energy) squid = dim(name='Squid', start=-1.5, stop=0.5, pt=201, scale=flux0) elem = handler(name='mag/phase', start=0, stop=10, pt=1) squid.flux0 = 2.07e-15 # Tm^2; flux0=h/(2*charging energy) elem.Z0 = 50 # R; Input impedance elem.Z1 = 50 # R; Impedance of transmission piece 1 elem.Z2 = 50 # R; Impedance of Coplanar Waveguide elem.Z3 = 50 elem.L1 = 0.44 elem.L2 = 0.04 elem.L3 = 0.01 elem.Z4 = 0.1 # Ohm; Wire bonds conductance to GND (-45dB isolation) # m/s; approx. velocity in a coaxial 2/3 * speed of light elem.v = 2.0e8 squid.Wb = 1000e-12 # Aprox: Lw (nH) = 5.08x10-3 * L * (ln(4*L/D) - 1) squid.LOOP = 1e-30 squid.ALP = 1.0
from lmfit import minimize, Parameters, report_fit # , Parameter import matplotlib.pyplot as plt import Gnuplot as gp plt.ion() # load hdf5 data # filename = 'S1_203.hdf5' filename = 'S1_945_S11_4p1_BPF7.hdf5' # filename = 'S1_945_S11_4p8_BPF7.hdf5' measdata = get_hdf5data(filename) flux0 = 2.07e-15 # Tm^2; Flux quanta: flux0 = h / (2*charging energy) squid = dim(name='Squid', start=-1.5, stop=0.5, pt=201, scale=flux0) elem = handler(name='mag/phase', start=0, stop=10, pt=1) squid.flux0 = 2.07e-15 # Tm^2; flux0=h/(2*charging energy) elem.Z0 = 50 # R; Input impedance elem.Z1 = 50 # R; Impedance of transmission piece 1 elem.Z2 = 50 # R; Impedance of Coplanar Waveguide elem.Z3 = 50 elem.L1 = 0.44 elem.L2 = 0.04 elem.L3 = 0.01 elem.Z4 = 0.1 # Ohm; Wire bonds conductance to GND (-45dB isolation) # m/s; approx. velocity in a coaxial 2/3 * speed of light elem.v = 2.0e8 squid.Wb = 1000e-12 # Aprox: Lw (nH) = 5.08x10-3 * L * (ln(4*L/D) - 1) squid.LOOP = 1e-30 squid.ALP = 1.0
R = 2.3e3 # Ohm
Cap = 450e-15 # 450.0e-15 # F
Y4 = 1/0.1 # 1/Ohm; Wire bonds conductance to GND (-45dB isolation)
magnet = dim(name = 'Flux (Phi0)',
start = -1,
stop = 1,
pt = 1001,
scale = flux0)
freq = dim(name = 'Frequency (GHz)',
start = 4,
stop = 8,
pt = 101,
scale = 1e9)
dim_3 = handler(name = 'mag/phase',
start = 0,
stop = 10,
pt = 11) #8 pts for S 4x2 values
dim_3._Z0 = 50
head1 = make_header(magnet, freq, dim_3, 'S11 S12 S21 S22 Z L')
dim_3.prepare_data_save(magnet, freq, dim_3)
for jj, f0 in enumerate(freq.lin):
for ii, flux in enumerate(magnet.lin):
# b = k = 2pi/wavelength; wavelength = velocity / frequency
b = 2.0*pi*f0/v
L = flux0 / (Ic*2.0*pi* abs(cos(pi*flux/flux0)))
Ysq = (1.0/R + 1.0/(i*2.0*pi*f0*L +i*1e-90) + i*2.0*pi*f0*Cap)
Zsq = 1.0/Ysq