radio/util/sinus.py

   1 #!/bin/env python
   2
   3 from __future__ import print_function
   4
   5 import math
   6
   7 samples = 1024
   8 amplitudes = (32000, 0, 4000, 0, 4000)
   9
  10 max = 0
  11 min = 0
  12
  13 for i in range(samples):
  14     sample = 0.0
  15     for harmonic, amplitude in enumerate(amplitudes):
  16         sample += math.sin(math.pi * 2 * i * (harmonic + 1) / samples) * amplitude / 2
  17     sample = int(sample)
  18     if sample > max:
  19         max = sample
  20     elif sample < min:
  21         min = sample
  22     print("%d," % sample, end=' ')
  23     if i % 10 == 9:
  24         print()
  25
  26 print()
  27 print('Range is:', min, max)
  28
  29 if max > 32767 or min < -32768:
  30     print("Invalid range!")