if domain == 'freq':
db_s11 = 10. * np.log10(np.abs(simulation.gainFrequency))
pha_s11 = np.angle(simulation.gainFrequency)
x = simulation.fAxis
elif domain == 'delay': #only allow delay kernel for simulations for now.
if kernel and meastype == 'simulation':
y = simulation.delay_kernel()
else:
y = simulation.gainDelay
db_s11 = 10. * np.log10(np.abs(y))
pha_s11 = np.angle(y)
x = simulation.tAxis
elif meastype == 'balun':
balunmeas = ABM()
balunmeas.read_files(prefix,
postfix,
prefixb,
postfixb,
filetype,
portmap[0],
portmap[1],
portmap[2],
fMin=fmin,
fMax=fmax,
changeZ=changez,
z0=zi,
z1=zf)
if domain == 'freq':
db_s11 = 10. * np.log10(np.abs(balunmeas.antenna_gain_frequency))
import numpy as np
from gainData import AntennaBalunMeasurement as ABM
from gainData import AntennaDiffMeasurement as ADM
from gainData import GainData as GD
import matplotlib.pyplot as plt
#load measurements
hybrid_coupler_A = ABM()
hybrid_coupler_B = ABM()
cambridge_balun = ABM()
no_balun = ADM()
hybrid_coupler_A.read_files('../Rooftop_Antenna_Measurements_August_10th/',
'_hybrid_A_port_A_antenna_',
'../hybrid_coupler_A_Sierra/A_', '_',
'ANRITSU_CSV', '1', '3', '4', 0.05, 0.25)
hybrid_coupler_B.read_files('../Rooftop_Antenna_Measurements_August_10th/',
'_hybrid_B_port_A_antenna_',
'../hybrid_coupler_B/1', '_', 'ANRITSU_CSV', '1',
'3', '4', 0.05, 0.25)
cambridge_balun.read_files('../Rooftop_Antenna_Measurements_August_10th/',
'_cambridge_balun_N_antenna_0_',
'../Cambridge_Balun_Measurements_N/', '_',
'ANRITSU_CSV', '1', '2', '3', 0.05, 0.25)
no_balun.read_files('../Rooftop_Antenna_Measurements_August_10th/',
'_no_balun_antenna_', 'ANRITSU_CSV', 0.05, 0.25)
simulation = GD()
simulation.read_files(
'../Rooftop_Antenna_Measurements_August_10th/Simulation/simulation_s11_rooftop',
simulation.read_files(
'../August30thSinuousFeedOverDish/Simulation/S11_0.80-30-175_dish-imp-100-band_no-skirt-1.2-0.3-backplane-50-0.99',
'CST_S11',
fMin=0.05,
fMax=0.25)
simulation.gainFrequency = simulation.gainFrequency
db_s11_sim = 10. * np.log10(np.abs(simulation.gainFrequency))
pha_s11_sim = np.angle(simulation.gainFrequency)
fig1 = plt.figure()
fig2 = plt.figure()
ax1 = fig1.add_axes([.1, .1, .8, .8])
ax2 = fig2.add_axes([.1, .1, .8, .8])
for fnum, fname in enumerate(measurement_names):
hybrid_coupler_B = ABM()
hybrid_coupler_B.read_files('../August30thSinuousFeedOverDish/', fname,
'../hybrid_coupler_B/1', '_', 'ANRITSU_CSV',
'1', '3', '4', 0.05, 0.25)
db_s11_B_corr = 10. * np.log10(
np.abs(hybrid_coupler_B.antenna_gain_frequency))
db_s11_B = 10. * np.log10(
np.abs(hybrid_coupler_B.antenna_raw.gainFrequency))
pha_s11_B_corr = np.angle(hybrid_coupler_B.antenna_gain_frequency)
pha_s11_B = np.angle(hybrid_coupler_B.antenna_raw.gainFrequency)
ax1.plot(hybrid_coupler_B.fAxis,
db_s11_B,
color=colors[fnum],
ls='--',
alpha=.75)