00001 /* -*- c++ -*- */ 00002 /* 00003 * Copyright 2006 Free Software Foundation, Inc. 00004 * 00005 * This file is part of GNU Radio 00006 * 00007 * GNU Radio is free software; you can redistribute it and/or modify 00008 * it under the terms of the GNU General Public License as published by 00009 * the Free Software Foundation; either version 2, or (at your option) 00010 * any later version. 00011 * 00012 * GNU Radio is distributed in the hope that it will be useful, 00013 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00014 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00015 * GNU General Public License for more details. 00016 * 00017 * You should have received a copy of the GNU General Public License 00018 * along with GNU Radio; see the file COPYING. If not, write to 00019 * the Free Software Foundation, Inc., 59 Temple Place - Suite 330, 00020 * Boston, MA 02111-1307, USA. 00021 */ 00022 00023 #ifndef _GRI_AGC_CC_H_ 00024 #define _GRI_AGC_CC_H_ 00025 00026 #include <math.h> 00027 00034 class gri_agc_cc { 00035 00036 public: 00037 gri_agc_cc (float rate = 1e-4, float reference = 1.0, float gain = 1.0) 00038 : _rate(rate), _reference(reference), _gain(gain) {}; 00039 00040 float rate () const { return _rate; } 00041 float reference () const { return _reference; } 00042 float gain () const { return _gain; } 00043 00044 void set_rate (float rate) { _rate = rate; } 00045 void set_reference (float reference) { _reference = reference; } 00046 void set_gain (float gain) { _gain = gain; } 00047 00048 gr_complex scale (gr_complex input){ 00049 gr_complex output = input * _gain; 00050 _gain += (_reference - sqrt(output.real()*output.real()+output.imag()*output.imag())) * _rate;//use abs or cabs to get approximation by absolute value, 00051 //note that abs is computationally more intensive then norm for a complex number 00052 return output; 00053 } 00054 00055 void scaleN (gr_complex output[], const gr_complex input[], unsigned n){ 00056 for (unsigned i = 0; i < n; i++) 00057 output[i] = scale (input[i]); 00058 } 00059 00060 protected: 00061 float _rate; // adjustment rate 00062 float _reference; // reference value 00063 float _gain; // current gain 00064 }; 00065 00066 #endif /* _GRI_AGC_CC_H_ */