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)