iqgen: added unit tests

Added unit tests. Fixed issues with PRN code handling and TCXO.

Author: valeri-atamaniouk

.gitignore

@@ -56,3 +56,8 @@ tests/test_data_old*
 .project
 .cproject
 .settings/
+
+# Coverage
+.coverage
+htmlcov/
+

peregrine/iqgen/bits/amplitude_base.py

@@ -206,7 +206,7 @@ def convertUnits2SNR(value, units, noiseParams):
       snrDb = 10 * numpy.log10(snr)
     elif units == AmplitudeBase.UNITS_SNR_DB:
       snrDb = value
-    else:
+    else:  # pragma: no cover
       assert False
     return snrDb
 
@@ -240,6 +240,6 @@ def convertUnits2Amp(value, units, noiseParams):
       snrDb = value
       snr = 10. ** (0.1 * snrDb)
       amp = numpy.sqrt(4. * snr / freqTimesTau) * noiseSigma
-    else:
+    else:  # pragma: no cover
       assert False
     return amp

peregrine/iqgen/bits/doppler_base.py

@@ -181,12 +181,8 @@ def computeBatch(self,
     -------
     signal : numpy.ndarray(n_samples, dtype=float)
       Generated samples
-    userTimeX_s : float
-      End of interval time in seconds
-    chipAll_idx : numpy.ndarray(n_samples, dtype=float)
-      Code chip phases for the samples
-    chips : numpy.ndarray(n_samples, dtype=int)
-      Code combined with data
+    dopplerAll_hz : numpy.ndarray(n_samples, dtype=float)
+      Doppler values in Hz if debug is enabled
     '''
 
     userTimeAll_s = self.applySignalDelays(userTimeAll_s, carrierSignal)

peregrine/iqgen/bits/doppler_sine.py

@@ -65,21 +65,6 @@ def __str__(self):
         format(self.distance0_m, self.tec_epm2, self.speed0_mps,
                self.amplutude_mps, self.period_s, self.codeDopplerIgnored)
 
-  def __repr__(self):
-    '''
-    Constructs python expression presentation of object.
-
-    Returns
-    -------
-    string
-      Python expression presentation of object
-    '''
-    return "Doppler({}, {}, {}, {}, {})".format(self.distance0_m,
-                                                self.tec_epm2,
-                                                self.speed0_mps,
-                                                self.amplutude_mps,
-                                                self.period_s)
-
   def computeDistanceM(self, svTime_s):
     '''
     Computes doppler shift in meters.

peregrine/iqgen/bits/prn_glo_l1l2.py

@@ -14,9 +14,10 @@
 
 """
 import numpy
-from peregrine.include.glo_ca_code import value as GLONASS_CA_Code
+from peregrine.include.glo_ca_code import value as GLONASS_CACode
 
-caCode = GLONASS_CA_Code[:]
+# Binary CA code (0/1)
+caCode = (GLONASS_CACode < 0).astype(numpy.uint8)
 
 
 class PrnCode(object):
@@ -36,11 +37,7 @@ def __init__(self, prnNo):
       SV identifier
     '''
     super(PrnCode, self).__init__()
-    self.caCode = caCode[:]
-    tmp = numpy.asarray(self.caCode, dtype=numpy.int8)
-    tmp -= 1
-    tmp /= -2
-    self.binCode = tmp
+    self.binCode = caCode
     self.prnNo = prnNo
     self.bitLookup = numpy.asarray([1, -1], dtype=numpy.int8)
 

peregrine/iqgen/bits/prn_gps_l1ca.py

@@ -16,9 +16,9 @@
 
 """
 
-import peregrine.include.generateCAcode
+from peregrine.include.generateCAcode import caCodes as L1CACodes
 
-caCodes = peregrine.include.generateCAcode.caCodes
+caCodes = (L1CACodes < 0).astype(numpy.uint8)
 
 
 class PrnCode(object):
@@ -38,11 +38,7 @@ def __init__(self, prnNo):
       SV identifier
     '''
     super(PrnCode, self).__init__()
-    self.caCode = caCodes[prnNo - 1][:]
-    tmp = numpy.asarray(self.caCode, dtype=numpy.int8)
-    tmp -= 1
-    tmp /= -2
-    self.binCode = tmp
+    self.binCode = caCodes[prnNo - 1]
     self.prnNo = prnNo
     self.bitLookup = numpy.asarray([1, -1], dtype=numpy.int8)
 

peregrine/iqgen/bits/prn_gps_l2c.py

@@ -19,6 +19,8 @@
 
 from peregrine.include.generateL2CMcode import L2CMCodes
 
+caCodes = (L2CMCodes < 0).astype(numpy.uint8)
+
 
 class PrnCode(object):
   '''
@@ -42,7 +44,7 @@ def __init__(self, prnNo):
         SV identifier
       '''
       super(PrnCode.CM_Code, self).__init__()
-      self.binCode = numpy.asarray(L2CMCodes[prnNo - 1], dtype=numpy.int8) < 0
+      self.binCode = caCodes[prnNo - 1]
 
     def getCodeBits(self):
       return self.binCode

peregrine/iqgen/bits/satellite_base.py

@@ -96,10 +96,10 @@ def getAmplitude(self):
     return self.amplitude
 
   def __str__(self):
-    return self.getSvName()
-
-  def __repr__(self):
-    return self.getSvName()
+    '''
+    Returns string representation of SV object
+    '''
+    return self.getName()
 
   def getBatchSignals(self,
                       userTimeAll_s,

peregrine/iqgen/bits/tcxo_poly.py

@@ -53,12 +53,6 @@ def __str__(self, *args, **kwargs):
     '''
     return "TCXOPoly: coeffs=%s" % str(self.coeffs)
 
-  def __repr__(self):
-    '''
-    Provides string representation of the object
-    '''
-    return "TCXOPoly(%s)" % repr(self.coeffs)
-
   def computeTcxoTime(self, fromSample, toSample, outputConfig):
     '''
     Method generates time vector for the given sample index range depending on

peregrine/iqgen/bits/tcxo_sine.py

@@ -43,7 +43,7 @@ def __init__(self, initial_ppm, amplitude_ppm, period_s):
     self.initial_ppm = initial_ppm
     self.amplitude_ppm = amplitude_ppm
     self.period_s = period_s
-    self.c0 = -2. * scipy.constants.pi * amplitude_ppm * 1e-6
+    self.c0 = -amplitude_ppm * 1e-6 * self.period_s / (2. * scipy.constants.pi)
     self.c1 = 2. * scipy.constants.pi / period_s
     self.c2 = initial_ppm * 1e-6
 
@@ -54,13 +54,6 @@ def __str__(self, *args, **kwargs):
     return "TCXOSine: initial_ppm=%f amplitude_ppm=%f period_s=%f" % \
            (self.initial_ppm, self.amplitude_ppm, self.period_s)
 
-  def __repr__(self):
-    '''
-    Provides string representation of the object
-    '''
-    return "TCXOSine(%f, %f, %f)" % \
-           (self.initial_ppm, self.amplitude_ppm, self.period_s)
-
   def computeTcxoTime(self, fromSample, toSample, outputConfig):
     '''
     Method generates time vector for the given sample index range depending on
@@ -86,13 +79,13 @@ def computeTcxoTime(self, fromSample, toSample, outputConfig):
     time0_s = fromSample / outputConfig.SAMPLE_RATE_HZ
     timeX_s = toSample / outputConfig.SAMPLE_RATE_HZ
 
-    timeAll_s = numpy.linspace(time0_s * c1,
-                               timeX_s * c1,
+    timeAll_s = numpy.linspace(time0_s,
+                               timeX_s,
                                toSample - fromSample,
                                endpoint=False,
                                dtype=numpy.float)
 
-    result = numpy.cos(timeAll_s)
+    result = numpy.cos(timeAll_s * c1)
     result += -1.
     result *= c0
     if c2: