class SignalIntegrityAppHeadless(object):
projectStack = ProjectStack()
def __init__(self):
# make absolutely sure the directory of this file is the first in the
# python path
thisFileDir = os.path.dirname(os.path.realpath(__file__))
sys.path = [thisFileDir] + sys.path
SignalIntegrity.App.Preferences = Preferences()
SignalIntegrity.App.InstallDir = os.path.dirname(
os.path.abspath(__file__))
self.Drawing = DrawingHeadless(self)
def NullCommand(self):
pass
def OpenProjectFile(self, filename):
if filename is None:
filename = ''
if isinstance(filename, tuple):
filename = ''
filename = str(filename)
if filename == '':
return False
try:
self.fileparts = FileParts(filename)
os.chdir(self.fileparts.AbsoluteFilePath())
self.fileparts = FileParts(filename)
SignalIntegrity.App.Project = ProjectFile().Read(
self.fileparts.FullFilePathExtension('.si'))
self.Drawing.InitFromProject()
except:
return False
self.Drawing.schematic.Consolidate()
for device in self.Drawing.schematic.deviceList:
device.selected = False
for wireProject in SignalIntegrity.App.Project[
'Drawing.Schematic.Wires']:
for vertexProject in wireProject['Vertices']:
vertexProject['Selected'] = False
return True
def SaveProjectToFile(self, filename):
self.fileparts = FileParts(filename)
os.chdir(self.fileparts.AbsoluteFilePath())
self.fileparts = FileParts(filename)
SignalIntegrity.App.Project.Write(self, filename)
def SaveProject(self):
if self.fileparts.filename == '':
return
filename = self.fileparts.AbsoluteFilePath(
) + '/' + self.fileparts.FileNameWithExtension(ext='.si')
self.SaveProjectToFile(filename)
def NetListText(self):
return self.Drawing.schematic.NetList().Text(
) + SignalIntegrity.App.Project['PostProcessing'].NetListLines()
def config(self, cursor=None):
pass
def CalculateSParameters(self):
import SignalIntegrity.Lib as si
if not hasattr(self.Drawing, 'canCalculate'):
self.Drawing.DrawSchematic()
if self.Drawing.canCalculateSParametersFromNetworkAnalyzerModel:
try:
sp = self.SimulateNetworkAnalyzerModel(SParameters=True)
except si.SignalIntegrityException as e:
return None
return (sp,
self.fileparts.FullFilePathExtension('s' + str(sp.m_P) +
'p'))
elif not self.Drawing.canCalculateSParameters:
return None
netListText = self.NetListText()
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
spnp = si.p.SystemSParametersNumericParser(
si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.
Project['CalculationProperties.EndFrequency'], SignalIntegrity.
App.Project['CalculationProperties.FrequencyPoints']),
cacheFileName=cacheFileName)
spnp.AddLines(netListText)
try:
sp = spnp.SParameters()
except si.SignalIntegrityException as e:
return None
return (sp,
self.fileparts.FullFilePathExtension('s' + str(sp.m_P) + 'p'))
def Simulate(self):
if not hasattr(self.Drawing, 'canCalculate'):
self.Drawing.DrawSchematic()
if self.Drawing.canSimulateNetworkAnalyzerModel:
return self.SimulateNetworkAnalyzerModel(SParameters=False)
elif not self.Drawing.canSimulate:
return None
netList = self.Drawing.schematic.NetList()
netListText = self.NetListText()
import SignalIntegrity.Lib as si
fd = si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints'])
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
snp = si.p.SimulatorNumericParser(fd, cacheFileName=cacheFileName)
snp.AddLines(netListText)
try:
transferMatrices = snp.TransferMatrices()
except si.SignalIntegrityException as e:
return None
outputWaveformLabels = netList.OutputNames()
try:
inputWaveformList = self.Drawing.schematic.InputWaveforms()
sourceNames = netList.SourceNames()
except si.SignalIntegrityException as e:
return None
transferMatricesProcessor = si.td.f.TransferMatricesProcessor(
transferMatrices)
si.td.wf.Waveform.adaptionStrategy = 'SinX' if SignalIntegrity.App.Preferences[
'Calculation.UseSinX'] else 'Linear'
try:
outputWaveformList = transferMatricesProcessor.ProcessWaveforms(
inputWaveformList)
except si.SignalIntegrityException as e:
return None
for outputWaveformIndex in range(len(outputWaveformList)):
outputWaveform = outputWaveformList[outputWaveformIndex]
outputWaveformLabel = outputWaveformLabels[outputWaveformIndex]
for device in self.Drawing.schematic.deviceList:
if device['partname'].GetValue() in [
'Output', 'DifferentialVoltageOutput', 'CurrentOutput'
]:
if device['ref'].GetValue() == outputWaveformLabel:
# probes may have different kinds of gain specified
gainProperty = device['gain']
gain = gainProperty.GetValue()
offset = device['offset'].GetValue()
delay = device['td'].GetValue()
if gain != 1.0 or offset != 0.0 or delay != 0.0:
outputWaveform = outputWaveform.DelayBy(
delay) * gain + offset
outputWaveformList[
outputWaveformIndex] = outputWaveform
break
userSampleRate = SignalIntegrity.App.Project[
'CalculationProperties.UserSampleRate']
outputWaveformList = [
wf.Adapt(
si.td.wf.TimeDescriptor(
wf.td.H, int(wf.td.K * userSampleRate / wf.td.Fs),
userSampleRate)) for wf in outputWaveformList
]
return (sourceNames, outputWaveformLabels, transferMatrices,
outputWaveformList)
def VirtualProbe(self):
netList = self.Drawing.schematic.NetList()
netListText = self.NetListText()
import SignalIntegrity.Lib as si
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
snp = si.p.VirtualProbeNumericParser(si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints']),
cacheFileName=cacheFileName)
snp.AddLines(netListText)
try:
transferMatrices = snp.TransferMatrices()
except si.SignalIntegrityException as e:
return None
transferMatricesProcessor = si.td.f.TransferMatricesProcessor(
transferMatrices)
si.td.wf.Waveform.adaptionStrategy = 'SinX' if SignalIntegrity.App.Preferences[
'Calculation.UseSinX'] else 'Linear'
try:
inputWaveformList = self.Drawing.schematic.InputWaveforms()
sourceNames = netList.MeasureNames()
except si.SignalIntegrityException as e:
return None
try:
outputWaveformList = transferMatricesProcessor.ProcessWaveforms(
inputWaveformList)
except si.SignalIntegrityException as e:
return None
outputWaveformLabels = netList.OutputNames()
for outputWaveformIndex in range(len(outputWaveformList)):
outputWaveform = outputWaveformList[outputWaveformIndex]
outputWaveformLabel = outputWaveformLabels[outputWaveformIndex]
for device in self.Drawing.schematic.deviceList:
if device['partname'].GetValue() in [
'Output', 'DifferentialVoltageOutput', 'CurrentOutput'
]:
if device['ref'].GetValue() == outputWaveformLabel:
# probes may have different kinds of gain specified
gainProperty = device['gain']
gain = gainProperty.GetValue()
offset = device['offset'].GetValue()
delay = device['td'].GetValue()
if gain != 1.0 or offset != 0.0 or delay != 0.0:
outputWaveform = outputWaveform.DelayBy(
delay) * gain + offset
outputWaveformList[
outputWaveformIndex] = outputWaveform
break
userSampleRate = SignalIntegrity.App.Project[
'CalculationProperties.UserSampleRate']
outputWaveformList = [
wf.Adapt(
si.td.wf.TimeDescriptor(
wf.td.H, int(wf.td.K * userSampleRate / wf.td.Fs),
userSampleRate)) for wf in outputWaveformList
]
return (sourceNames, outputWaveformLabels, transferMatrices,
outputWaveformList)
def Deembed(self):
netListText = self.NetListText()
import SignalIntegrity.Lib as si
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
dnp = si.p.DeembedderNumericParser(si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints']),
cacheFileName=cacheFileName)
dnp.AddLines(netListText)
try:
sp = dnp.Deembed()
except si.SignalIntegrityException as e:
return None
unknownNames = dnp.m_sd.UnknownNames()
if len(unknownNames) == 1:
sp = [sp]
return (unknownNames, sp)
filename = []
for u in range(len(unknownNames)):
extension = '.s' + str(sp[u].m_P) + 'p'
filename = unknownNames[u] + extension
if self.fileparts.filename != '':
filename.append(self.fileparts.filename + '_' + filename)
return (unknownNames, sp, filename)
def CalculateErrorTerms(self):
if not hasattr(self.Drawing, 'canCalculate'):
self.Drawing.DrawSchematic()
if not self.Drawing.canCalculateErrorTerms:
return None
netList = self.Drawing.schematic.NetList()
netListText = self.NetListText()
import SignalIntegrity.Lib as si
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
etnp = si.p.CalibrationNumericParser(si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints']),
cacheFileName=cacheFileName)
etnp.AddLines(netListText)
try:
cal = etnp.CalculateCalibration()
except si.SignalIntegrityException as e:
return None
return (cal, self.fileparts.FullFilePathExtension('cal'))
def Device(self, ref):
"""
accesses a device by it's reference string
@param ref string reference designator
@return device if found otherwise None
Some examples of how to use this (proj is a project name)
gain=proj.Device('U1')['gain']['Value']
proj.Device('U1')['gain']['Value']=gain
"""
devices = self.Drawing.schematic.deviceList
for device in devices:
if device['ref']['Value'] == ref:
return device
return None
def SimulateNetworkAnalyzerModel(self, SParameters=False):
netList = self.Drawing.schematic.NetList().Text()
import SignalIntegrity.Lib as si
fd = si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints'])
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle() + '_DUTSParameters'
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
spnp = si.p.DUTSParametersNumericParser(fd,
cacheFileName=cacheFileName)
spnp.AddLines(netList)
try:
(DUTSp, NetworkAnalyzerProjectFile) = spnp.SParameters()
except si.SignalIntegrityException as e:
return None
netListText = None
if NetworkAnalyzerProjectFile != None:
level = SignalIntegrityAppHeadless.projectStack.Push()
try:
app = SignalIntegrityAppHeadless()
if app.OpenProjectFile(
os.path.realpath(NetworkAnalyzerProjectFile)):
app.Drawing.DrawSchematic()
netList = app.Drawing.schematic.NetList()
netListText = netList.Text()
else:
pass
except:
pass
finally:
SignalIntegrityAppHeadless.projectStack.Pull(level)
else:
netList = self.Drawing.schematic.NetList()
netListText = self.NetListText()
if netListText == None:
return None
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle(
) + '_TransferMatrices'
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
snp = si.p.NetworkAnalyzerSimulationNumericParser(
fd,
DUTSp,
spnp.NetworkAnalyzerPortConnectionList,
cacheFileName=cacheFileName)
snp.AddLines(netListText)
level = SignalIntegrityAppHeadless.projectStack.Push()
try:
os.chdir(
FileParts(os.path.abspath(
NetworkAnalyzerProjectFile)).AbsoluteFilePath())
transferMatrices = snp.TransferMatrices()
except si.SignalIntegrityException as e:
return None
finally:
SignalIntegrityAppHeadless.projectStack.Pull(level)
sourceNames = snp.m_sd.SourceVector()
gdoDict = {}
if NetworkAnalyzerProjectFile != None:
level = SignalIntegrityAppHeadless.projectStack.Push()
try:
app = SignalIntegrityAppHeadless()
if app.OpenProjectFile(
os.path.realpath(NetworkAnalyzerProjectFile)):
app.Drawing.DrawSchematic()
# get output gain, offset, delay
for name in [rdn[2] for rdn in snp.m_sd.pOutputList]:
gdoDict[name] = {
'gain': float(app.Device(name)['gain']['Value']),
'offset':
float(app.Device(name)['offset']['Value']),
'delay': float(app.Device(name)['td']['Value'])
}
stateList = [
app.Device(sourceNames[port])['state']['Value']
for port in range(snp.simulationNumPorts)
]
self.wflist = []
for driven in range(snp.simulationNumPorts):
thiswflist = []
for port in range(snp.simulationNumPorts):
app.Device(sourceNames[port])['state'][
'Value'] = 'on' if port == driven else 'off'
for wfIndex in range(len(sourceNames)):
thiswflist.append(
app.Device(sourceNames[wfIndex]).Waveform())
self.wflist.append(thiswflist)
for port in range(snp.simulationNumPorts):
app.Device(sourceNames[port]
)['state']['Value'] = stateList[port]
else:
pass
except:
pass
finally:
SignalIntegrityAppHeadless.projectStack.Pull(level)
else:
stateList = [
app.Device(sourceNames[port])['state']['Value']
for port in range(snp.simulationNumPorts)
]
self.wflist = []
for name in [rdn[2] for rdn in snp.m_sd.pOutputList]:
gdoDict[name] = {
'gain': float(app.Device()[name]['gain']['Value']),
'offset': float(app.Device()[name]['offset']['Value']),
'delay': float(app.Device()[name]['td']['Value'])
}
for driven in range(snp.simulationNumPorts):
thiswflist = []
for port in range(snp.simulationNumPorts):
app.Device(
sourceNames[port]
)['state']['Value'] = 'on' if port == driven else 'off'
for wfIndex in range(len(sourceNames)):
thiswflist.append(
app.Device(sourceNames[wfIndex]).Waveform())
self.wflist.append(thiswflist)
for port in range(snp.simulationNumPorts):
app.Device(
sourceNames[port])['state']['Value'] = stateList[port]
self.transferMatriceProcessor = si.td.f.TransferMatricesProcessor(
transferMatrices)
si.td.wf.Waveform.adaptionStrategy = 'SinX' if SignalIntegrity.App.Preferences[
'Calculation.UseSinX'] else 'Linear'
try:
outputwflist = []
for port in range(len(self.wflist)):
outputwflist.append(
self.transferMatriceProcessor.ProcessWaveforms(
self.wflist[port], adaptToLargest=True))
except si.SignalIntegrityException as e:
return None
#
# The list of list of input waveforms have been processed processed, generating a list of list of output waveforms in
# self.outputwflist. The names of the output waveforms are in snp.m_sd.pOutputList.
#
outputwflist = [[
wf.Adapt(
si.td.wf.TimeDescriptor(wf.td[wf.td.IndexOfTime(-5e-9)],
fd.TimeDescriptor().K, wf.td.Fs))
for wf in driven
] for driven in outputwflist]
# The port connection list, which is a list of True or False for each port on the network analyzer, is
# converted to a list of network port indices corresponding to the driven ports.
#
portConnections = []
for pci in range(len(snp.PortConnectionList)):
if snp.PortConnectionList[pci]: portConnections.append(pci)
outputWaveformList = []
outputWaveformLabels = []
for r in range(len(outputwflist)):
wflist = outputwflist[r]
for c in range(len(wflist)):
wf = wflist[c]
wfName = snp.m_sd.pOutputList[c][2]
gain = gdoDict[wfName]['gain']
offset = gdoDict[wfName]['offset']
delay = gdoDict[wfName]['delay']
if gain != 1.0 or offset != 0.0 or delay != 0.0:
wf = wf.DelayBy(delay) * gain + offset
outputWaveformList.append(wf)
outputWaveformLabels.append(snp.m_sd.pOutputList[c][2] +
str(portConnections[r] + 1))
userSampleRate = SignalIntegrity.App.Project[
'CalculationProperties.UserSampleRate']
outputWaveformList = [
wf.Adapt(
si.td.wf.TimeDescriptor(
wf.td.H, int(wf.td.K * userSampleRate / wf.td.Fs),
userSampleRate)) for wf in outputWaveformList
]
td = si.td.wf.TimeDescriptor(
-5e-9,
SignalIntegrity.App.Project['CalculationProperties.TimePoints'],
SignalIntegrity.App.Project['CalculationProperties.BaseSampleRate']
)
if snp.simulationType != 'CW':
# note this matrix is transposed from what is normally expected
Vmat = [[
outputWaveformList[outputWaveformLabels.index(
'V' + str(portConnections[r] + 1) +
str(portConnections[c] + 1))]
for r in range(len(portConnections))
] for c in range(len(portConnections))]
for vli in range(len(Vmat)):
tdr = si.m.tdr.TDRWaveformToSParameterConverter(
WindowForwardHalfWidthTime=200e-12,
WindowReverseHalfWidthTime=200e-12,
WindowRaisedCosineDuration=50e-12,
Step=(snp.simulationType == 'TDRStep'),
Length=0,
Denoise=True,
DenoisePercent=20.,
Inverted=False,
fd=td.FrequencyList())
tdr.Convert(Vmat[vli], vli)
for r in range(len(portConnections)):
outputWaveformList.append(tdr.IncidentWaveform if r ==
vli else si.td.wf.Waveform(td))
outputWaveformLabels.append('A' +
str(portConnections[r] + 1) +
str(portConnections[vli] + 1))
for r in range(len(portConnections)):
outputWaveformList.append(tdr.ReflectWaveforms[r])
outputWaveformLabels.append('B' +
str(portConnections[r] + 1) +
str(portConnections[vli] + 1))
if not SParameters:
return (sourceNames, outputWaveformLabels, transferMatrices,
outputWaveformList)
else:
# waveforms are adapted this way to give the horizontal offset that it already has closest to
#-5 ns, with the correct number of points without resampling the waveform in any way.
frequencyContentList = []
for wf in outputWaveformList:
td = si.td.wf.TimeDescriptor(
-5e-9, SignalIntegrity.App.
Project['CalculationProperties.TimePoints'],
SignalIntegrity.App.
Project['CalculationProperties.BaseSampleRate'])
td.H = wf.TimeDescriptor()[wf.TimeDescriptor().IndexOfTime(
td.H)]
fc = wf.Adapt(td).FrequencyContent()
frequencyContentList.append(fc)
Afc = [[
frequencyContentList[outputWaveformLabels.index(
'A' + str(portConnections[r] + 1) +
str(portConnections[c] + 1))]
for c in range(len(portConnections))
] for r in range(len(portConnections))]
Bfc = [[
frequencyContentList[outputWaveformLabels.index(
'B' + str(portConnections[r] + 1) +
str(portConnections[c] + 1))]
for c in range(len(portConnections))
] for r in range(len(portConnections))]
from numpy import array
from numpy.linalg import inv
frequencyList = td.FrequencyList()
data = [None for _ in range(len(frequencyList))]
for n in range(len(frequencyList)):
B = [[Bfc[r][c][n] for c in range(snp.simulationNumPorts)]
for r in range(snp.simulationNumPorts)]
A = [[Afc[r][c][n] for c in range(snp.simulationNumPorts)]
for r in range(snp.simulationNumPorts)]
data[n] = (array(B).dot(inv(array(A)))).tolist()
sp = si.sp.SParameters(frequencyList, data)
return sp
class SignalIntegrityAppHeadless(object):
def __init__(self):
# make absolutely sure the directory of this file is the first in the
# python path
thisFileDir = os.path.dirname(os.path.realpath(__file__))
sys.path = [thisFileDir] + sys.path
SignalIntegrity.App.Preferences = Preferences()
SignalIntegrity.App.InstallDir = os.path.dirname(
os.path.abspath(__file__))
self.Drawing = DrawingHeadless(self)
def NullCommand(self):
pass
def OpenProjectFile(self, filename):
if filename is None:
filename = ''
if isinstance(filename, tuple):
filename = ''
filename = str(filename)
if filename == '':
return False
try:
self.fileparts = FileParts(filename)
os.chdir(self.fileparts.AbsoluteFilePath())
self.fileparts = FileParts(filename)
SignalIntegrity.App.Project = ProjectFile().Read(
self.fileparts.FullFilePathExtension('.si'))
self.Drawing.InitFromProject()
except:
return False
self.Drawing.schematic.Consolidate()
for device in self.Drawing.schematic.deviceList:
device.selected = False
for wireProject in SignalIntegrity.App.Project[
'Drawing.Schematic.Wires']:
for vertexProject in wireProject['Vertices']:
vertexProject['Selected'] = False
return True
def SaveProjectToFile(self, filename):
self.fileparts = FileParts(filename)
os.chdir(self.fileparts.AbsoluteFilePath())
self.fileparts = FileParts(filename)
SignalIntegrity.App.Project.Write(self, filename)
def SaveProject(self):
if self.fileparts.filename == '':
return
filename = self.fileparts.AbsoluteFilePath(
) + '/' + self.fileparts.FileNameWithExtension(ext='.si')
self.SaveProjectToFile(filename)
def config(self, cursor=None):
pass
def CalculateSParameters(self):
netList = self.Drawing.schematic.NetList().Text()
import SignalIntegrity.Lib as si
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
spnp = si.p.SystemSParametersNumericParser(
si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.
Project['CalculationProperties.EndFrequency'], SignalIntegrity.
App.Project['CalculationProperties.FrequencyPoints']),
cacheFileName=cacheFileName)
spnp.AddLines(netList)
try:
sp = spnp.SParameters()
except si.SignalIntegrityException as e:
return None
return (sp,
self.fileparts.FullFilePathExtension('s' + str(sp.m_P) + 'p'))
def Simulate(self):
netList = self.Drawing.schematic.NetList()
netListText = netList.Text()
import SignalIntegrity.Lib as si
fd = si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints'])
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
snp = si.p.SimulatorNumericParser(fd, cacheFileName=cacheFileName)
snp.AddLines(netListText)
try:
transferMatrices = snp.TransferMatrices()
except si.SignalIntegrityException as e:
return None
outputWaveformLabels = netList.OutputNames()
try:
inputWaveformList = self.Drawing.schematic.InputWaveforms()
sourceNames = netList.SourceNames()
except si.SignalIntegrityException as e:
return None
transferMatricesProcessor = si.td.f.TransferMatricesProcessor(
transferMatrices)
si.td.wf.Waveform.adaptionStrategy = 'SinX' if SignalIntegrity.App.Preferences[
'Calculation.UseSinX'] else 'Linear'
try:
outputWaveformList = transferMatricesProcessor.ProcessWaveforms(
inputWaveformList)
except si.SignalIntegrityException as e:
return None
for outputWaveformIndex in range(len(outputWaveformList)):
outputWaveform = outputWaveformList[outputWaveformIndex]
outputWaveformLabel = outputWaveformLabels[outputWaveformIndex]
for device in self.Drawing.schematic.deviceList:
if device['partname'].GetValue() in [
'Output', 'DifferentialVoltageOutput', 'CurrentOutput'
]:
if device['ref'].GetValue() == outputWaveformLabel:
# probes may have different kinds of gain specified
gainProperty = device['gain']
gain = gainProperty.GetValue()
offset = device['offset'].GetValue()
delay = device['td'].GetValue()
if gain != 1.0 or offset != 0.0 or delay != 0.0:
outputWaveform = outputWaveform.DelayBy(
delay) * gain + offset
outputWaveformList[
outputWaveformIndex] = outputWaveform
break
userSampleRate = SignalIntegrity.App.Project[
'CalculationProperties.UserSampleRate']
outputWaveformList = [
wf.Adapt(
si.td.wf.TimeDescriptor(
wf.td.H, int(wf.td.K * userSampleRate / wf.td.Fs),
userSampleRate)) for wf in outputWaveformList
]
return (sourceNames, outputWaveformLabels, transferMatrices,
outputWaveformList)
def VirtualProbe(self):
netList = self.Drawing.schematic.NetList()
netListText = netList.Text()
import SignalIntegrity.Lib as si
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
snp = si.p.VirtualProbeNumericParser(si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints']),
cacheFileName=cacheFileName)
snp.AddLines(netListText)
try:
transferMatrices = snp.TransferMatrices()
except si.SignalIntegrityException as e:
return None
transferMatricesProcessor = si.td.f.TransferMatricesProcessor(
transferMatrices)
si.td.wf.Waveform.adaptionStrategy = 'SinX' if SignalIntegrity.App.Preferences[
'Calculation.UseSinX'] else 'Linear'
try:
inputWaveformList = self.Drawing.schematic.InputWaveforms()
sourceNames = netList.MeasureNames()
except si.SignalIntegrityException as e:
return None
try:
outputWaveformList = transferMatricesProcessor.ProcessWaveforms(
inputWaveformList)
except si.SignalIntegrityException as e:
return None
outputWaveformLabels = netList.OutputNames()
for outputWaveformIndex in range(len(outputWaveformList)):
outputWaveform = outputWaveformList[outputWaveformIndex]
outputWaveformLabel = outputWaveformLabels[outputWaveformIndex]
for device in self.Drawing.schematic.deviceList:
if device['partname'].GetValue() in [
'Output', 'DifferentialVoltageOutput', 'CurrentOutput'
]:
if device['ref'].GetValue() == outputWaveformLabel:
# probes may have different kinds of gain specified
gainProperty = device['gain']
gain = gainProperty.GetValue()
offset = device['offset'].GetValue()
delay = device['td'].GetValue()
if gain != 1.0 or offset != 0.0 or delay != 0.0:
outputWaveform = outputWaveform.DelayBy(
delay) * gain + offset
outputWaveformList[
outputWaveformIndex] = outputWaveform
break
userSampleRate = SignalIntegrity.App.Project[
'CalculationProperties.UserSampleRate']
outputWaveformList = [
wf.Adapt(
si.td.wf.TimeDescriptor(
wf.td.H, int(wf.td.K * userSampleRate / wf.td.Fs),
userSampleRate)) for wf in outputWaveformList
]
return (sourceNames, outputWaveformLabels, transferMatrices,
outputWaveformList)
def Deembed(self):
netList = self.Drawing.schematic.NetList().Text()
import SignalIntegrity.Lib as si
cacheFileName = None
if SignalIntegrity.App.Preferences['Cache.CacheResults']:
cacheFileName = self.fileparts.FileNameTitle()
si.sd.Numeric.trySVD = SignalIntegrity.App.Preferences[
'Calculation.TrySVD']
dnp = si.p.DeembedderNumericParser(si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.
Project['CalculationProperties.FrequencyPoints']),
cacheFileName=cacheFileName)
dnp.AddLines(netList)
try:
sp = dnp.Deembed()
except si.SignalIntegrityException as e:
return None
unknownNames = dnp.m_sd.UnknownNames()
if len(unknownNames) == 1:
sp = [sp]
return (unknownNames, sp)
filename = []
for u in range(len(unknownNames)):
extension = '.s' + str(sp[u].m_P) + 'p'
filename = unknownNames[u] + extension
if self.fileparts.filename != '':
filename.append(self.fileparts.filename + '_' + filename)
return (unknownNames, sp, filename)
class SParametersDialog(tk.Toplevel):
def __init__(self, parent,sp,filename=None,title=None,buttonLabels=None):
tk.Toplevel.__init__(self, parent)
self.parent=parent
self.withdraw()
self.fileparts=FileParts(filename)
if title is None:
if self.fileparts.filename =='':
self.title('S-parameters')
else:
self.title('S-parameters: '+self.fileparts.FileNameTitle())
else:
if filename is None:
self.title(title)
else:
self.title(title+': '+self.fileparts.FileNameTitle())
img = tk.PhotoImage(file=SignalIntegrity.App.IconsBaseDir+'AppIcon2.gif')
self.tk.call('wm', 'iconphoto', self._w, img)
self.protocol("WM_DELETE_WINDOW", self.onClosing)
self.variableLineWidth = tk.BooleanVar()
self.showPassivityViolations = tk.BooleanVar()
self.showCausalityViolations = tk.BooleanVar()
self.showImpedance = tk.BooleanVar()
self.logScale = tk.BooleanVar()
# the Doers - the holder of the commands, menu elements, toolbar elements, and key bindings
self.ReadSParametersFromFileDoer = Doer(self.onReadSParametersFromFile).AddKeyBindElement(self,'<Control-o>').AddHelpElement('Control-Help:Open-S-parameter-File')
self.WriteSParametersToFileDoer = Doer(self.onWriteSParametersToFile).AddKeyBindElement(self,'<Control-s>').AddHelpElement('Control-Help:Save-S-parameter-File')
self.Matplotlib2tikzDoer = Doer(self.onMatplotlib2TikZ)
# ------
self.CalculationPropertiesDoer = Doer(self.onCalculationProperties).AddHelpElement('Control-Help:Calculation-Properties')
self.ResampleDoer = Doer(self.onResample).AddHelpElement('Control-Help:Resample')
self.EnforcePassivityDoer = Doer(self.onEnforcePassivity).AddHelpElement('Control-Help:Enforce-Passivity')
self.EnforceCausalityDoer = Doer(self.onEnforceCausality).AddHelpElement('Control-Help:Enforce-Causality')
self.WaveletDenoiseDoer = Doer(self.onWaveletDenoise).AddHelpElement('Control-Help:Wavelet-Denoise')
# ------
self.HelpDoer = Doer(self.onHelp).AddHelpElement('Control-Help:Open-Help-File')
self.ControlHelpDoer = Doer(self.onControlHelp).AddHelpElement('Control-Help:Control-Help')
# ------
self.VariableLineWidthDoer = Doer(self.PlotSParameter).AddHelpElement('Control-Help:Variable-Line-Width')
self.ShowPassivityViolationsDoer = Doer(self.PlotSParameter).AddHelpElement('Control-Help:Show-Passivity-Violations')
self.ShowCausalityViolationsDoer = Doer(self.PlotSParameter).AddHelpElement('Control-Help:Show-Causality-Violations')
self.ShowImpedanceDoer = Doer(self.PlotSParameter).AddHelpElement('Control-Help:Show-Impedance')
self.LogScaleDoer = Doer(self.PlotSParameter).AddHelpElement('Control-Help:Log-Scale')
# ------
self.EscapeDoer = Doer(self.onEscape).AddKeyBindElement(self,'<Escape>').DisableHelp()
# The menu system
TheMenu=tk.Menu(self)
self.config(menu=TheMenu)
# ------
FileMenu=tk.Menu(self)
TheMenu.add_cascade(label='File',menu=FileMenu,underline=0)
self.WriteSParametersToFileDoer.AddMenuElement(FileMenu,label="Save",accelerator='Ctrl+S',underline=0)
self.ReadSParametersFromFileDoer.AddMenuElement(FileMenu,label="Open File",accelerator='Ctrl+O',underline=0)
FileMenu.add_separator()
self.Matplotlib2tikzDoer.AddMenuElement(FileMenu,label='Output to LaTeX (TikZ)',underline=10)
# ------
CalcMenu=tk.Menu(self)
TheMenu.add_cascade(label='Calculate',menu=CalcMenu,underline=0)
self.CalculationPropertiesDoer.AddMenuElement(CalcMenu,label='Calculation Properties',underline=0)
#CalcMenu.add_separator()
self.ResampleDoer.AddMenuElement(CalcMenu,label='Resample',underline=0)
#CalcMenu.add_separator()
self.EnforcePassivityDoer.AddMenuElement(CalcMenu,label='Enforce Passivity',underline=8)
self.EnforceCausalityDoer.AddMenuElement(CalcMenu,label='Enforce Causality',underline=9)
self.WaveletDenoiseDoer.AddMenuElement(CalcMenu,label='Wavelet Denoise',underline=0)
# ------
ViewMenu=tk.Menu(self)
TheMenu.add_cascade(label='View',menu=ViewMenu,underline=0)
self.VariableLineWidthDoer.AddCheckButtonMenuElement(ViewMenu,label='Variable Line Width',underline=9,onvalue=True,offvalue=False,variable=self.variableLineWidth)
self.ShowPassivityViolationsDoer.AddCheckButtonMenuElement(ViewMenu,label='Show Passivity Violations',underline=5,onvalue=True,offvalue=False,variable=self.showPassivityViolations)
self.ShowCausalityViolationsDoer.AddCheckButtonMenuElement(ViewMenu,label='Show Causality Violations',underline=6,onvalue=True,offvalue=False,variable=self.showCausalityViolations)
self.ShowImpedanceDoer.AddCheckButtonMenuElement(ViewMenu,label='Show Impedance',underline=5,onvalue=True,offvalue=False,variable=self.showImpedance)
self.LogScaleDoer.AddCheckButtonMenuElement(ViewMenu,label='Log Scale',underline=4,onvalue=True,offvalue=False,variable=self.logScale)
# ------
HelpMenu=tk.Menu(self)
TheMenu.add_cascade(label='Help',menu=HelpMenu,underline=0)
self.HelpDoer.AddMenuElement(HelpMenu,label='Open Help File',underline=0)
self.ControlHelpDoer.AddMenuElement(HelpMenu,label='Control Help',underline=0)
# The Toolbar
ToolBarFrame = tk.Frame(self)
ToolBarFrame.pack(side=tk.TOP,fill=tk.X,expand=tk.NO)
self.ReadSParametersFromFileDoer.AddToolBarElement(ToolBarFrame,iconfile=SignalIntegrity.App.IconsDir+'document-open-2.gif').Pack(side=tk.LEFT,fill=tk.NONE,expand=tk.NO)
self.WriteSParametersToFileDoer.AddToolBarElement(ToolBarFrame,iconfile=SignalIntegrity.App.IconsDir+'document-save-2.gif').Pack(side=tk.LEFT,fill=tk.NONE,expand=tk.NO)
tk.Frame(self,bd=2,relief=tk.SUNKEN).pack(side=tk.LEFT,fill=tk.X,padx=5,pady=5)
self.CalculationPropertiesDoer.AddToolBarElement(ToolBarFrame,iconfile=SignalIntegrity.App.IconsDir+'tooloptions.gif').Pack(side=tk.LEFT,fill=tk.NONE,expand=tk.NO)
tk.Frame(ToolBarFrame,height=2,bd=2,relief=tk.RAISED).pack(side=tk.LEFT,fill=tk.X,padx=5,pady=5)
self.HelpDoer.AddToolBarElement(ToolBarFrame,iconfile=SignalIntegrity.App.IconsDir+'help-contents-5.gif').Pack(side=tk.LEFT,fill=tk.NONE,expand=tk.NO)
self.ControlHelpDoer.AddToolBarElement(ToolBarFrame,iconfile=SignalIntegrity.App.IconsDir+'help-3.gif').Pack(side=tk.LEFT,fill=tk.NONE,expand=tk.NO)
topFrame=tk.Frame(self)
topFrame.pack(side=tk.TOP,fill=tk.BOTH,expand=tk.YES)
bottomFrame=tk.Frame(self)
bottomFrame.pack(side=tk.TOP,fill=tk.BOTH,expand=tk.YES)
topLeftFrame=tk.Frame(topFrame)
topLeftFrame.pack(side=tk.LEFT,fill=tk.BOTH,expand=tk.YES)
topRightFrame=tk.Frame(topFrame)
topRightFrame.pack(side=tk.LEFT,fill=tk.BOTH,expand=tk.YES)
bottomLeftFrame=tk.Frame(bottomFrame)
bottomLeftFrame.pack(side=tk.LEFT,fill=tk.BOTH,expand=tk.YES)
bottomRightFrame=tk.Frame(bottomFrame)
bottomRightFrame.pack(side=tk.LEFT,fill=tk.BOTH,expand=tk.YES)
self.topLeftFigure=Figure(figsize=(5,2), dpi=100)
self.topLeftPlot=self.topLeftFigure.add_subplot(111)
self.topLeftCanvas=FigureCanvasTkAgg(self.topLeftFigure, master=topLeftFrame)
self.topLeftCanvas.get_tk_widget().pack(side=tk.TOP, fill=tk.X, expand=1)
self.topLeftToolbar = NavigationToolbar2Tk( self.topLeftCanvas, topLeftFrame )
self.topLeftToolbar.update()
self.topLeftCanvas._tkcanvas.pack(side=tk.TOP, fill=tk.BOTH, expand=1)
self.topRightFigure=Figure(figsize=(5,2), dpi=100)
self.topRightPlot=self.topRightFigure.add_subplot(111)
self.topRightCanvas=FigureCanvasTkAgg(self.topRightFigure, master=topRightFrame)
self.topRightCanvas.get_tk_widget().pack(side=tk.TOP, fill=tk.X, expand=1)
self.topRightToolbar = NavigationToolbar2Tk( self.topRightCanvas, topRightFrame )
self.topRightToolbar.update()
self.topRightCanvas._tkcanvas.pack(side=tk.TOP, fill=tk.BOTH, expand=1)
self.topRightCanvasControlsFrame=tk.Frame(topRightFrame)
self.topRightCanvasControlsFrame.pack(side=tk.TOP, fill=tk.X, expand=tk.NO)
tk.Button(self.topRightCanvasControlsFrame,text='unwrap',command=self.onUnwrap).pack(side=tk.LEFT,expand=tk.NO,fill=tk.NONE)
self.delay=PartPropertyDelay(0.)
self.delayViewerProperty=ViewerProperty(self.topRightCanvasControlsFrame,self.delay,self.onDelayEntered)
self.bottomLeftFigure=Figure(figsize=(5,2), dpi=100)
self.bottomLeftPlot=self.bottomLeftFigure.add_subplot(111)
self.bottomLeftCanvas=FigureCanvasTkAgg(self.bottomLeftFigure, master=bottomLeftFrame)
self.bottomLeftCanvas.get_tk_widget().pack(side=tk.TOP, fill=tk.X, expand=1)
self.bottomLeftToolbar = NavigationToolbar2Tk( self.bottomLeftCanvas, bottomLeftFrame )
self.bottomLeftToolbar.update()
self.bottomLeftCanvas._tkcanvas.pack(side=tk.TOP, fill=tk.BOTH, expand=1)
self.bottomRightFigure=Figure(figsize=(5,2), dpi=100)
self.bottomRightPlot=self.bottomRightFigure.add_subplot(111)
self.bottomRightCanvas=FigureCanvasTkAgg(self.bottomRightFigure, master=bottomRightFrame)
self.bottomRightCanvas.get_tk_widget().pack(side=tk.TOP, fill=tk.X, expand=1)
self.bottomRightToolbar = NavigationToolbar2Tk( self.bottomRightCanvas, bottomRightFrame )
self.bottomRightToolbar.update()
self.bottomRightCanvas._tkcanvas.pack(side=tk.TOP, fill=tk.BOTH, expand=1)
controlsFrame = tk.Frame(self)
controlsFrame.pack(side=tk.TOP,fill=tk.X,expand=tk.NO)
self.sButtonsFrame = tk.Frame(controlsFrame, bd=1, relief=tk.SUNKEN)
self.sButtonsFrame.pack(side=tk.LEFT,expand=tk.NO,fill=tk.NONE)
self.resampleButton=tk.Button(controlsFrame,text='resample',command=self.onResample)
self.resampleButton.pack(side=tk.LEFT,expand=tk.NO,fill=tk.NONE)
self.sp=sp
if buttonLabels is None:
numPorts=self.sp.m_P
buttonLabels=[['s'+str(toP+1)+str(fromP+1) for fromP in range(numPorts)] for toP in range(numPorts)]
self.referenceImpedance=PartPropertyReferenceImpedance(self.sp.m_Z0)
self.referenceImpedanceProperty=ViewerProperty(controlsFrame,self.referenceImpedance,self.onReferenceImpedanceEntered)
else:
# button labels are a proxy for transfer parameters (until I do something better)
self.showPassivityViolations.set(False)
self.ShowPassivityViolationsDoer.Activate(False)
self.ShowCausalityViolationsDoer.Activate(False)
self.ShowImpedanceDoer.Activate(False)
#self.LogScaleDoer.Activate(False)
self.EnforcePassivityDoer.Activate(False)
self.EnforceCausalityDoer.Activate(False)
self.WaveletDenoiseDoer.Activate(False)
self.ReadSParametersFromFileDoer.Activate(False)
self.buttonLabels=buttonLabels
self.buttons=[]
for toP in range(len(buttonLabels)):
buttonrow=[]
rowFrame=tk.Frame(self.sButtonsFrame)
rowFrame.pack(side=tk.TOP,expand=tk.NO,fill=tk.NONE)
for fromP in range(len(buttonLabels[0])):
thisButton=tk.Button(rowFrame,text=buttonLabels[toP][fromP],width=len(buttonLabels[toP][fromP]),command=lambda x=toP+1,y=fromP+1: self.onSelectSParameter(x,y))
thisButton.pack(side=tk.LEFT,fill=tk.NONE,expand=tk.NO)
buttonrow.append(thisButton)
self.buttons.append(buttonrow)
self.fromPort = 1
self.toPort = 1
try:
from matplotlib2tikz import save as tikz_save
except:
self.Matplotlib2tikzDoer.Activate(False)
self.buttons[self.toPort-1][self.fromPort-1].config(relief=tk.SUNKEN)
self.PlotSParameter()
self.deiconify()
# self.geometry("%+d%+d" % (self.parent.root.winfo_x()+self.parent.root.winfo_width()/2-self.winfo_width()/2,
# self.parent.root.winfo_y()+self.parent.root.winfo_height()/2-self.winfo_height()/2))
def onClosing(self):
self.withdraw()
self.destroy()
def destroy(self):
tk.Toplevel.withdraw(self)
tk.Toplevel.destroy(self)
def PlotSParameter(self):
import SignalIntegrity.Lib as si
self.topLeftPlot.cla()
self.topRightPlot.cla()
self.bottomLeftPlot.cla()
self.bottomRightPlot.cla()
if not SignalIntegrity.App.Preferences['Appearance.PlotCursorValues']:
self.topLeftPlot.format_coord = lambda x, y: ''
self.topRightPlot.format_coord = lambda x, y: ''
self.bottomLeftPlot.format_coord = lambda x, y: ''
self.bottomRightPlot.format_coord = lambda x, y: ''
fr=self.sp.FrequencyResponse(self.toPort,self.fromPort)
ir=fr.ImpulseResponse()
y=fr.Response('dB')
self.freqLabel=ToSI(fr.Frequencies()[-1],'Hz')[-3:]
freqLabelDivisor=FromSI('1. '+self.freqLabel,'Hz')
x=fr.Frequencies(freqLabelDivisor)
if self.showPassivityViolations.get():
self.passivityViolations=[]
s=self.sp._LargestSingularValues()
for n in range(len(s)):
if s[n]-1 > 1e-15:
dotsize=max(min(20.,math.log10(s[0])/math.log10(1.01)*20.),1e-15)
self.passivityViolations.append([x[n],y[n],dotsize])
lw=[min(1.,math.sqrt(w))*1.5 for w in fr.Response('mag')]
fastway=True
# this works - properly displays on log plot, only it is just too slow!
# if self.logScale.get():
# fastway=False
if self.variableLineWidth.get():
if fastway:
segments = [[[x[i],y[i]],[x[i+1],y[i+1]]] for i in range(len(x)-1)]
slw=lw[:-1]
lc = LineCollection(segments, linewidths=slw,color='blue')
self.topLeftPlot.add_collection(lc)
else:
for i in range(len(x)-1):
if self.logScale.get():
self.topLeftPlot.semilogx(x[i:i+2],y[i:i+2],linewidth=lw[i],color='blue')
else:
self.topLeftPlot.plot(x[i:i+2],y[i:i+2],linewidth=lw[i],color='blue')
else:
if self.logScale.get():
self.topLeftPlot.semilogx(x,y)
else:
self.topLeftPlot.plot(x,y)
if self.showPassivityViolations.get():
self.topLeftPlot.scatter(
[c[0] for c in self.passivityViolations],
[c[1] for c in self.passivityViolations],
s=[c[2] for c in self.passivityViolations],
color='red')
self.topLeftPlot.set_xlim(xmin=min(x))
self.topLeftPlot.set_xlim(xmax=max(x))
self.topLeftPlot.set_ylim(ymin=max(min(y)-1.,-60.0))
self.topLeftPlot.set_ylim(ymax=max(y)+1.)
self.topLeftPlot.set_ylabel('magnitude (dB)',fontsize=10)
self.topLeftPlot.set_xlabel('frequency ('+self.freqLabel+')',fontsize=10)
y=fr.Response('deg')
x=fr.Frequencies(freqLabelDivisor)
if self.variableLineWidth.get():
if fastway:
segments = [[[x[i],y[i]],[x[i+1],y[i+1]]] for i in range(len(x)-1)]
slw=lw[:-1]
lc = LineCollection(segments, linewidths=slw,color='blue')
self.topRightPlot.add_collection(lc)
else:
for i in range(len(x)-1):
if self.logScale.get():
self.topRightPlot.semilogx(x[i:i+2],y[i:i+2],linewidth=lw[i],color='blue')
else:
self.topRightPlot.plot(x[i:i+2],y[i:i+2],linewidth=lw[i],color='blue')
else:
if self.logScale.get():
self.topRightPlot.semilogx(x,y)
else:
self.topRightPlot.plot(x,y)
self.topRightPlot.set_xlim(xmin=min(x))
self.topRightPlot.set_xlim(xmax=max(x))
self.topRightPlot.set_ylim(ymin=min(y)-1)
self.topRightPlot.set_ylim(ymax=max(y)+1)
self.topRightPlot.set_ylabel('phase (degrees)',fontsize=10)
self.topRightPlot.set_xlabel('frequency ('+self.freqLabel+')',fontsize=10)
if ir is not None:
if self.buttonLabels[self.toPort-1][self.fromPort-1][:2]=='i/' or self.buttonLabels[self.toPort-1][self.fromPort-1][:3]=='di/':
print('Integrate')
ir=si.td.wf.ImpulseResponse(ir.Integral(addPoint=False))
if self.buttonLabels[self.toPort-1][self.fromPort-1][:3]=='di/':
print('Integrate')
ir=si.td.wf.ImpulseResponse(ir.Integral(addPoint=False)*ir.td.Fs)
y=ir.Values()
timeLabel=ToSI(ir.Times()[-1],'s')[-2:]
timeLabelDivisor=FromSI('1. '+timeLabel,'s')
x=ir.Times(timeLabelDivisor)
self.bottomLeftPlot.plot(x,y)
if self.showCausalityViolations.get():
self.causalityViolations=[]
Ts=1./ir.td.Fs/1e-9
for k in range(len(x)):
if x[k]<=-Ts and abs(y[k])>0:
dotsize=max(min(20.,abs(y[k])/0.1*20.),1e-15)
self.causalityViolations.append([x[k],y[k],dotsize])
self.bottomLeftPlot.scatter(
[c[0] for c in self.causalityViolations],
[c[1] for c in self.causalityViolations],
s=[c[2] for c in self.causalityViolations],
color='red')
self.bottomLeftPlot.set_ylim(ymin=min(min(y)*1.05,-0.1))
self.bottomLeftPlot.set_ylim(ymax=max(max(y)*1.05,0.1))
self.bottomLeftPlot.set_xlim(xmin=min(x))
self.bottomLeftPlot.set_xlim(xmax=max(x))
self.bottomLeftPlot.set_ylabel('amplitude',fontsize=10)
self.bottomLeftPlot.set_xlabel('time ('+timeLabel+')',fontsize=10)
firFilter=ir.FirFilter()
stepWaveformTimeDescriptor=ir.td/firFilter.FilterDescriptor()
stepWaveform=si.td.wf.StepWaveform(stepWaveformTimeDescriptor)
stepResponse=stepWaveform*firFilter
y=stepResponse.Values()
x=stepResponse.Times(timeLabelDivisor)
if self.showImpedance.get() and (self.fromPort == self.toPort):
Z0=self.referenceImpedance.GetValue()
y=[3000. if (1-yv)<=.000001 else min(Z0*(1+yv)/(1-yv),3000) for yv in y]
x=[xv/2 for xv in x]
self.bottomRightPlot.set_ylabel('impedance (Ohms)',fontsize=10)
self.bottomRightPlot.set_xlabel('length ('+timeLabel+')',fontsize=10)
self.bottomRightPlot.set_ylim(ymin=min(min(y)*1.05,Z0-1))
else:
self.bottomRightPlot.set_ylabel('amplitude',fontsize=10)
self.bottomRightPlot.set_xlabel('time ('+timeLabel+')',fontsize=10)
self.bottomRightPlot.set_ylim(ymin=min(min(y)*1.05,-0.1))
self.bottomRightPlot.plot(x,y)
self.bottomRightPlot.set_ylim(ymax=max(max(y)*1.05,0.1))
self.bottomRightPlot.set_xlim(xmin=min(x))
self.bottomRightPlot.set_xlim(xmax=max(x))
self.topLeftCanvas.draw()
self.topRightCanvas.draw()
self.bottomLeftCanvas.draw()
self.bottomRightCanvas.draw()
def onSelectSParameter(self,toP,fromP):
self.buttons[self.toPort-1][self.fromPort-1].config(relief=tk.RAISED)
self.toPort = toP
self.fromPort = fromP
self.buttons[self.toPort-1][self.fromPort-1].config(relief=tk.SUNKEN)
self.delay.SetValueFromString(str(0))
self.delayViewerProperty.onUntouched(None)
self.PlotSParameter()
def onAutoscale(self):
self.plt.autoscale(True)
self.f.canvas.draw()
def onUnwrap(self):
fr=self.sp.FrequencyResponse(self.toPort,self.fromPort)
ir=fr.ImpulseResponse()
if ir is not None:
idx = ir.Values('abs').index(max(ir.Values('abs')))
TD = ir.Times()[idx] # the time of the main peak
else:
TD=0.
self.delay.SetValueFromString(str(TD))
self.delayViewerProperty.onUntouched(None)
def onDelayEntered(self):
self.topRightPlot.cla()
fr=self.sp.FrequencyResponse(self.toPort,self.fromPort)
TD = self.delay.GetValue()
fr=fr._DelayBy(-TD)
lw=[min(1.,math.sqrt(w))*1.5 for w in fr.Response('mag')]
y=fr.Response('deg')
x=fr.Frequencies('GHz')
fastway=True
if self.variableLineWidth.get():
if fastway:
segments = [[[x[i],y[i]],[x[i+1],y[i+1]]] for i in range(len(x)-1)]
slw=lw[:-1]
lc = LineCollection(segments, linewidths=slw,color='blue')
self.topRightPlot.add_collection(lc)
else:
for i in range(len(x)-1):
if self.logScale.get():
self.topRightPlot.semilogx(x[i:i+2],y[i:i+2],linewidth=lw[i],color='blue')
else:
self.topRightPlot.plot(x[i:i+2],y[i:i+2],linewidth=lw[i],color='blue')
else:
if self.logScale.get():
self.topRightPlot.semilogx(x,y)
else:
self.topRightPlot.plot(x,y)
self.topRightPlot.set_xlim(xmin=min(x))
self.topRightPlot.set_xlim(xmax=max(x))
self.topRightPlot.set_ylim(ymin=min(y)-1)
self.topRightPlot.set_ylim(ymax=max(y)+1)
self.topRightPlot.set_ylabel('phase (degrees)',fontsize=10)
self.topRightPlot.set_xlabel('frequency ('+self.freqLabel+')',fontsize=10)
self.topRightCanvas.draw()
def onReferenceImpedanceEntered(self):
self.sp.SetReferenceImpedance(self.referenceImpedance.GetValue())
self.PlotSParameter()
def onReadSParametersFromFile(self):
import SignalIntegrity.Lib as si
filename=AskOpenFileName(filetypes=[('s-parameter files', ('*.s*p'))],
initialdir=self.fileparts.AbsoluteFilePath(),
parent=self)
if filename is None:
return
self.fileparts=FileParts(filename)
if self.fileparts.fileext=='':
return
self.title('S-parameters: '+self.fileparts.FileNameTitle())
self.sp=si.sp.SParameterFile(filename)
self.referenceImpedance.SetValueFromString(str(self.sp.m_Z0))
self.referenceImpedanceProperty.propertyString.set(self.referenceImpedance.PropertyString(stype='entry'))
for widget in self.sButtonsFrame.winfo_children():
widget.destroy()
numPorts=self.sp.m_P
self.buttonLabels=[['s'+str(toP+1)+str(fromP+1) for fromP in range(numPorts)] for toP in range(numPorts)]
self.buttons=[]
for toP in range(numPorts):
buttonrow=[]
rowFrame=tk.Frame(self.sButtonsFrame)
rowFrame.pack(side=tk.TOP,expand=tk.NO,fill=tk.NONE)
for fromP in range(numPorts):
thisButton=tk.Button(rowFrame,text=self.buttonLabels[toP][fromP],command=lambda x=toP+1,y=fromP+1: self.onSelectSParameter(x,y))
thisButton.pack(side=tk.LEFT,fill=tk.NONE,expand=tk.NO)
buttonrow.append(thisButton)
self.buttons.append(buttonrow)
self.fromPort = 1
self.toPort = 1
self.buttons[self.toPort-1][self.fromPort-1].config(relief=tk.SUNKEN)
self.PlotSParameter()
def onWriteSParametersToFile(self):
ports=self.sp.m_P
extension='.s'+str(ports)+'p'
filename=AskSaveAsFilename(filetypes=[('s-parameters', extension)],
defaultextension=extension,
initialdir=self.fileparts.AbsoluteFilePath(),
initialfile=self.fileparts.FileNameWithExtension(extension),
parent=self)
if filename is None:
return
self.fileparts=FileParts(filename)
self.sp.WriteToFile(filename,'R '+str(self.sp.m_Z0))
def onResample(self):
import SignalIntegrity.Lib as si
self.sp=self.sp.Resample(si.fd.EvenlySpacedFrequencyList(
SignalIntegrity.App.Project['CalculationProperties.EndFrequency'],
SignalIntegrity.App.Project['CalculationProperties.FrequencyPoints']))
self.PlotSParameter()
def onCalculationProperties(self):
self.parent.onCalculationProperties()
def onMatplotlib2TikZ(self):
filename=AskSaveAsFilename(parent=self,filetypes=[('tex', '.tex')],
defaultextension='.tex',
initialdir=self.fileparts.AbsoluteFilePath(),
initialfile=self.fileparts.filename+'Magnitude.tex')
if filename is None:
return
try:
PlotTikZ(filename,self.topLeftFigure)
except:
messagebox.showerror('Export LaTeX','LaTeX could not be generated or written ')
fp=FileParts(filename.replace('Magnitude.tex', ''))
filename=fp.filename
filename=AskSaveAsFilename(parent=self,filetypes=[('tex', '.tex')],
defaultextension='.tex',
initialdir=self.fileparts.AbsoluteFilePath(),
initialfile=filename+'Phase.tex')
if filename is None:
return
try:
PlotTikZ(filename,self.topRightFigure)
except:
messagebox.showerror('Export LaTeX','LaTeX could not be generated or written ')
fp=FileParts(filename.replace('Phase.tex', ''))
filename=fp.filename
filename=AskSaveAsFilename(parent=self,filetypes=[('tex', '.tex')],
defaultextension='.tex',
initialdir=self.fileparts.AbsoluteFilePath(),
initialfile=filename+'ImpulseResponse.tex')
if filename is None:
return
try:
PlotTikZ(filename,self.bottomLeftFigure)
except:
messagebox.showerror('Export LaTeX','LaTeX could not be generated or written ')
fp=FileParts(filename.replace('ImpulseResponse.tex', ''))
filename=fp.filename
filename=AskSaveAsFilename(parent=self,filetypes=[('tex', '.tex')],
defaultextension='.tex',
initialdir=self.fileparts.AbsoluteFilePath(),
initialfile=filename+'StepResponse.tex')
if filename is None:
return
try:
PlotTikZ(filename,self.bottomRightFigure)
except:
messagebox.showerror('Export LaTeX','LaTeX could not be generated or written ')
def onHelp(self):
if Doer.helpKeys is None:
messagebox.showerror('Help System','Cannot find or open this help element')
return
Doer.helpKeys.Open('sec:S-parameter-Viewer')
def onControlHelp(self):
Doer.inHelp=True
self.config(cursor='question_arrow')
def onEscape(self):
Doer.inHelp=False
self.config(cursor='left_ptr')
def onEnforcePassivity(self):
self.sp.EnforcePassivity()
self.PlotSParameter()
def onEnforceCausality(self):
self.sp.EnforceCausality()
self.PlotSParameter()
def onWaveletDenoise(self):
self.sp.WaveletDenoise()
self.PlotSParameter()
