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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
/* -*- c++ -*- */
/*
* Copyright 2011 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.
*/
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#include <digital_mpsk_snr_est_cc.h>
#include <gr_io_signature.h>
#include <cstdio>
digital_mpsk_snr_est_cc_sptr
digital_make_mpsk_snr_est_cc(snr_est_type_t type,
double alpha)
{
return gnuradio::get_initial_sptr(new digital_mpsk_snr_est_cc(type, alpha));
}
digital_mpsk_snr_est_cc::digital_mpsk_snr_est_cc(snr_est_type_t type,
double alpha)
: gr_sync_block ("mpsk_snr_est_cc",
gr_make_io_signature(1, 1, sizeof(gr_complex)),
gr_make_io_signature(1, 1, sizeof(gr_complex)))
{
d_snr_est = NULL;
d_type = type;
set_alpha(alpha);
set_type(type);
// at least 1 estimator has to look back
set_history(2);
}
digital_mpsk_snr_est_cc::~digital_mpsk_snr_est_cc()
{
if(d_snr_est)
delete d_snr_est;
}
int
digital_mpsk_snr_est_cc::work(int noutput_items,
gr_vector_const_void_star &input_items,
gr_vector_void_star &output_items)
{
// This is a pass-through block; copy input to output
memcpy(output_items[0], input_items[0],
noutput_items * sizeof(gr_complex));
// Update the SNR estimate registers from the current
return d_snr_est->update(noutput_items, input_items);
}
double
digital_mpsk_snr_est_cc::snr()
{
if(d_snr_est)
return d_snr_est->snr();
else
throw std::runtime_error("digital_mpsk_snr_est_cc:: No SNR estimator defined.\n");
}
snr_est_type_t
digital_mpsk_snr_est_cc::type() const
{
return d_type;
}
double
digital_mpsk_snr_est_cc::alpha() const
{
return d_alpha;
}
void
digital_mpsk_snr_est_cc::set_type(snr_est_type_t t)
{
d_type = t;
if(d_snr_est)
delete d_snr_est;
switch (d_type) {
case(SNR_EST_SIMPLE):
d_snr_est = new digital_impl_mpsk_snr_est_simple(d_alpha);
break;
case(SNR_EST_SKEW):
d_snr_est = new digital_impl_mpsk_snr_est_skew(d_alpha);
break;
case(SNR_EST_M2M4):
d_snr_est = new digital_impl_mpsk_snr_est_m2m4(d_alpha);
break;
case(SNR_EST_SVR):
d_snr_est = new digital_impl_mpsk_snr_est_svr(d_alpha);
break;
default:
throw std::invalid_argument("digital_mpsk_snr_est_cc: unknown type specified.\n");
}
}
void
digital_mpsk_snr_est_cc::set_alpha(double alpha)
{
d_alpha = alpha;
if(d_snr_est)
d_snr_est->set_alpha(d_alpha);
}
|
