def calc_and_massage_a_th_Bode(f, th_bode, a_bode, \
**kwargs):
a_th_bode = a_bode/th_bode
a_th_bode.seedfreq = 10.0
a_th_bode.seedphase = 0.0
a_th_bode.PhaseMassage(f)
a_th_bode = SFLR_TMM.a_massage(a_th_bode, f, **kwargs)
return a_th_bode
def calc_and_massage_Bodes(f, params, func=mod.Bodes, \
**kwargs):
s = 2.0j*pi*f
th_comp, a_comp = func(s, params)
a_bode = rwkbode.rwkbode(compin=a_comp)
th_bode = rwkbode.rwkbode(compin=th_comp)
a_bode.seedfreq = 2.0
a_bode.seedphase = 30.0
a_bode.PhaseMassage(f)
a_bode = SFLR_TMM.a_massage(a_bode, f, **kwargs)
th_bode.seedfreq = 1.2
th_bode.seedphase = -90.0
th_bode.PhaseMassage(f)
return th_bode, a_bode
