def formula_comparison(pl=None, **kwargs):
a = IDT()
frq = linspace(3e9, 7e9, 10000)
print a._get_coupling(frq)
print a.max_coupling
pl, pf = line(frq,
a._get_coupling(frq) / a.max_coupling,
plotter=pl,
linewidth=1.0,
**kwargs)
line(frq,
a._get_Lamb_shift(frq) / a.max_coupling,
plotter=pl,
color="red",
linewidth=1.0)
#line(frq, coup, color="purple", plotter=pl, linewidth=0.3)
line(frq,
a._get_full_coupling(frq) / a.max_coupling,
plotter=pl,
color="green",
linewidth=0.3)
line(frq,
a._get_full_Lamb_shift(frq) / a.max_coupling,
plotter=pl,
color="black",
linewidth=0.3)
return pl
def formula_comparison(pl=None, **kwargs):
a=IDT()
frq=linspace(3e9, 7e9, 10000)
print a._get_coupling(frq)
print a.max_coupling
pl, pf=line(frq, a._get_coupling(frq)/a.max_coupling, plotter=pl, linewidth=1.0, **kwargs)
line(frq, a._get_Lamb_shift(frq)/a.max_coupling, plotter=pl, color="red", linewidth=1.0)
#line(frq, coup, color="purple", plotter=pl, linewidth=0.3)
line(frq, a._get_full_coupling(frq)/a.max_coupling, plotter=pl, color="green", linewidth=0.3)
line(frq, a._get_full_Lamb_shift(frq)/a.max_coupling, plotter=pl, color="black", linewidth=0.3)
return pl
X=a._get_X(f=frq)
Np=a.Np
f0=a.f0
#coup=(sqrt(2)*cos(pi*frq/(4*f0))*(1.0/Np)*sin(X)/sin(X/Np))**2
coup=(sin(X)/X)**2
#coup=(1.0/Np*sin(X)/sin(X/Np))**2
element_factor_plot(a)
pl=Plotter()
line(frq, a._get_coupling(frq)/a.max_coupling, plotter=pl, linewidth=1.0)
line(frq, a._get_Lamb_shift(frq)/a.max_coupling, plotter=pl, color="red", linewidth=1.0)
line(frq, coup, color="purple", plotter=pl, linewidth=0.3)
line(frq, a._get_full_coupling(frq)/a.max_coupling/1.247**2/2, plotter=pl, color="green", linewidth=0.3)
line(frq, a._get_full_Lamb_shift(frq)/a.max_coupling/1.247**2/2, plotter=pl, color="black", linewidth=0.3)
#hb=hilbert(coup) #a._get_coupling(frq))
#line(frq, real(hb), plotter=pl, color="green", linewidth=0.3)
#line(frq, imag(hb), plotter=pl, color="black", linewidth=0.3)
#Baa= (1+cos(X/(4*Np)))*(1.0/Np)**2*2*(Np*sin(2*X/Np)-sin(2*X))/(2*(1-cos(2*X/Np)))
Baa=-(sin(2.0*X)-2.0*X)/(2.0*X**2)
#Baa= (1.0/Np)**2*2*(Np*sin(2*X/Np)-sin(2*X))/(2*(1-cos(2*X/Np)))
line(frq, Baa, plotter=pl, color="cyan", linewidth=0.3)
pl.show()
b=IDT(ft="single")
a.ft_mult=5
print a.mu_mult, a.ft_mult, b.mu_mult, b.ft_mult
#coup=(sqrt(2)*cos(pi*frq/(4*f0))*(1.0/Np)*sin(X)/sin(X/Np))**2
coup = (sin(X) / X)**2
#coup=(1.0/Np*sin(X)/sin(X/Np))**2
element_factor_plot(a)
pl = Plotter()
line(frq, a._get_coupling(frq) / a.max_coupling, plotter=pl, linewidth=1.0)
line(frq,
a._get_Lamb_shift(frq) / a.max_coupling,
plotter=pl,
color="red",
linewidth=1.0)
line(frq, coup, color="purple", plotter=pl, linewidth=0.3)
line(frq,
a._get_full_coupling(frq) / a.max_coupling / 1.247**2 / 2,
plotter=pl,
color="green",
linewidth=0.3)
line(frq,
a._get_full_Lamb_shift(frq) / a.max_coupling / 1.247**2 / 2,
plotter=pl,
color="black",
linewidth=0.3)
#hb=hilbert(coup) #a._get_coupling(frq))
#line(frq, real(hb), plotter=pl, color="green", linewidth=0.3)
#line(frq, imag(hb), plotter=pl, color="black", linewidth=0.3)
#Baa= (1+cos(X/(4*Np)))*(1.0/Np)**2*2*(Np*sin(2*X/Np)-sin(2*X))/(2*(1-cos(2*X/Np)))
Baa = -(sin(2.0 * X) - 2.0 * X) / (2.0 * X**2)
#Baa= (1.0/Np)**2*2*(Np*sin(2*X/Np)-sin(2*X))/(2*(1-cos(2*X/Np)))