23 #ifndef INCLUDED_DIGITAL_KURTOTIC_EQUALIZER_CC_H
24 #define INCLUDED_DIGITAL_KURTOTIC_EQUALIZER_CC_H
52 float d_alpha_p, d_alpha_q, d_alpha_m;
60 float re = (float)(x.real() >= 0.0f);
61 float im = (float)(x.imag() >= 0.0f);
75 float nrm = norm(out);
77 float epsilon_f = 1e-12;
81 d_p = (1-d_alpha_p)*d_p + (d_alpha_p)*nrm + epsilon_f;
82 d_q = (1-d_alpha_q)*d_q + (d_alpha_q)*out*out + epsilon_c;
83 d_m = (1-d_alpha_m)*d_m + (d_alpha_m)*nrm*nrm + epsilon_f;
84 d_u = d_m - 2.0f*(d_p*d_p) - d_q*d_q;
87 (sign(d_u) * (nrm*cnj - 2.0f*d_p*cnj - conj(d_q)*out) -
94 float re =
gr_clip(F.real(), 1.0);
95 float im =
gr_clip(F.imag(), 1.0);
108 throw std::out_of_range(
"digital_kurtotic_equalizer::set_gain: Gain value must be >= 0");