Advertisement
acclivity

pyEquivSoundLevel

Mar 13th, 2021
193
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 8.68 KB | None | 0 0
  1.  
  2. # Compute Equivalent Sound Level (SPL)
  3. # Mike Kerry - 13-March-2020
  4.  
  5. import math
  6.  
  7. Lp = [67.32336646966162,
  8.  67.32336646966162,
  9.  69.9935171407445,
  10.  65.36501192192063,
  11.  63.49667762408437,
  12.  64.84263622147162,
  13.  63.64175067957247,
  14.  63.10697994124109,
  15.  63.23307543161612,
  16.  68.12800787111614,
  17.  67.63668598721365,
  18.  71.89437038368278,
  19.  70.40953147305076,
  20.  69.75656233762591,
  21.  69.0814815187725,
  22.  70.07610095905218,
  23.  70.8479041574523,
  24.  70.8404325602515,
  25.  70.47081462154567,
  26.  62.79695976154486,
  27.  56.165021815541614,
  28.  55.31001880674193,
  29.  62.93374301271796,
  30.  59.082944726472455,
  31.  61.72420310873101,
  32.  55.881092084692696,
  33.  62.30552747106572,
  34.  61.526031137288,
  35.  67.88645952632314,
  36.  65.3676689883366,
  37.  68.47341182066452,
  38.  64.89517285778184,
  39.  56.66148204140534,
  40.  57.30401284887815,
  41.  64.0120958532131,
  42.  65.40200906840025,
  43.  64.79262166350605,
  44.  69.68729780299549,
  45.  70.20518626086832,
  46.  69.79132202164664,
  47.  62.88110142172912,
  48.  65.73666581828397,
  49.  63.3769177710237,
  50.  60.86098554449683,
  51.  66.15340237411513,
  52.  66.42446742255484,
  53.  64.46709000513836,
  54.  68.1895199934836,
  55.  64.98426618376253,
  56.  65.54634705942196,
  57.  63.48012174724825,
  58.  63.02781756177382,
  59.  65.41152792815667,
  60.  65.40085905046551,
  61.  56.632426479316734,
  62.  59.29155731607538,
  63.  64.8824601146881,
  64.  64.76190390102522,
  65.  69.40690384406328,
  66.  68.30757442007692,
  67.  70.17333228777665,
  68.  68.50757522754799,
  69.  63.39930397412935,
  70.  65.65784119950544,
  71.  65.70544791692254,
  72.  63.06806189573864,
  73.  62.11888565194215,
  74.  62.84950187150918,
  75.  63.36953520266584,
  76.  68.58765847909596,
  77.  66.1642753541416,
  78.  66.11228963673442,
  79.  63.31567972469273,
  80.  58.70144213426909,
  81.  61.46959467874335,
  82.  61.94598868862434,
  83.  58.617801600202554,
  84.  62.179308842114025,
  85.  64.69684417979369,
  86.  64.62264177151313,
  87.  65.98924693502121,
  88.  61.85060900200202,
  89.  64.82719968238375,
  90.  69.83424417444755,
  91.  68.22189166626248,
  92.  67.84510679436752,
  93.  68.50853531956085,
  94.  67.85270623259387,
  95.  63.46641381340176,
  96.  55.62306410900811,
  97.  68.23159827973078,
  98.  65.57855594322733,
  99.  62.80686044489778,
  100.  59.69530951319552,
  101.  59.23301382738342,
  102.  62.466274555716446,
  103.  62.01519365725744,
  104.  61.49089128733455,
  105.  58.72314917640075,
  106.  70.26448551398036,
  107.  68.20394683040813,
  108.  70.4560554005297,
  109.  68.67097090693561,
  110.  68.49228856871692,
  111.  69.31423290341216,
  112.  68.94936575242964,
  113.  62.74657292194825,
  114.  58.619091321969506,
  115.  56.69279320884756,
  116.  56.62451888074863,
  117.  54.99024947793882,
  118.  58.58112192460169,
  119.  65.28647461432915,
  120.  65.54212915787434,
  121.  59.579214370350954,
  122.  65.05575615679689,
  123.  66.07515580463593,
  124.  64.64027104642534,
  125.  61.3466886002372,
  126.  61.88103805406414,
  127.  60.20248061566963,
  128.  55.85836879989958,
  129.  57.05920152636161,
  130.  56.11177233800258,
  131.  63.70908521976456,
  132.  65.57894028996395,
  133.  65.15596114414907,
  134.  65.56921450216043,
  135.  64.96999288658019,
  136.  50.76644369083102,
  137.  67.07661409032553,
  138.  64.93133719563896,
  139.  62.4539514017717,
  140.  67.70428031139994,
  141.  62.91311923435906,
  142.  61.80258750153895,
  143.  62.030321964338114,
  144.  68.23220239410021,
  145.  63.30054308830876,
  146.  67.94219171067066,
  147.  66.81718517462494,
  148.  61.458971328963294,
  149.  67.50253944183788,
  150.  57.31289497917814,
  151.  70.60477818510792,
  152.  65.55294076565086,
  153.  57.24597049485472,
  154.  68.82582307711594,
  155.  64.58162668532566,
  156.  62.872701907113196,
  157.  70.06537131421248,
  158.  66.32143646210767,
  159.  65.4886863912912,
  160.  64.97054189628867,
  161.  66.87656461841945,
  162.  66.16584130582459,
  163.  62.63168652021886,
  164.  72.2412316764658,
  165.  67.29878732863389,
  166.  64.75751145431248,
  167.  56.406947195484946,
  168.  71.97026260035784,
  169.  68.33770274572201,
  170.  57.883463570918245,
  171.  67.8308901755753,
  172.  69.19360002820311,
  173.  53.56379409059188,
  174.  66.73371931801356,
  175.  70.9798786602943,
  176.  52.47694811252207,
  177.  45.514560185053654,
  178.  47.72825133962793,
  179.  45.67585334494089,
  180.  46.22760332701799,
  181.  56.61233708999946,
  182.  46.88324190604364,
  183.  45.193173965071196,
  184.  53.48348797933126,
  185.  67.70217609099717,
  186.  67.25500981467864,
  187.  62.43797773195063,
  188.  65.54452144176976]
  189.  
  190. t = [0,
  191.      0.12498866213151928,
  192.      0.24997732426303856,
  193.      0.3749659863945578,
  194.      0.49995464852607713,
  195.      0.6249433106575963,
  196.      0.7499319727891156,
  197.      0.874920634920635,
  198.      0.9999092970521543,
  199.      1.1248979591836734,
  200.      1.2498866213151927,
  201.      1.374875283446712,
  202.      1.4998639455782312,
  203.      1.6248526077097507,
  204.      1.74984126984127,
  205.      1.8748299319727892,
  206.      1.9998185941043085,
  207.      2.124807256235828,
  208.      2.249795918367347,
  209.      2.3747845804988663,
  210.      2.4997732426303854,
  211.      2.624761904761905,
  212.      2.749750566893424,
  213.      2.8747392290249434,
  214.      2.9997278911564624,
  215.      3.124716553287982,
  216.      3.2497052154195014,
  217.      3.3746938775510205,
  218.      3.49968253968254,
  219.      3.624671201814059,
  220.      3.7496598639455785,
  221.      3.8746485260770975,
  222.      3.999637188208617,
  223.      4.124625850340136,
  224.      4.249614512471656,
  225.      4.374603174603175,
  226.      4.499591836734694,
  227.      4.624580498866213,
  228.      4.749569160997733,
  229.      4.874557823129252,
  230.      4.999546485260771,
  231.      5.12453514739229,
  232.      5.24952380952381,
  233.      5.374512471655329,
  234.      5.499501133786848,
  235.      5.624489795918367,
  236.      5.749478458049887,
  237.      5.874467120181406,
  238.      5.999455782312925,
  239.      6.124444444444444,
  240.      6.249433106575964,
  241.      6.374421768707483,
  242.      6.499410430839003,
  243.      6.624399092970521,
  244.      6.749387755102041,
  245.      6.87437641723356,
  246.      6.99936507936508,
  247.      7.1243537414965985,
  248.      7.249342403628118,
  249.      7.3743310657596375,
  250.      7.499319727891157,
  251.      7.624308390022676,
  252.      7.749297052154195,
  253.      7.8742857142857146,
  254.      7.999274376417234,
  255.      8.124263038548753,
  256.      8.249251700680272,
  257.      8.374240362811792,
  258.      8.499229024943311,
  259.      8.62421768707483,
  260.      8.74920634920635,
  261.      8.874195011337868,
  262.      8.999183673469387,
  263.      9.124172335600907,
  264.      9.249160997732426,
  265.      9.374149659863946,
  266.      9.499138321995465,
  267.      9.624126984126985,
  268.      9.749115646258504,
  269.      9.874104308390024,
  270.      9.999092970521541,
  271.      10.124081632653061,
  272.      10.24907029478458,
  273.      10.3740589569161,
  274.      10.49904761904762,
  275.      10.624036281179139,
  276.      10.749024943310658,
  277.      10.874013605442178,
  278.      10.999002267573696,
  279.      11.123990929705215,
  280.      11.248979591836735,
  281.      11.373968253968254,
  282.      11.498956916099774,
  283.      11.623945578231293,
  284.      11.748934240362813,
  285.      11.873922902494332,
  286.      11.99891156462585,
  287.      12.12390022675737,
  288.      12.248888888888889,
  289.      12.373877551020408,
  290.      12.498866213151928,
  291.      12.623854875283447,
  292.      12.748843537414967,
  293.      12.873832199546486,
  294.      12.998820861678006,
  295.      13.123809523809523,
  296.      13.248798185941043,
  297.      13.373786848072562,
  298.      13.498775510204082,
  299.      13.623764172335601,
  300.      13.74875283446712,
  301.      13.87374149659864,
  302.      13.99873015873016,
  303.      14.123718820861678,
  304.      14.248707482993197,
  305.      14.373696145124716,
  306.      14.498684807256236,
  307.      14.623673469387755,
  308.      14.748662131519275,
  309.      14.873650793650794,
  310.      14.998639455782314,
  311.      15.123628117913832,
  312.      15.248616780045351,
  313.      15.37360544217687,
  314.      15.49859410430839,
  315.      15.62358276643991,
  316.      15.748571428571429,
  317.      15.873560090702949,
  318.      15.998548752834468,
  319.      16.123537414965988,
  320.      16.248526077097505,
  321.      16.373514739229027,
  322.      16.498503401360544,
  323.      16.623492063492066,
  324.      16.748480725623583,
  325.      16.8734693877551,
  326.      16.998458049886622,
  327.      17.12344671201814,
  328.      17.24843537414966,
  329.      17.37342403628118,
  330.      17.4984126984127,
  331.      17.623401360544218,
  332.      17.748390022675736,
  333.      17.873378684807257,
  334.      17.998367346938775,
  335.      18.123356009070296,
  336.      18.248344671201814,
  337.      18.373333333333335,
  338.      18.498321995464853,
  339.      18.623310657596374,
  340.      18.74829931972789,
  341.      18.87328798185941,
  342.      18.99827664399093,
  343.      19.12326530612245,
  344.      19.24825396825397,
  345.      19.373242630385487,
  346.      19.49823129251701,
  347.      19.623219954648526,
  348.      19.748208616780047,
  349.      19.873197278911565,
  350.      19.998185941043083,
  351.      20.123174603174604,
  352.      20.248163265306122,
  353.      20.373151927437643,
  354.      20.49814058956916,
  355.      20.623129251700682,
  356.      20.7481179138322,
  357.      20.873106575963718,
  358.      20.99809523809524,
  359.      21.123083900226757,
  360.      21.248072562358278,
  361.      21.373061224489796,
  362.      21.498049886621317,
  363.      21.623038548752834,
  364.      21.748027210884356,
  365.      21.873015873015873,
  366.      21.99800453514739,
  367.      22.122993197278912,
  368.      22.24798185941043,
  369.      22.37297052154195,
  370.      22.49795918367347,
  371.      22.62294784580499]
  372.  
  373. sumpowers = 0.0
  374. output_list = []
  375.  
  376. for x in range(len(Lp) -1):
  377.     sublist = Lp[x: x+1].copy()
  378.     for spl in sublist:
  379.         sumpowers += 10 ** (spl / 10.0)
  380.  
  381.     res = 10.0 * math.log(sumpowers, 10)
  382.     output_list.append(res)
  383.  
  384. print("Output list: ")
  385. x = 0
  386. for item in output_list:
  387.     print("Time: ", t[x], item, " dB")
  388.     x += 1
  389.  
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement