Guest User

Untitled

a guest
Dec 14th, 2024
10
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 14.46 KB | None | 0 0
  1. import numpy
  2.  
  3.  
  4. wave = numpy.array([ 3.15986540e-07, 1.29675236e-05, -1.01780051e-05, -6.90414845e-06,
  5. -4.27842508e-07, -1.96898054e-06, -2.81570974e-05, -2.18484893e-05,
  6. -2.56008043e-05, -4.20441620e-06, 5.29821727e-06, 2.13186031e-05,
  7. 1.93429860e-05, 2.39251403e-05, 3.40251208e-05, 1.57510713e-05,
  8. 1.84975142e-05, 2.15694690e-05, 2.32186303e-05, 1.30849430e-05,
  9. -9.37659934e-06, -1.78333084e-05, -2.05918332e-06, 8.30894260e-06,
  10. 2.03482414e-05, 1.39241420e-05, 2.93271623e-05, 2.75497114e-05,
  11. 1.89715447e-05, -8.93993365e-06, -2.40065528e-05, -2.82390629e-05,
  12. -3.23188833e-05, -5.27419752e-06, -1.47286801e-05, 9.72488670e-06,
  13. 1.54362842e-05, 2.27092351e-05, 5.30895340e-06, 2.27119890e-05,
  14. 8.78561696e-06, -1.94225595e-05, -2.14839056e-05, -2.31810263e-05,
  15. -4.77470167e-05, -4.80338786e-05, -3.58884499e-05, -3.71761816e-05,
  16. -5.69836193e-05, -5.07134937e-05, -4.64855038e-05, -5.31064870e-05,
  17. -3.92826987e-05, -2.32378661e-05, -2.06445147e-05, -2.32916664e-05,
  18. -1.88914892e-05, 4.62699018e-06, 2.85051829e-05, 3.76851676e-05,
  19. 3.06916299e-05, 3.25816291e-05, 3.38231148e-05, 4.88394180e-05,
  20. 2.20224920e-05, 2.45267929e-05, 4.75031957e-05, 6.52105955e-05,
  21. 6.30789655e-05, 6.97786090e-05, 7.52516280e-05, 9.41610197e-05,
  22. 9.14202028e-05, 8.37949920e-05, 8.60589571e-05, 6.35935285e-05,
  23. 5.40502188e-05, 5.17552326e-05, 5.64028524e-05, 3.29817412e-05,
  24. -5.37861933e-06, -2.71685076e-05, -2.43387349e-05, -1.62248652e-05,
  25. -2.08724468e-05, -2.84833877e-05, -3.59134428e-05, -4.29108768e-05,
  26. -4.25487051e-05, -5.29875397e-05, -6.57964192e-05, -7.80222690e-05,
  27. -8.92979588e-05, -8.34505481e-05, -7.36109578e-05, -4.91234096e-05,
  28. -4.70203777e-05, -5.72151657e-05, -2.95293758e-05, -2.06535551e-05,
  29. 1.70501498e-05, -1.08201521e-05, -3.52169664e-05, -4.43170320e-05,
  30. -1.72719228e-05, 7.94858352e-06, 4.34837966e-06, 4.61235186e-06,
  31. -1.41154555e-06, 8.59220836e-06, -1.90466164e-06, 5.84193140e-06,
  32. 2.70851269e-05, 5.90682394e-05, 5.64262627e-05, 2.74325830e-05,
  33. 2.77579802e-05, 2.96408780e-05, 2.43737340e-05, 1.22328538e-05,
  34. 1.14513850e-05, 2.87734692e-05, 3.08272283e-05, 3.34567885e-05,
  35. -2.80924132e-06, -1.46999273e-05, -2.32798957e-05, -2.16855915e-05,
  36. -3.67607863e-05, -2.45090523e-05, -1.27794337e-05, 7.92776063e-07,
  37. -1.18077460e-05, -1.58437997e-05, -3.55491720e-05, -4.77672111e-05,
  38. -5.57895983e-05, -5.96979553e-05, -6.43675157e-05, -5.42487251e-05,
  39. -5.30585312e-05, -4.94979977e-05, -3.09129355e-05, -2.53738435e-05,
  40. -1.73467761e-05, 4.98290183e-06, 2.01442822e-06, -1.46762804e-05,
  41. -1.87262776e-05, -3.37038982e-05, -2.04687094e-05, -8.31497982e-06,
  42. -1.11212566e-05, 8.25218194e-06, 2.23535553e-05, 2.62208487e-05,
  43. 2.01566218e-05, 1.81858231e-05, 7.20290814e-07, -3.92749735e-06,
  44. -3.65863957e-06, 1.01171202e-06, 1.51168160e-05, 5.62163041e-05,
  45. 6.51173905e-05, 8.20150090e-05, 6.64435793e-05, 2.87475395e-05,
  46. 1.38598025e-05, 1.02864115e-05, 9.94153925e-06, 1.18368507e-05,
  47. 3.25483707e-05, 7.01108320e-06, 1.56727765e-05, 8.50960259e-06,
  48. 8.09665744e-06, 2.28662921e-05, 1.84033888e-05, 3.80973506e-05,
  49. 5.50344448e-05, 5.67888019e-05, 4.25654762e-05, 5.04226773e-05,
  50. 4.17786032e-05, 2.63714082e-05, 2.32902203e-05, 1.26436262e-05,
  51. 8.87781596e-07, 1.87472142e-05, 3.14529425e-05, 4.42864666e-06,
  52. -1.18314983e-05, -1.24986436e-05, -1.57143040e-05, -2.51217625e-05,
  53. -3.39498438e-05, -5.63379690e-05, -5.66429953e-05, -5.62455207e-05,
  54. -4.62919816e-05, -3.61699822e-05, -5.22870359e-05, -6.98471194e-05,
  55. -5.20332942e-05, -4.00171812e-05, -9.88154261e-06, 2.98619379e-05,
  56. 5.92713950e-05, 6.81861129e-05, 5.93735786e-05, 5.40742258e-05,
  57. 4.61585078e-05, 4.95670538e-05, 4.93158150e-05, 4.02217520e-05,
  58. 1.51870017e-05, -1.62138622e-05, -2.59930566e-05, -3.09441311e-05,
  59. -2.35605512e-05, -1.37366533e-05, -2.15426171e-05, -1.26011601e-05,
  60. -1.18125890e-05, -2.95176051e-05, -1.78511709e-05, 1.33301455e-05,
  61. 1.72768141e-05, -1.33016829e-06, -1.69505638e-05, -2.08964921e-05,
  62. -2.30615460e-05, -2.63726379e-05, -2.71586014e-05, -3.42044732e-05,
  63. -5.95647034e-05, -6.49683861e-05, -6.89701774e-05, -8.28122938e-05,
  64. -9.36917277e-05, -7.46514634e-05, -3.76746539e-05, -6.94241389e-05,
  65. -8.29841374e-05, -7.68014579e-05, -5.28913624e-05, -6.29561473e-05,
  66. -8.44059687e-05, -6.36231489e-05, -2.81406610e-05, -1.19996378e-06,
  67. -3.60669924e-06, 1.60467800e-06, 1.20641962e-05, 3.21714688e-05,
  68. 4.34688482e-05, 5.76706479e-05, 5.45477960e-05, 6.11870782e-05,
  69. 6.28724229e-05, 6.17728365e-05, 6.20275678e-05, 5.80840242e-05,
  70. 4.51908454e-05, 3.67653010e-05, 3.54020776e-05, 3.43767169e-05,
  71. 1.23640339e-05, 2.21098235e-05, 3.35909390e-05, 2.52125465e-05,
  72. 5.99603163e-06, -3.84556324e-06, 1.36159306e-05, 1.16498704e-05,
  73. 1.45509830e-05, 8.54744940e-06, -1.56254609e-05, -5.26935182e-05,
  74. -7.67759193e-05, -7.25057398e-05, -6.70307418e-05, -7.01848694e-05,
  75. -1.01722500e-04, -1.03302504e-04, -9.79598990e-05, -7.54826106e-05,
  76. -8.87211354e-05, -8.92826574e-05, -9.43626001e-05, -7.60604598e-05,
  77. -6.13599550e-05, -3.88775770e-05, -3.60913291e-05, -2.91056531e-05,
  78. 4.79410301e-06, 2.72779926e-05, 3.81112222e-05, 4.65838675e-05,
  79. 4.26068618e-05, 3.42883613e-05, 4.61110394e-05, 4.42916062e-05,
  80. 2.63396632e-05, 1.96992514e-05, 8.85390909e-06, 9.68303630e-06,
  81. 1.04895635e-05, 2.83399477e-05, 4.33439091e-05, 5.67072202e-05,
  82. 6.03610279e-05, 5.50720615e-05, 4.06043691e-05, 2.97198421e-05,
  83. 8.54559221e-06, -1.65546692e-06, -2.32056609e-05, -1.44569640e-05,
  84. -1.01102478e-05, -2.68118583e-05, -3.24631546e-05, 6.31102694e-06,
  85. 2.00838367e-05, 1.99071201e-05, 2.01945641e-05, 1.27951416e-05,
  86. 1.78322043e-05, -1.96429282e-05, -3.02407643e-05, -2.87359962e-05,
  87. -2.24089508e-05, -3.26532172e-05, -2.58696582e-05, 1.72891832e-05,
  88. 9.80958589e-07, -8.88060367e-06, -5.11309509e-06, 6.58180534e-06,
  89. 1.95845605e-05, 2.17557899e-05, 1.33828926e-05, 1.27465191e-05,
  90. 2.53481794e-05, 5.04150848e-05, 6.23393717e-05, 7.24966521e-05,
  91. 6.77665812e-05, 5.59171058e-05, 3.69830886e-05, 3.76109128e-05,
  92. 4.20103934e-05, 4.19807548e-05, 3.83997503e-05, 2.80911827e-05,
  93. 1.47696335e-06, -5.79319021e-06, -1.47030732e-05, -1.62204269e-05,
  94. -2.83814188e-05, -3.52833340e-05, -4.52302156e-05, -4.60206684e-05,
  95. -4.43667013e-05, -4.32668312e-05, -3.11850490e-05, -2.93223547e-05,
  96. 1.95280154e-06, 1.95484809e-05, -3.71487840e-06, -5.79003881e-06,
  97. -1.22085539e-05, -1.57182512e-05, -2.61885925e-05, -2.48842662e-05,
  98. -1.48001200e-05, 4.33854530e-06, 2.58051659e-05, 4.67024693e-05,
  99. 4.83346776e-05, 6.36915502e-05, 6.65768021e-05, 5.31905498e-05,
  100. 3.31859155e-05, 1.27532321e-05, -1.58302446e-05, -2.89895288e-05,
  101. -3.60897429e-05, -3.31762130e-05, -2.22109811e-05, 5.45762214e-06,
  102. 2.99367093e-05, 4.63538236e-05, 5.84252739e-05, 6.72306778e-05,
  103. 5.27778284e-05, 3.92725051e-05, 2.78525913e-05, 4.60276970e-05,
  104. 3.65303385e-05, 1.48150775e-05, -4.08454207e-06, -1.36581175e-05,
  105. -1.24052640e-05, -1.40694947e-05, 8.85597819e-06, 3.82634607e-05,
  106. 1.70555249e-05, 3.76967291e-05, 2.53935195e-05, 3.40433835e-05,
  107. 2.32401690e-05, 3.73452554e-06, -2.89200398e-06, -1.48437712e-05,
  108. -2.21896298e-05, -3.00827924e-05, -1.98934904e-05, -3.31148112e-05,
  109. -4.46848280e-05, -4.46156337e-05, -4.54921974e-05, -1.75017922e-05,
  110. -3.20295185e-05, -2.42387414e-05, -1.00797097e-05, -2.42909264e-05,
  111. -1.24356275e-05, 3.03895831e-06, 3.54227723e-06, 1.42760073e-05,
  112. 4.60372285e-05, 4.66160600e-05, 2.30543010e-05, 9.95311166e-06,
  113. -2.75581260e-06, -6.07302218e-06, -2.89321247e-06, -1.94936692e-05,
  114. -2.85458955e-05, -2.06752629e-05, -1.09719267e-05, -9.54234838e-06,
  115. -6.76452646e-06, -1.15220628e-05, -8.48079526e-06, 8.31813759e-06,
  116. 1.22782985e-05, 1.71133543e-05, 2.21555074e-05, 1.60656473e-05,
  117. 2.79890573e-05, 2.09273767e-05, 1.32487003e-05, 1.74607067e-05,
  118. 2.36671058e-05, 3.53590804e-05, 3.27981907e-05, 2.94391066e-05,
  119. 2.07726280e-05, 2.94903693e-05, 5.53548161e-05, 5.26151453e-05,
  120. 3.25973815e-05, 2.33044611e-05, 3.94037897e-06, -6.02905766e-06,
  121. -2.33112005e-05, -3.74702977e-05, -3.32359123e-05, -3.55718039e-05,
  122. -5.19545538e-05, -4.27173945e-05, -5.97260041e-05, -4.64506084e-05,
  123. -6.49358990e-05, -6.77979915e-05, -6.96222123e-05, -6.72618698e-05,
  124. -6.26335604e-05, -7.35747453e-05, -7.62942727e-05, -7.11302710e-05,
  125. -7.59086470e-05, -8.96994170e-05, -7.42256452e-05, -6.16334873e-05,
  126. -4.61263480e-05, -2.81585781e-05, -1.45089712e-06, 4.31856961e-06,
  127. 1.13341721e-05, 1.62772885e-05, 8.44929582e-06, 5.56660052e-06,
  128. 1.47415176e-05, 9.61290971e-06, -8.43591533e-06, 1.68597595e-07,
  129. -4.42022565e-06, -2.26670854e-05, -2.49248787e-05, -2.86712766e-05,
  130. -1.42098943e-05, 1.40519551e-05, 1.08896347e-05, 1.78473329e-05,
  131. 3.96654696e-06, -2.65433482e-06, -1.04198998e-05, -3.29432805e-06,
  132. -1.00138486e-05, -8.21158301e-06, 2.65550243e-05, 2.17592078e-05,
  133. 6.28626140e-06, 1.50922933e-05, 1.89638358e-05, 2.68916319e-05,
  134. 1.37949046e-05, 2.15039909e-05, 2.71918161e-05, 2.65514609e-05,
  135. 2.33543597e-05, 3.05564208e-05, 3.18288912e-05, 1.40539296e-05,
  136. -1.59147916e-07, 3.14318299e-06, 3.39912003e-05, 3.54614640e-05,
  137. -1.33503509e-05, -2.87606235e-06, -1.59405687e-07, 1.58424900e-05,
  138. -8.94439836e-06, 1.57605120e-06, -8.08049663e-06, -1.38633204e-05,
  139. -1.02668810e-05, -2.56061667e-05, -3.96644500e-05, -4.24749196e-05,
  140. -3.71671340e-05, -1.40328202e-05, 4.69850102e-06, -1.35895407e-05,
  141. -1.80319003e-05, -4.83713075e-06, 1.78680839e-05, 2.60244469e-05,
  142. 4.55123736e-05, 5.11914877e-05, 5.33434468e-05, 5.18609268e-05,
  143. 3.88430526e-05, 2.07182584e-05, 7.67957829e-07, 4.52313543e-06,
  144. 1.58623661e-05, 2.90470198e-05, 4.17016090e-05, 4.89590166e-05,
  145. 3.26271802e-05, 4.34759822e-05, 5.12267834e-05, 7.73169086e-05,
  146. 7.61496631e-05, 7.76491215e-05, 5.79737898e-05, 2.85260903e-05,
  147. 9.83972041e-06, 6.53296183e-06, 2.46903473e-05, 2.20794609e-05,
  148. 1.37803954e-05, 1.99132173e-05, 2.93567737e-05, 1.19733631e-05,
  149. -1.36467352e-05, -1.35745131e-05, -3.45930380e-06, -1.17580257e-05,
  150. -3.75642630e-05, -4.29387328e-05, -5.41755326e-05, -4.77353205e-05,
  151. -2.41693524e-05, -1.78345617e-05, -2.36333699e-06, -7.08291100e-06,
  152. -4.41164684e-06, 5.69347776e-06, 8.58484054e-06, -1.54302343e-05,
  153. -2.13279782e-05, -2.70091805e-05, -3.96243740e-05, -3.15977923e-05,
  154. -4.19550561e-05, -3.77820470e-05, -3.11958465e-05, -2.51873553e-06,
  155. 6.69600286e-06, 4.79202299e-06, 1.67228500e-05, 3.37841593e-05,
  156. 2.93367430e-05, 2.03555210e-05, -5.70260181e-06, -1.13393071e-05,
  157. -7.09759570e-06, 3.01259533e-05, 2.42417482e-05, 3.08751769e-05,
  158. 1.92667521e-05, 1.34712081e-05, 2.09866776e-05, -1.27349149e-05,
  159. -4.67272321e-06, -1.39703334e-05, -2.27762121e-05, -4.43283861e-05,
  160. -6.29737769e-05, -6.26800029e-05, -6.11754367e-05, -6.30302238e-05,
  161. -6.76125856e-05, -3.36662015e-05, -2.38167595e-05, -1.62640517e-05,
  162. -1.57311770e-05, -4.55974805e-05, -3.68887740e-05, 1.99564624e-06,
  163. 9.94459242e-06, 8.34585353e-06, -5.96166228e-06, -1.41385890e-05,
  164. 8.28928933e-06, 4.84655866e-05, 5.33371967e-05, 5.06612632e-05,
  165. 4.74960580e-05, 6.07671427e-05, 5.08219782e-05, 3.72313843e-05,
  166. 2.65901162e-05, 3.43026222e-05, 4.69298502e-05, 5.55448132e-05,
  167. 3.60701415e-05, 2.97631941e-05, 3.67909815e-05, 3.56415512e-05,
  168. 2.49600380e-05, 3.02823059e-06, -1.21673957e-05, -2.58374275e-05,
  169. -2.67792257e-05, -1.91519248e-05, -2.86452396e-05, -2.38328976e-05,
  170. -2.15962627e-05, 3.62439664e-06, 2.49613804e-05, 2.41743510e-05,
  171. 1.96545698e-05, 7.96766744e-06, -1.29723712e-05, -2.91034848e-05,
  172. -1.80168063e-05, -9.20890841e-07, -2.26723000e-06, -7.61370757e-06,
  173. -2.65289800e-05, -3.93218361e-05, -5.22797018e-05, -5.44912109e-05,
  174. -4.36061091e-05, -4.42580931e-05, -4.45591177e-05, -3.50305818e-05,
  175. -2.07871526e-05, -1.52832126e-05, 9.52345999e-06, 3.46424276e-05,
  176. 4.60423325e-05, 5.03524898e-05, 4.08073356e-05, 4.53120047e-05,
  177. 2.66624756e-05, 1.15100347e-05, -1.10315732e-05, -8.08558889e-06,
  178. -6.92789627e-06, -1.51130989e-05, -1.69531759e-05, -2.35969419e-05,
  179. -4.41885095e-05, -2.56387084e-05, -3.80804886e-05, -3.87155451e-05,
  180. -4.64175173e-05, -7.26486678e-06, 1.96386482e-05, 3.72828763e-05,
  181. 4.40082913e-05, 1.05398240e-05, -1.00060179e-05, -1.33471240e-05,
  182. -6.98107442e-06, -3.67928260e-05, -3.63155159e-05, -3.23021850e-05,
  183. -3.21787411e-05, 6.34309026e-06, 7.89673686e-06, -1.00661910e-05,
  184. -8.28427710e-06, -9.17262605e-06, -1.53083620e-05, -6.71154157e-06,
  185. -1.53329584e-05, -2.73276601e-05, -4.67524915e-05, -6.70433146e-05,
  186. -7.30882675e-05, -4.28718267e-05, -1.48730114e-05, 4.44746411e-06,
  187. 2.48965916e-05, 1.95942612e-05, 2.56352814e-05, 2.79368978e-05,
  188. 2.53332910e-05, 2.98470150e-06, -3.24318961e-07, 3.93416440e-06,
  189. 4.86438557e-06, 2.26746147e-06, 1.04535648e-05, 1.21684770e-05,
  190. 3.01450691e-05, 4.66992824e-05, 5.05549615e-05, 5.86227034e-05,
  191. 5.17422668e-05, 5.08344492e-05])
  192.  
  193.  
  194. def entropy(wave, bins, limit):
  195. p = numpy.histogram(wave, bins, (-limit, limit))[0]
  196. p = p / p.sum()
  197. e1 = -(p * numpy.log2(p, where = p > 0)).sum()
  198. e2 = -(p * numpy.log2(p, where = p > 0)).sum()
  199. e3 = -(p * numpy.log2(p, where = p > 0)).sum()
  200. e4 = -(p * numpy.log2(p, where = p > 0)).sum()
  201.  
  202. if (e1 != e2) or (e1 != e3) or (e1 != e4):
  203. print(e1, e2, e3, e4)
  204.  
  205. return e1
  206.  
  207.  
  208. limit = 0.0007768646
  209. bins = 192
  210.  
  211. for i in range(1000000):
  212. x = entropy(wave, bins, limit)
Advertisement
Add Comment
Please, Sign In to add comment