pl88 = scatter(d0527.frequency[::nskip] / 1e9,
20.0 * log10(magabs88[::nskip]) -
bg_A4(d0527.frequency[::nskip]),
facecolor="blue",
edgecolor="blue",
pl=pl88,
nrows=2,
ncols=3,
nplot=1,
fig_width=fig_width,
fig_height=fig_height,
marker_size=10)
(S11, S12, S13, S21, S22, S23, S31, S32,
S33) = idt88._get_simple_S(f=d0527.frequency)
S13xS31 = S13 * S31
print idt88.Np, idt88.K2
print idt88.f0
print d0527.comment
print -d0527.fridge_atten + d0527.fridge_gain - d0527.rt_atten + d0527.rt_gain - 10
line(
d0527.frequency / 1e9,
20 * log10(absolute(S13xS31)) - 4,
color="red",
pl=pl88,
auto_ylim=False,
y_min=-40,
y_max=0,
auto_xlim=False,
if __name__ == "__main__":
pl = line(a.frequency, 10 * log10(absolute(a.MagcomData[:, 0]))) #.show()
#pl=line(a.frequency, 10*log10(absolute(a.MagcomData[:, 1])))#.show()
magfilt = MagcomFilt(a)
magabs = absolute(magfilt)
pl1 = line(a.frequency, magabs)
pl = line(a.frequency, 10.0 * log10(magabs), color="red", pl=pl)
#idt.Np=56
#idt.f0=4.46e9 #4.452
idt.K2 = 0.039
#idt.f0=4.437e9
(S11, S12, S13, S21, S22, S23, S31, S32,
S33) = idt._get_simple_S(f=a.frequency)
print idt.f0
print a.comment
print -a.fridge_atten + a.fridge_gain - a.rt_atten + a.rt_gain - 10
line(a.frequency,
10 * log10(absolute(S33)) - 10,
color="green",
pl=pl,
auto_ylim=False,
y_min=-30,
y_max=-8,
auto_xlim=False,
x_min=4e9,
x_max=5e9,
if __name__=="__main__":
nskip=50
d0527.read_data()
magfilt88=MagcomFilt88(d0527)
magabs88=absolute(magfilt88)
#line(d0527.frequency/1e9, magabs88)
pl88=scatter(d0527.frequency[::nskip]/1e9, 20.0*log10(magabs88[::nskip])-bg_A4(d0527.frequency[::nskip]),
facecolor="blue", edgecolor="blue", pl=pl88, nrows=2, ncols=3, nplot=1,
fig_width=fig_width, fig_height=fig_height, marker_size=10)
(S11, S12, S13,
S21, S22, S23,
S31, S32, S33)=idt88._get_simple_S(f=d0527.frequency)
S13xS31=S13*S31
print idt88.Np, idt88.K2
print idt88.f0
print d0527.comment
print -d0527.fridge_atten+d0527.fridge_gain-d0527.rt_atten+d0527.rt_gain-10
line(d0527.frequency/1e9, 20*log10(absolute(S13xS31))-4, color="red", pl=pl88,
auto_ylim=False, y_min=-40, y_max=0,
auto_xlim=False, x_min=4.2, x_max=4.7, xlabel="Frequency (GHz)",
ylabel="Transmission (dB)", linewidth=1.0,
)#.show()
pl88.axes.set_xticks(linspace(4.2, 4.7, 3))
pl88.axes.set_yticks(linspace(-40, 0, 5))
magfilt=MagcomFilt(a)
magabs=absolute(magfilt)
line(a.frequency, magabs)
nskip=50
pl=scatter(a.frequency[::nskip], 20.0*log10(magabs[::nskip])-bg(a.frequency[::nskip]), facecolor="red", edgecolor="red", pl=pl)
#idt.Np=56
#idt.f0=4.46e9 #4.452
#idt.K2=0.032
idt.Np=36.5
#idt.f0=4.452e9
(S11, S12, S13,
S21, S22, S23,
S31, S32, S33)=idt._get_simple_S(f=a.frequency)
print idt.Np, idt.K2
print idt.f0
print a.comment
print -a.fridge_atten+a.fridge_gain-a.rt_atten+a.rt_gain-10
line(a.frequency, 20*log10(absolute(S13*S31))-4*1, color="green", pl=pl,
auto_ylim=False, y_min=-40, y_max=-10*0,
auto_xlim=False, x_min=4e9, x_max=5e9, xlabel="Frequency (Hz)", ylabel="Transmission (dB)",
title="Pickup IDT (37 fingers)")#.show()
line(a.frequency, angle(magfilt))
X=idt.Np*pi*(a.frequency-idt.f0)/idt.f0
#line(a.frequency, 20*log10(0.5*(sin(X)/(X))**2), pl=pl)
facecolor="red", edgecolor="red", pl=pl, nrows=2, ncols=2, nplot=1,
fig_width=3.5, fig_height=3.5)
plbg=scatter(d0527.frequency[::nskip]/1e9, 20.0*log10(magabs[::nskip])-bg_A4(d0527.frequency[::nskip]),
facecolor="red", edgecolor="red", pl=plbg)
#idt.Np=56
#idt.f0=4.46e9 #4.452
#idt.K2=0.032
#idt.Np=36.5
#idt.f0=4.452e9
(S11, S12, S13,
S21, S22, S23,
S31, S32, S33)=idt._get_simple_S(f=d0527.frequency)
print idt.Np, idt.K2
print idt.f0
print d0527.comment
print -d0527.fridge_atten+d0527.fridge_gain-d0527.rt_atten+d0527.rt_gain-10
line(d0527.frequency/1e9, 20*log10(absolute(S13*S31))-4, color="green", pl=pl,
auto_ylim=False, y_min=-40, y_max=0,
auto_xlim=False, x_min=4.2, x_max=4.7, xlabel="Frequency (GHz)",
ylabel="Transmission (dB)",
)#.show()
pl.axes.set_xticks(linspace(4.2, 4.7, 3))
pl.axes.set_yticks(linspace(-40, 0, 5))
line(d0527.frequency/1e9, 20*log10(absolute(S13*S31))-4, color="green", pl=plbg,