def maxandmin_intime(): plot(time*1e6, dB(S11c[:, 148]), label="Max refl, yoko=3.15 V".format(yoko[315])) plot(time*1e6, dB(S11c[:, 200]), label="Min refl, yoko=2.21 V".format(yoko[221])) xlabel("Time (us)") ylabel("Reflection (dB)") title("Time domain reflection at 4.46 GHz") ylim(-65, -30) legend()
def maxandmin_intime(): plot(time * 1e6, dB(S11c[:, 148]), label="Max refl, yoko=3.15 V".format(yoko[315])) plot(time * 1e6, dB(S11c[:, 200]), label="Min refl, yoko=2.21 V".format(yoko[221])) xlabel("Time (us)") ylabel("Reflection (dB)") title("Time domain reflection at 4.46 GHz") ylim(-65, -30) legend()
def VNA_twotone_colormesh(): pcolormesh(yoko2, pwr2[:,0]-fridge_att, dB(Magcom[0, :, :])) ylim(-187, -87) xlim(5.0, 6.0) xlabel("Flux (V)") ylabel("Control power (dBm)") title("Reflection (dB) Probe at 4.403 GHz, -137 dBm,\n Control at 4.285 GHz")
def plotmapdBtime(): pcolormesh(yoko, time * 1e6, dB(S11c), vmin=-65, vmax=-30) title("Reflection (dB) vs flux \n and time (1 us pulse) at 4.46 GHz") xlabel("Flux (V)") ylabel("Time (us)") #ylim(0, 1.5) colorbar()
def plotmapdBtime(): pcolormesh(yoko, time*1e6, dB(S11c), vmin=-65, vmax=-30) title("Reflection (dB) vs flux \n and time (1 us pulse) at 4.46 GHz") xlabel("Flux (V)") ylabel("Time (us)") #ylim(0, 1.5) colorbar()
def VNA_twotonesat(): powsat = mean(dB(Magcom[0, :, 19:31]), axis=1) print shape(pwr2) print pwr2[:, 0] print yoko2[0, :] print shape(yoko2) fridge_att = 87.0 + 15.0 * 0.0 pwrlin = 0.001 * 10.0**((pwr2[:, 0] - fridge_att) / 10.0) powsat_theory = lorentzian(flux_parabola(flux_rescale(yoko2[0, :])), 4.285e9, pwrlin, 50.0e6) print shape(powsat_theory) plot(pwr2[:, 0] - fridge_att, normalize(absolute(powsat_theory[:, 56])), label="Theory 50 MHz coupling") plot(pwr2[:, 0] - fridge_att, normalize(powsat), label="Data") #mean(dB(Magcom[0, :, 19:31]), axis=1)) xlabel("Control power (dBm)") ylabel("Reflection (normalized)") title( "Reflection (normalized) Probe at 4.403 GHz, -137 dBm, \n Control at 4.285 GHz" ) legend(loc=3)
def cs_dB_meanyoko_refl(): plot(freq, dB(mean(Magcom[:, :], axis=1))) #plot(freq, dB(Magcom[:, :])) xlabel("Frequency (Hz)") ylabel("Reflection (dB)") title("Reflection") ylim(-50, -20)
def plotdB_colormap(pwi): pcolormesh(yoko, freq, dB(Magcom[:, pwi, :])) xlabel("Flux (V)") ylabel("Frequency (Hz)") caption="{} dBm".format(pwr[pwi]-fridge_attn) title("Reflection fluxmap at {}".format(caption)) ylim(amin(freq), amax(freq)) colorbar() plotdB_colormap.caption=caption
def VNA_twotone_colormesh(): pcolormesh(yoko2, pwr2[:, 0] - fridge_att, dB(Magcom[0, :, :])) ylim(-187, -87) xlim(5.0, 6.0) xlabel("Flux (V)") ylabel("Control power (dBm)") title( "Reflection (dB) Probe at 4.403 GHz, -137 dBm,\n Control at 4.285 GHz" )
def plotdB_colormap_refl(): pcolormesh(yoko, freq, dB(Magcom[:, :])) xlabel("Flux (V)") ylabel("Frequency (Hz)") caption = "{} dBm".format(-30 - fridge_attn) title("Reflection fluxmap at {}".format(caption)) ylim(amin(freq), amax(freq)) colorbar() plotdB_colormap.caption = caption
def time_cuts(): plot(yoko, dB(S11c[35, :]), label="20 ns".format(time[35]-time[30])) plot(yoko, dB(S11c[46, :]), label="64 ns".format(time[46]-time[30])) plot(yoko, dB(S11c[83, :]), label="212 ns".format(time[83]-time[30])) plot(yoko, dB(S11c[90, :]), label="240 ns".format(time[90]-time[30])) plot(yoko, dB(S11c[100, :]), label="280 ns".format(time[100]-time[30])) plot(yoko, dB(S11c[110, :]), label="320 ns".format(time[110]-time[30])) plot(yoko, dB(S11c[153, :]), label="492 ns".format(time[153]-time[30])) plot(yoko, dB(S11c[230, :]), label="800 ns".format(time[230]-time[30])) plot(yoko, dB(S11c[294, :]), label="1056 ns".format(time[294]-time[30])) title("Flux modulation at various times at 4.46 GHz") xlabel("Flux (V)") ylabel("Reflection (dB)") legend()
def time_cuts(): plot(yoko, dB(S11c[35, :]), label="20 ns".format(time[35] - time[30])) plot(yoko, dB(S11c[46, :]), label="64 ns".format(time[46] - time[30])) plot(yoko, dB(S11c[83, :]), label="212 ns".format(time[83] - time[30])) plot(yoko, dB(S11c[90, :]), label="240 ns".format(time[90] - time[30])) plot(yoko, dB(S11c[100, :]), label="280 ns".format(time[100] - time[30])) plot(yoko, dB(S11c[110, :]), label="320 ns".format(time[110] - time[30])) plot(yoko, dB(S11c[153, :]), label="492 ns".format(time[153] - time[30])) plot(yoko, dB(S11c[230, :]), label="800 ns".format(time[230] - time[30])) plot(yoko, dB(S11c[294, :]), label="1056 ns".format(time[294] - time[30])) title("Flux modulation at various times at 4.46 GHz") xlabel("Flux (V)") ylabel("Reflection (dB)") legend()
def VNA_twotonesat(): powsat=mean(dB(Magcom[0, :, 19:31]), axis=1) print shape(pwr2) print pwr2[:,0] print yoko2[0,:] print shape(yoko2) fridge_att=87.0+15.0*0.0 pwrlin=0.001*10.0**((pwr2[:,0]-fridge_att)/10.0) powsat_theory=lorentzian(flux_parabola(flux_rescale(yoko2[0,:])), 4.285e9, pwrlin, 50.0e6) print shape(powsat_theory) plot(pwr2[:,0]-fridge_att, normalize(absolute(powsat_theory[:,56])), label="Theory 50 MHz coupling") plot(pwr2[:,0]-fridge_att, normalize(powsat), label="Data")#mean(dB(Magcom[0, :, 19:31]), axis=1)) xlabel("Control power (dBm)") ylabel("Reflection (normalized)") title("Reflection (normalized) Probe at 4.403 GHz, -137 dBm, \n Control at 4.285 GHz") legend(loc=3)
def cs_dB(fqi, pwi): plot(yoko, dB(Magcom[fqi, pwi, :]))
print shape(Magcom) Magcom=squeeze(Magcom) # print shape(Magcom) # print shape(yoko) # print yoko #print freq #Magcom=mean(Magcom, axis=1) powind=4 print pwr[powind] Magabs=Magcom[:, :, :]-mean(Magcom[:, 197:200, :], axis=1, keepdims=True) fridge_att=87.0+20.0+5.0 pwrlin=0*0.001*10.0**((pwr[powind]-fridge_att)/10.0) if 0: pcolormesh(dB(Magcom[:, :, powind])) show() def cs_refl_lowfrq(fqi, pwi): plot(yoko, normalize(absolute(Magabs[fqi, :, pwi]))) RR=lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.02)), freq[fqi], [pwrlin], g) #RR=1/(1-1j*detuning(yoko*0.195)/(2.0*pi*10.0e6)) plot(yoko, normalize(absolute(RR)), label="50 MHz {}".format(g)) title("Cross section reflection (normalized) \n at {frq:.3f} GHz and {pwr} dBm".format(frq=freq[fqi]/1.0e9, pwr=pwr[pwi])) xlabel("Flux (V)") ylabel("Reflection (normalized)") def cs_refl_lowfrq_dB(fqi, pwi): fqi=103 plot(yoko, normalize(dB(Magabs[fqi, :, pwi]))) RR=lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.07)), freq[fqi], [0*pwrlin], g)
def gatecolormesh(pwi): pcolormesh(yoko, freq, dB(Magcom[:, 0, :])) xlabel("Flux (V)") ylabel("Frequency (Hz)") title("Reflection fluxmap at -117 dBm") colorbar()
def gatecolormesh_noindex(pwi): pcolormesh(dB(Magcom[:, pwi, :]))
def cs_refl_lowfrq_dB(fqi, pwi): fqi=103 plot(yoko, normalize(dB(Magabs[fqi, :, pwi]))) RR=lorentzian(flux_parabola(flux_rescale(yoko, offset=-0.07)), freq[fqi], [0*pwrlin], g) #RR=1/(1-1j*detuning(yoko*0.195)/(2.0*pi*10.0e6)) plot(yoko, normalize(absolute(RR)), label="50 MHz {}".format(g))
def cm_refl_lowfrq_dB(pwi): pcolormesh(yoko, freq, dB(Magcom[:, :, pwi])) ylim(amin(freq), amax(freq))
def cs_dB_meanyoko_trans(): plot(freq, dB(mean(Magcom[:, :], axis=1))) xlabel("Frequency (Hz)") ylabel("Transmission (dB)") title("Transmission") ylim(-60, -20)
def plotdBnoindex(): pcolormesh( dB(Magcom[:, 0, :]))
def plotdBnoindex(): pcolormesh(dB(Magcom[:, :]))