# Output results errors = [] deltaInput = [] deltaSnap = [] deltaConverge = [] # Simulate the transmission simulationData = sim.simulateTransmission(inputFile, logDir, predictionInterval, samplingInterval, heartbeat, drThreshold, delay, jitter, packetLoss) inputData = simulationData[0] deltaInput = sim.findDistanceBetweenSamples(inputData, jumpThreshold, spacing) for reconThreshold in reconThresholds: convergeTxData = sim.simulateLinearConvergence(simulationData[4], logDir, "_ex1", samplingInterval, interpolationType, reconThreshold) snapRxData = sim.simulateSnapRecon(simulationData[4], logDir, "_ex1", samplingInterval) errors.append( sim.findDistanceError(inputData, convergeTxData[0]) ) deltaSnap.append( sim.findDistanceBetweenSamples(snapRxData, jumpThreshold, spacing) ) deltaConverge.append( sim.findDistanceBetweenSamples(convergeTxData, jumpThreshold, spacing) ) # Prepare the data for plotting temp = [] for error in errors: temp.append( scipy.stats.mean(error) ) errors = temp #pylab.figure(1) #pylab.plot(reconThresholds, errors) pylab.figure(2) n, bins, patches = pylab.hist(deltaInput, 50)
# Simulate the transmission and part of the reconstruction
startTimeTotal = time.time()
testData = []
staticJumps = [ [ [], [], [] ],
[ [], [], [] ],
[ [], [], [] ],
[ [], [], [] ] ]
for inputFile in inputFiles:
simulationData = sim.simulateTransmission(inputFile, logDir, predictionInterval,
samplingInterval, heartbeat, drThreshold,
delay, jitter, packetLoss)
simNumber = inputFile.split('/')[-2][-1]
snapReconstruction = sim.simulateSnapRecon(simulationData[4], logDir, simNumber, samplingInterval)[0]
snapLimitReconstruction = sim.simulateSnapLimitRecon(simulationData[4], logDir, simNumber, samplingInterval,
interpolationType, closeThreshold, snapLimit)[0]
convergedReconstruction = sim.simulateLinearConvergence(simulationData[4], logDir, simNumber, samplingInterval,
interpolationType, closeThreshold)[0]
testData.append([simulationData[0], simulationData[4], snapReconstruction,
snapLimitReconstruction, convergedReconstruction])
# Find intra-sample jump
threshold = 0
spacing = 1
jumpDataSets = []
jumpDataSets.append( sim.findDistanceBetweenSamples(simulationData[0], threshold, spacing) )
jumpDataSets.append( sim.findDistanceBetweenSamples(snapReconstruction, threshold, spacing) )
jumpDataSets.append( sim.findDistanceBetweenSamples(snapLimitReconstruction, threshold, spacing) )
jumpDataSets.append( sim.findDistanceBetweenSamples(convergedReconstruction, threshold, spacing) )
for i, data in enumerate(jumpDataSets):
mean, std, median = sim.calculateStats(data)
staticJumps[i][0].append(mean)
staticJumps[i][1].append(std)
jitters = range(5, 101, 5)
jitters.insert(0, 1)
packetLoss = 0
varJitterResults = []
for jitter in jitters:
for i, inputFile in enumerate(inputFiles):
# Simulate the transmission
simulationData = sim.simulateTransmission(inputFile, logDir, predictionInterval,
samplingInterval, heartbeat, drThreshold,
delay, jitter, packetLoss)
# Simulate the reconstruction
simNumber = inputFile.split('/')[-2][-1]
reconstructionSnap = sim.simulateSnapRecon(simulationData[4], logDir, simNumber, samplingInterval)[0]
reconstructionInt = sim.simulateLinearConvergence(simulationData[4], logDir, simNumber, samplingInterval,
interpolationType, reconThreshold)[0]
reconstructionFilter = sim.simulateFilterRecon(reconstructionSnap, logDir, simNumber, samplingInterval, coefficients)[0]
sim.amplifyData(reconstructionFilter)
s =napError = 1
intError = 1
filterError = 1
varJitterResults.append([snapError, intError, filterError])
# Second -> variable delay, fixed jitter, fixed packet loss
delays = range(10,401,10)
delays.insert(0,1)
jitter = 20
packetLoss = 0
varDelayResults = []
