def updateLabels(self, s11data: List[Datapoint], s21data: List[Datapoint]):
from NanoVNASaver.Chart import PhaseChart
from NanoVNASaver.NanoVNASaver import NanoVNASaver
if self.location != -1:
im50, re50, vswr = NanoVNASaver.vswr(s11data[self.location])
if im50 < 0:
im50str = " - j" + str(round(-1 * im50, 3))
else:
im50str = " + j" + str(round(im50, 3))
self.frequency_label.setText(
NanoVNASaver.formatFrequency(s11data[self.location].freq))
self.impedance_label.setText(str(round(re50, 3)) + im50str)
self.returnloss_label.setText(
str(round(20 * math.log10((vswr - 1) / (vswr + 1)), 3)) +
" dB")
reactance = NanoVNASaver.reactanceEquivalent(
im50, s11data[self.location].freq)
self.reactance_label.setText(reactance)
self.vswr_label.setText(str(round(vswr, 3)))
if len(s21data) == len(s11data):
_, _, vswr = NanoVNASaver.vswr(s21data[self.location])
self.gain_label.setText(
str(round(20 * math.log10((vswr - 1) / (vswr + 1)), 3)) +
" dB")
self.phase_label.setText(
str(round(PhaseChart.angle(s21data[self.location]), 2)) +
"\N{DEGREE SIGN}")
def updateLabels(self, s11data: List[Datapoint], s21data: List[Datapoint]):
from NanoVNASaver.Chart import PhaseChart
from NanoVNASaver.NanoVNASaver import NanoVNASaver
if self.location != -1:
im50, re50, vswr = NanoVNASaver.vswr(s11data[self.location])
rp = (re50**2 + im50**2) / re50
xp = (re50**2 + im50**2) / im50
re50 = round(re50, 4 - max(0, math.floor(math.log10(abs(re50)))))
rp = round(rp, 4 - max(0, math.floor(math.log10(abs(rp)))))
im50 = round(im50, 4 - max(0, math.floor(math.log10(abs(im50)))))
xp = round(xp, 4 - max(0, math.floor(math.log10(abs(xp)))))
if im50 < 0:
im50str = " -j" + str(-1 * im50)
else:
im50str = " +j" + str(im50)
im50str += "\N{OHM SIGN}"
if xp < 0:
xpstr = NanoVNASaver.capacitanceEquivalent(
xp, s11data[self.location].freq)
else:
xpstr = NanoVNASaver.inductanceEquivalent(
xp, s11data[self.location].freq)
self.frequency_label.setText(
NanoVNASaver.formatFrequency(s11data[self.location].freq))
self.impedance_label.setText(str(re50) + im50str)
self.parallel_r_label.setText(str(rp) + "\N{OHM SIGN}")
self.parallel_x_label.setText(xpstr)
self.returnloss_label.setText(
str(round(20 * math.log10((vswr - 1) / (vswr + 1)), 3)) +
" dB")
capacitance = NanoVNASaver.capacitanceEquivalent(
im50, s11data[self.location].freq)
inductance = NanoVNASaver.inductanceEquivalent(
im50, s11data[self.location].freq)
self.inductance_label.setText(inductance)
self.capacitance_label.setText(capacitance)
vswr = round(vswr, 3)
if vswr < 0:
vswr = "-"
self.vswr_label.setText(str(vswr))
self.quality_factor_label.setText(
str(
round(NanoVNASaver.qualifyFactor(s11data[self.location]),
1)))
self.s11_phase_label.setText(
str(round(-PhaseChart.angle(s11data[self.location]), 2)) +
"\N{DEGREE SIGN}")
if len(s21data) == len(s11data):
_, _, vswr = NanoVNASaver.vswr(s21data[self.location])
self.gain_label.setText(
str(round(20 * math.log10((vswr - 1) / (vswr + 1)), 3)) +
" dB")
self.s21_phase_label.setText(
str(round(-PhaseChart.angle(s21data[self.location]), 2)) +
"\N{DEGREE SIGN}")
def runAnalysis(self):
from NanoVNASaver.NanoVNASaver import NanoVNASaver
self.reset()
pass_band_location = self.app.markers[0].location
logger.debug("Pass band location: %d", pass_band_location)
if len(self.app.data21) == 0:
logger.debug("No data to analyse")
self.result_label.setText("No data to analyse.")
return
if pass_band_location < 0:
logger.debug("No location for %s", self.app.markers[0].name)
self.result_label.setText("Please place " + self.app.markers[0].name + " in the passband.")
return
pass_band_db = NanoVNASaver.gain(self.app.data21[pass_band_location])
logger.debug("Initial passband gain: %d", pass_band_db)
initial_cutoff_location = -1
for i in range(pass_band_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found a cutoff location
initial_cutoff_location = i
break
if initial_cutoff_location < 0:
self.result_label.setText("Cutoff location not found.")
return
initial_cutoff_frequency = self.app.data21[initial_cutoff_location].freq
logger.debug("Found initial cutoff frequency at %d", initial_cutoff_frequency)
peak_location = -1
peak_db = NanoVNASaver.gain(self.app.data21[initial_cutoff_location])
for i in range(0, initial_cutoff_location):
db = NanoVNASaver.gain(self.app.data21[i])
if db > peak_db:
peak_db = db
peak_location = i
logger.debug("Found peak of %f at %d", peak_db, self.app.data[peak_location].freq)
self.app.markers[0].setFrequency(str(self.app.data21[peak_location].freq))
self.app.markers[0].frequencyInput.setText(str(self.app.data21[peak_location].freq))
cutoff_location = -1
pass_band_db = peak_db
for i in range(peak_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found the cutoff location
cutoff_location = i
break
cutoff_frequency = self.app.data21[cutoff_location].freq
cutoff_gain = NanoVNASaver.gain(self.app.data21[cutoff_location]) - pass_band_db
if cutoff_gain < -4:
logger.debug("Cutoff frequency found at %f dB - insufficient data points for true -3 dB point.",
cutoff_gain)
logger.debug("Found true cutoff frequency at %d", cutoff_frequency)
self.cutoff_label.setText(NanoVNASaver.formatFrequency(cutoff_frequency) +
" (" + str(round(cutoff_gain, 1)) + " dB)")
self.app.markers[1].setFrequency(str(cutoff_frequency))
self.app.markers[1].frequencyInput.setText(str(cutoff_frequency))
six_db_location = -1
for i in range(cutoff_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 6:
# We found 6dB location
six_db_location = i
break
if six_db_location < 0:
self.result_label.setText("6 dB location not found.")
return
six_db_cutoff_frequency = self.app.data21[six_db_location].freq
self.six_db_label.setText(NanoVNASaver.formatFrequency(six_db_cutoff_frequency))
ten_db_location = -1
for i in range(cutoff_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 10:
# We found 6dB location
ten_db_location = i
break
twenty_db_location = -1
for i in range(cutoff_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 20:
# We found 6dB location
twenty_db_location = i
break
sixty_db_location = -1
for i in range(six_db_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 60:
# We found 60dB location! Wow.
sixty_db_location = i
break
if sixty_db_location > 0:
sixty_db_cutoff_frequency = self.app.data21[sixty_db_location].freq
self.sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency))
else:
# # We derive 60 dB instead
# factor = 10 * (-54 / decade_attenuation)
# sixty_db_cutoff_frequency = round(six_db_cutoff_frequency + six_db_cutoff_frequency * factor)
# self.sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency) + " (derived)")
self.sixty_db_label.setText("Not calculated")
if ten_db_location > 0 and twenty_db_location > 0:
octave_attenuation, decade_attenuation = self.calculateRolloff(ten_db_location, twenty_db_location)
self.db_per_octave_label.setText(str(round(octave_attenuation, 3)) + " dB / octave")
self.db_per_decade_label.setText(str(round(decade_attenuation, 3)) + " dB / decade")
else:
self.db_per_octave_label.setText("Not calculated")
self.db_per_decade_label.setText("Not calculated")
self.result_label.setText("Analysis complete (" + str(len(self.app.data)) + " points)")
def runAnalysis(self):
from NanoVNASaver.NanoVNASaver import NanoVNASaver
self.reset()
if len(self.app.data21) == 0:
logger.debug("No data to analyse")
self.result_label.setText("No data to analyse.")
return
peak_location = -1
peak_db = NanoVNASaver.gain(self.app.data21[0])
for i in range(len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if db > peak_db:
peak_db = db
peak_location = i
logger.debug("Found peak of %f at %d", peak_db, self.app.data[peak_location].freq)
lower_cutoff_location = -1
pass_band_db = peak_db
for i in range(len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found the cutoff location
lower_cutoff_location = i
break
lower_cutoff_frequency = self.app.data21[lower_cutoff_location].freq
lower_cutoff_gain = NanoVNASaver.gain(self.app.data21[lower_cutoff_location]) - pass_band_db
if lower_cutoff_gain < -4:
logger.debug("Lower cutoff frequency found at %f dB - insufficient data points for true -3 dB point.",
lower_cutoff_gain)
logger.debug("Found true lower cutoff frequency at %d", lower_cutoff_frequency)
self.lower_cutoff_label.setText(NanoVNASaver.formatFrequency(lower_cutoff_frequency) +
" (" + str(round(lower_cutoff_gain, 1)) + " dB)")
self.app.markers[1].setFrequency(str(lower_cutoff_frequency))
self.app.markers[1].frequencyInput.setText(str(lower_cutoff_frequency))
upper_cutoff_location = -1
for i in range(len(self.app.data21)-1, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found the cutoff location
upper_cutoff_location = i
break
upper_cutoff_frequency = self.app.data21[upper_cutoff_location].freq
upper_cutoff_gain = NanoVNASaver.gain(self.app.data21[upper_cutoff_location]) - pass_band_db
if upper_cutoff_gain < -4:
logger.debug("Upper cutoff frequency found at %f dB - insufficient data points for true -3 dB point.",
upper_cutoff_gain)
logger.debug("Found true upper cutoff frequency at %d", upper_cutoff_frequency)
self.upper_cutoff_label.setText(NanoVNASaver.formatFrequency(upper_cutoff_frequency) +
" (" + str(round(upper_cutoff_gain, 1)) + " dB)")
self.app.markers[2].setFrequency(str(upper_cutoff_frequency))
self.app.markers[2].frequencyInput.setText(str(upper_cutoff_frequency))
span = upper_cutoff_frequency - lower_cutoff_frequency
center_frequency = math.sqrt(lower_cutoff_frequency * upper_cutoff_frequency)
q = center_frequency / span
self.span_label.setText(NanoVNASaver.formatFrequency(span))
self.center_frequency_label.setText(NanoVNASaver.formatFrequency(center_frequency))
self.quality_label.setText(str(round(q, 2)))
self.app.markers[0].setFrequency(str(round(center_frequency)))
self.app.markers[0].frequencyInput.setText(str(round(center_frequency)))
# Lower roll-off
lower_six_db_location = -1
for i in range(lower_cutoff_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 6:
# We found 6dB location
lower_six_db_location = i
break
if lower_six_db_location < 0:
self.result_label.setText("Lower 6 dB location not found.")
return
lower_six_db_cutoff_frequency = self.app.data21[lower_six_db_location].freq
self.lower_six_db_label.setText(NanoVNASaver.formatFrequency(lower_six_db_cutoff_frequency))
ten_db_location = -1
for i in range(lower_cutoff_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 10:
# We found 6dB location
ten_db_location = i
break
twenty_db_location = -1
for i in range(lower_cutoff_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 20:
# We found 6dB location
twenty_db_location = i
break
sixty_db_location = -1
for i in range(lower_six_db_location, len(self.app.data21)):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 60:
# We found 60dB location! Wow.
sixty_db_location = i
break
if sixty_db_location > 0:
sixty_db_cutoff_frequency = self.app.data21[sixty_db_location].freq
self.lower_sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency))
else:
# # We derive 60 dB instead
# factor = 10 * (-54 / decade_attenuation)
# sixty_db_cutoff_frequency = round(six_db_cutoff_frequency + six_db_cutoff_frequency * factor)
# self.sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency) + " (derived)")
self.lower_sixty_db_label.setText("Not calculated")
if ten_db_location > 0 and twenty_db_location > 0:
octave_attenuation, decade_attenuation = self.calculateRolloff(ten_db_location, twenty_db_location)
self.lower_db_per_octave_label.setText(str(round(octave_attenuation, 3)) + " dB / octave")
self.lower_db_per_decade_label.setText(str(round(decade_attenuation, 3)) + " dB / decade")
else:
self.lower_db_per_octave_label.setText("Not calculated")
self.lower_db_per_decade_label.setText("Not calculated")
# Upper roll-off
upper_six_db_location = -1
for i in range(upper_cutoff_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 6:
# We found 6dB location
upper_six_db_location = i
break
if upper_six_db_location < 0:
self.result_label.setText("Upper 6 dB location not found.")
return
upper_six_db_cutoff_frequency = self.app.data21[upper_six_db_location].freq
self.upper_six_db_label.setText(NanoVNASaver.formatFrequency(upper_six_db_cutoff_frequency))
six_db_span = upper_six_db_cutoff_frequency - lower_six_db_cutoff_frequency
self.six_db_span_label.setText(NanoVNASaver.formatFrequency(six_db_span))
ten_db_location = -1
for i in range(upper_cutoff_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 10:
# We found 6dB location
ten_db_location = i
break
twenty_db_location = -1
for i in range(upper_cutoff_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 20:
# We found 6dB location
twenty_db_location = i
break
sixty_db_location = -1
for i in range(upper_six_db_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 60:
# We found 60dB location! Wow.
sixty_db_location = i
break
if sixty_db_location > 0:
sixty_db_cutoff_frequency = self.app.data21[sixty_db_location].freq
self.upper_sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency))
else:
# # We derive 60 dB instead
# factor = 10 * (-54 / decade_attenuation)
# sixty_db_cutoff_frequency = round(six_db_cutoff_frequency + six_db_cutoff_frequency * factor)
# self.sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency) + " (derived)")
self.upper_sixty_db_label.setText("Not calculated")
if ten_db_location > 0 and twenty_db_location > 0:
octave_attenuation, decade_attenuation = self.calculateRolloff(ten_db_location, twenty_db_location)
self.upper_db_per_octave_label.setText(str(round(octave_attenuation, 3)) + " dB / octave")
self.upper_db_per_decade_label.setText(str(round(decade_attenuation, 3)) + " dB / decade")
else:
self.upper_db_per_octave_label.setText("Not calculated")
self.upper_db_per_decade_label.setText("Not calculated")
if upper_cutoff_gain < -4 or lower_cutoff_gain < -4:
self.result_label.setText("Analysis complete (" + str(len(self.app.data)) + " points)\n" +
"Insufficient data for analysis. Increase segment count.")
else:
self.result_label.setText("Analysis complete (" + str(len(self.app.data)) + " points)")
def runAnalysis(self):
from NanoVNASaver.NanoVNASaver import NanoVNASaver
self.reset()
pass_band_location = self.app.markers[0].location
logger.debug("Pass band location: %d", pass_band_location)
if len(self.app.data21) == 0:
logger.debug("No data to analyse")
self.result_label.setText("No data to analyse.")
return
if pass_band_location < 0:
logger.debug("No location for %s", self.app.markers[0].name)
self.result_label.setText("Please place " +
self.app.markers[0].name +
" in the passband.")
return
pass_band_db = NanoVNASaver.gain(self.app.data21[pass_band_location])
logger.debug("Initial passband gain: %d", pass_band_db)
initial_lower_cutoff_location = -1
for i in range(pass_band_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found a cutoff location
initial_lower_cutoff_location = i
break
if initial_lower_cutoff_location < 0:
self.result_label.setText("Lower cutoff location not found.")
return
initial_lower_cutoff_frequency = self.app.data21[
initial_lower_cutoff_location].freq
logger.debug("Found initial lower cutoff frequency at %d",
initial_lower_cutoff_frequency)
initial_upper_cutoff_location = -1
for i in range(pass_band_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found a cutoff location
initial_upper_cutoff_location = i
break
if initial_upper_cutoff_location < 0:
self.result_label.setText("Upper cutoff location not found.")
return
initial_upper_cutoff_frequency = self.app.data21[
initial_upper_cutoff_location].freq
logger.debug("Found initial upper cutoff frequency at %d",
initial_upper_cutoff_frequency)
peak_location = -1
peak_db = NanoVNASaver.gain(
self.app.data21[initial_lower_cutoff_location])
for i in range(initial_lower_cutoff_location,
initial_upper_cutoff_location, 1):
db = NanoVNASaver.gain(self.app.data21[i])
if db > peak_db:
peak_db = db
peak_location = i
logger.debug("Found peak of %f at %d", peak_db,
self.app.data[peak_location].freq)
lower_cutoff_location = -1
pass_band_db = peak_db
for i in range(peak_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found the cutoff location
lower_cutoff_location = i
break
lower_cutoff_frequency = self.app.data21[lower_cutoff_location].freq
lower_cutoff_gain = NanoVNASaver.gain(
self.app.data21[lower_cutoff_location]) - pass_band_db
if lower_cutoff_gain < -4:
logger.debug(
"Lower cutoff frequency found at %f dB - insufficient data points for true -3 dB point.",
lower_cutoff_gain)
logger.debug("Found true lower cutoff frequency at %d",
lower_cutoff_frequency)
self.lower_cutoff_label.setText(
NanoVNASaver.formatFrequency(lower_cutoff_frequency) + " (" +
str(round(lower_cutoff_gain, 1)) + " dB)")
self.app.markers[1].setFrequency(str(lower_cutoff_frequency))
self.app.markers[1].frequencyInput.setText(str(lower_cutoff_frequency))
upper_cutoff_location = -1
pass_band_db = peak_db
for i in range(peak_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found the cutoff location
upper_cutoff_location = i
break
upper_cutoff_frequency = self.app.data21[upper_cutoff_location].freq
upper_cutoff_gain = NanoVNASaver.gain(
self.app.data21[upper_cutoff_location]) - pass_band_db
if upper_cutoff_gain < -4:
logger.debug(
"Upper cutoff frequency found at %f dB - insufficient data points for true -3 dB point.",
upper_cutoff_gain)
logger.debug("Found true upper cutoff frequency at %d",
upper_cutoff_frequency)
self.upper_cutoff_label.setText(
NanoVNASaver.formatFrequency(upper_cutoff_frequency) + " (" +
str(round(upper_cutoff_gain, 1)) + " dB)")
self.app.markers[2].setFrequency(str(upper_cutoff_frequency))
self.app.markers[2].frequencyInput.setText(str(upper_cutoff_frequency))
span = upper_cutoff_frequency - lower_cutoff_frequency
center_frequency = math.sqrt(lower_cutoff_frequency *
upper_cutoff_frequency)
q = center_frequency / span
self.span_label.setText(NanoVNASaver.formatFrequency(span))
self.center_frequency_label.setText(
NanoVNASaver.formatFrequency(center_frequency))
self.quality_label.setText(str(round(q, 2)))
self.app.markers[0].setFrequency(str(round(center_frequency)))
self.app.markers[0].frequencyInput.setText(str(
round(center_frequency)))
# Lower roll-off
lower_six_db_location = -1
for i in range(lower_cutoff_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 6:
# We found 6dB location
lower_six_db_location = i
break
if lower_six_db_location < 0:
self.result_label.setText("Lower 6 dB location not found.")
return
lower_six_db_cutoff_frequency = self.app.data21[
lower_six_db_location].freq
self.lower_six_db_label.setText(
NanoVNASaver.formatFrequency(lower_six_db_cutoff_frequency))
ten_db_location = -1
for i in range(lower_cutoff_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 10:
# We found 6dB location
ten_db_location = i
break
twenty_db_location = -1
for i in range(lower_cutoff_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 20:
# We found 6dB location
twenty_db_location = i
break
sixty_db_location = -1
for i in range(lower_six_db_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 60:
# We found 60dB location! Wow.
sixty_db_location = i
break
if sixty_db_location > 0:
if sixty_db_location > 0:
sixty_db_cutoff_frequency = self.app.data21[
sixty_db_location].freq
self.lower_sixty_db_label.setText(
NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency))
elif ten_db_location != -1 and twenty_db_location != -1:
ten = self.app.data21[ten_db_location].freq
twenty = self.app.data21[twenty_db_location].freq
sixty_db_frequency = ten * 10**(
5 * (math.log10(twenty) - math.log10(ten)))
self.lower_sixty_db_label.setText(
NanoVNASaver.formatFrequency(sixty_db_frequency) +
" (derived)")
else:
self.lower_sixty_db_label.setText("Not calculated")
if ten_db_location > 0 and twenty_db_location > 0 and ten_db_location != twenty_db_location:
octave_attenuation, decade_attenuation = self.calculateRolloff(
ten_db_location, twenty_db_location)
self.lower_db_per_octave_label.setText(
str(round(octave_attenuation, 3)) + " dB / octave")
self.lower_db_per_decade_label.setText(
str(round(decade_attenuation, 3)) + " dB / decade")
else:
self.lower_db_per_octave_label.setText("Not calculated")
self.lower_db_per_decade_label.setText("Not calculated")
# Upper roll-off
upper_six_db_location = -1
for i in range(upper_cutoff_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 6:
# We found 6dB location
upper_six_db_location = i
break
if upper_six_db_location < 0:
self.result_label.setText("Upper 6 dB location not found.")
return
upper_six_db_cutoff_frequency = self.app.data21[
upper_six_db_location].freq
self.upper_six_db_label.setText(
NanoVNASaver.formatFrequency(upper_six_db_cutoff_frequency))
six_db_span = upper_six_db_cutoff_frequency - lower_six_db_cutoff_frequency
self.six_db_span_label.setText(
NanoVNASaver.formatFrequency(six_db_span))
ten_db_location = -1
for i in range(upper_cutoff_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 10:
# We found 6dB location
ten_db_location = i
break
twenty_db_location = -1
for i in range(upper_cutoff_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 20:
# We found 6dB location
twenty_db_location = i
break
sixty_db_location = -1
for i in range(upper_six_db_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 60:
# We found 60dB location! Wow.
sixty_db_location = i
break
if sixty_db_location > 0:
sixty_db_cutoff_frequency = self.app.data21[sixty_db_location].freq
self.upper_sixty_db_label.setText(
NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency))
elif ten_db_location != -1 and twenty_db_location != -1:
ten = self.app.data21[ten_db_location].freq
twenty = self.app.data21[twenty_db_location].freq
sixty_db_frequency = ten * 10**(
5 * (math.log10(twenty) - math.log10(ten)))
self.upper_sixty_db_label.setText(
NanoVNASaver.formatFrequency(sixty_db_frequency) +
" (derived)")
else:
self.upper_sixty_db_label.setText("Not calculated")
if ten_db_location > 0 and twenty_db_location > 0 and ten_db_location != twenty_db_location:
octave_attenuation, decade_attenuation = self.calculateRolloff(
ten_db_location, twenty_db_location)
self.upper_db_per_octave_label.setText(
str(round(octave_attenuation, 3)) + " dB / octave")
self.upper_db_per_decade_label.setText(
str(round(decade_attenuation, 3)) + " dB / decade")
else:
self.upper_db_per_octave_label.setText("Not calculated")
self.upper_db_per_decade_label.setText("Not calculated")
if upper_cutoff_gain < -4 or lower_cutoff_gain < -4:
self.result_label.setText(
"Analysis complete (" + str(len(self.app.data)) +
" points)\n" +
"Insufficient data for analysis. Increase segment count.")
else:
self.result_label.setText("Analysis complete (" +
str(len(self.app.data)) + " points)")
def runAnalysis(self):
from NanoVNASaver.NanoVNASaver import NanoVNASaver
self.reset()
pass_band_location = self.app.markers[0].location
logger.debug("Pass band location: %d", pass_band_location)
if len(self.app.data21) == 0:
logger.debug("No data to analyse")
self.result_label.setText("No data to analyse.")
return
if pass_band_location < 0:
logger.debug("No location for %s", self.app.markers[0].name)
self.result_label.setText("Please place " +
self.app.markers[0].name +
" in the passband.")
return
pass_band_db = NanoVNASaver.gain(self.app.data21[pass_band_location])
logger.debug("Initial passband gain: %d", pass_band_db)
initial_lower_cutoff_location = -1
for i in range(pass_band_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found a cutoff location
initial_lower_cutoff_location = i
break
if initial_lower_cutoff_location < 0:
self.result_label.setText("Lower cutoff location not found.")
return
initial_lower_cutoff_frequency = self.app.data21[
initial_lower_cutoff_location].freq
logger.debug("Found initial lower cutoff frequency at %d",
initial_lower_cutoff_frequency)
initial_upper_cutoff_location = -1
for i in range(pass_band_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found a cutoff location
initial_upper_cutoff_location = i
break
if initial_upper_cutoff_location < 0:
self.result_label.setText("Upper cutoff location not found.")
return
initial_upper_cutoff_frequency = self.app.data21[
initial_upper_cutoff_location].freq
logger.debug("Found initial upper cutoff frequency at %d",
initial_upper_cutoff_frequency)
peak_location = -1
peak_db = NanoVNASaver.gain(
self.app.data21[initial_lower_cutoff_location])
for i in range(initial_lower_cutoff_location,
initial_upper_cutoff_location, 1):
db = NanoVNASaver.gain(self.app.data21[i])
if db > peak_db:
peak_db = db
peak_location = i
logger.debug("Found peak of %f at %d", peak_db,
self.app.data[peak_location].freq)
lower_cutoff_location = -1
pass_band_db = peak_db
for i in range(peak_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found the cutoff location
lower_cutoff_location = i
break
lower_cutoff_frequency = self.app.data21[lower_cutoff_location].freq
lower_cutoff_gain = NanoVNASaver.gain(
self.app.data21[lower_cutoff_location]) - pass_band_db
if lower_cutoff_gain < -4:
logger.debug(
"Lower cutoff frequency found at %f dB - insufficient data points for true -3 dB point.",
lower_cutoff_gain)
logger.debug("Found true lower cutoff frequency at %d",
lower_cutoff_frequency)
self.lower_cutoff_label.setText(
NanoVNASaver.formatFrequency(lower_cutoff_frequency) + " (" +
str(round(lower_cutoff_gain, 1)) + " dB)")
self.app.markers[1].setFrequency(str(lower_cutoff_frequency))
self.app.markers[1].frequencyInput.setText(str(lower_cutoff_frequency))
upper_cutoff_location = -1
pass_band_db = peak_db
for i in range(peak_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 3:
# We found the cutoff location
upper_cutoff_location = i
break
upper_cutoff_frequency = self.app.data21[upper_cutoff_location].freq
upper_cutoff_gain = NanoVNASaver.gain(
self.app.data21[upper_cutoff_location]) - pass_band_db
if upper_cutoff_gain < -4:
logger.debug(
"Upper cutoff frequency found at %f dB - insufficient data points for true -3 dB point.",
upper_cutoff_gain)
logger.debug("Found true upper cutoff frequency at %d",
upper_cutoff_frequency)
self.upper_cutoff_label.setText(
NanoVNASaver.formatFrequency(upper_cutoff_frequency) + " (" +
str(round(upper_cutoff_gain, 1)) + " dB)")
self.app.markers[2].setFrequency(str(upper_cutoff_frequency))
self.app.markers[2].frequencyInput.setText(str(upper_cutoff_frequency))
span = upper_cutoff_frequency - lower_cutoff_frequency
center_frequency = math.sqrt(lower_cutoff_frequency *
upper_cutoff_frequency)
q = center_frequency / span
self.span_label.setText(NanoVNASaver.formatFrequency(span))
self.center_frequency_label.setText(
NanoVNASaver.formatFrequency(center_frequency))
self.quality_label.setText(str(round(q, 2)))
self.app.markers[0].setFrequency(str(round(center_frequency)))
self.app.markers[0].frequencyInput.setText(str(
round(center_frequency)))
# Lower roll-off
lower_six_db_location = -1
for i in range(lower_cutoff_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 6:
# We found 6dB location
lower_six_db_location = i
break
if lower_six_db_location < 0:
self.result_label.setText("Lower 6 dB location not found.")
return
lower_six_db_cutoff_frequency = self.app.data21[
lower_six_db_location].freq
self.lower_six_db_label.setText(
NanoVNASaver.formatFrequency(lower_six_db_cutoff_frequency))
lower_six_db_attenuation = NanoVNASaver.gain(
self.app.data21[lower_six_db_location])
lower_max_attenuation = NanoVNASaver.gain(self.app.data21[0])
frequency_factor = self.app.data21[
0].freq / lower_six_db_cutoff_frequency
lower_attenuation = (lower_max_attenuation - lower_six_db_attenuation)
logger.debug("Measured points: %d Hz and %d Hz",
lower_six_db_cutoff_frequency, self.app.data21[0].freq)
logger.debug("%d dB over %f factor", lower_attenuation,
frequency_factor)
octave_attenuation = lower_attenuation / (
math.log10(frequency_factor) / math.log10(2))
self.lower_db_per_octave_label.setText(
str(round(octave_attenuation, 3)) + " dB / octave")
decade_attenuation = lower_attenuation / math.log10(frequency_factor)
self.lower_db_per_decade_label.setText(
str(round(decade_attenuation, 3)) + " dB / decade")
lower_sixty_db_location = -1
for i in range(lower_six_db_location, -1, -1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 60:
# We found 60dB location! Wow.
lower_sixty_db_location = i
break
if lower_sixty_db_location < 0:
# # We derive 60 dB instead
# factor = 10 * (-54 / decade_attenuation)
# sixty_db_cutoff_frequency = round(six_db_cutoff_frequency + six_db_cutoff_frequency * factor)
# self.upper_sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency) + " (derived)")
self.lower_sixty_db_label.setText("Not calculated")
else:
lower_sixty_db_cutoff_frequency = self.app.data21[
lower_sixty_db_location].freq
self.lower_sixty_db_label.setText(
NanoVNASaver.formatFrequency(lower_sixty_db_cutoff_frequency))
# Upper roll-off
upper_six_db_location = -1
for i in range(upper_cutoff_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 6:
# We found 6dB location
upper_six_db_location = i
break
if upper_six_db_location < 0:
self.result_label.setText("Upper 6 dB location not found.")
return
upper_six_db_cutoff_frequency = self.app.data21[
upper_six_db_location].freq
self.upper_six_db_label.setText(
NanoVNASaver.formatFrequency(upper_six_db_cutoff_frequency))
upper_six_db_attenuation = NanoVNASaver.gain(
self.app.data21[upper_six_db_location])
upper_max_attenuation = NanoVNASaver.gain(self.app.data21[0])
frequency_factor = self.app.data21[
0].freq / upper_six_db_cutoff_frequency
upper_attenuation = (upper_max_attenuation - upper_six_db_attenuation)
logger.debug("Measured points: %d Hz and %d Hz",
upper_six_db_cutoff_frequency, self.app.data21[0].freq)
logger.debug("%d dB over %f factor", upper_attenuation,
frequency_factor)
octave_attenuation = upper_attenuation / (
math.log10(frequency_factor) / math.log10(2))
self.upper_db_per_octave_label.setText(
str(round(octave_attenuation, 3)) + " dB / octave")
decade_attenuation = upper_attenuation / math.log10(frequency_factor)
self.upper_db_per_decade_label.setText(
str(round(decade_attenuation, 3)) + " dB / decade")
upper_sixty_db_location = -1
for i in range(upper_six_db_location, len(self.app.data21), 1):
db = NanoVNASaver.gain(self.app.data21[i])
if (pass_band_db - db) > 60:
# We found 60dB location! Wow.
upper_sixty_db_location = i
break
if upper_sixty_db_location < 0:
# # We derive 60 dB instead
# factor = 10 * (-54 / decade_attenuation)
# sixty_db_cutoff_frequency = round(six_db_cutoff_frequency + six_db_cutoff_frequency * factor)
# self.upper_sixty_db_label.setText(NanoVNASaver.formatFrequency(sixty_db_cutoff_frequency) + " (derived)")
self.upper_sixty_db_label.setText("Not calculated")
else:
upper_sixty_db_cutoff_frequency = self.app.data21[
upper_sixty_db_location].freq
self.upper_sixty_db_label.setText(
NanoVNASaver.formatFrequency(upper_sixty_db_cutoff_frequency))
if upper_cutoff_gain < -4 or lower_cutoff_gain < -4:
self.result_label.setText(
"Analysis complete (" + str(len(self.app.data)) +
" points)\n" +
"Insufficient data for analysis. Increase segment count.")
else:
self.result_label.setText("Analysis complete (" +
str(len(self.app.data)) + " points)")