1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
/* -*- c++ -*- */
/*
* Copyright 2004 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 3, 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.
*/
#ifndef INCLUDED_DAB_OFDM_FFE_ALL_IN_ONE_IMPL_H
#define INCLUDED_DAB_OFDM_FFE_ALL_IN_ONE_IMPL_H
#include <grdab/ofdm_ffe_all_in_one.h>
namespace gr {
namespace dab {
/*!
* \brief calculates fine frequency error estimation and averages it
* \ingroup DAB
* \param symbol_length number of samples in an OFDM symbol
* \param fft_length number of samples in an OFDM symbol without the cyclic prefix
* \param num_symbols number of symbols to use for averaging (more symbols is better, but symbols towards the end of the frame tend to have larger time offsets and worse values)
* \param alpha how fast should we adapt to new FFS error values (1=immediately)
* \param sample_rate sampling rate - needed to calculate the offset estimation in Hz
*
* input: port 0: complex - actual data; port 1: byte - trigger signal indicating the start of a frame
* output: float fine frequency offset estimation (in radian per sample)
*
* this is an all in one version of ffe in ofdm_sync_dab.py, because the flow graph does not allow to only calculate the estimation when its needed
*/
class ofdm_ffe_all_in_one_impl : public ofdm_ffe_all_in_one
{
private:
float calc_ffe_estimate(const gr_complex *iptr);
unsigned int d_symbol_length; // length of a symbol in samples
unsigned int d_fft_length; // length of a symbol without cyclic prefix in samples
unsigned int d_num_symbols; // number of symbols per frame to average over
float d_alpha; // adjustment speed factor
unsigned int d_sample_rate; // sample rate -- only needed to print the ffs error in Hz
unsigned int d_cur_symbol; // which symbol in the frame is currently under observation?
unsigned int d_cur_sample; // which sample in the symbol is currently under observation?
float d_ffs_error_sum; // sum of error samples in current frame
float d_estimated_error; // total estimated error
float d_estimated_error_per_sample; // total estimated error / fft_length
public:
ofdm_ffe_all_in_one_impl (unsigned int symbol_length, unsigned int fft_length, unsigned int num_symbols, float alpha, unsigned int sample_rate);
/*! \return fine frequency error estimate in Hz */
float ffe_estimate() { return d_estimated_error_per_sample*d_sample_rate/(2*M_PI); }
int work (int noutput_items,
gr_vector_const_void_star &input_items,
gr_vector_void_star &output_items);
};
}
}
#endif /* INCLUDED_DAB_OFDM_FFE_ALL_IN_ONE_H */
|