~ubuntu-branches/ubuntu/trusty/gnuradio/trusty

« back to all changes in this revision

Viewing changes to gnuradio-core/src/python/gnuradio/blks2impl/logpwrfft.py

  • Committer: Bazaar Package Importer
  • Author(s): Kamal Mostafa
  • Date: 2010-03-13 07:46:01 UTC
  • mfrom: (2.1.2 sid)
  • Revision ID: james.westby@ubuntu.com-20100313074601-zjsa893a87bozyh7
Tags: 3.2.2.dfsg-1ubuntu1
* Fix build for Ubuntu lucid (LP: #260406)
  - add binary package dep for libusrp0, libusrp2-0: adduser
  - debian/rules clean: remove pre-built Qt moc files

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
#
 
2
# Copyright 2008 Free Software Foundation, Inc.
 
3
#
 
4
# This file is part of GNU Radio
 
5
#
 
6
# GNU Radio is free software; you can redistribute it and/or modify
 
7
# it under the terms of the GNU General Public License as published by
 
8
# the Free Software Foundation; either version 3, or (at your option)
 
9
# any later version.
 
10
#
 
11
# GNU Radio is distributed in the hope that it will be useful,
 
12
# but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
14
# GNU General Public License for more details.
 
15
#
 
16
# You should have received a copy of the GNU General Public License
 
17
# along with GNU Radio; see the file COPYING.  If not, write to
 
18
# the Free Software Foundation, Inc., 51 Franklin Street,
 
19
# Boston, MA 02110-1301, USA.
 
20
#
 
21
 
 
22
from gnuradio import gr, window
 
23
from stream_to_vector_decimator import stream_to_vector_decimator
 
24
import math
 
25
 
 
26
class _logpwrfft_base(gr.hier_block2):
 
27
    """
 
28
    Create a log10(abs(fft)) stream chain, with real or complex input.
 
29
    """
 
30
 
 
31
    def __init__(self, sample_rate, fft_size, ref_scale, frame_rate, avg_alpha, average):
 
32
        """
 
33
        Create an log10(abs(fft)) stream chain.
 
34
        Provide access to the setting the filter and sample rate.
 
35
        @param sample_rate        Incoming stream sample rate
 
36
        @param fft_size                Number of FFT bins
 
37
        @param ref_scale        Sets 0 dB value input amplitude
 
38
        @param frame_rate        Output frame rate
 
39
        @param avg_alpha        FFT averaging (over time) constant [0.0-1.0]
 
40
        @param average                Whether to average [True, False]
 
41
        """
 
42
        gr.hier_block2.__init__(self, self._name,
 
43
                                gr.io_signature(1, 1, self._item_size),          # Input signature
 
44
                                gr.io_signature(1, 1, gr.sizeof_float*fft_size)) # Output signature
 
45
 
 
46
        self._sd = stream_to_vector_decimator(item_size=self._item_size,
 
47
                                              sample_rate=sample_rate,
 
48
                                              vec_rate=frame_rate,
 
49
                                              vec_len=fft_size)
 
50
 
 
51
        fft_window = window.blackmanharris(fft_size)
 
52
        fft = self._fft_block[0](fft_size, True, fft_window)
 
53
        window_power = sum(map(lambda x: x*x, fft_window))
 
54
 
 
55
        c2mag = gr.complex_to_mag(fft_size)
 
56
        self._avg = gr.single_pole_iir_filter_ff(1.0, fft_size)
 
57
        self._log = gr.nlog10_ff(20, fft_size,
 
58
                                 -10*math.log10(fft_size)              # Adjust for number of bins
 
59
                                 -10*math.log10(window_power/fft_size) # Adjust for windowing loss
 
60
                                 -20*math.log10(ref_scale/2))          # Adjust for reference scale
 
61
        self.connect(self, self._sd, fft, c2mag, self._avg, self._log, self)
 
62
 
 
63
        self._average = average
 
64
        self._avg_alpha = avg_alpha
 
65
        self.set_avg_alpha(avg_alpha)
 
66
        self.set_average(average)
 
67
 
 
68
    def set_decimation(self, decim):
 
69
        """
 
70
        Set the decimation on stream decimator.
 
71
        @param decim the new decimation
 
72
        """
 
73
        self._sd.set_decimation(decim)
 
74
 
 
75
    def set_vec_rate(self, vec_rate):
 
76
        """
 
77
        Set the vector rate on stream decimator.
 
78
        @param vec_rate the new vector rate
 
79
        """
 
80
        self._sd.set_vec_rate(vec_rate)
 
81
 
 
82
    def set_sample_rate(self, sample_rate):
 
83
        """
 
84
        Set the new sampling rate
 
85
        @param sample_rate the new rate
 
86
        """
 
87
        self._sd.set_sample_rate(sample_rate)
 
88
 
 
89
    def set_average(self, average):
 
90
        """
 
91
        Set the averaging filter on/off.
 
92
        @param average true to set averaging on
 
93
        """
 
94
        self._average = average
 
95
        if self._average:
 
96
            self._avg.set_taps(self._avg_alpha)
 
97
        else:
 
98
            self._avg.set_taps(1.0)
 
99
 
 
100
    def set_avg_alpha(self, avg_alpha):
 
101
        """
 
102
        Set the average alpha and set the taps if average was on.
 
103
        @param avg_alpha the new iir filter tap
 
104
        """
 
105
        self._avg_alpha = avg_alpha
 
106
        self.set_average(self._average)
 
107
 
 
108
    def sample_rate(self):
 
109
        """
 
110
        Return the current sample rate.
 
111
        """
 
112
        return self._sd.sample_rate()
 
113
 
 
114
    def decimation(self):
 
115
        """
 
116
        Return the current decimation.
 
117
        """
 
118
        return self._sd.decimation()
 
119
 
 
120
    def frame_rate(self):
 
121
        """
 
122
        Return the current frame rate.
 
123
        """
 
124
        return self._sd.frame_rate()
 
125
 
 
126
    def average(self):
 
127
        """
 
128
        Return whether or not averaging is being performed.
 
129
        """
 
130
        return self._average
 
131
 
 
132
    def avg_alpha(self):
 
133
        """
 
134
        Return averaging filter constant.
 
135
        """
 
136
        return self._avg_alpha
 
137
 
 
138
class logpwrfft_f(_logpwrfft_base):
 
139
        """
 
140
        Create an fft block chain, with real input.
 
141
        """
 
142
        _name = "logpwrfft_f"
 
143
        _item_size = gr.sizeof_float
 
144
        _fft_block = (gr.fft_vfc, )
 
145
 
 
146
class logpwrfft_c(_logpwrfft_base):
 
147
        """
 
148
        Create an fft block chain, with complex input.
 
149
        """
 
150
        _name = "logpwrfft_c"
 
151
        _item_size = gr.sizeof_gr_complex
 
152
        _fft_block = (gr.fft_vcc, )
 
153