예제 #1
0
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
예제 #2
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
예제 #3
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