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
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
|
//
// Copyright 2010 Ettus Research LLC
//
// This program 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 of the License, or
// (at your option) any later version.
//
// This program 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 this program. If not, see <http://www.gnu.org/licenses/>.
//
#ifndef INCLUDED_LIBUHD_USRP_DSP_UTILS_HPP
#define INCLUDED_LIBUHD_USRP_DSP_UTILS_HPP
#include <uhd/config.hpp>
#include <uhd/utils/assert.hpp>
#include <boost/cstdint.hpp>
#include <boost/math/special_functions/round.hpp>
namespace uhd{ namespace usrp{
namespace dsp_type1{
template <class T> T ceil_log2(T num){
return std::ceil(std::log(num)/std::log(T(2)));
}
/*!
* Calculate the rx mux word from properties.
* \param is_quadrature true if the subdev is complex
* \param is_iq_swapped true if the i and q are reversed
* \param the 32-bit rx mux control word
*/
static inline boost::uint32_t calc_rx_mux_word(
bool is_quadrature,
bool is_iq_swapped
){
boost::uint32_t rx_mux = 0;
if (is_quadrature){
rx_mux = (0x01 << 2) | (0x00 << 0); //Q=ADC1, I=ADC0
}else{
rx_mux = (0x11 << 2) | (0x00 << 0); //Q=ZERO, I=ADC0
}
if (is_iq_swapped){
rx_mux = (rx_mux << 2) | (rx_mux >> 2);
}
return rx_mux;
}
/*!
* Calculate the tx mux word from properties.
* \param is_iq_swapped true if the i and q are reversed
* \param the 32-bit tx mux control word
*/
static inline boost::uint32_t calc_tx_mux_word(bool is_iq_swapped){
boost::uint32_t tx_mux = 0x10;
if (is_iq_swapped){
tx_mux = (tx_mux << 4) | (tx_mux >> 4);
}
return tx_mux;
}
/*!
* Calculate the cordic word from the frequency and clock rate.
* The frequency will be set to the actual (possible) frequency.
*
* \param freq the requested frequency in Hz
* \param codec_rate the dsp codec rate in Hz
* \param the 32-bit cordic control word
*/
static inline boost::uint32_t calc_cordic_word_and_update(
double &freq,
double codec_rate
){
UHD_ASSERT_THROW(std::abs(freq) < codec_rate/2.0);
static const double scale_factor = std::pow(2.0, 32);
//calculate the freq register word
boost::uint32_t freq_word = boost::math::iround((freq / codec_rate) * scale_factor);
//update the actual frequency
freq = (double(freq_word) / scale_factor) * codec_rate;
return freq_word;
}
/*!
* Calculate the CIC filter word from the rate.
* Check if requested decim/interp rate is:
* multiple of 4, enable two halfband filters
* multiple of 2, enable one halfband filter
* handle remainder in CIC
*
* \param rate the requested rate in Sps
* \return the 32-bit cic filter control word
*/
template <typename dsp_rate_type>
static inline boost::uint32_t calc_cic_filter_word(dsp_rate_type rate){
int hb0 = 0, hb1 = 0;
if (not (rate & 0x1)){
hb0 = 1;
rate /= 2;
}
if (not (rate & 0x1)){
hb1 = 1;
rate /= 2;
}
return (hb1 << 9) | (hb0 << 8) | (rate & 0xff);
}
/*!
* Calculate the IQ scale factor word from I and Q components.
* \param i the I component of the scalar
* \param q the Q component of the scalar
* \return the 32-bit scale factor control word
*/
static inline boost::uint32_t calc_iq_scale_word(
boost::int16_t i, boost::int16_t q
){
return (boost::uint32_t(i) << 16) | (boost::uint32_t(q) << 0);
}
/*!
* Calculate the IQ scale factor word from the rate.
* \param rate the requested rate in Sps
* \return the 32-bit scale factor control word
*/
template <typename dsp_rate_type>
static inline boost::uint32_t calc_iq_scale_word(dsp_rate_type rate){
// Calculate CIC interpolation (i.e., without halfband interpolators)
dsp_rate_type tmp_rate = calc_cic_filter_word(rate) & 0xff;
// Calculate closest multiplier constant to reverse gain absent scale multipliers
double rate_cubed = std::pow(double(tmp_rate), 3);
boost::int16_t scale = boost::math::iround((4096*std::pow(2, ceil_log2(rate_cubed)))/(1.65*rate_cubed));
return calc_iq_scale_word(scale, scale);
}
} //namespace dsp_type1
}} //namespace
#endif /* INCLUDED_LIBUHD_USRP_DSP_UTILS_HPP */
|
