testdsp4.py 1.4 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071
  1. import cmsisdsp as dsp
  2. import numpy as np
  3. from scipy import signal
  4. from scipy.fftpack import dct
  5. import cmsisdsp.fixedpoint as f
  6. from pyquaternion import Quaternion
  7. import colorama
  8. from colorama import init,Fore, Back, Style
  9. import statsmodels.tsa.stattools
  10. import scipy.spatial
  11. init()
  12. def printTitle(s):
  13. print("\n" + Fore.GREEN + Style.BRIGHT + s + Style.RESET_ALL)
  14. def printSubTitle(s):
  15. print("\n" + Style.BRIGHT + s + Style.RESET_ALL)
  16. printTitle("Max and AbsMax")
  17. a=np.array([1.,-3.,4.,0.,-10.,8.])
  18. i=dsp.arm_absmax_f32(a)
  19. printSubTitle("Fixed point tests")
  20. # Normalize for fixed point tests
  21. a = a / i[0]
  22. a7 = f.toQ7(a)
  23. a15 = f.toQ15(a)
  24. a31 = f.toQ31(a)
  25. print(a31)
  26. print(dsp.arm_absmax_no_idx_q31(a31))
  27. print(dsp.arm_max_no_idx_q31(a31))
  28. print(a15)
  29. print(dsp.arm_absmax_no_idx_q15(a15))
  30. print(dsp.arm_max_no_idx_q15(a15))
  31. print(a7)
  32. print(dsp.arm_absmax_no_idx_q7(a7))
  33. print(dsp.arm_max_no_idx_q7(a7))
  34. printTitle("Min and AbsMin")
  35. a=np.array([1.,-3.,4.,0.5,-10.,8.])
  36. i=dsp.arm_absmax_f32(a)
  37. printSubTitle("Fixed point tests")
  38. # Normalize for fixed point tests
  39. a = a / i[0]
  40. a7 = f.toQ7(a)
  41. a15 = f.toQ15(a)
  42. a31 = f.toQ31(a)
  43. print(a31)
  44. print(dsp.arm_absmin_no_idx_q31(a31))
  45. print(dsp.arm_min_no_idx_q31(a31))
  46. print(a15)
  47. print(dsp.arm_absmin_no_idx_q15(a15))
  48. print(dsp.arm_min_no_idx_q15(a15))
  49. print(a7)
  50. print(dsp.arm_absmin_no_idx_q7(a7))
  51. print(dsp.arm_min_no_idx_q7(a7))