THD calculation seems off by a lot

[ftab]

I'm having trouble getting a good THD / THD+N calculation from either THD() or THDN().

If you have a moment, could you take a look at this sample code/data and tell me why?

I'm trying to get a distortion figure on a 1khz sine wave on a device where an HP 8920b told me the THD was about 0.056%

However I am fresh out of 8920s to share and need to concoct as good a result as I can using...an 8 bit USB scope and a Raspberry Pi with Python!

Output:

Frequency: 1000.729392 Hz
fundamental amplitude: 300737.188
36118.519 117.290 24.098 128.721 52.875 156.427 66.371 25.883 140.883 73.748 129.345 87.715 131.846 
THD: 12.387467%
0.9079618095454111

import matplotlib.pyplot as plt
import numpy as np
from waveform_analysis import thd

t = np.arange(0.0, 0.005, 1/1000000)
s = [1, 1, 2, 3, 3, 4, 4, 6, 7, 7, 9, 9, 10, 11, 12, 11, 12, 13, 14, 15, 16, 16, 17, 18, 18, 19, 20, 20, 22, 23, 24, 24, 25, 25, 25, 27, 28, 29, 29, 29, 30, 32, 31, 33, 33, 34, 35, 35, 36, 37, 38, 40, 39, 40, 41, 41, 42, 43, 43, 44, 45, 45, 46, 46, 48, 49, 49, 49, 50, 51, 52, 52, 53, 54, 55, 56, 56, 56, 57, 57, 59, 58, 58, 60, 61, 62, 62, 63, 64, 65, 65, 65, 66, 67, 68, 69, 68, 69, 71, 71, 72, 71, 72, 73, 73, 74, 74, 75, 77, 77, 78, 78, 79, 80, 80, 80, 80, 81, 82, 83, 83, 82, 83, 84, 86, 85, 85, 86, 87, 88, 88, 89, 89, 90, 91, 91, 91, 92, 94, 94, 93, 92, 94, 96, 94, 95, 95, 96, 97, 97, 98, 97, 99, 100, 100, 100, 101, 102, 103, 102, 101, 102, 103, 104, 104, 103, 104, 105, 106, 105, 105, 105, 106, 108, 107, 106, 108, 109, 110, 109, 110, 110, 111, 112, 111, 110, 111, 111, 112, 111, 112, 112, 112, 113, 112, 112, 113, 114, 115, 114, 114, 116, 116, 117, 116, 116, 117, 116, 117, 116, 116, 117, 117, 118, 115, 116, 118, 118, 118, 117, 117, 119, 120, 120, 119, 119, 119, 120, 120, 118, 119, 120, 120, 119, 119, 119, 120, 119, 120, 119, 120, 120, 121, 120, 120, 121, 121, 121, 121, 121, 121, 121, 121, 120, 120, 120, 121, 119, 119, 119, 119, 120, 120, 119, 120, 120, 121, 120, 120, 120, 120, 121, 120, 119, 119, 119, 120, 118, 118, 118, 118, 118, 118, 116, 116, 118, 118, 116, 116, 117, 117, 117, 116, 115, 116, 116, 117, 114, 114, 115, 115, 114, 112, 112, 113, 113, 113, 111, 111, 112, 112, 112, 110, 110, 111, 112, 111, 109, 110, 110, 110, 109, 108, 108, 107, 108, 106, 105, 105, 106, 106, 105, 103, 104, 105, 104, 104, 103, 103, 103, 102, 101, 101, 101, 101, 100, 99, 97, 98, 97, 96, 96, 95, 94, 95, 96, 93, 94, 94, 94, 93, 92, 92, 92, 92, 90, 89, 89, 89, 88, 87, 86, 86, 85, 86, 84, 83, 84, 84, 84, 82, 81, 81, 82, 82, 80, 79, 79, 78, 78, 76, 76, 76, 75, 74, 72, 71, 72, 72, 71, 70, 70, 70, 69, 69, 68, 67, 68, 67, 65, 65, 64, 63, 62, 62, 60, 60, 59, 60, 57, 57, 57, 57, 56, 55, 54, 54, 54, 53, 52, 51, 51, 50, 49, 48, 47, 47, 46, 46, 44, 43, 43, 42, 42, 41, 40, 40, 39, 39, 37, 37, 37, 36, 35, 34, 33, 33, 32, 31, 30, 29, 29, 28, 27, 26, 26, 26, 24, 24, 22, 22, 22, 21, 20, 20, 19, 18, 17, 16, 14, 14, 15, 13, 12, 12, 11, 10, 10, 9, 8, 7, 7, 5, 5, 5, 3, 3, 1, 0, 0, -1, -2, -2, -2, -3, -3, -4, -5, -6, -6, -8, -8, -9, -10, -11, -12, -13, -13, -14, -15, -16, -17, -18, -18, -19, -20, -20, -22, -20, -23, -23, -24, -23, -25, -26, -26, -28, -29, -28, -29, -31, -31, -31, -32, -34, -35, -34, -35, -35, -36, -37, -38, -39, -39, -40, -42, -42, -43, -42, -44, -46, -46, -47, -47, -48, -49, -49, -50, -50, -51, -52, -52, -52, -53, -55, -55, -56, -56, -56, -56, -59, -59, -59, -59, -61, -62, -62, -62, -64, -65, -65, -66, -66, -67, -67, -69, -68, -68, -69, -70, -72, -71, -72, -73, -73, -74, -75, -74, -76, -77, -78, -77, -78, -79, -80, -81, -80, -80, -82, -82, -83, -82, -83, -84, -85, -86, -84, -86, -87, -87, -89, -88, -88, -90, -90, -91, -91, -92, -92, -94, -94, -93, -94, -95, -95, -96, -95, -96, -97, -97, -98, -97, -98, -100, -100, -99, -100, -101, -102, -103, -103, -101, -103, -105, -104, -104, -103, -104, -105, -105, -105, -105, -106, -108, -107, -107, -108, -108, -110, -110, -110, -109, -111, -111, -110, -110, -110, -111, -112, -111, -111, -112, -112, -113, -112, -112, -113, -114, -115, -114, -114, -115, -117, -117, -116, -115, -116, -117, -118, -115, -115, -116, -116, -118, -116, -116, -117, -118, -119, -118, -118, -119, -120, -120, -119, -119, -120, -120, -121, -119, -119, -120, -120, -119, -118, -118, -119, -120, -119, -118, -119, -121, -120, -121, -120, -120, -122, -121, -121, -119, -121, -121, -120, -121, -120, -119, -121, -120, -119, -119, -119, -120, -120, -120, -119, -119, -120, -120, -120, -119, -120, -120, -120, -119, -119, -119, -119, -119, -118, -116, -118, -118, -118, -116, -116, -117, -118, -116, -117, -116, -117, -118, -116, -116, -116, -116, -117, -115, -114, -114, -114, -115, -114, -112, -113, -113, -113, -112, -111, -112, -112, -113, -111, -110, -112, -111, -111, -110, -110, -110, -110, -109, -108, -108, -108, -107, -108, -105, -105, -106, -105, -106, -104, -104, -105, -104, -104, -103, -103, -103, -103, -102, -100, -100, -101, -100, -99, -97, -98, -98, -98, -96, -95, -95, -96, -95, -94, -94, -95, -94, -94, -93, -91, -92, -92, -91, -89, -88, -88, -89, -88, -87, -87, -86, -86, -84, -84, -83, -83, -84, -82, -82, -81, -81, -81, -79, -78, -78, -79, -78, -76, -75, -75, -74, -75, -73, -72, -72, -72, -72, -70, -69, -69, -70, -69, -67, -66, -66, -66, -65, -64, -63, -64, -63, -63, -60, -60, -59, -58, -58, -57, -56, -56, -55, -56, -55, -53, -54, -53, -52, -51, -51, -50, -49, -49, -47, -47, -46, -45, -45, -43, -42, -42, -42, -42, -40, -40, -39, -39, -38, -36, -36, -36, -35, -34, -33, -33, -32, -31, -31, -29, -28, -28, -28, -26, -26, -25, -24, -24, -23, -22, -21, -22, -20, -19, -18, -17, -17, -16, -15, -15, -14, -13, -12, -12, -11, -11, -10, -9, -9, -7, -7, -6, -4, -5, -3, -3, -2, -1, 1, 1, 1, 3, 2, 2, 4, 5, 6, 6, 7, 7, 7, 8, 10, 10, 11, 12, 14, 15, 15, 16, 17, 17, 18, 19, 20, 21, 22, 22, 22, 22, 23, 24, 25, 25, 26, 27, 27, 28, 29, 30, 32, 32, 32, 33, 35, 36, 36, 36, 36, 37, 38, 39, 39, 40, 41, 41, 43, 42, 44, 45, 45, 46, 46, 48, 49, 50, 51, 50, 51, 52, 53, 54, 53, 54, 55, 55, 56, 56, 57, 58, 59, 60, 60, 62, 63, 63, 64, 63, 65, 65, 66, 66, 67, 67, 69, 69, 69, 69, 70, 71, 71, 72, 72, 73, 75, 76, 75, 75, 77, 79, 78, 78, 78, 80, 80, 80, 80, 81, 82, 83, 82, 82, 83, 85, 86, 86, 86, 87, 88, 90, 89, 89, 90, 89, 92, 91, 91, 92, 93, 93, 93, 93, 94, 94, 96, 95, 95, 97, 97, 98, 99, 98, 99, 100, 100, 100, 100, 102, 102, 102, 102, 102, 103, 104, 104, 103, 104, 105, 106, 105, 105, 106, 107, 108, 108, 108, 108, 109, 110, 110, 109, 110, 111, 111, 111, 110, 110, 112, 112, 112, 112, 112, 113, 113, 113, 113, 114, 116, 116, 115, 115, 116, 117, 117, 116, 116, 116, 117, 117, 116, 116, 117, 118, 116, 116, 117, 117, 119, 118, 118, 118, 119, 121, 120, 120, 119, 120, 120, 120, 119, 118, 119, 120, 119, 119, 119, 120, 120, 119, 119, 120, 121, 122, 120, 120, 121, 121, 122, 120, 119, 120, 121, 121, 119, 119, 120, 120, 120, 120, 119, 120, 121, 120, 120, 120, 121, 121, 120, 120, 119, 120, 120, 119, 118, 118, 119, 119, 118, 117, 118, 118, 118, 118, 118, 117, 117, 118, 117, 117, 118, 117, 117, 116, 115, 116, 116, 115, 115, 113, 115, 114, 114, 113, 113, 113, 114, 112, 113, 111, 112, 113, 112, 111, 110, 111, 111, 110, 109, 109, 109, 109, 108, 107, 108, 107, 108, 106, 105, 105, 105, 105, 104, 104, 104, 104, 104, 102, 102, 103, 102, 103, 100, 99, 100, 100, 99, 97, 97, 98, 97, 96, 96, 95, 96, 95, 96, 93, 93, 94, 94, 94, 92, 90, 92, 91, 90, 88, 88, 88, 88, 86, 85, 86, 86, 86, 84, 83, 83, 84, 83, 82, 81, 81, 81, 80, 78, 78, 78, 77, 77, 76, 74, 75, 75, 75, 72, 72, 72, 72, 71, 70, 70, 70, 69, 68, 67, 67, 66, 65, 65, 63, 63, 63, 62, 60, 59, 60, 60, 59, 57, 57, 57, 56, 56, 55, 53, 54, 53, 52, 52, 50, 50, 49, 48, 47, 46, 46, 46, 45, 43, 43, 43, 42, 42, 40, 40, 40, 39, 39, 37, 36, 36, 35, 35, 33, 32, 31, 31, 30, 29, 29, 28, 28, 27, 26, 25, 25, 24, 22, 22, 22, 22, 20, 20, 18, 17, 17, 17, 16, 15, 14, 13, 13, 12, 11, 11, 9, 9, 8, 8, 7, 6, 6, 5, 3, 2, 1, 1, 0, -1, -1, -2, -3, -4, -5, -6, -6, -6, -7, -7, -9, -10, -10, -10, -11, -13, -14, -14, -15, -15, -17, -18, -18, -18, -19, -20, -21, -22, -22, -23, -23, -24, -24, -25, -26, -28, -29, -28, -29, -30, -31, -32, -33, -33, -34, -35, -36, -36, -36, -38, -38, -39, -39, -39, -41, -41, -43, -43, -43, -44, -45, -46, -47, -47, -48, -50, -50, -50, -50, -51, -53, -53, -53, -53, -55, -55, -56, -56, -56, -58, -59, -60, -60, -61, -62, -62, -64, -63, -64, -65, -66, -66, -66, -67, -68, -70, -69, -68, -69, -71, -72, -72, -72, -73, -74, -75, -76, -76, -77, -78, -78, -79, -79, -79, -80, -81, -81, -80, -82, -82, -83, -83, -83, -84, -86, -86, -86, -87, -88, -88, -88, -89, -89, -91, -91, -90, -91, -91, -92, -94, -94, -93, -93, -94, -95, -95, -96, -97, -97, -98, -98, -97, -100, -100, -101, -100, -100, -102, -102, -103, -101, -102, -103, -104, -104, -103, -103, -105, -105, -106, -105, -105, -106, -107, -108, -107, -108, -110, -110, -110, -109, -109, -111, -111, -111, -109, -110, -111, -111, -111, -110, -112, -114, -113, -113, -113, -114, -115, -116, -116, -115, -115, -116, -116, -116, -115, -115, -117, -116, -115, -116, -117, -117, -117, -116, -116, -118, -119, -119, -119, -119, -119, -121, -120, -120, -119, -120, -120, -119, -120, -119, -118, -120, -119, -119, -119, -119, -120, -120, -120, -120, -121, -122, -121, -120, -121, -121, -121, -120, -120, -120, -120, -120, -119, -119, -119, -119, -121, -120, -119, -120, -120, -122, -120, -119, -121, -121, -121, -119, -119, -119, -119, -120, -119, -119, -118, -119, -119, -117, -117, -119, -118, -118, -117, -117, -119, -118, -117, -116, -116, -117, -118, -116, -115, -115, -115, -115, -114, -114, -114, -114, -114, -113, -113, -113, -114, -113, -112, -111, -111, -112, -112, -111, -109, -111, -111, -110, -109, -108, -109, -110, -108, -107, -107, -107, -107, -106, -106, -105, -106, -106, -105, -104, -104, -104, -105, -103, -102, -102, -103, -102, -100, -100, -100, -99, -99, -98, -97, -97, -97, -97, -95, -95, -95, -95, -96, -94, -93, -94, -94, -94, -92, -90, -91, -90, -90, -88, -87, -88, -87, -87, -85, -84, -85, -84, -85, -84, -82, -84, -84, -82, -81, -80, -81, -80, -80, -78, -77, -77, -76, -76, -74, -74, -74, -74, -74, -72, -71, -72, -71, -71, -69, -69, -69, -68, -67, -66, -66, -66, -65, -63, -61, -61, -61, -60, -59, -58, -59, -60, -58, -56, -56, -56, -56, -55, -54, -53, -53, -53, -52, -50, -50, -50, -48, -47, -46, -45, -45, -45, -44, -42, -42, -42, -42, -41, -40, -39, -39, -39, -36, -36, -36, -35, -34, -33, -32, -32, -31, -31, -29, -28, -28, -27, -26, -26, -25, -25, -23, -23, -22, -22, -21, -20, -20, -18, -17, -17, -16, -16, -15, -14, -13, -12, -11, -11, -10, -10, -9, -9, -7, -7, -6, -5, -4, -4, -2, -2, -1, 0, 1, 1, 1, 3, 3, 4, 6, 6, 6, 7, 8, 8, 9, 10, 10, 12, 12, 12, 14, 15, 16, 16, 17, 18, 18, 19, 21, 21, 22, 22, 22, 23, 25, 24, 25, 27, 28, 28, 29, 30, 30, 31, 32, 32, 33, 34, 35, 36, 36, 36, 37, 38, 39, 40, 40, 40, 41, 42, 42, 44, 45, 45, 47, 48, 47, 48, 49, 50, 50, 50, 51, 52, 53, 53, 54, 54, 55, 56, 56, 56, 58, 59, 61, 60, 61, 62, 63, 63, 64, 64, 65, 65, 66, 66, 66, 67, 69, 70, 70, 70, 72, 72, 73, 72, 74, 75, 75, 76, 75, 77, 78, 78, 78, 78, 79, 80, 80, 81, 80, 81, 82, 83, 84, 83, 84, 86, 86, 87, 86, 87, 89, 90, 89, 89, 90, 91, 92, 91, 91, 92, 94, 93, 93, 93, 95, 95, 96, 96, 96, 97, 98, 100, 98, 99, 100, 101, 101, 101, 101, 102, 102, 101, 102, 102, 103, 104, 103, 104, 104, 105, 106, 106, 106, 107, 107, 109, 108, 108, 109, 109, 111, 109, 109, 110, 111, 111, 110, 110, 110, 111, 112, 112, 112, 113, 113, 115, 114, 114, 115, 116, 116, 115, 115, 116, 116, 117, 115, 115, 116, 116, 116, 116, 116, 117, 117, 118, 117, 118, 119, 119, 119, 119, 119, 120, 120, 120, 119, 119, 119, 120, 118, 118, 118, 119, 119, 119, 118, 119, 121, 121, 120, 119, 120, 122, 121, 121, 120, 121, 121, 121, 119, 118, 120, 121, 120, 120, 119, 120, 121, 120, 120, 119, 120, 122, 121, 120, 119, 120, 121, 119, 118, 119, 120, 119, 119, 118, 118, 118, 118, 118, 118, 118, 118, 118, 117, 117, 119, 117, 119, 117, 117, 117, 117, 117, 117, 115, 116, 115, 115, 114, 114, 115, 114, 114, 112, 112, 113, 113, 112, 112, 112, 112, 113, 112, 111, 110, 111, 111, 110, 109, 109, 110, 109, 107, 106, 107, 108, 107, 105, 105, 105, 106, 105, 103, 103, 104, 105, 104, 102, 101, 102, 102, 102, 100, 99, 100, 100, 98, 97, 97, 97, 97, 97, 95, 95, 96, 96, 95, 94, 94, 94, 94, 93, 92, 92, 91, 90, 89, 87, 87, 86, 88, 86, 85, 85, 85, 85, 84, 84, 84, 83, 83, 81, 81, 82, 80, 79, 78, 77, 77, 77, 76, 74, 74, 74, 73, 73, 72, 72, 72, 72, 71, 69, 69, 70, 68, 67, 65, 66, 65, 64, 64, 62, 62, 62, 60, 61, 59, 58, 59, 58, 57, 57, 57, 56, 56, 53, 54, 53, 53, 51, 50, 49, 48, 48, 48, 46, 46, 46, 45, 45, 43, 43, 43, 41, 41, 40, 39, 40, 39, 38, 36, 35, 35, 35, 34, 33, 32, 31, 31, 30, 28, 28, 28, 27, 26, 24, 25, 24, 23, 22, 21, 21, 20, 20, 20, 18, 18, 17, 16, 15, 14, 14, 13, 12, 11, 10, 10, 9, 9, 8, 6, 6, 5, 5, 4, 3, 2, 1, 0, -1, -1, -2, -2, -3, -3, -4, -6, -6, -6, -7, -7, -9, -9, -10, -11, -11, -13, -14, -15, -15, -16, -18, -18, -18, -20, -20, -20, -22, -22, -22, -24, -24, -24, -26, -25, -26, -28, -28, -30, -30, -32, -33, -32, -33, -34, -34, -36, -36, -36, -38, -39, -39, -39, -40, -40, -40, -42, -43, -44, -44, -45, -46, -47, -47, -47, -48, -50, -50, -51, -51, -52, -53, -53, -52, -54, -55, -56, -56, -56, -57, -59, -60, -60, -60, -62, -63, -64, -64, -64, -64, -66, -67, -66, -66, -67, -68, -69, -68, -69, -70, -71, -73, -72, -72, -74, -74, -77, -75, -76, -77, -79, -79, -79, -78, -80, -80, -80, -80, -79, -81, -83, -83, -83, -84, -86, -87, -87, -86, -87, -89, -90, -90, -89, -90, -91, -91, -91, -91, -92, -93, -93, -94, -93, -94, -96, -96, -96, -96, -97, -99, -99, -98, -98, -99, -100, -100, -101, -100, -101, -103, -102, -101, -102, -102, -104, -104, -104, -103, -105, -106, -106, -106, -107, -108, -108, -107, -107, -108, -109, -110, -108, -108, -110, -110, -111, -109, -110, -110, -111, -112, -111, -111, -114, -114, -114, -114, -113, -115, -115, -115, -114, -114, -116, -116, -117, -114, -114, -116, -116, -117, -115, -116, -118, -117, -117, -117, -117, -118, -119, -119, -118, -118, -119, -119, -119, -118, -118, -119, -119, -119, -119, -118, -120, -119, -120, -119, -119, -121, -120, -120, -120, -120, -121, -121, -121, -119, -119, -121, -120, -120, -120, -120, -121, -120, -119, -119, -119, -120, -120, -119, -119, -119, -121, -120, -119, -119, -120, -121, -120, -119, -119, -120, -119, -119, -118, -118, -118, -119, -118, -118, -117, -117, -118, -117, -116, -117, -118, -118, -117, -116, -116, -116, -117, -115, -114, -116, -115, -115, -113, -114, -114, -114, -114, -113, -112, -114, -112, -113, -111, -111, -113, -112, -112, -111, -110, -111, -110, -110, -108, -108, -109, -108, -108, -106, -106, -107, -107, -106, -105, -106, -107, -106, -105, -104, -104, -104, -103, -102, -102, -102, -103, -100, -99, -98, -99, -100, -98, -97, -96, -98, -98, -96, -96, -95, -95, -96, -94, -94, -93, -94, -93, -91, -90, -89, -89, -90, -89, -87, -88, -87, -88, -87, -86, -84, -85, -85, -84, -83, -83, -82, -83, -81, -80, -79, -79, -79, -78, -76, -77, -76, -77, -74, -74, -75, -74, -74, -73, -72, -73, -72, -71, -70, -68, -68, -68, -67, -65, -64, -65, -65, -64, -62, -62, -62, -61, -60, -58, -59, -59, -58, -58, -56, -56, -56, -55, -54, -53, -53, -53, -51, -50, -49, -49, -49, -47, -46, -45, -46, -45, -45, -44, -42, -42, -42, -41, -40, -39, -38, -38, -38, -37, -35, -35, -35, -32, -32, -31, -31, -31, -30, -29, -28, -28, -28, -27, -26, -25, -24, -23, -23, -22, -21, -21, -20, -18, -18, -16, -16, -16, -15, -14, -14, -12, -12, -11, -10, -10, -9, -9, -8, -7, -7, -6, -4, -3, -3, -1, -1, 0, 1, 1, 3, 3, 4, 4, 5, 6, 6, 7, 7, 8, 9, 10, 10, 11, 12, 13, 14, 15, 16, 17, 17, 18, 19, 20, 20, 21, 22, 23, 22, 23, 24, 25, 25, 26, 27, 27, 29, 29, 30, 31, 32, 34, 33, 33, 35, 36, 36, 37, 37, 38, 39, 39, 40, 41, 41, 42, 44, 44, 44, 46, 46, 47, 48, 48, 49, 49, 50, 52, 51, 53, 53, 54, 54, 54, 55, 56, 56, 57, 57, 60, 60, 60, 61, 62, 62, 65, 65, 64, 65, 65, 66, 66, 66, 67, 68, 69, 69, 70, 71, 72, 73, 73, 73, 75, 76, 76, 76, 76, 78, 79, 78, 78, 79, 79, 80, 81, 80, 81, 82, 83, 84, 83, 84, 86, 86, 87, 86, 87, 89, 90, 90, 89, 88, 91, 91, 91, 91, 92, 94, 94, 93, 93, 94, 96, 95, 97, 97, 97, 98, 99, 99, 98, 99, 100, 101, 100, 100, 101, 102, 101, 102, 101, 103, 104, 104, 104, 104, 106, 107, 106, 107, 106, 107, 109, 109, 108, 108, 109, 110, 110, 108, 110, 110, 111, 110, 111, 110, 112, 113, 112, 112, 113, 115, 114, 114, 114, 114, 116, 116, 115, 115, 114, 115, 116, 115, 116, 115, 116, 117, 115, 115, 118, 118, 118, 118, 118, 120, 119, 120, 118, 119, 118, 120, 119, 118, 119, 119, 119, 118, 118, 118, 120, 121, 120, 119, 119, 121, 121, 121, 120, 120, 122, 121, 121, 120, 120, 121, 120, 120, 119, 120, 121, 119, 120, 120, 120, 121, 121, 120, 120, 121, 121, 120, 119, 120, 121, 121, 120, 118, 119, 119, 120, 118, 118, 118, 118, 119, 118, 118, 118, 119, 119, 118, 118, 118, 118, 118, 116, 117, 116, 117, 117, 115, 114, 115, 115, 115, 113, 113, 114, 114, 114, 113, 113, 113, 114, 114, 112, 113, 113, 113, 113, 110, 110, 111, 110, 109, 109, 108, 108, 108, 108, 107, 106, 108, 107, 106, 105, 106, 107, 106, 105, 105, 104, 104, 103, 103, 102, 102, 102, 101, 99, 99, 100, 100, 98, 97, 97, 98, 98, 97, 97, 96, 96, 95, 95, 94, 94, 94, 94, 92, 90, 89, 90, 90, 88, 88, 88, 87, 89, 86, 86, 87, 86, 86, 84, 84, 84, 84, 83, 81, 81, 81, 80, 79, 77, 77, 78, 77, 76, 74, 75, 75, 74, 74, 72, 73, 72, 72, 71, 69, 69, 69, 69, 68, 66, 66, 65, 65, 63, 63, 62, 62, 61, 61, 60, 59, 59, 59, 57, 56, 57, 56, 55, 54, 54, 53, 53, 52, 50, 50, 48, 49, 48, 47, 46, 46, 46, 45, 44, 43, 43, 42, 42, 40, 40, 40, 38, 38, 37, 36, 35, 35, 33, 32, 32, 32, 31, 31, 29, 28, 28, 26, 26, 25, 25, 24, 24, 23, 22, 22, 21, 20, 20, 17, 17, 17, 16, 15, 14, 14, 13, 12, 12, 10, 10, 10, 9, 8, 8, 6, 6, 4, 4, 2, 2, 1, 0, -1, -2, -2, -2, -4, -4, -4, -5, -6, -7, -7, -8, -9, -10, -10, -10, -12, -13, -14, -15, -15, -16, -17, -18, -18, -19, -19, -21, -20, -23, -22, -23, -24, -24, -25, -26, -27, -28, -28, -30, -30, -30, -32, -33, -33, -33, -34, -36, -36, -36, -36, -38, -39, -40, -40, -40, -42, -43, -43, -44, -44, -45, -46, -46, -47, -47, -49, -49, -49, -50, -51, -52, -52, -52, -53, -54, -55, -56, -56, -57, -58, -59, -60, -61, -60, -61, -61, -64, -63, -63, -64, -65, -67, -66, -66, -68, -68, -69, -69, -70, -70, -72, -72, -72, -73, -74, -75, -75, -76, -76, -78, -78, -78, -78, -78, -80, -80, -80, -81, -81, -83, -83, -84, -83, -83, -86, -86, -87, -86, -88, -88, -88, -90, -88, -89, -91, -91, -91, -91, -91, -93, -93, -92, -92, -93, -96, -96, -97, -97, -96, -98, -98, -98, -98, -100, -100, -100, -100, -100, -101, -103, -102, -102, -102, -102, -104, -104, -103, -104, -105, -107, -107, -107, -106, -108, -109, -108, -108, -108, -109, -110, -109, -108, -108, -109, -110, -109, -110, -110, -111, -113, -111, -113, -113, -114, -115, -114, -113, -115, -114, -115, -114, -114, -114, -115, -115, -114, -114, -115, -115, -117, -115, -116, -117, -118, -119, -117, -117, -119, -119, -119, -117, -118, -119, -119, -119, -118, -118, -118, -119, -119, -118, -119, -119, -120, -120, -119, -120, -120, -121, -121, -120, -120, -121, -121, -120, -119, -119, -120, -120, -119, -119, -119, -120, -120, -119, -119, -120, -120, -120, -120, -120, -120, -122, -120, -120, -119, -120, -121, -119, -119, -119, -118, -119, -119, -118, -118, -118, -120, -118, -118, -118, -117, -118, -118, -118, -118, -117, -118, -116, -116, -117, -116, -116, -115, -114, -115, -114, -115, -113, -114, -114, -114, -114, -112, -113, -114, -113, -113, -112, -111, -112, -112, -112, -110, -109, -111, -110, -109, -108, -108, -108, -108, -108, -107, -106, -108, -107, -107, -106, -106, -106, -106, -106, -103, -103, -104, -103, -103, -101, -101, -101, -101, -100, -99, -98, -100, -99, -98, -98, -98, -98, -97, -95, -95, -95, -96, -95, -94, -93, -93, -93, -92, -90, -90, -89, -91, -89, -88, -87, -88, -88, -87, -85, -85, -84, -86, -84, -83, -82, -83, -82, -80, -80, -80, -80, -80, -78, -78, -78, -77, -77, -74, -74, -74, -74, -74, -73, -72, -73, -71, -71, -68, -69, -68, -68, -68, -66, -65, -65, -64, -64, -63, -62, -63, -61, -61, -59, -60, -59, -60, -58, -56, -56, -56, -55, -55, -52, -52, -52, -50, -51, -50, -48, -48, -47, -47, -46, -46, -45, -45, -44, -42, -43, -42, -41, -40, -39, -38, -38, -36, -36, -35, -34, -34, -34, -33, -32, -32, -31, -31, -30, -29, -28, -28, -26, -25, -25, -24, -24, -23, -22, -20, -20, -19, -18, -17, -17, -16, -16, -15, -14, -14, -13, -12, -12, -11, -10, -10, -8, -8, -7, -6, -5, -5, -3, -2, -2, -2, 0, 0, 2, 2, 2, 3, 5, 4, 5, 6, 6, 7, 8, 9, 10, 10, 11, 12, 14, 14, 14, 15, 16, 17, 18, 18, 18, 19, 21, 21, 22, 22, 23, 23, 25, 25, 26, 28, 28, 29, 29, 30, 30, 32, 32, 33, 33, 34, 35, 36, 37, 36, 38, 38, 40, 40, 40, 42, 43, 44, 44, 44, 46, 46, 47, 47, 48, 48, 49, 50, 50, 50, 52, 52, 54, 53, 54, 56, 56, 57, 57, 57, 59, 61, 61, 61, 62, 61, 64, 64, 63, 64, 65, 66, 67, 67, 67, 69, 69, 70, 69, 70, 72, 73, 73, 73, 74, 75, 76, 76, 76, 77, 78, 78, 78, 78, 79, 80, 81, 81, 81, 83, 84, 84, 84, 84, 86, 86, 87, 87, 87, 88, 89, 89, 89, 88, 90, 91, 92, 91, 92, 94, 94, 94, 94, 94, 96, 97, 97, 96, 97, 98, 99, 98, 98, 99, 100, 101, 101, 100, 100, 102, 103, 102, 102, 103, 104, 105, 104, 105, 106, 107, 107, 106, 107, 107, 109, 108, 108, 107, 109, 109, 109, 108, 109, 110, 110, 111, 111, 110, 111, 113, 113, 112, 113, 114, 115, 114, 113, 114, 116, 115, 114, 114, 114, 115, 116, 114, 114, 116, 117, 118, 117, 117, 118, 118, 118, 118, 118, 120, 120, 120, 118, 118, 118, 119, 119, 118, 117, 119, 119, 118, 118, 118, 120, 120, 121, 120, 120, 121, 121, 121, 120, 120, 121, 120, 120, 120, 118, 120, 119, 119, 119, 120, 121, 121, 121, 120, 121, 122, 122, 121, 120, 121, 121, 120, 120, 119, 119, 120, 119, 118, 119, 118, 120, 119, 118, 118, 119, 120, 120, 118, 118, 119, 120, 118, 118, 117, 118, 118, 116, 116, 116, 116, 117, 116, 115, 116, 116, 115, 115, 113, 115, 116, 115, 113, 114, 114, 114, 114, 112, 111, 112, 113, 111, 110, 109, 110, 110, 110, 109, 108, 109, 109, 108, 107, 107, 108, 107, 106, 106, 105, 106, 106, 105, 105, 104, 105, 104, 103, 102, 101, 101, 101, 100, 99, 99, 100, 100, 98, 97, 98, 98, 97, 97, 96, 96, 96, 96, 94, 93, 93, 93, 92, 91, 91, 90, 91, 89, 88, 88, 87, 88, 87, 86, 87, 86, 86, 84, 84, 84, 83, 83, 81, 80, 80, 80, 80, 78, 77, 77, 78, 77, 76, 75, 75, 74, 74, 73, 72, 73, 73, 71, 69, 68, 69, 68, 66, 65, 65, 65, 64, 64, 62, 62, 62, 62, 62, 60, 60, 60, 58, 59, 57, 57, 56, 55, 55, 53, 52, 52, 51, 50, 49, 48, 48, 48, 47, 46, 46, 46, 46, 44, 44, 42, 41, 42, 40, 39, 39, 38, 37, 36, 35, 35, 35, 33, 33, 33, 32, 31, 30, 30, 29, 29, 27, 26, 26, 25, 25, 24, 22, 22, 21, 20, 20, 19, 17, 17, 17, 16, 15, 14, 14, 14, 12, 12, 11, 10, 10, 9, 8, 7, 6, 6, 5, 4, 3, 2, 1, 0, 0, -1, -1, -2, -3, -4, -4, -5, -6, -6, -7, -8, -9, -9, -10, -11, -12, -13, -13, -14, -15, -16, -17, -17, -18, -19, -20, -21, -20, -22, -22, -23, -25, -25, -25, -26, -26, -27, -28, -29, -29, -30, -32, -32, -33, -33, -34, -35, -35, -37, -36, -38, -39, -39, -40, -40, -42, -43, -43, -44, -44, -44, -46, -47, -47, -47, -48, -50, -49, -51, -50, -51, -53, -53, -53, -54, -55, -56, -56, -57, -58, -58, -60, -60, -60, -60, -61, -62, -63, -63, -63, -65, -66, -66, -66, -66, -68, -70, -69, -70, -71, -71, -73, -72, -73, -73, -74, -76, -74, -76, -76, -78, -78, -78, -78, -79, -80, -81, -81, -80, -82, -83, -84, -84, -84, -86, -86, -87, -86, -86, -88, -88, -89, -88, -88, -90, -90, -92, -90, -91, -92, -93, -94, -94, -94, -95, -96, -96, -96, -97, -98, -99, -98, -98, -99, -100, -100, -100, -99, -100, -102, -101, -101, -102, -102, -104, -104, -104, -104, -105, -106, -107, -107, -107, -108, -108, -107, -107, -107, -109, -109, -109, -109, -109, -110, -111, -110, -111, -111, -112, -112, -113, -112, -113, -113, -115, -114, -114, -114, -114, -114, -115, -114, -114, -114, -116, -115, -114, -116, -116, -118, -116, -117, -118, -118, -119, -118, -118, -118, -118, -119, -118, -118, -118, -118, -119, -118, -117, -118, -119, -120, -118, -118, -120, -120, -121, -120, -120, -120, -122, -122, -120, -121, -121, -120, -120, -119, -119, -120, -120, -120, -119, -119, -120, -120, -121, -120, -120, -122, -121, -121, -120, -120, -121, -121, -120, -119, -120, -120, -118, -119, -118, -118, -120, -118, -119, -118, -118, -120, -118, -119, -118, -118, -119, -119, -117, -117, -117, -118, -118, -116, -115, -116, -116, -114, -115, -114, -114, -115, -115, -114, -114, -114, -114, -114, -114, -113, -114, -114, -112, -111, -112, -112, -112, -111, -109, -110, -109, -110, -109, -108, -108, -109, -110, -108, -107, -106, -108, -108, -106, -106, -106, -106, -106, -104, -103, -104, -103, -102, -101, -100, -101, -101, -100, -99, -99, -100, -100, -99, -97, -98, -98, -96, -97, -95, -95, -96, -95, -94, -92, -92, -93, -92, -91, -90, -90, -90, -89, -89, -88, -88, -89, -87, -87, -86, -85, -86, -84, -84, -82, -83, -83, -82, -81, -79, -80, -80, -78, -77, -76, -76, -77, -76, -74, -75, -74, -75, -73, -72, -71, -72, -70, -69, -69, -68, -67, -68, -66, -65, -65, -63, -64, -63, -62, -62, -62, -60, -60, -59, -58, -58, -58, -56, -55, -56, -55, -54, -54, -52, -51, -51, -50, -50, -49, -49, -48, -48, -46, -45, -45, -45, -45, -42, -42, -42, -41, -41, -39, -38, -38, -38, -37, -35, -34, -34, -33, -33, -31, -31, -30, -30, -30, -28, -28, -27, -26, -26, -24, -24, -24, -23, -22, -20, -20, -20, -19, -17, -17, -17, -15, -15, -14, -14, -13, -13, -12, -11, -10, -10, -10, -8, -8, -6, -5, -5, -4, -2, -2, -2, 0]
thd.THD(s, 1000000)
thdn = thd.THDN(s, 1000000)
print(thdn)
fig, ax = plt.subplots()
ax.plot(t, s)
ax.set(xlabel='time (s)', ylabel='counts', title='Channel 1')
ax.grid()
plt.show()