Exemplo n.º 1
0
    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,
Exemplo n.º 2
0
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,
Exemplo n.º 3
0

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))
Exemplo n.º 4
0
    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)
Exemplo n.º 5
0
               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,