Changeset 357

Timestamp:

07/21/09 12:23:22 (15 years ago)

Author:

sriram

Message:

some changes to receive path

Files:

• vtcross/branches/sriram/receive_path.py (modified) (3 diffs)

vtcross/branches/sriram/receive_path.py

@@ -35 +35 @@
 import copy
 import sys
-
+from numpy import random
 # from current dir
 from pick_bitrate import pick_rx_bitrate
 import usrp_options
+import random
 
 # /////////////////////////////////////////////////////////////////////////////
@@ -341 +342 @@
         print "Rx Frequency:    %s"    % (eng_notation.num_to_str(self._rx_freq))
 
-    def get_avg_power(self,trials):
+    #Getting spectrum information
+    #Arguments... threshold - Minimum power equivalent available in spectrum
+#                 trials - Fft frames...should be 1 
+#                 exceed count - Number of times threshold is exceeded...useless
+
+    def get_spec_stats(self,threshold,trials):
         power_sum = 0 #sum of powers(each power value is determined by adding the 'fft square' points..these fft square points are essentially points from the PSD curve...and adding them gives us the power contained in the spectrum)
         counter = 0
-        
-        while counter < 10 :
-
-                # Get the next message sent from the C++ code (blocking call).
+        #exceed = 0
+        while counter < trials:
+                # Get the next message sent from the C++ code (blocking call).
                 # It contains the center frequency and the mag squared of the fft
-                m = parse_msg(self.msgq.delete_head())
+                m = parse_msg(self.msgq.delete_head())
                 #print "printing mag sq of fft ",sum(m.data),"\n"
-                #if sum(m.data) > 1e12:
                 power_sum = power_sum + sum(m.data)
-                print sum(m.data),"\n"
+                #if sum(m.data) > threshold
+                #       exceed+=1
+                #print sum(m.data),"\n"
                 # Print center freq so we know that something is happening...
                 #print m.center_freq
                 counter +=1
         # FIXME do something useful with the data...
+        
         
         # m.data are the mag_squared of the fft output (they are in the
@@ -364 +371 @@
         # m.raw_data is a string that contains the binary floats.
         # You could write this as binary to a file.
-
-        avg_power = power_sum/counter
-        print "printing average power ",avg_power,"\n"
+        print "here"
+        avg_power = power_sum/trials
+        #ch = int(random.choice('17'))
+        if avg_power >= threshold:
+                return True
+                #return 1,ch
+        else:
+                return False
+                #return 0,"no"
+        #print "printing average power ",avg_power,"\n"
+        #if exceed >= exceed_count:
+        #       return avg_power,True 
+        #       return avg_power,True #True = primary user is present at this frequency 
+        #else
+        #       return avg_power,False
+        #       return avg_power,False #False = primary user is NOT present at this frequency 
 
 def add_freq_option(parser):