[Commit-gnuradio] r4135 - gnuradio/branches/developers/n4hy/ofdm/gnuradio-examples/python/ofdm

[trondeau]

Author: trondeau
Date: 2006-12-17 23:34:40 -0700 (Sun, 17 Dec 2006)
New Revision: 4135

Added:
   
gnuradio/branches/developers/n4hy/ofdm/gnuradio-examples/python/ofdm/ofdm_test_nogui.py
Log:
adding ofdm test routine without the gui

Added: 
gnuradio/branches/developers/n4hy/ofdm/gnuradio-examples/python/ofdm/ofdm_test_nogui.py
===================================================================
--- 
gnuradio/branches/developers/n4hy/ofdm/gnuradio-examples/python/ofdm/ofdm_test_nogui.py
                             (rev 0)
+++ 
gnuradio/branches/developers/n4hy/ofdm/gnuradio-examples/python/ofdm/ofdm_test_nogui.py
     2006-12-18 06:34:40 UTC (rev 4135)
@@ -0,0 +1,126 @@
+#!/usr/bin/env python
+#
+# Copyright 2004,2005 Free Software Foundation, Inc.
+# 
+# This file is part of GNU Radio
+# 
+# GNU Radio is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2, or (at your option)
+# any later version.
+# 
+# GNU Radio is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU General Public License for more details.
+# 
+# You should have received a copy of the GNU General Public License
+# along with GNU Radio; see the file COPYING.  If not, write to
+# the Free Software Foundation, Inc., 51 Franklin Street,
+# Boston, MA 02110-1301, USA.
+# 
+
+from gnuradio import gr, gru
+from gnuradio import eng_notation
+from gnuradio.eng_option import eng_option
+from optparse import OptionParser
+import os, sys, random, math
+
+from ofdm_receiver import ofdm_receiver
+
+class my_flow_graph(gr.flow_graph):
+    def __init__(self):
+        gr.flow_graph.__init__(self)
+                
+        parser = OptionParser(option_class=eng_option)
+        parser.add_option("-W", "--waterfall", action="store_true", 
default=False,
+                          help="Enable waterfall display")
+        parser.add_option("-S", "--oscilloscope", action="store_true", 
default=False,
+                          help="Enable oscilloscope display")
+        parser.add_option("-I", "--audio-input", type="string", default="",
+                          help="pcm input device name.  E.g., hw:0,0 or 
/dev/dsp")
+        parser.add_option("-r", "--sample-rate", type="eng_float", 
default=48000,
+                          help="set sample rate to RATE (48000)")
+
+        (options, args) = parser.parse_args()   
+       sample_rate = int(options.sample_rate)
+        
+       if len(args) != 0:
+            parser.print_help()
+            sys.exit(1)
+
+        self.show_debug_info = True
+        
+        # parameters
+        occupied_tones = 70
+        fftsize = 128
+        cpsize = 32
+        symbol_length = fftsize + cpsize
+
+        SNR_db = 100.0
+        SNR = 10.0**(SNR_db/10.0)
+        frequency_offset = -0.0
+
+        power_in_signal = occupied_tones
+        noise_power_in_channel = power_in_signal/SNR
+        noise_power_required = noise_power_in_channel * fftsize / 
occupied_tones
+        noise_voltage = math.sqrt(noise_power_required)
+        rho = SNR / (SNR + 1.0)
+
+        print rho
+        
+        win = [1 for i in range(fftsize)]
+        #data = [1,0,0,1]* 5000
+        data = [random.randint(0,1) for i in range(100000)]
+       self.src = gr.vector_source_b(data, True)
+
+        # Modulator
+        self.ofdm = gr.ofdm_bpsk_mapper(occupied_tones, fftsize)
+        self.ifft = gr.fft_vcc(fftsize, False, win)
+        self.cp_adder = gr.ofdm_cyclic_prefixer(fftsize,symbol_length)
+
+        self.connect(self.src, self.ofdm, self.ifft, self.cp_adder)
+
+        # Channel
+        fdiff = frequency_offset * sample_rate / fftsize
+        self.noise_adder = gr.add_cc()
+        self.noise = gr.noise_source_c(gr.GR_GAUSSIAN,noise_voltage,1)
+        self.offset = 
gr.sig_source_c((sample_rate*1.0),gr.GR_SIN_WAVE,fdiff,1.0,0.0)
+        self.mixer_offset = gr.multiply_cc()
+        self.throttle = gr.throttle(gr.sizeof_gr_complex, 5e3)
+
+        self.connect(self.offset,(self.mixer_offset,0))
+        self.connect(self.noise,(self.noise_adder,0))
+        self.connect(self.mixer_offset, self.throttle, (self.noise_adder,1))
+
+        # Put modulated signal through channel
+        #self.scope = fftsink.fft_sink_c (self, panel, fft_size=1024, 
sample_rate=sample_rate, fft_rate=30)
+        self.connect(self.cp_adder, (self.mixer_offset,1))
+
+
+        # ML Sync
+        self.ofdm_rx = ofdm_receiver(self, fftsize, symbol_length, rho)
+        
+        # OFDM Demod
+        self.fftdemod = gr.fft_vcc(fftsize,True,win)
+        self.ofdm_demod = gr.ofdm_bpsk_demapper(occupied_tones, fftsize)
+        self.unpack = gr.unpack_k_bits_bb(1)
+
+        self.connect(self.noise_adder, self.ofdm_rx, self.fftdemod, 
self.ofdm_demod, self.unpack)
+        self.connect(self.ofdm_demod, gr.file_sink(gr.sizeof_char, "src"))
+             
+        # save output
+        self.outdat = gr.file_sink(gr.sizeof_char,"output_data")
+        self.connect(self.unpack, self.outdat)
+        
+        
+def main ():
+    fg = my_flow_graph()
+    fg.start()
+    fg.wait()
+    
+if __name__ == '__main__':
+    try:
+        main ()
+    except KeyboardInterrupt:
+        pass


Property changes on: 
gnuradio/branches/developers/n4hy/ofdm/gnuradio-examples/python/ofdm/ofdm_test_nogui.py
___________________________________________________________________
Name: svn:executable
   + *