# C.append(Cev)
C.append(Ceb)
T_InputNames = sigExcList
T_OutputNames = sigFbList
# Compute Gain, Phase, Crit Distance
gain_mag = []
phase_deg = []
rCritNom_mag = []
rCrit_mag = []
sNom = []
for iSeg in range(0, len(oDataSegs)):
gainElem_mag, _, _ = FreqTrans.DistCritCirc(T[iSeg], pCrit = 0+0j, typeNorm = 'RSS')
gain_mag.append(gainElem_mag)
phase_deg.append(FreqTrans.Phase(T[iSeg], phaseUnit = 'deg', unwrap = True))
nom_mag, _, _= FreqTrans.DistCritCirc(T[iSeg], pCrit = -1+0j, typeNorm = 'RSS')
rCritNom_mag.append(nom_mag)
sNom_seg, _ = FreqTrans.Sigma(T[iSeg]) # Singular Value Decomp, U @ S @ Vh == T[...,i]
sNom.append(sNom_seg)
#%% Sigma Plot
fig = None
T_InputNames = sigExcList
T_OutputNames = sigFbList
# Compute Gain, Phase, Crit Distance
gain_mag = []
gainUnc_mag = []
phase_deg = []
rCritNom_mag = []
rCritUnc_mag = []
rCrit_mag = []
sNom = []
sUnc = []
for iSeg in range(0, len(oDataSegs)):
gainElem_mag, gainElemUnc_mag, gainElemDiff_mag = FreqTrans.DistCritCirc(
T[iSeg], TUnc[iSeg], pCrit=0 + 0j, typeNorm='RSS')
gain_mag.append(gainElem_mag)
gainUnc_mag.append(gainElemUnc_mag)
phase_deg.append(FreqTrans.Phase(T[iSeg], phaseUnit='deg', unwrap=True))
nom_mag, unc_mag, diff_mag = FreqTrans.DistCritCirc(T[iSeg],
TUnc[iSeg],
pCrit=-1 + 0j,
typeNorm='RSS')
rCritNom_mag.append(nom_mag)
rCritUnc_mag.append(unc_mag)
rCrit_mag.append(diff_mag)
# Compute Gain, Phase, Crit Distance
gainLiEstNomList_mag = []
phaseLiEstNomList_deg = []
rCritLiEstNomList_mag = []
for iSeg in range(0, len(oDataSegs)):
gain_mag, phase_deg = FreqTrans.GainPhase(LiEstNomList[iSeg],
magUnit='mag',
phaseUnit='deg',
unwrap=True)
gainLiEstNomList_mag.append(gain_mag)
phaseLiEstNomList_deg.append(phase_deg)
# rCritLiEstNom_mag, _, _ = FreqTrans.DistCrit(LiEstNomList[iSeg], typeUnc = 'ellipse')
rCritLiEstNom_mag, _, _ = FreqTrans.DistCritCirc(LiEstNomList[iSeg])
rCritLiEstNomList_mag.append(rCritLiEstNom_mag)
#%% Sigma Plot
fig = None
for iSeg in range(0, len(oDataSegs)):
Cmin = np.min(np.min(LiEstCohList[iSeg], axis=0), axis=0)
sNomMin = np.min(svLiEstNomList[iSeg], axis=0)
fig = FreqTrans.PlotSigma(freq_hz[0],
svLiEstNomList[iSeg],
coher_nd=Cmin,
fig=fig,
color=rtsmSegList[iSeg]['color'],
linestyle='-',
