comicsansg/CubicSDR
comicsansg/CubicSDR
1
0
Fork 0
mirror of https://github.com/cjcliffe/CubicSDR.git synced 2024-11-23 12:18:37 -05:00
Code Issues Projects Releases Wiki Activity
98c7c30aee
CubicSDR/external/liquid-dsp/include/liquid/liquid.h

6891 lines
319 KiB
C++
Raw Blame History

/*
* Copyright (c) 2007 - 2016 Joseph Gaeddert
*
* Permission is hereby granted, free of charge, to any person obtaining a copy
* of this software and associated documentation files (the "Software"), to deal
* in the Software without restriction, including without limitation the rights
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
* copies of the Software, and to permit persons to whom the Software is
* furnished to do so, subject to the following conditions:
*
* The above copyright notice and this permission notice shall be included in
* all copies or substantial portions of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
* THE SOFTWARE.
*/
#ifndef __LIQUID_H__
#define __LIQUID_H__
#ifdef _MSC_VER
#pragma warning( push )
#pragma warning( disable : 4003 )
#endif
#ifdef __cplusplus
extern "C" {
# define LIQUID_USE_COMPLEX_H 0
#else
# define LIQUID_USE_COMPLEX_H 1
#endif // __cplusplus
// common headers
#include <inttypes.h>
//
// Make sure the version and version number macros weren't defined by
// some prevoiusly included header file.
//
#ifdef LIQUID_VERSION
# undef LIQUID_VERSION
#endif
#ifdef LIQUID_VERSION_NUMBER
# undef LIQUID_VERSION_NUMBER
#endif
//
// Compile-time version numbers
//
// LIQUID_VERSION = "X.Y.Z"
// LIQUID_VERSION_NUMBER = (X*1000000 + Y*1000 + Z)
//
#define LIQUID_VERSION "1.2.0"
#define LIQUID_VERSION_NUMBER 1002000
//
// Run-time library version numbers
//
extern const char liquid_version[];
const char * liquid_libversion(void);
int liquid_libversion_number(void);
// run-time library validation
#define LIQUID_VALIDATE_LIBVERSION \
if (LIQUID_VERSION_NUMBER != liquid_libversion_number()) { \
fprintf(stderr,"%s:%u: ", __FILE__,__LINE__); \
fprintf(stderr,"error: invalid liquid runtime library\n"); \
exit(1); \
} \
#define LIQUID_CONCAT(prefix, name) prefix ## name
#define LIQUID_VALIDATE_INPUT
/*
* Compile-time complex data type definitions
*
* Default: use the C99 complex data type, otherwise
* define complex type compatible with the C++ complex standard,
* otherwise resort to defining binary compatible array.
*/
#if LIQUID_USE_COMPLEX_H==1
# include <complex.h>
# define LIQUID_DEFINE_COMPLEX(R,C) typedef R _Complex C
#elif defined _GLIBCXX_COMPLEX || defined _LIBCPP_COMPLEX
# define LIQUID_DEFINE_COMPLEX(R,C) typedef std::complex<R> C
#else
# define LIQUID_DEFINE_COMPLEX(R,C) typedef struct {R real; R imag;} C;
#endif
//# define LIQUID_DEFINE_COMPLEX(R,C) typedef R C[2]
LIQUID_DEFINE_COMPLEX(float, liquid_float_complex);
LIQUID_DEFINE_COMPLEX(double, liquid_double_complex);
//
// MODULE : agc (automatic gain control)
//
#define AGC_MANGLE_CRCF(name) LIQUID_CONCAT(agc_crcf, name)
#define AGC_MANGLE_RRRF(name) LIQUID_CONCAT(agc_rrrf, name)
// large macro
// AGC : name-mangling macro
// T : primitive data type
// TC : input/output data type
#define LIQUID_AGC_DEFINE_API(AGC,T,TC) \
typedef struct AGC(_s) * AGC(); \
\
/* create automatic gain control object */ \
AGC() AGC(_create)(void); \
\
/* destroy object, freeing all internally-allocated memory */ \
void AGC(_destroy)(AGC() _q); \
\
/* print object properties to stdout */ \
void AGC(_print)(AGC() _q); \
\
/* reset object's internal state */ \
void AGC(_reset)(AGC() _q); \
\
/* execute automatic gain control on an single input sample */ \
/* _q : automatic gain control object */ \
/* _x : input sample */ \
/* _y : output sample */ \
void AGC(_execute)(AGC() _q, \
TC _x, \
TC * _y); \
\
/* execute automatic gain control on block of samples */ \
/* _q : automatic gain control object */ \
/* _x : input data array, [size: _n x 1] */ \
/* _n : number of input, output samples */ \
/* _y : output data array, [size: _n x 1] */ \
void AGC(_execute_block)(AGC() _q, \
TC * _x, \
unsigned int _n, \
TC * _y); \
\
/* lock/unlock gain control */ \
void AGC(_lock)( AGC() _q); \
void AGC(_unlock)(AGC() _q); \
\
/* get/set loop filter bandwidth; attack/release time */ \
float AGC(_get_bandwidth)(AGC() _q); \
void AGC(_set_bandwidth)(AGC() _q, float _bt); \
\
/* get/set signal level (linear) relative to unity energy */ \
float AGC(_get_signal_level)(AGC() _q); \
void AGC(_set_signal_level)(AGC() _q, float _signal_level); \
\
/* get/set signal level (dB) relative to unity energy */ \
float AGC(_get_rssi)(AGC() _q); \
void AGC(_set_rssi)(AGC() _q, float _rssi); \
\
/* get/set gain value (linear) relative to unity energy */ \
float AGC(_get_gain)(AGC() _q); \
void AGC(_set_gain)(AGC() _q, float _gain); \
\
/* initialize internal gain on input array */ \
/* _q : automatic gain control object */ \
/* _x : input data array, [size: _n x 1] */ \
/* _n : number of input, output samples */ \
void AGC(_init)(AGC() _q, \
TC * _x, \
unsigned int _n); \
// Define agc APIs
LIQUID_AGC_DEFINE_API(AGC_MANGLE_CRCF, float, liquid_float_complex)
LIQUID_AGC_DEFINE_API(AGC_MANGLE_RRRF, float, float)
//
// MODULE : audio
//
// CVSD: continuously variable slope delta
typedef struct cvsd_s * cvsd;
// create cvsd object
// _num_bits : number of adjacent bits to observe (4 recommended)
// _zeta : slope adjustment multiplier (1.5 recommended)
// _alpha : pre-/post-emphasis filter coefficient (0.9 recommended)
// NOTE: _alpha must be in [0,1]
cvsd cvsd_create(unsigned int _num_bits,
float _zeta,
float _alpha);
// destroy cvsd object
void cvsd_destroy(cvsd _q);
// print cvsd object parameters
void cvsd_print(cvsd _q);
// encode/decode single sample
unsigned char cvsd_encode(cvsd _q, float _audio_sample);
float cvsd_decode(cvsd _q, unsigned char _bit);
// encode/decode 8 samples at a time
void cvsd_encode8(cvsd _q, float * _audio, unsigned char * _data);
void cvsd_decode8(cvsd _q, unsigned char _data, float * _audio);
//
// MODULE : buffer
//
// circular buffer
#define CBUFFER_MANGLE_FLOAT(name) LIQUID_CONCAT(cbufferf, name)
#define CBUFFER_MANGLE_CFLOAT(name) LIQUID_CONCAT(cbuffercf, name)
// large macro
// CBUFFER : name-mangling macro
// T : data type
#define LIQUID_CBUFFER_DEFINE_API(CBUFFER,T) \
typedef struct CBUFFER(_s) * CBUFFER(); \
\
/* create circular buffer object of a particular size */ \
CBUFFER() CBUFFER(_create)(unsigned int _max_size); \
\
/* create circular buffer object of a particular size and */ \
/* specify the maximum number of elements that can be read */ \
/* at any given time. */ \
CBUFFER() CBUFFER(_create_max)(unsigned int _max_size, \
unsigned int _max_read); \
\
/* destroy cbuffer object, freeing all internal memory */ \
void CBUFFER(_destroy)(CBUFFER() _q); \
\
/* print cbuffer object properties */ \
void CBUFFER(_print)(CBUFFER() _q); \
\
/* print cbuffer object properties and internal state */ \
void CBUFFER(_debug_print)(CBUFFER() _q); \
\
/* clear internal buffer */ \
void CBUFFER(_reset)(CBUFFER() _q); \
\
/* get the number of elements currently in the buffer */ \
unsigned int CBUFFER(_size)(CBUFFER() _q); \
\
/* get the maximum number of elements the buffer can hold */ \
unsigned int CBUFFER(_max_size)(CBUFFER() _q); \
\
/* get the maximum number of elements you may read at once */ \
unsigned int CBUFFER(_max_read)(CBUFFER() _q); \
\
/* get the number of available slots (max_size - size) */ \
unsigned int CBUFFER(_space_available)(CBUFFER() _q); \
\
/* is buffer full? */ \
int CBUFFER(_is_full)(CBUFFER() _q); \
\
/* write a single sample into the buffer */ \
/* _q : circular buffer object */ \
/* _v : input sample */ \
void CBUFFER(_push)(CBUFFER() _q, \
T _v); \
\
/* write samples to the buffer */ \
/* _q : circular buffer object */ \
/* _v : output array */ \
/* _n : number of samples to write */ \
void CBUFFER(_write)(CBUFFER() _q, \
T * _v, \
unsigned int _n); \
\
/* remove and return a single element from the buffer */ \
/* _q : circular buffer object */ \
/* _v : pointer to sample output */ \
void CBUFFER(_pop)(CBUFFER() _q, \
T * _v); \
\
/* read buffer contents */ \
/* _q : circular buffer object */ \
/* _num_requested : number of elements requested */ \
/* _v : output pointer */ \
/* _nr : number of elements referenced by _v */ \
void CBUFFER(_read)(CBUFFER() _q, \
unsigned int _num_requested, \
T ** _v, \
unsigned int * _num_read); \
\
/* release _n samples from the buffer */ \
void CBUFFER(_release)(CBUFFER() _q, \
unsigned int _n); \
// Define buffer APIs
LIQUID_CBUFFER_DEFINE_API(CBUFFER_MANGLE_FLOAT, float)
LIQUID_CBUFFER_DEFINE_API(CBUFFER_MANGLE_CFLOAT, liquid_float_complex)
// Windowing functions
#define WINDOW_MANGLE_FLOAT(name) LIQUID_CONCAT(windowf, name)
#define WINDOW_MANGLE_CFLOAT(name) LIQUID_CONCAT(windowcf, name)
// large macro
// WINDOW : name-mangling macro
// T : data type
#define LIQUID_WINDOW_DEFINE_API(WINDOW,T) \
\
typedef struct WINDOW(_s) * WINDOW(); \
\
/* create window buffer object of length _n */ \
WINDOW() WINDOW(_create)(unsigned int _n); \
\
/* recreate window buffer object with new length */ \
/* _q : old window object */ \
/* _n : new window length */ \
WINDOW() WINDOW(_recreate)(WINDOW() _q, unsigned int _n); \
\
/* destroy window object, freeing all internally memory */ \
void WINDOW(_destroy)(WINDOW() _q); \
\
/* print window object to stdout */ \
void WINDOW(_print)(WINDOW() _q); \
\
/* print window object to stdout (with extra information) */ \
void WINDOW(_debug_print)(WINDOW() _q); \
\
/* reset window object (initialize to zeros) */ \
void WINDOW(_reset)(WINDOW() _q); \
\
/* read window buffer contents */ \
/* _q : window object */ \
/* _v : output pointer (set to internal array) */ \
void WINDOW(_read)(WINDOW() _q, T ** _v); \
\
/* index single element in buffer at a particular index */ \
/* _q : window object */ \
/* _i : index of element to read */ \
/* _v : output value pointer */ \
void WINDOW(_index)(WINDOW() _q, \
unsigned int _i, \
T * _v); \
\
/* push single element onto window buffer */ \
/* _q : window object */ \
/* _v : single input element */ \
void WINDOW(_push)(WINDOW() _q, \
T _v); \
\
/* write array of elements onto window buffer */ \
/* _q : window object */ \
/* _v : input array of values to write */ \
/* _n : number of input values to write */ \
void WINDOW(_write)(WINDOW() _q, \
T * _v, \
unsigned int _n); \
// Define window APIs
LIQUID_WINDOW_DEFINE_API(WINDOW_MANGLE_FLOAT, float)
LIQUID_WINDOW_DEFINE_API(WINDOW_MANGLE_CFLOAT, liquid_float_complex)
//LIQUID_WINDOW_DEFINE_API(WINDOW_MANGLE_UINT, unsigned int)
// wdelay functions : windowed-delay
// Implements an efficient z^-k delay with minimal memory
#define WDELAY_MANGLE_FLOAT(name) LIQUID_CONCAT(wdelayf, name)
#define WDELAY_MANGLE_CFLOAT(name) LIQUID_CONCAT(wdelaycf, name)
#define WDELAY_MANGLE_UINT(name) LIQUID_CONCAT(wdelayui, name)
// large macro
// WDELAY : name-mangling macro
// T : data type
#define LIQUID_WDELAY_DEFINE_API(WDELAY,T) \
typedef struct WDELAY(_s) * WDELAY(); \
\
/* create delay buffer object with '_delay' samples */ \
WDELAY() WDELAY(_create)(unsigned int _delay); \
\
/* re-create delay buffer object with '_delay' samples */ \
/* _q : old delay buffer object */ \
/* _delay : delay for new object */ \
WDELAY() WDELAY(_recreate)(WDELAY() _q, \
unsigned int _delay); \
\
/* destroy delay buffer object, freeing internal memory */ \
void WDELAY(_destroy)(WDELAY() _q); \
\
/* print delay buffer object's state to stdout */ \
void WDELAY(_print)(WDELAY() _q); \
\
/* clear/reset state of object */ \
void WDELAY(_reset)(WDELAY() _q); \
\
/* read delayed sample from delay buffer object */ \
/* _q : delay buffer object */ \
/* _v : value of delayed element */ \
void WDELAY(_read)(WDELAY() _q, \
T * _v); \
\
/* push new sample into delay buffer object */ \
/* _q : delay buffer object */ \
/* _v : new value to be added to buffer */ \
void WDELAY(_push)(WDELAY() _q, \
T _v); \
// Define wdelay APIs
LIQUID_WDELAY_DEFINE_API(WDELAY_MANGLE_FLOAT, float)
LIQUID_WDELAY_DEFINE_API(WDELAY_MANGLE_CFLOAT, liquid_float_complex)
//LIQUID_WDELAY_DEFINE_API(WDELAY_MANGLE_UINT, unsigned int)
//
// MODULE : channel
//
#define CHANNEL_MANGLE_CCCF(name) LIQUID_CONCAT(channel_cccf,name)
// large macro
// CHANNEL : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_CHANNEL_DEFINE_API(CHANNEL,TO,TC,TI) \
\
typedef struct CHANNEL(_s) * CHANNEL(); \
\
/* create channel object with default parameters */ \
CHANNEL() CHANNEL(_create)(void); \
\
/* destroy channel object, freeing all internal memory */ \
void CHANNEL(_destroy)(CHANNEL() _q); \
\
/* print channel object internals to standard output */ \
void CHANNEL(_print)(CHANNEL() _q); \
\
/* apply additive white Gausss noise impairment */ \
/* _q : channel object */ \
/* _noise_floor_dB : noise floor power spectral density */ \
/* _SNR_dB : signal-to-noise ratio [dB] */ \
void CHANNEL(_add_awgn)(CHANNEL() _q, \
float _noise_floor_dB, \
float _SNRdB); \
\
/* apply carrier offset impairment */ \
/* _q : channel object */ \
/* _frequency : carrier frequency offse [radians/sample */ \
/* _phase : carrier phase offset [radians] */ \
void CHANNEL(_add_carrier_offset)(CHANNEL() _q, \
float _frequency, \
float _phase); \
\
/* apply multi-path channel impairment */ \
/* _q : channel object */ \
/* _h : channel coefficients (NULL for random) */ \
/* _h_len : number of channel coefficients */ \
void CHANNEL(_add_multipath)(CHANNEL() _q, \
TC * _h, \
unsigned int _h_len); \
\
/* apply slowly-varying shadowing impairment */ \
/* _q : channel object */ \
/* _sigma : std. deviation for log-normal shadowing */ \
/* _fd : Doppler frequency, _fd in (0,0.5) */ \
void CHANNEL(_add_shadowing)(CHANNEL() _q, \
float _sigma, \
float _fd); \
\
/* apply channel impairments on single input sample */ \
/* _q : channel object */ \
/* _x : input sample */ \
/* _y : pointer to output sample */ \
void CHANNEL(_execute)(CHANNEL() _q, \
TI _x, \
TO * _y); \
\
/* apply channel impairments on block of samples */ \
/* _q : channel object */ \
/* _x : input array [size: _n x 1] */ \
/* _n : input array length */ \
/* _y : output array [size: _n x 1] */ \
void CHANNEL(_execute_block)(CHANNEL() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
LIQUID_CHANNEL_DEFINE_API(CHANNEL_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// time-varying multi-path channel
//
#define TVMPCH_MANGLE_CCCF(name) LIQUID_CONCAT(tvmpch_cccf,name)
// large macro
// TVMPCH : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_TVMPCH_DEFINE_API(TVMPCH,TO,TC,TI) \
\
typedef struct TVMPCH(_s) * TVMPCH(); \
\
/* create channel object with default parameters */ \
/* create time-varying multi-path channel emulator object */ \
/* _n : number of coefficients, _n > 0 */ \
/* _std : standard deviation */ \
/* _tau : coherence time */ \
TVMPCH() TVMPCH(_create)(unsigned int _n, \
float _std, \
float _tau); \
\
/* destroy channel object, freeing all internal memory */ \
void TVMPCH(_destroy)(TVMPCH() _q); \
\
/* reset object */ \
void TVMPCH(_reset)(TVMPCH() _q); \
\
/* print channel object internals to standard output */ \
void TVMPCH(_print)(TVMPCH() _q); \
\
/* push sample into emulator */ \
/* _q : channel object */ \
/* _x : input sample */ \
void TVMPCH(_push)(TVMPCH() _q, \
TI _x); \
\
/* compute output sample */ \
/* _q : channel object */ \
/* _y : output sample */ \
void TVMPCH(_execute)(TVMPCH() _q, \
TO * _y); \
\
/* apply channel impairments on a block of samples */ \
/* _q : channel object */ \
/* _x : input array [size: _nx x 1] */ \
/* _nx : input array length */ \
/* _y : output array */ \
void TVMPCH(_execute_block)(TVMPCH() _q, \
TI * _x, \
unsigned int _nx, \
TO * _y); \
LIQUID_TVMPCH_DEFINE_API(TVMPCH_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// MODULE : dotprod (vector dot product)
//
#define DOTPROD_MANGLE_RRRF(name) LIQUID_CONCAT(dotprod_rrrf,name)
#define DOTPROD_MANGLE_CCCF(name) LIQUID_CONCAT(dotprod_cccf,name)
#define DOTPROD_MANGLE_CRCF(name) LIQUID_CONCAT(dotprod_crcf,name)
// large macro
// DOTPROD : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_DOTPROD_DEFINE_API(DOTPROD,TO,TC,TI) \
\
/* run dot product without creating object [unrolled loop] */ \
/* _v : coefficients array [size: _n x 1] */ \
/* _x : input array [size: _n x 1] */ \
/* _n : dotprod length, _n > 0 */ \
/* _y : output sample pointer */ \
void DOTPROD(_run)( TC *_v, TI *_x, unsigned int _n, TO *_y); \
void DOTPROD(_run4)(TC *_v, TI *_x, unsigned int _n, TO *_y); \
\
typedef struct DOTPROD(_s) * DOTPROD(); \
\
/* create dot product object */ \
/* _v : coefficients array [size: _n x 1] */ \
/* _n : dotprod length, _n > 0 */ \
DOTPROD() DOTPROD(_create)(TC * _v, \
unsigned int _n); \
\
/* re-create dot product object */ \
/* _q : old dotprod object */ \
/* _v : coefficients array [size: _n x 1] */ \
/* _n : dotprod length, _n > 0 */ \
DOTPROD() DOTPROD(_recreate)(DOTPROD() _q, \
TC * _v, \
unsigned int _n); \
\
/* destroy dotprod object, freeing all internal memory */ \
void DOTPROD(_destroy)(DOTPROD() _q); \
\
/* print dotprod object internals to standard output */ \
void DOTPROD(_print)(DOTPROD() _q); \
\
/* execute dot product */ \
/* _q : dotprod object */ \
/* _x : input array [size: _n x 1] */ \
/* _y : output sample pointer */ \
void DOTPROD(_execute)(DOTPROD() _q, \
TI * _x, \
TO * _y); \
LIQUID_DOTPROD_DEFINE_API(DOTPROD_MANGLE_RRRF,
float,
float,
float)
LIQUID_DOTPROD_DEFINE_API(DOTPROD_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
LIQUID_DOTPROD_DEFINE_API(DOTPROD_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
//
// sum squared methods
//
float liquid_sumsqf(float * _v,
unsigned int _n);
float liquid_sumsqcf(liquid_float_complex * _v,
unsigned int _n);
//
// MODULE : equalization
//
// least mean-squares (LMS)
#define EQLMS_MANGLE_RRRF(name) LIQUID_CONCAT(eqlms_rrrf,name)
#define EQLMS_MANGLE_CCCF(name) LIQUID_CONCAT(eqlms_cccf,name)
// large macro
// EQLMS : name-mangling macro
// T : data type
#define LIQUID_EQLMS_DEFINE_API(EQLMS,T) \
typedef struct EQLMS(_s) * EQLMS(); \
\
/* create LMS EQ initialized with external coefficients */ \
/* _h : filter coefficients (NULL for {1,0,0...}) */ \
/* _h_len : filter length */ \
EQLMS() EQLMS(_create)(T * _h, \
unsigned int _h_len); \
\
/* create LMS EQ initialized with square-root Nyquist */ \
/* _type : filter type (e.g. LIQUID_FIRFILT_RRC) */ \
/* _k : samples/symbol */ \
/* _m : filter delay (symbols) */ \
/* _beta : rolloff factor (0 < beta <= 1) */ \
/* _dt : fractional sample delay */ \
EQLMS() EQLMS(_create_rnyquist)(int _type, \
unsigned int _k, \
unsigned int _m, \
float _beta, \
float _dt); \
\
/* create LMS EQ initialized with low-pass filter */ \
/* _h_len : filter length */ \
/* _fc : filter cut-off, _fc in (0,0.5] */ \
EQLMS() EQLMS(_create_lowpass)(unsigned int _h_len, \
float _fc); \
\
/* re-create EQ initialized with external coefficients */ \
/* _q : equalizer object */ \
/* _h : filter coefficients (NULL for {1,0,0...}) */ \
/* _h_len : filter length */ \
EQLMS() EQLMS(_recreate)(EQLMS() _q, \
T * _h, \
unsigned int _h_len); \
\
/* destroy equalizer object, freeing all internal memory */ \
void EQLMS(_destroy)(EQLMS() _q); \
\
/* reset equalizer object, clearing internal state */ \
void EQLMS(_reset)(EQLMS() _q); \
\
/* print equalizer internal state */ \
void EQLMS(_print)(EQLMS() _q); \
\
/* get/set equalizer learning rate */ \
float EQLMS(_get_bw)(EQLMS() _q); \
void EQLMS(_set_bw)(EQLMS() _q, \
float _lambda); \
\
/* push sample into equalizer internal buffer */ \
void EQLMS(_push)(EQLMS() _q, \
T _x); \
\
/* push sample into equalizer internal buffer as block */ \
/* _q : equalizer object */ \
/* _x : input sample array */ \
/* _n : input sample array length */ \
void EQLMS(_push_block)(EQLMS() _q, \
T * _x, \
unsigned int _n); \
\
/* execute internal dot product and return result */ \
/* _q : equalizer object */ \
/* _y : output sample */ \
void EQLMS(_execute)(EQLMS() _q, \
T * _y); \
\
/* execute equalizer with block of samples using constant */ \
/* modulus algorithm, operating on a decimation rate of _k */ \
/* samples. */ \
/* _q : equalizer object */ \
/* _k : down-sampling rate */ \
/* _x : input sample array [size: _n x 1] */ \
/* _n : input sample array length */ \
/* _y : output sample array [size: _n x 1] */ \
void EQLMS(_execute_block)(EQLMS() _q, \
unsigned int _k, \
T * _x, \
unsigned int _n, \
T * _y); \
\
/* step through one cycle of equalizer training */ \
/* _q : equalizer object */ \
/* _d : desired output */ \
/* _d_hat : actual output */ \
void EQLMS(_step)(EQLMS() _q, \
T _d, \
T _d_hat); \
\
/* step through one cycle of equalizer training (blind) */ \
/* _q : equalizer object */ \
/* _d_hat : actual output */ \
void EQLMS(_step_blind)(EQLMS() _q, \
T _d_hat); \
\
/* get equalizer's internal coefficients */ \
/* _q : equalizer object */ \
/* _w : weights [size: _p x 1] */ \
void EQLMS(_get_weights)(EQLMS() _q, \
T * _w); \
\
/* train equalizer object on group of samples */ \
/* _q : equalizer object */ \
/* _w : input/output weights [size: _p x 1] */ \
/* _x : received sample vector [size: _n x 1] */ \
/* _d : desired output vector [size: _n x 1] */ \
/* _n : input, output vector length */ \
void EQLMS(_train)(EQLMS() _q, \
T * _w, \
T * _x, \
T * _d, \
unsigned int _n); \
LIQUID_EQLMS_DEFINE_API(EQLMS_MANGLE_RRRF, float);
LIQUID_EQLMS_DEFINE_API(EQLMS_MANGLE_CCCF, liquid_float_complex);
// recursive least-squares (RLS)
#define EQRLS_MANGLE_RRRF(name) LIQUID_CONCAT(eqrls_rrrf,name)
#define EQRLS_MANGLE_CCCF(name) LIQUID_CONCAT(eqrls_cccf,name)
// large macro
// EQRLS : name-mangling macro
// T : data type
#define LIQUID_EQRLS_DEFINE_API(EQRLS,T) \
typedef struct EQRLS(_s) * EQRLS(); \
\
/* create RLS EQ initialized with external coefficients */ \
/* _h : filter coefficients (NULL for {1,0,0...}) */ \
/* _p : filter length */ \
EQRLS() EQRLS(_create)(T * _h, \
unsigned int _p); \
\
/* re-create RLS EQ initialized with external coefficients */ \
/* _q : initial equalizer object */ \
/* _h : filter coefficients (NULL for {1,0,0...}) */ \
/* _p : filter length */ \
EQRLS() EQRLS(_recreate)(EQRLS() _q, \
T * _h, \
unsigned int _p); \
\
/* destroy equalizer object, freeing all internal memory */ \
void EQRLS(_destroy)(EQRLS() _q); \
\
/* print equalizer internal state */ \
void EQRLS(_print)(EQRLS() _q); \
\
/* reset equalizer object, clearing internal state */ \
void EQRLS(_reset)(EQRLS() _q); \
\
/* get/set equalizer learning rate */ \
float EQRLS(_get_bw)(EQRLS() _q); \
void EQRLS(_set_bw)(EQRLS() _q, \
float _mu); \
\
/* push sample into equalizer internal buffer */ \
void EQRLS(_push)(EQRLS() _q, T _x); \
\
/* execute internal dot product and return result */ \
/* _q : equalizer object */ \
/* _y : output sample */ \
void EQRLS(_execute)(EQRLS() _q, T * _y); \
\
/* step through one cycle of equalizer training */ \
/* _q : equalizer object */ \
/* _d : desired output */ \
/* _d_hat : actual output */ \
void EQRLS(_step)(EQRLS() _q, T _d, T _d_hat); \
\
/* retrieve internal filter coefficients */ \
/* _q : equalizer object */ \
/* _w : weights [size: _p x 1] */ \
void EQRLS(_get_weights)(EQRLS() _q, \
T * _w); \
\
/* train equalizer object on group of samples */ \
/* _q : equalizer object */ \
/* _w : input/output weights [size: _p x 1] */ \
/* _x : received sample vector [size: _n x 1] */ \
/* _d : desired output vector [size: _n x 1] */ \
/* _n : input, output vector length */ \
void EQRLS(_train)(EQRLS() _q, \
T * _w, \
T * _x, \
T * _d, \
unsigned int _n);
LIQUID_EQRLS_DEFINE_API(EQRLS_MANGLE_RRRF, float);
LIQUID_EQRLS_DEFINE_API(EQRLS_MANGLE_CCCF, liquid_float_complex);
//
// MODULE : fec (forward error-correction)
//
// soft bit values
#define LIQUID_SOFTBIT_0 (0)
#define LIQUID_SOFTBIT_1 (255)
#define LIQUID_SOFTBIT_ERASURE (127)
// available CRC schemes
#define LIQUID_CRC_NUM_SCHEMES 7
typedef enum {
LIQUID_CRC_UNKNOWN=0, // unknown/unavailable CRC scheme
LIQUID_CRC_NONE, // no error-detection
LIQUID_CRC_CHECKSUM, // 8-bit checksum
LIQUID_CRC_8, // 8-bit CRC
LIQUID_CRC_16, // 16-bit CRC
LIQUID_CRC_24, // 24-bit CRC
LIQUID_CRC_32 // 32-bit CRC
} crc_scheme;
// pretty names for crc schemes
extern const char * crc_scheme_str[LIQUID_CRC_NUM_SCHEMES][2];
// Print compact list of existing and available CRC schemes
void liquid_print_crc_schemes();
// returns crc_scheme based on input string
crc_scheme liquid_getopt_str2crc(const char * _str);
// get length of CRC (bytes)
unsigned int crc_get_length(crc_scheme _scheme);
// generate error-detection key
// _scheme : error-detection scheme
// _msg : input data message, [size: _n x 1]
// _n : input data message size
unsigned int crc_generate_key(crc_scheme _scheme,
unsigned char * _msg,
unsigned int _n);
// generate error-detection key and append to end of message
// _scheme : error-detection scheme (resulting in 'p' bytes)
// _msg : input data message, [size: _n+p x 1]
// _n : input data message size (excluding key at end)
void crc_append_key(crc_scheme _scheme,
unsigned char * _msg,
unsigned int _n);
// validate message using error-detection key
// _scheme : error-detection scheme
// _msg : input data message, [size: _n x 1]
// _n : input data message size
// _key : error-detection key
int crc_validate_message(crc_scheme _scheme,
unsigned char * _msg,
unsigned int _n,
unsigned int _key);
// check message with key appended to end of array
// _scheme : error-detection scheme (resulting in 'p' bytes)
// _msg : input data message, [size: _n+p x 1]
// _n : input data message size (excluding key at end)
int crc_check_key(crc_scheme _scheme,
unsigned char * _msg,
unsigned int _n);
// get size of key (bytes)
unsigned int crc_sizeof_key(crc_scheme _scheme);
// available FEC schemes
#define LIQUID_FEC_NUM_SCHEMES 28
typedef enum {
LIQUID_FEC_UNKNOWN=0, // unknown/unsupported scheme
LIQUID_FEC_NONE, // no error-correction
LIQUID_FEC_REP3, // simple repeat code, r1/3
LIQUID_FEC_REP5, // simple repeat code, r1/5
LIQUID_FEC_HAMMING74, // Hamming (7,4) block code, r1/2 (really 4/7)
LIQUID_FEC_HAMMING84, // Hamming (7,4) with extra parity bit, r1/2
LIQUID_FEC_HAMMING128, // Hamming (12,8) block code, r2/3
LIQUID_FEC_GOLAY2412, // Golay (24,12) block code, r1/2
LIQUID_FEC_SECDED2216, // SEC-DED (22,16) block code, r8/11
LIQUID_FEC_SECDED3932, // SEC-DED (39,32) block code
LIQUID_FEC_SECDED7264, // SEC-DED (72,64) block code, r8/9
// codecs not defined internally (see http://www.ka9q.net/code/fec/)
LIQUID_FEC_CONV_V27, // r1/2, K=7, dfree=10
LIQUID_FEC_CONV_V29, // r1/2, K=9, dfree=12
LIQUID_FEC_CONV_V39, // r1/3, K=9, dfree=18
LIQUID_FEC_CONV_V615, // r1/6, K=15, dfree<=57 (Heller 1968)
// punctured (perforated) codes
LIQUID_FEC_CONV_V27P23, // r2/3, K=7, dfree=6
LIQUID_FEC_CONV_V27P34, // r3/4, K=7, dfree=5
LIQUID_FEC_CONV_V27P45, // r4/5, K=7, dfree=4
LIQUID_FEC_CONV_V27P56, // r5/6, K=7, dfree=4
LIQUID_FEC_CONV_V27P67, // r6/7, K=7, dfree=3
LIQUID_FEC_CONV_V27P78, // r7/8, K=7, dfree=3
LIQUID_FEC_CONV_V29P23, // r2/3, K=9, dfree=7
LIQUID_FEC_CONV_V29P34, // r3/4, K=9, dfree=6
LIQUID_FEC_CONV_V29P45, // r4/5, K=9, dfree=5
LIQUID_FEC_CONV_V29P56, // r5/6, K=9, dfree=5
LIQUID_FEC_CONV_V29P67, // r6/7, K=9, dfree=4
LIQUID_FEC_CONV_V29P78, // r7/8, K=9, dfree=4
// Reed-Solomon codes
LIQUID_FEC_RS_M8 // m=8, n=255, k=223
} fec_scheme;
// pretty names for fec schemes
extern const char * fec_scheme_str[LIQUID_FEC_NUM_SCHEMES][2];
// Print compact list of existing and available FEC schemes
void liquid_print_fec_schemes();
// returns fec_scheme based on input string
fec_scheme liquid_getopt_str2fec(const char * _str);
// fec object (pointer to fec structure)
typedef struct fec_s * fec;
// return the encoded message length using a particular error-
// correction scheme (object-independent method)
// _scheme : forward error-correction scheme
// _msg_len : raw, uncoded message length
unsigned int fec_get_enc_msg_length(fec_scheme _scheme,
unsigned int _msg_len);
// get the theoretical rate of a particular forward error-
// correction scheme (object-independent method)
float fec_get_rate(fec_scheme _scheme);
// create a fec object of a particular scheme
// _scheme : error-correction scheme
// _opts : (ignored)
fec fec_create(fec_scheme _scheme,
void *_opts);
// recreate fec object
// _q : old fec object
// _scheme : new error-correction scheme
// _opts : (ignored)
fec fec_recreate(fec _q,
fec_scheme _scheme,
void *_opts);
// destroy fec object
void fec_destroy(fec _q);
// print fec object internals
void fec_print(fec _q);
// encode a block of data using a fec scheme
// _q : fec object
// _dec_msg_len : decoded message length
// _msg_dec : decoded message
// _msg_enc : encoded message
void fec_encode(fec _q,
unsigned int _dec_msg_len,
unsigned char * _msg_dec,
unsigned char * _msg_enc);
// decode a block of data using a fec scheme
// _q : fec object
// _dec_msg_len : decoded message length
// _msg_enc : encoded message
// _msg_dec : decoded message
void fec_decode(fec _q,
unsigned int _dec_msg_len,
unsigned char * _msg_enc,
unsigned char * _msg_dec);
// decode a block of data using a fec scheme (soft decision)
// _q : fec object
// _dec_msg_len : decoded message length
// _msg_enc : encoded message (soft bits)
// _msg_dec : decoded message
void fec_decode_soft(fec _q,
unsigned int _dec_msg_len,
unsigned char * _msg_enc,
unsigned char * _msg_dec);
//
// Packetizer
//
// computes the number of encoded bytes after packetizing
//
// _n : number of uncoded input bytes
// _crc : error-detecting scheme
// _fec0 : inner forward error-correction code
// _fec1 : outer forward error-correction code
unsigned int packetizer_compute_enc_msg_len(unsigned int _n,
int _crc,
int _fec0,
int _fec1);
// computes the number of decoded bytes before packetizing
//
// _k : number of encoded bytes
// _crc : error-detecting scheme
// _fec0 : inner forward error-correction code
// _fec1 : outer forward error-correction code
unsigned int packetizer_compute_dec_msg_len(unsigned int _k,
int _crc,
int _fec0,
int _fec1);
typedef struct packetizer_s * packetizer;
// create packetizer object
//
// _n : number of uncoded input bytes
// _crc : error-detecting scheme
// _fec0 : inner forward error-correction code
// _fec1 : outer forward error-correction code
packetizer packetizer_create(unsigned int _dec_msg_len,
int _crc,
int _fec0,
int _fec1);
// re-create packetizer object
//
// _p : initialz packetizer object
// _n : number of uncoded input bytes
// _crc : error-detecting scheme
// _fec0 : inner forward error-correction code
// _fec1 : outer forward error-correction code
packetizer packetizer_recreate(packetizer _p,
unsigned int _dec_msg_len,
int _crc,
int _fec0,
int _fec1);
// destroy packetizer object
void packetizer_destroy(packetizer _p);
// print packetizer object internals
void packetizer_print(packetizer _p);
// access methods
unsigned int packetizer_get_dec_msg_len(packetizer _p);
unsigned int packetizer_get_enc_msg_len(packetizer _p);
crc_scheme packetizer_get_crc (packetizer _p);
fec_scheme packetizer_get_fec0 (packetizer _p);
fec_scheme packetizer_get_fec1 (packetizer _p);
// Execute the packetizer on an input message
//
// _p : packetizer object
// _msg : input message (uncoded bytes)
// _pkt : encoded output message
void packetizer_encode(packetizer _p,
const unsigned char * _msg,
unsigned char * _pkt);
// Execute the packetizer to decode an input message, return validity
// check of resulting data
//
// _p : packetizer object
// _pkt : input message (coded bytes)
// _msg : decoded output message
int packetizer_decode(packetizer _p,
const unsigned char * _pkt,
unsigned char * _msg);
// Execute the packetizer to decode an input message, return validity
// check of resulting data
//
// _p : packetizer object
// _pkt : input message (coded soft bits)
// _msg : decoded output message
int packetizer_decode_soft(packetizer _p,
const unsigned char * _pkt,
unsigned char * _msg);
//
// interleaver
//
typedef struct interleaver_s * interleaver;
// create interleaver
// _n : number of bytes
interleaver interleaver_create(unsigned int _n);
// destroy interleaver object
void interleaver_destroy(interleaver _q);
// print interleaver object internals
void interleaver_print(interleaver _q);
// set depth (number of internal iterations)
// _q : interleaver object
// _depth : depth
void interleaver_set_depth(interleaver _q,
unsigned int _depth);
// execute forward interleaver (encoder)
// _q : interleaver object
// _msg_dec : decoded (un-interleaved) message
// _msg_enc : encoded (interleaved) message
void interleaver_encode(interleaver _q,
unsigned char * _msg_dec,
unsigned char * _msg_enc);
// execute forward interleaver (encoder) on soft bits
// _q : interleaver object
// _msg_dec : decoded (un-interleaved) message
// _msg_enc : encoded (interleaved) message
void interleaver_encode_soft(interleaver _q,
unsigned char * _msg_dec,
unsigned char * _msg_enc);
// execute reverse interleaver (decoder)
// _q : interleaver object
// _msg_enc : encoded (interleaved) message
// _msg_dec : decoded (un-interleaved) message
void interleaver_decode(interleaver _q,
unsigned char * _msg_enc,
unsigned char * _msg_dec);
// execute reverse interleaver (decoder) on soft bits
// _q : interleaver object
// _msg_enc : encoded (interleaved) message
// _msg_dec : decoded (un-interleaved) message
void interleaver_decode_soft(interleaver _q,
unsigned char * _msg_enc,
unsigned char * _msg_dec);
//
// MODULE : fft (fast Fourier transform)
//
// type of transform
typedef enum {
LIQUID_FFT_UNKNOWN = 0, // unknown transform type
// regular complex one-dimensional transforms
LIQUID_FFT_FORWARD = +1, // complex one-dimensional FFT
LIQUID_FFT_BACKWARD = -1, // complex one-dimensional inverse FFT
// discrete cosine transforms
LIQUID_FFT_REDFT00 = 10, // real one-dimensional DCT-I
LIQUID_FFT_REDFT10 = 11, // real one-dimensional DCT-II
LIQUID_FFT_REDFT01 = 12, // real one-dimensional DCT-III
LIQUID_FFT_REDFT11 = 13, // real one-dimensional DCT-IV
// discrete sine transforms
LIQUID_FFT_RODFT00 = 20, // real one-dimensional DST-I
LIQUID_FFT_RODFT10 = 21, // real one-dimensional DST-II
LIQUID_FFT_RODFT01 = 22, // real one-dimensional DST-III
LIQUID_FFT_RODFT11 = 23, // real one-dimensional DST-IV
// modified discrete cosine transform
LIQUID_FFT_MDCT = 30, // MDCT
LIQUID_FFT_IMDCT = 31, // IMDCT
} liquid_fft_type;
#define LIQUID_FFT_MANGLE_FLOAT(name) LIQUID_CONCAT(fft,name)
// Macro : FFT
// FFT : name-mangling macro
// T : primitive data type
// TC : primitive data type (complex)
#define LIQUID_FFT_DEFINE_API(FFT,T,TC) \
\
typedef struct FFT(plan_s) * FFT(plan); \
\
/* create regular complex one-dimensional transform */ \
/* _n : transform size */ \
/* _x : pointer to input array [size: _n x 1] */ \
/* _y : pointer to output array [size: _n x 1] */ \
/* _dir : direction (e.g. LIQUID_FFT_FORWARD) */ \
/* _flags : options, optimization */ \
FFT(plan) FFT(_create_plan)(unsigned int _n, \
TC * _x, \
TC * _y, \
int _dir, \
int _flags); \
\
/* create real-to-real transform */ \
/* _n : transform size */ \
/* _x : pointer to input array [size: _n x 1] */ \
/* _y : pointer to output array [size: _n x 1] */ \
/* _type : transform type (e.g. LIQUID_FFT_REDFT00) */ \
/* _flags : options, optimization */ \
FFT(plan) FFT(_create_plan_r2r_1d)(unsigned int _n, \
T * _x, \
T * _y, \
int _type, \
int _flags); \
\
/* destroy transform */ \
void FFT(_destroy_plan)(FFT(plan) _p); \
\
/* print transform plan and internal strategy */ \
void FFT(_print_plan)(FFT(plan) _p); \
\
/* run the transform */ \
void FFT(_execute)(FFT(plan) _p); \
\
/* object-independent methods */ \
\
/* perform n-point FFT allocating plan internally */ \
/* _nfft : fft size */ \
/* _x : input array [size: _nfft x 1] */ \
/* _y : output array [size: _nfft x 1] */ \
/* _dir : fft direction: LIQUID_FFT_{FORWARD,BACKWARD} */ \
/* _flags : fft flags */ \
void FFT(_run)(unsigned int _n, \
TC * _x, \
TC * _y, \
int _dir, \
int _flags); \
\
/* perform n-point real FFT allocating plan internally */ \
/* _nfft : fft size */ \
/* _x : input array [size: _nfft x 1] */ \
/* _y : output array [size: _nfft x 1] */ \
/* _type : fft type, e.g. LIQUID_FFT_REDFT10 */ \
/* _flags : fft flags */ \
void FFT(_r2r_1d_run)(unsigned int _n, \
T * _x, \
T * _y, \
int _type, \
int _flags); \
\
/* perform _n-point fft shift */ \
void FFT(_shift)(TC * _x, \
unsigned int _n); \
LIQUID_FFT_DEFINE_API(LIQUID_FFT_MANGLE_FLOAT,float,liquid_float_complex)
// antiquated fft methods
// FFT(plan) FFT(_create_plan_mdct)(unsigned int _n,
// T * _x,
// T * _y,
// int _kind,
// int _flags);
//
// spectral periodogram
//
#define LIQUID_SPGRAM_MANGLE_CFLOAT(name) LIQUID_CONCAT(spgramcf,name)
#define LIQUID_SPGRAM_MANGLE_FLOAT(name) LIQUID_CONCAT(spgramf, name)
// Macro : SPGRAM
// SPGRAM : name-mangling macro
// T : primitive data type
// TC : primitive data type (complex)
// TI : primitive data type (input)
#define LIQUID_SPGRAM_DEFINE_API(SPGRAM,T,TC,TI) \
\
typedef struct SPGRAM(_s) * SPGRAM(); \
\
/* create spgram object */ \
/* _nfft : FFT size */ \
/* _window : window [size: _window_len x 1] */ \
/* _window_len : window length, _window_len in [1,_nfft] */ \
/* _delay : delay between transforms, _delay > 0 */ \
SPGRAM() SPGRAM(_create)(unsigned int _nfft, \
int _wtype, \
unsigned int _window_len, \
unsigned int _delay); \
\
/* create default spgram object (Kaiser-Bessel window) */ \
SPGRAM() SPGRAM(_create_default)(unsigned int _nfft); \
\
/* destroy spgram object */ \
void SPGRAM(_destroy)(SPGRAM() _q); \
\
/* resets the internal state of the spgram object */ \
void SPGRAM(_reset)(SPGRAM() _q); \
\
/* print internal state of the spgram object */ \
void SPGRAM(_print)(SPGRAM() _q); \
\
/* set methods */ \
int SPGRAM(_set_alpha)(SPGRAM() _q, float _alpha); \
\
/* access methods */ \
unsigned int SPGRAM(_get_nfft) (SPGRAM() _q); \
unsigned int SPGRAM(_get_window_len) (SPGRAM() _q); \
unsigned int SPGRAM(_get_delay) (SPGRAM() _q); \
uint64_t SPGRAM(_get_num_samples) (SPGRAM() _q); \
uint64_t SPGRAM(_get_num_samples_total) (SPGRAM() _q); \
uint64_t SPGRAM(_get_num_transforms) (SPGRAM() _q); \
uint64_t SPGRAM(_get_num_transforms_total)(SPGRAM() _q); \
float SPGRAM(_get_alpha) (SPGRAM() _q); \
\
/* push a single sample into the spgram object */ \
/* _q : spgram object */ \
/* _x : input sample */ \
void SPGRAM(_push)(SPGRAM() _q, \
TI _x); \
\
/* write a block of samples to the spgram object */ \
/* _q : spgram object */ \
/* _x : input buffer [size: _n x 1] */ \
/* _n : input buffer length */ \
void SPGRAM(_write)(SPGRAM() _q, \
TI * _x, \
unsigned int _n); \
\
/* compute spectral periodogram output (fft-shifted values */ \
/* in dB) from current buffer contents */ \
/* _q : spgram object */ \
/* _X : output spectrum (dB) [size: _nfft x 1] */ \
void SPGRAM(_get_psd)(SPGRAM() _q, \
T * _X); \
\
/* export gnuplot file */ \
/* _q : spgram object */ \
/* _filename : input buffer [size: _n x 1] */ \
int SPGRAM(_export_gnuplot)(SPGRAM() _q, \
const char * _filename); \
\
/* object-independent methods */ \
\
/* estimate spectrum on input signal */ \
/* _nfft : FFT size */ \
/* _x : input signal [size: _n x 1] */ \
/* _n : input signal length */ \
/* _psd : output spectrum, [size: _nfft x 1] */ \
void SPGRAM(_estimate_psd)(unsigned int _nfft, \
TI * _x, \
unsigned int _n, \
T * _psd); \
LIQUID_SPGRAM_DEFINE_API(LIQUID_SPGRAM_MANGLE_CFLOAT,
float,
liquid_float_complex,
liquid_float_complex)
LIQUID_SPGRAM_DEFINE_API(LIQUID_SPGRAM_MANGLE_FLOAT,
float,
liquid_float_complex,
float)
//
// asgram : ascii spectral periodogram
//
#define LIQUID_ASGRAM_MANGLE_CFLOAT(name) LIQUID_CONCAT(asgramcf,name)
#define LIQUID_ASGRAM_MANGLE_FLOAT(name) LIQUID_CONCAT(asgramf, name)
// Macro : ASGRAM
// ASGRAM : name-mangling macro
// T : primitive data type
// TC : primitive data type (complex)
// TI : primitive data type (input)
#define LIQUID_ASGRAM_DEFINE_API(ASGRAM,T,TC,TI) \
\
typedef struct ASGRAM(_s) * ASGRAM(); \
\
/* create asgram object with size _nfft */ \
ASGRAM() ASGRAM(_create)(unsigned int _nfft); \
\
/* destroy asgram object */ \
void ASGRAM(_destroy)(ASGRAM() _q); \
\
/* resets the internal state of the asgram object */ \
void ASGRAM(_reset)(ASGRAM() _q); \
\
/* set scale and offset for spectrogram */ \
/* _q : asgram object */ \
/* _offset : signal offset level [dB] */ \
/* _scale : signal scale [dB] */ \
void ASGRAM(_set_scale)(ASGRAM() _q, \
float _offset, \
float _scale); \
\
/* push a single sample into the asgram object */ \
/* _q : asgram object */ \
/* _x : input sample */ \
void ASGRAM(_push)(ASGRAM() _q, \
TI _x); \
\
/* write a block of samples to the asgram object */ \
/* _q : asgram object */ \
/* _x : input buffer [size: _n x 1] */ \
/* _n : input buffer length */ \
void ASGRAM(_write)(ASGRAM() _q, \
TI * _x, \
unsigned int _n); \
\
/* compute spectral periodogram output from current buffer */ \
/* contents */ \
/* _q : spgram object */ \
/* _ascii : output ASCII string [size: _nfft x 1] */ \
/* _peakval : peak power spectral density value [dB] */ \
/* _peakfreq : peak power spectral density frequency */ \
void ASGRAM(_execute)(ASGRAM() _q, \
char * _ascii, \
float * _peakval, \
float * _peakfreq); \
\
/* compute spectral periodogram output from current buffer */ \
/* contents and print standard format to stdout */ \
void ASGRAM(_print)(ASGRAM() _q); \
LIQUID_ASGRAM_DEFINE_API(LIQUID_ASGRAM_MANGLE_CFLOAT,
float,
liquid_float_complex,
liquid_float_complex)
LIQUID_ASGRAM_DEFINE_API(LIQUID_ASGRAM_MANGLE_FLOAT,
float,
liquid_float_complex,
float)
//
// MODULE : filter
//
//
// firdes: finite impulse response filter design
//
// prototypes
typedef enum {
LIQUID_FIRFILT_UNKNOWN=0, // unknown filter type
// Nyquist filter prototypes
LIQUID_FIRFILT_KAISER, // Nyquist Kaiser filter
LIQUID_FIRFILT_PM, // Parks-McClellan filter
LIQUID_FIRFILT_RCOS, // raised-cosine filter
LIQUID_FIRFILT_FEXP, // flipped exponential
LIQUID_FIRFILT_FSECH, // flipped hyperbolic secant
LIQUID_FIRFILT_FARCSECH, // flipped arc-hyperbolic secant
// root-Nyquist filter prototypes
LIQUID_FIRFILT_ARKAISER, // root-Nyquist Kaiser (approximate optimum)
LIQUID_FIRFILT_RKAISER, // root-Nyquist Kaiser (true optimum)
LIQUID_FIRFILT_RRC, // root raised-cosine
LIQUID_FIRFILT_hM3, // harris-Moerder-3 filter
LIQUID_FIRFILT_GMSKTX, // GMSK transmit filter
LIQUID_FIRFILT_GMSKRX, // GMSK receive filter
LIQUID_FIRFILT_RFEXP, // flipped exponential
LIQUID_FIRFILT_RFSECH, // flipped hyperbolic secant
LIQUID_FIRFILT_RFARCSECH, // flipped arc-hyperbolic secant
} liquid_firfilt_type;
// Design (root-)Nyquist filter from prototype
// _type : filter type (e.g. LIQUID_FIRFILT_RRC)
// _k : samples/symbol, _k > 1
// _m : symbol delay, _m > 0
// _beta : excess bandwidth factor, _beta in [0,1)
// _dt : fractional sample delay, _dt in [-1,1]
// _h : output coefficient buffer (length: 2*_k*_m+1)
void liquid_firdes_prototype(liquid_firfilt_type _type,
unsigned int _k,
unsigned int _m,
float _beta,
float _dt,
float * _h);
// returns filter type based on input string
int liquid_getopt_str2firfilt(const char * _str);
// estimate required filter length given
// _df : transition bandwidth (0 < _b < 0.5)
// _As : stop-band attenuation [dB], _As > 0
unsigned int estimate_req_filter_len(float _df,
float _As);
// estimate filter stop-band attenuation given
// _df : transition bandwidth (0 < _b < 0.5)
// _N : filter length
float estimate_req_filter_As(float _df,
unsigned int _N);
// estimate filter transition bandwidth given
// _As : stop-band attenuation [dB], _As > 0
// _N : filter length
float estimate_req_filter_df(float _As,
unsigned int _N);
// returns the Kaiser window beta factor give the filter's target
// stop-band attenuation (As) [Vaidyanathan:1993]
// _As : target filter's stop-band attenuation [dB], _As > 0
float kaiser_beta_As(float _As);
// Design FIR filter using Parks-McClellan algorithm
// band type specifier
typedef enum {
LIQUID_FIRDESPM_BANDPASS=0, // regular band-pass filter
LIQUID_FIRDESPM_DIFFERENTIATOR, // differentiating filter
LIQUID_FIRDESPM_HILBERT // Hilbert transform
} liquid_firdespm_btype;
// weighting type specifier
typedef enum {
LIQUID_FIRDESPM_FLATWEIGHT=0, // flat weighting
LIQUID_FIRDESPM_EXPWEIGHT, // exponential weighting
LIQUID_FIRDESPM_LINWEIGHT, // linear weighting
} liquid_firdespm_wtype;
// run filter design (full life cycle of object)
// _h_len : length of filter (number of taps)
// _num_bands : number of frequency bands
// _bands : band edges, f in [0,0.5], [size: _num_bands x 2]
// _des : desired response [size: _num_bands x 1]
// _weights : response weighting [size: _num_bands x 1]
// _wtype : weight types (e.g. LIQUID_FIRDESPM_FLATWEIGHT) [size: _num_bands x 1]
// _btype : band type (e.g. LIQUID_FIRDESPM_BANDPASS)
// _h : output coefficients array [size: _h_len x 1]
void firdespm_run(unsigned int _h_len,
unsigned int _num_bands,
float * _bands,
float * _des,
float * _weights,
liquid_firdespm_wtype * _wtype,
liquid_firdespm_btype _btype,
float * _h);
// structured object
typedef struct firdespm_s * firdespm;
// create firdespm object
// _h_len : length of filter (number of taps)
// _num_bands : number of frequency bands
// _bands : band edges, f in [0,0.5], [size: _num_bands x 2]
// _des : desired response [size: _num_bands x 1]
// _weights : response weighting [size: _num_bands x 1]
// _wtype : weight types (e.g. LIQUID_FIRDESPM_FLATWEIGHT) [size: _num_bands x 1]
// _btype : band type (e.g. LIQUID_FIRDESPM_BANDPASS)
firdespm firdespm_create(unsigned int _h_len,
unsigned int _num_bands,
float * _bands,
float * _des,
float * _weights,
liquid_firdespm_wtype * _wtype,
liquid_firdespm_btype _btype);
// destroy firdespm object
void firdespm_destroy(firdespm _q);
// print firdespm object internals
void firdespm_print(firdespm _q);
// execute filter design, storing result in _h
void firdespm_execute(firdespm _q, float * _h);
// Design FIR using kaiser window
// _n : filter length, _n > 0
// _fc : cutoff frequency, 0 < _fc < 0.5
// _As : stop-band attenuation [dB], _As > 0
// _mu : fractional sample offset, -0.5 < _mu < 0.5
// _h : output coefficient buffer, [size: _n x 1]
void liquid_firdes_kaiser(unsigned int _n,
float _fc,
float _As,
float _mu,
float *_h);
// Design FIR doppler filter
// _n : filter length
// _fd : normalized doppler frequency (0 < _fd < 0.5)
// _K : Rice fading factor (K >= 0)
// _theta : LoS component angle of arrival
// _h : output coefficient buffer
void liquid_firdes_doppler(unsigned int _n,
float _fd,
float _K,
float _theta,
float * _h);
// Design Nyquist raised-cosine filter
// _k : samples/symbol
// _m : symbol delay
// _beta : rolloff factor (0 < beta <= 1)
// _dt : fractional sample delay
// _h : output coefficient buffer (length: 2*k*m+1)
void liquid_firdes_rcos(unsigned int _k,
unsigned int _m,
float _beta,
float _dt,
float * _h);
// Design root-Nyquist raised-cosine filter
void liquid_firdes_rrcos(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Design root-Nyquist Kaiser filter
void liquid_firdes_rkaiser(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Design (approximate) root-Nyquist Kaiser filter
void liquid_firdes_arkaiser(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Design root-Nyquist harris-Moerder filter
void liquid_firdes_hM3(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Design GMSK transmit and receive filters
void liquid_firdes_gmsktx(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
void liquid_firdes_gmskrx(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Design flipped exponential Nyquist/root-Nyquist filters
void liquid_firdes_fexp( unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
void liquid_firdes_rfexp(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Design flipped hyperbolic secand Nyquist/root-Nyquist filters
void liquid_firdes_fsech( unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
void liquid_firdes_rfsech(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Design flipped arc-hyperbolic secand Nyquist/root-Nyquist filters
void liquid_firdes_farcsech( unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
void liquid_firdes_rfarcsech(unsigned int _k, unsigned int _m, float _beta, float _dt, float * _h);
// Compute group delay for an FIR filter
// _h : filter coefficients array
// _n : filter length
// _fc : frequency at which delay is evaluated (-0.5 < _fc < 0.5)
float fir_group_delay(float * _h,
unsigned int _n,
float _fc);
// Compute group delay for an IIR filter
// _b : filter numerator coefficients
// _nb : filter numerator length
// _a : filter denominator coefficients
// _na : filter denominator length
// _fc : frequency at which delay is evaluated (-0.5 < _fc < 0.5)
float iir_group_delay(float * _b,
unsigned int _nb,
float * _a,
unsigned int _na,
float _fc);
// liquid_filter_autocorr()
//
// Compute auto-correlation of filter at a specific lag.
//
// _h : filter coefficients [size: _h_len x 1]
// _h_len : filter length
// _lag : auto-correlation lag (samples)
float liquid_filter_autocorr(float * _h,
unsigned int _h_len,
int _lag);
// liquid_filter_crosscorr()
//
// Compute cross-correlation of two filters at a specific lag.
//
// _h : filter coefficients [size: _h_len]
// _h_len : filter length
// _g : filter coefficients [size: _g_len]
// _g_len : filter length
// _lag : cross-correlation lag (samples)
float liquid_filter_crosscorr(float * _h,
unsigned int _h_len,
float * _g,
unsigned int _g_len,
int _lag);
// liquid_filter_isi()
//
// Compute inter-symbol interference (ISI)--both RMS and
// maximum--for the filter _h.
//
// _h : filter coefficients [size: 2*_k*_m+1 x 1]
// _k : filter over-sampling rate (samples/symbol)
// _m : filter delay (symbols)
// _rms : output root mean-squared ISI
// _max : maximum ISI
void liquid_filter_isi(float * _h,
unsigned int _k,
unsigned int _m,
float * _rms,
float * _max);
// Compute relative out-of-band energy
//
// _h : filter coefficients [size: _h_len x 1]
// _h_len : filter length
// _fc : analysis cut-off frequency
// _nfft : fft size
float liquid_filter_energy(float * _h,
unsigned int _h_len,
float _fc,
unsigned int _nfft);
//
// IIR filter design
//
// IIR filter design filter type
typedef enum {
LIQUID_IIRDES_BUTTER=0,
LIQUID_IIRDES_CHEBY1,
LIQUID_IIRDES_CHEBY2,
LIQUID_IIRDES_ELLIP,
LIQUID_IIRDES_BESSEL
} liquid_iirdes_filtertype;
// IIR filter design band type
typedef enum {
LIQUID_IIRDES_LOWPASS=0,
LIQUID_IIRDES_HIGHPASS,
LIQUID_IIRDES_BANDPASS,
LIQUID_IIRDES_BANDSTOP
} liquid_iirdes_bandtype;
// IIR filter design coefficients format
typedef enum {
LIQUID_IIRDES_SOS=0,
LIQUID_IIRDES_TF
} liquid_iirdes_format;
// IIR filter design template
// _ftype : filter type (e.g. LIQUID_IIRDES_BUTTER)
// _btype : band type (e.g. LIQUID_IIRDES_BANDPASS)
// _format : coefficients format (e.g. LIQUID_IIRDES_SOS)
// _n : filter order
// _fc : low-pass prototype cut-off frequency
// _f0 : center frequency (band-pass, band-stop)
// _Ap : pass-band ripple in dB
// _As : stop-band ripple in dB
// _B : numerator
// _A : denominator
void liquid_iirdes(liquid_iirdes_filtertype _ftype,
liquid_iirdes_bandtype _btype,
liquid_iirdes_format _format,
unsigned int _n,
float _fc,
float _f0,
float _Ap,
float _As,
float * _B,
float * _A);
// compute analog zeros, poles, gain for specific filter types
void butter_azpkf(unsigned int _n,
liquid_float_complex * _za,
liquid_float_complex * _pa,
liquid_float_complex * _ka);
void cheby1_azpkf(unsigned int _n,
float _ep,
liquid_float_complex * _z,
liquid_float_complex * _p,
liquid_float_complex * _k);
void cheby2_azpkf(unsigned int _n,
float _es,
liquid_float_complex * _z,
liquid_float_complex * _p,
liquid_float_complex * _k);
void ellip_azpkf(unsigned int _n,
float _ep,
float _es,
liquid_float_complex * _z,
liquid_float_complex * _p,
liquid_float_complex * _k);
void bessel_azpkf(unsigned int _n,
liquid_float_complex * _z,
liquid_float_complex * _p,
liquid_float_complex * _k);
// compute frequency pre-warping factor
float iirdes_freqprewarp(liquid_iirdes_bandtype _btype,
float _fc,
float _f0);
// convert analog z/p/k form to discrete z/p/k form (bilinear z-transform)
// _za : analog zeros [length: _nza]
// _nza : number of analog zeros
// _pa : analog poles [length: _npa]
// _npa : number of analog poles
// _m : frequency pre-warping factor
// _zd : output digital zeros [length: _npa]
// _pd : output digital poles [length: _npa]
// _kd : output digital gain (should actually be real-valued)
void bilinear_zpkf(liquid_float_complex * _za,
unsigned int _nza,
liquid_float_complex * _pa,
unsigned int _npa,
liquid_float_complex _ka,
float _m,
liquid_float_complex * _zd,
liquid_float_complex * _pd,
liquid_float_complex * _kd);
// digital z/p/k low-pass to high-pass
// _zd : digital zeros (low-pass prototype), [length: _n]
// _pd : digital poles (low-pass prototype), [length: _n]
// _n : low-pass filter order
// _zdt : output digital zeros transformed [length: _n]
// _pdt : output digital poles transformed [length: _n]
void iirdes_dzpk_lp2hp(liquid_float_complex * _zd,
liquid_float_complex * _pd,
unsigned int _n,
liquid_float_complex * _zdt,
liquid_float_complex * _pdt);
// digital z/p/k low-pass to band-pass
// _zd : digital zeros (low-pass prototype), [length: _n]
// _pd : digital poles (low-pass prototype), [length: _n]
// _n : low-pass filter order
// _f0 : center frequency
// _zdt : output digital zeros transformed [length: 2*_n]
// _pdt : output digital poles transformed [length: 2*_n]
void iirdes_dzpk_lp2bp(liquid_float_complex * _zd,
liquid_float_complex * _pd,
unsigned int _n,
float _f0,
liquid_float_complex * _zdt,
liquid_float_complex * _pdt);
// convert discrete z/p/k form to transfer function
// _zd : digital zeros [length: _n]
// _pd : digital poles [length: _n]
// _n : filter order
// _kd : digital gain
// _b : output numerator [length: _n+1]
// _a : output denominator [length: _n+1]
void iirdes_dzpk2tff(liquid_float_complex * _zd,
liquid_float_complex * _pd,
unsigned int _n,
liquid_float_complex _kd,
float * _b,
float * _a);
// convert discrete z/p/k form to second-order sections
// _zd : digital zeros [length: _n]
// _pd : digital poles [length: _n]
// _n : filter order
// _kd : digital gain
// _B : output numerator [size: 3 x L+r]
// _A : output denominator [size: 3 x L+r]
// where r = _n%2, L = (_n-r)/2
void iirdes_dzpk2sosf(liquid_float_complex * _zd,
liquid_float_complex * _pd,
unsigned int _n,
liquid_float_complex _kd,
float * _B,
float * _A);
// additional IIR filter design templates
// design 2nd-order IIR filter (active lag)
// 1 + t2 * s
// F(s) = ------------
// 1 + t1 * s
//
// _w : filter bandwidth
// _zeta : damping factor (1/sqrt(2) suggested)
// _K : loop gain (1000 suggested)
// _b : output feed-forward coefficients [size: 3 x 1]
// _a : output feed-back coefficients [size: 3 x 1]
void iirdes_pll_active_lag(float _w,
float _zeta,
float _K,
float * _b,
float * _a);
// design 2nd-order IIR filter (active PI)
// 1 + t2 * s
// F(s) = ------------
// t1 * s
//
// _w : filter bandwidth
// _zeta : damping factor (1/sqrt(2) suggested)
// _K : loop gain (1000 suggested)
// _b : output feed-forward coefficients [size: 3 x 1]
// _a : output feed-back coefficients [size: 3 x 1]
void iirdes_pll_active_PI(float _w,
float _zeta,
float _K,
float * _b,
float * _a);
// checks stability of iir filter
// _b : feed-forward coefficients [size: _n x 1]
// _a : feed-back coefficients [size: _n x 1]
// _n : number of coefficients
int iirdes_isstable(float * _b,
float * _a,
unsigned int _n);
//
// linear prediction
//
// compute the linear prediction coefficients for an input signal _x
// _x : input signal [size: _n x 1]
// _n : input signal length
// _p : prediction filter order
// _a : prediction filter [size: _p+1 x 1]
// _e : prediction error variance [size: _p+1 x 1]
void liquid_lpc(float * _x,
unsigned int _n,
unsigned int _p,
float * _a,
float * _g);
// solve the Yule-Walker equations using Levinson-Durbin recursion
// for _symmetric_ autocorrelation
// _r : autocorrelation array [size: _p+1 x 1]
// _p : filter order
// _a : output coefficients [size: _p+1 x 1]
// _e : error variance [size: _p+1 x 1]
//
// NOTES:
// By definition _a[0] = 1.0
void liquid_levinson(float * _r,
unsigned int _p,
float * _a,
float * _e);
//
// auto-correlator (delay cross-correlation)
//
#define AUTOCORR_MANGLE_CCCF(name) LIQUID_CONCAT(autocorr_cccf,name)
#define AUTOCORR_MANGLE_RRRF(name) LIQUID_CONCAT(autocorr_rrrf,name)
// Macro:
// AUTOCORR : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_AUTOCORR_DEFINE_API(AUTOCORR,TO,TC,TI) \
\
typedef struct AUTOCORR(_s) * AUTOCORR(); \
\
/* create auto-correlator object */ \
/* _window_size : size of the correlator window */ \
/* _delay : correlator delay [samples] */ \
AUTOCORR() AUTOCORR(_create)(unsigned int _window_size, \
unsigned int _delay); \
\
/* destroy auto-correlator object, freeing internal memory */ \
void AUTOCORR(_destroy)(AUTOCORR() _q); \
\
/* reset auto-correlator object's internals */ \
void AUTOCORR(_reset)(AUTOCORR() _q); \
\
/* print auto-correlator parameters to stdout */ \
void AUTOCORR(_print)(AUTOCORR() _q); \
\
/* push sample into auto-correlator object */ \
void AUTOCORR(_push)(AUTOCORR() _q, \
TI _x); \
\
/* compute single auto-correlation output */ \
void AUTOCORR(_execute)(AUTOCORR() _q, \
TO * _rxx); \
\
/* compute auto-correlation on block of samples; the input */ \
/* and output arrays may have the same pointer */ \
/* _q : auto-correlation object */ \
/* _x : input array [size: _n x 1] */ \
/* _n : number of input, output samples */ \
/* _rxx : input array [size: _n x 1] */ \
void AUTOCORR(_execute_block)(AUTOCORR() _q, \
TI * _x, \
unsigned int _n, \
TO * _rxx); \
\
/* return sum of squares of buffered samples */ \
float AUTOCORR(_get_energy)(AUTOCORR() _q); \
LIQUID_AUTOCORR_DEFINE_API(AUTOCORR_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
LIQUID_AUTOCORR_DEFINE_API(AUTOCORR_MANGLE_RRRF,
float,
float,
float)
//
// Finite impulse response filter
//
#define FIRFILT_MANGLE_RRRF(name) LIQUID_CONCAT(firfilt_rrrf,name)
#define FIRFILT_MANGLE_CRCF(name) LIQUID_CONCAT(firfilt_crcf,name)
#define FIRFILT_MANGLE_CCCF(name) LIQUID_CONCAT(firfilt_cccf,name)
// Macro:
// FIRFILT : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_FIRFILT_DEFINE_API(FIRFILT,TO,TC,TI) \
typedef struct FIRFILT(_s) * FIRFILT(); \
\
FIRFILT() FIRFILT(_create)(TC * _h, unsigned int _n); \
\
/* create using Kaiser-Bessel windowed sinc method */ \
/* _n : filter length, _n > 0 */ \
/* _fc : filter cut-off frequency 0 < _fc < 0.5 */ \
/* _As : filter stop-band attenuation [dB], _As > 0 */ \
/* _mu : fractional sample offset, -0.5 < _mu < 0.5 */ \
FIRFILT() FIRFILT(_create_kaiser)(unsigned int _n, \
float _fc, \
float _As, \
float _mu); \
\
/* create from square-root Nyquist prototype */ \
/* _type : filter type (e.g. LIQUID_FIRFILT_RRC) */ \
/* _k : nominal samples/symbol, _k > 1 */ \
/* _m : filter delay [symbols], _m > 0 */ \
/* _beta : rolloff factor, 0 < beta <= 1 */ \
/* _mu : fractional sample offset,-0.5 < _mu < 0.5 */ \
FIRFILT() FIRFILT(_create_rnyquist)(int _type, \
unsigned int _k, \
unsigned int _m, \
float _beta, \
float _mu); \
\
/* create rectangular filter prototype */ \
FIRFILT() FIRFILT(_create_rect)(unsigned int _n); \
\
/* re-create filter */ \
/* _q : original filter object */ \
/* _h : pointer to filter coefficients [size: _n x 1] */ \
/* _n : filter length, _n > 0 */ \
FIRFILT() FIRFILT(_recreate)(FIRFILT() _q, \
TC * _h, \
unsigned int _n); \
\
/* destroy filter object and free all internal memory */ \
void FIRFILT(_destroy)(FIRFILT() _q); \
\
/* reset filter object's internal buffer */ \
void FIRFILT(_reset)(FIRFILT() _q); \
\
/* print filter object information */ \
void FIRFILT(_print)(FIRFILT() _q); \
\
/* set output scaling for filter */ \
void FIRFILT(_set_scale)(FIRFILT() _q, \
TC _scale); \
\
/* push sample into filter object's internal buffer */ \
/* _q : filter object */ \
/* _x : single input sample */ \
void FIRFILT(_push)(FIRFILT() _q, \
TI _x); \
\
/* execute the filter on internal buffer and coefficients */ \
/* _q : filter object */ \
/* _y : pointer to single output sample */ \
void FIRFILT(_execute)(FIRFILT() _q, \
TO * _y); \
\
/* execute the filter on a block of input samples; the */ \
/* input and output buffers may be the same */ \
/* _q : filter object */ \
/* _x : pointer to input array [size: _n x 1] */ \
/* _n : number of input, output samples */ \
/* _y : pointer to output array [size: _n x 1] */ \
void FIRFILT(_execute_block)(FIRFILT() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
\
/* return length of filter object */ \
unsigned int FIRFILT(_get_length)(FIRFILT() _q); \
\
/* compute complex frequency response of filter object */ \
/* _q : filter object */ \
/* _fc : frequency to evaluate */ \
/* _H : pointer to output complex frequency response */ \
void FIRFILT(_freqresponse)(FIRFILT() _q, \
float _fc, \
liquid_float_complex * _H); \
\
/* compute and return group delay of filter object */ \
/* _q : filter object */ \
/* _fc : frequency to evaluate */ \
float FIRFILT(_groupdelay)(FIRFILT() _q, \
float _fc); \
LIQUID_FIRFILT_DEFINE_API(FIRFILT_MANGLE_RRRF,
float,
float,
float)
LIQUID_FIRFILT_DEFINE_API(FIRFILT_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_FIRFILT_DEFINE_API(FIRFILT_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// FIR Hilbert transform
// 2:1 real-to-complex decimator
// 1:2 complex-to-real interpolator
//
#define FIRHILB_MANGLE_FLOAT(name) LIQUID_CONCAT(firhilbf, name)
//#define FIRHILB_MANGLE_DOUBLE(name) LIQUID_CONCAT(firhilb, name)
// NOTES:
// Although firhilb is a placeholder for both decimation and
// interpolation, separate objects should be used for each task.
#define LIQUID_FIRHILB_DEFINE_API(FIRHILB,T,TC) \
typedef struct FIRHILB(_s) * FIRHILB(); \
\
/* create finite impulse reponse Hilbert transform */ \
/* _m : filter semi-length, delay is 2*m+1 */ \
/* _As : filter stop-band attenuation [dB] */ \
FIRHILB() FIRHILB(_create)(unsigned int _m, \
float _As); \
\
/* destroy finite impulse reponse Hilbert transform */ \
void FIRHILB(_destroy)(FIRHILB() _q); \
\
/* print firhilb object internals to stdout */ \
void FIRHILB(_print)(FIRHILB() _q); \
\
/* reset firhilb object internal state */ \
void FIRHILB(_reset)(FIRHILB() _q); \
\
/* execute Hilbert transform (real to complex) */ \
/* _q : Hilbert transform object */ \
/* _x : real-valued input sample */ \
/* _y : complex-valued output sample */ \
void FIRHILB(_r2c_execute)(FIRHILB() _q, \
T _x, \
TC * _y); \
\
/* execute Hilbert transform (complex to real) */ \
/* _q : Hilbert transform object */ \
/* _x : complex-valued input sample */ \
/* _y : real-valued output sample */ \
void FIRHILB(_c2r_execute)(FIRHILB() _q, \
TC _x, \
T * _y); \
\
/* execute Hilbert transform decimator (real to complex) */ \
/* _q : Hilbert transform object */ \
/* _x : real-valued input array [size: 2 x 1] */ \
/* _y : complex-valued output sample */ \
void FIRHILB(_decim_execute)(FIRHILB() _q, \
T * _x, \
TC * _y); \
\
/* execute Hilbert transform decimator (real to complex) on */ \
/* a block of samples */ \
/* _q : Hilbert transform object */ \
/* _x : real-valued input array [size: 2*_n x 1] */ \
/* _n : number of *output* samples */ \
/* _y : complex-valued output array [size: _n x 1] */ \
void FIRHILB(_decim_execute_block)(FIRHILB() _q, \
T * _x, \
unsigned int _n, \
TC * _y); \
\
/* execute Hilbert transform interpolator (real to complex) */ \
/* _q : Hilbert transform object */ \
/* _x : complex-valued input sample */ \
/* _y : real-valued output array [size: 2 x 1] */ \
void FIRHILB(_interp_execute)(FIRHILB() _q, \
TC _x, \
T * _y); \
\
/* execute Hilbert transform interpolator (complex to real) */ \
/* on a block of samples */ \
/* _q : Hilbert transform object */ \
/* _x : complex-valued input array [size: _n x 1] */ \
/* _n : number of *input* samples */ \
/* _y : real-valued output array [size: 2*_n x 1] */ \
void FIRHILB(_interp_execute_block)(FIRHILB() _q, \
TC * _x, \
unsigned int _n, \
T * _y); \
LIQUID_FIRHILB_DEFINE_API(FIRHILB_MANGLE_FLOAT, float, liquid_float_complex)
//LIQUID_FIRHILB_DEFINE_API(FIRHILB_MANGLE_DOUBLE, double, liquid_double_complex)
//
// FFT-based finite impulse response filter
//
#define FFTFILT_MANGLE_RRRF(name) LIQUID_CONCAT(fftfilt_rrrf,name)
#define FFTFILT_MANGLE_CRCF(name) LIQUID_CONCAT(fftfilt_crcf,name)
#define FFTFILT_MANGLE_CCCF(name) LIQUID_CONCAT(fftfilt_cccf,name)
// Macro:
// FFTFILT : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_FFTFILT_DEFINE_API(FFTFILT,TO,TC,TI) \
typedef struct FFTFILT(_s) * FFTFILT(); \
\
/* create FFT-based FIR filter using external coefficients */ \
/* _h : filter coefficients [size: _h_len x 1] */ \
/* _h_len : filter length, _h_len > 0 */ \
/* _n : block size = nfft/2, at least _h_len-1 */ \
FFTFILT() FFTFILT(_create)(TC * _h, \
unsigned int _h_len, \
unsigned int _n); \
\
/* destroy filter object and free all internal memory */ \
void FFTFILT(_destroy)(FFTFILT() _q); \
\
/* reset filter object's internal buffer */ \
void FFTFILT(_reset)(FFTFILT() _q); \
\
/* print filter object information */ \
void FFTFILT(_print)(FFTFILT() _q); \
\
/* set output scaling for filter */ \
void FFTFILT(_set_scale)(FFTFILT() _q, \
TC _scale); \
\
/* execute the filter on internal buffer and coefficients */ \
/* _q : filter object */ \
/* _x : pointer to input data array [size: _n x 1] */ \
/* _y : pointer to output data array [size: _n x 1] */ \
void FFTFILT(_execute)(FFTFILT() _q, \
TI * _x, \
TO * _y); \
\
/* return length of filter object's internal coefficients */ \
unsigned int FFTFILT(_get_length)(FFTFILT() _q); \
LIQUID_FFTFILT_DEFINE_API(FFTFILT_MANGLE_RRRF,
float,
float,
float)
LIQUID_FFTFILT_DEFINE_API(FFTFILT_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_FFTFILT_DEFINE_API(FFTFILT_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Infinite impulse response filter
//
#define IIRFILT_MANGLE_RRRF(name) LIQUID_CONCAT(iirfilt_rrrf,name)
#define IIRFILT_MANGLE_CRCF(name) LIQUID_CONCAT(iirfilt_crcf,name)
#define IIRFILT_MANGLE_CCCF(name) LIQUID_CONCAT(iirfilt_cccf,name)
// Macro:
// IIRFILT : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_IIRFILT_DEFINE_API(IIRFILT,TO,TC,TI) \
\
typedef struct IIRFILT(_s) * IIRFILT(); \
\
/* create infinite impulse reponse filter */ \
/* _b : feed-forward coefficients [size: _nb x 1] */ \
/* _nb : number of feed-forward coefficients */ \
/* _a : feed-back coefficients [size: _na x 1] */ \
/* _na : number of feed-back coefficients */ \
IIRFILT() IIRFILT(_create)(TC * _b, \
unsigned int _nb, \
TC * _a, \
unsigned int _na); \
\
/* create IIR filter using 2nd-order secitons */ \
/* _B : feed-forward coefficients [size: _nsos x 3] */ \
/* _A : feed-back coefficients [size: _nsos x 3] */ \
IIRFILT() IIRFILT(_create_sos)(TC * _B, \
TC * _A, \
unsigned int _nsos); \
\
/* create IIR filter from design template */ \
/* _ftype : filter type (e.g. LIQUID_IIRDES_BUTTER) */ \
/* _btype : band type (e.g. LIQUID_IIRDES_BANDPASS) */ \
/* _format : coefficients format (e.g. LIQUID_IIRDES_SOS) */ \
/* _n : filter order */ \
/* _fc : low-pass prototype cut-off frequency */ \
/* _f0 : center frequency (band-pass, band-stop) */ \
/* _Ap : pass-band ripple in dB */ \
/* _As : stop-band ripple in dB */ \
IIRFILT() IIRFILT(_create_prototype)( \
liquid_iirdes_filtertype _ftype, \
liquid_iirdes_bandtype _btype, \
liquid_iirdes_format _format, \
unsigned int _order, \
float _fc, \
float _f0, \
float _Ap, \
float _As); \
\
/* create simplified low-pass Butterworth IIR filter */ \
/* _n : filter order */ \
/* _fc : low-pass prototype cut-off frequency */ \
IIRFILT() IIRFILT(_create_lowpass)( \
unsigned int _order, \
float _fc); \
\
/* create 8th-order integrator filter */ \
IIRFILT() IIRFILT(_create_integrator)(); \
\
/* create 8th-order differentiator filter */ \
IIRFILT() IIRFILT(_create_differentiator)(); \
\
/* create simple DC-blocking filter */ \
IIRFILT() IIRFILT(_create_dc_blocker)(float _alpha); \
\
/* create phase-locked loop iirfilt object */ \
/* _w : filter bandwidth */ \
/* _zeta : damping factor (1/sqrt(2) suggested) */ \
/* _K : loop gain (1000 suggested) */ \
IIRFILT() IIRFILT(_create_pll)(float _w, \
float _zeta, \
float _K); \
\
/* destroy iirfilt object, freeing all internal memory */ \
void IIRFILT(_destroy)(IIRFILT() _q); \
\
/* print iirfilt object properties to stdout */ \
void IIRFILT(_print)(IIRFILT() _q); \
\
/* clear/reset iirfilt object internals */ \
void IIRFILT(_reset)(IIRFILT() _q); \
\
/* compute filter output */ \
/* _q : iirfilt object */ \
/* _x : input sample */ \
/* _y : output sample pointer */ \
void IIRFILT(_execute)(IIRFILT() _q, \
TI _x, \
TO * _y); \
\
/* execute the filter on a block of input samples; the */ \
/* input and output buffers may be the same */ \
/* _q : filter object */ \
/* _x : pointer to input array [size: _n x 1] */ \
/* _n : number of input, output samples */ \
/* _y : pointer to output array [size: _n x 1] */ \
void IIRFILT(_execute_block)(IIRFILT() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
\
/* return iirfilt object's filter length (order + 1) */ \
unsigned int IIRFILT(_get_length)(IIRFILT() _q); \
\
/* compute complex frequency response of filter object */ \
/* _q : filter object */ \
/* _fc : frequency to evaluate */ \
/* _H : pointer to output complex frequency response */ \
void IIRFILT(_freqresponse)(IIRFILT() _q, \
float _fc, \
liquid_float_complex * _H); \
\
/* compute and return group delay of filter object */ \
/* _q : filter object */ \
/* _fc : frequency to evaluate */ \
float IIRFILT(_groupdelay)(IIRFILT() _q, float _fc); \
LIQUID_IIRFILT_DEFINE_API(IIRFILT_MANGLE_RRRF,
float,
float,
float)
LIQUID_IIRFILT_DEFINE_API(IIRFILT_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_IIRFILT_DEFINE_API(IIRFILT_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// FIR Polyphase filter bank
//
#define FIRPFB_MANGLE_RRRF(name) LIQUID_CONCAT(firpfb_rrrf,name)
#define FIRPFB_MANGLE_CRCF(name) LIQUID_CONCAT(firpfb_crcf,name)
#define FIRPFB_MANGLE_CCCF(name) LIQUID_CONCAT(firpfb_cccf,name)
// Macro:
// FIRPFB : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_FIRPFB_DEFINE_API(FIRPFB,TO,TC,TI) \
\
typedef struct FIRPFB(_s) * FIRPFB(); \
\
/* create firpfb from external coefficients */ \
/* _M : number of filters in the bank */ \
/* _h : coefficients [size: _M*_h_len x 1] */ \
/* _h_len : filter delay (symbols) */ \
FIRPFB() FIRPFB(_create)(unsigned int _M, \
TC * _h, \
unsigned int _h_len); \
\
/* create firpfb from external coefficients */ \
/* _M : number of filters in the bank */ \
/* _m : filter semi-length [samples] */ \
/* _fc : filter cut-off frequency 0 < _fc < 0.5 */ \
/* _As : filter stop-band suppression [dB] */ \
FIRPFB() FIRPFB(_create_kaiser)(unsigned int _M, \
unsigned int _m, \
float _fc, \
float _As); \
\
/* create firpfb from square-root Nyquist prototype */ \
/* _type : filter type (e.g. LIQUID_FIRFILT_RRC) */ \
/* _npfb : number of filters in the bank */ \
/* _k : nominal samples/symbol */ \
/* _m : filter delay (symbols) */ \
/* _beta : rolloff factor (0 < beta <= 1) */ \
FIRPFB() FIRPFB(_create_rnyquist)(int _type, \
unsigned int _npfb, \
unsigned int _k, \
unsigned int _m, \
float _beta); \
\
/* create from square-root derivative Nyquist prototype */ \
/* _type : filter type (e.g. LIQUID_FIRFILT_RRC) */ \
/* _npfb : number of filters in the bank */ \
/* _k : nominal samples/symbol */ \
/* _m : filter delay (symbols) */ \
/* _beta : rolloff factor (0 < beta <= 1) */ \
FIRPFB() FIRPFB(_create_drnyquist)(int _type, \
unsigned int _npfb, \
unsigned int _k, \
unsigned int _m, \
float _beta); \
\
/* re-create filterbank object */ \
/* _q : original firpfb object */ \
/* _M : number of filters in the bank */ \
/* _h : coefficients [size: _M x _h_len] */ \
/* _h_len : length of each filter */ \
FIRPFB() FIRPFB(_recreate)(FIRPFB() _q, \
unsigned int _M, \
TC * _h, \
unsigned int _h_len); \
\
/* destroy firpfb object, freeing all internal memory */ \
void FIRPFB(_destroy)(FIRPFB() _q); \
\
/* print firpfb object's parameters */ \
void FIRPFB(_print)(FIRPFB() _q); \
\
/* set output scaling for filter */ \
void FIRPFB(_set_scale)(FIRPFB() _q, \
TC _g); \
\
/* clear/reset firpfb object internal state */ \
void FIRPFB(_reset)(FIRPFB() _q); \
\
/* push sample into firpfb internal buffer */ \
void FIRPFB(_push)(FIRPFB() _q, TI _x); \
\
/* execute the filter on internal buffer and coefficients */ \
/* _q : firpfb object */ \
/* _i : index of filter to use */ \
/* _y : pointer to output sample */ \
void FIRPFB(_execute)(FIRPFB() _q, \
unsigned int _i, \
TO * _y); \
\
/* execute the filter on a block of input samples; the */ \
/* input and output buffers may be the same */ \
/* _q : firpfb object */ \
/* _i : index of filter to use */ \
/* _x : pointer to input array [size: _n x 1] */ \
/* _n : number of input, output samples */ \
/* _y : pointer to output array [size: _n x 1] */ \
void FIRPFB(_execute_block)(FIRPFB() _q, \
unsigned int _i, \
TI * _x, \
unsigned int _n, \
TO * _y); \
LIQUID_FIRPFB_DEFINE_API(FIRPFB_MANGLE_RRRF,
float,
float,
float)
LIQUID_FIRPFB_DEFINE_API(FIRPFB_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_FIRPFB_DEFINE_API(FIRPFB_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Interpolators
//
// firinterp : finite impulse response interpolator
#define FIRINTERP_MANGLE_RRRF(name) LIQUID_CONCAT(firinterp_rrrf,name)
#define FIRINTERP_MANGLE_CRCF(name) LIQUID_CONCAT(firinterp_crcf,name)
#define FIRINTERP_MANGLE_CCCF(name) LIQUID_CONCAT(firinterp_cccf,name)
#define LIQUID_FIRINTERP_DEFINE_API(FIRINTERP,TO,TC,TI) \
\
typedef struct FIRINTERP(_s) * FIRINTERP(); \
\
/* create interpolator from external coefficients */ \
/* _M : interpolation factor */ \
/* _h : filter coefficients [size: _h_len x 1] */ \
/* _h_len : filter length */ \
FIRINTERP() FIRINTERP(_create)(unsigned int _M, \
TC * _h, \
unsigned int _h_len); \
\
/* create interpolator from Kaiser prototype */ \
/* _M : interpolation factor */ \
/* _m : filter delay (symbols) */ \
/* _As : stop-band attenuation [dB] */ \
FIRINTERP() FIRINTERP(_create_kaiser)(unsigned int _M, \
unsigned int _m, \
float _As); \
\
/* create prorotype (root-)Nyquist interpolator */ \
/* _type : filter type (e.g. LIQUID_FIRFILT_RCOS) */ \
/* _k : samples/symbol, _k > 1 */ \
/* _m : filter delay (symbols), _m > 0 */ \
/* _beta : excess bandwidth factor, _beta < 1 */ \
/* _dt : fractional sample delay, _dt in (-1, 1) */ \
FIRINTERP() FIRINTERP(_create_prototype)(int _type, \
unsigned int _k, \
unsigned int _m, \
float _beta, \
float _dt); \
\
/* destroy firinterp object, freeing all internal memory */ \
void FIRINTERP(_destroy)(FIRINTERP() _q); \
\
/* print firinterp object's internal properties to stdout */ \
void FIRINTERP(_print)(FIRINTERP() _q); \
\
/* reset internal state */ \
void FIRINTERP(_reset)(FIRINTERP() _q); \
\
/* execute interpolation on single input sample */ \
/* _q : firinterp object */ \
/* _x : input sample */ \
/* _y : output sample array [size: _M x 1] */ \
void FIRINTERP(_execute)(FIRINTERP() _q, \
TI _x, \
TO * _y); \
\
/* execute interpolation on block of input samples */ \
/* _q : firinterp object */ \
/* _x : input array [size: _n x 1] */ \
/* _n : size of input array */ \
/* _y : output sample array [size: _M*_n x 1] */ \
void FIRINTERP(_execute_block)(FIRINTERP() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
LIQUID_FIRINTERP_DEFINE_API(FIRINTERP_MANGLE_RRRF,
float,
float,
float)
LIQUID_FIRINTERP_DEFINE_API(FIRINTERP_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_FIRINTERP_DEFINE_API(FIRINTERP_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
// iirinterp : infinite impulse response interpolator
#define IIRINTERP_MANGLE_RRRF(name) LIQUID_CONCAT(iirinterp_rrrf,name)
#define IIRINTERP_MANGLE_CRCF(name) LIQUID_CONCAT(iirinterp_crcf,name)
#define IIRINTERP_MANGLE_CCCF(name) LIQUID_CONCAT(iirinterp_cccf,name)
#define LIQUID_IIRINTERP_DEFINE_API(IIRINTERP,TO,TC,TI) \
typedef struct IIRINTERP(_s) * IIRINTERP(); \
\
/* create interpolator from external coefficients */ \
/* _M : interpolation factor */ \
/* _b : feed-back coefficients [size: _nb x 1] */ \
/* _nb : feed-back coefficients length */ \
/* _a : feed-forward coefficients [size: _na x 1] */ \
/* _na : feed-forward coefficients length */ \
IIRINTERP() IIRINTERP(_create)(unsigned int _M, \
TC * _b, \
unsigned int _nb, \
TC * _a, \
unsigned int _na); \
\
/* create decimator with default Butterworth prototype */ \
/* _M : decimation factor */ \
/* _order : filter order */ \
IIRINTERP() IIRINTERP(_create_default)(unsigned int _M, \
unsigned int _order); \
\
/* create interpolator from prototype */ \
/* _M : interpolation factor */ \
IIRINTERP() IIRINTERP(_create_prototype)( \
unsigned int _M, \
liquid_iirdes_filtertype _ftype, \
liquid_iirdes_bandtype _btype, \
liquid_iirdes_format _format, \
unsigned int _order, \
float _fc, \
float _f0, \
float _Ap, \
float _As); \
\
/* destroy interpolator object and free internal memory */ \
void IIRINTERP(_destroy)(IIRINTERP() _q); \
\
/* print interpolator object internals */ \
void IIRINTERP(_print)(IIRINTERP() _q); \
\
/* reset interpolator object */ \
void IIRINTERP(_reset)(IIRINTERP() _q); \
\
/* execute interpolation on single input sample */ \
/* _q : iirinterp object */ \
/* _x : input sample */ \
/* _y : output sample array [size: _M x 1] */ \
void IIRINTERP(_execute)(IIRINTERP() _q, \
TI _x, \
TO * _y); \
\
/* execute interpolation on block of input samples */ \
/* _q : iirinterp object */ \
/* _x : input array [size: _n x 1] */ \
/* _n : size of input array */ \
/* _y : output sample array [size: _M*_n x 1] */ \
void IIRINTERP(_execute_block)(IIRINTERP() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
\
/* get system group delay at frequency _fc */ \
float IIRINTERP(_groupdelay)(IIRINTERP() _q, float _fc); \
LIQUID_IIRINTERP_DEFINE_API(IIRINTERP_MANGLE_RRRF,
float,
float,
float)
LIQUID_IIRINTERP_DEFINE_API(IIRINTERP_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_IIRINTERP_DEFINE_API(IIRINTERP_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Decimators
//
// firdecim : finite impulse response decimator
#define FIRDECIM_MANGLE_RRRF(name) LIQUID_CONCAT(firdecim_rrrf,name)
#define FIRDECIM_MANGLE_CRCF(name) LIQUID_CONCAT(firdecim_crcf,name)
#define FIRDECIM_MANGLE_CCCF(name) LIQUID_CONCAT(firdecim_cccf,name)
#define LIQUID_FIRDECIM_DEFINE_API(FIRDECIM,TO,TC,TI) \
typedef struct FIRDECIM(_s) * FIRDECIM(); \
\
/* create decimator from external coefficients */ \
/* _M : decimation factor */ \
/* _h : filter coefficients [size: _h_len x 1] */ \
/* _h_len : filter coefficients length */ \
FIRDECIM() FIRDECIM(_create)(unsigned int _M, \
TC * _h, \
unsigned int _h_len); \
\
/* create decimator from Kaiser prototype */ \
/* _M : decimation factor */ \
/* _m : filter delay (symbols) */ \
/* _As : stop-band attenuation [dB] */ \
FIRDECIM() FIRDECIM(_create_kaiser)(unsigned int _M, \
unsigned int _m, \
float _As); \
\
/* create square-root Nyquist decimator */ \
/* _type : filter type (e.g. LIQUID_FIRFILT_RRC) */ \
/* _M : samples/symbol (decimation factor) */ \
/* _m : filter delay (symbols) */ \
/* _beta : rolloff factor (0 < beta <= 1) */ \
/* _dt : fractional sample delay */ \
FIRDECIM() FIRDECIM(_create_prototype)(int _type, \
unsigned int _M, \
unsigned int _m, \
float _beta, \
float _dt); \
\
/* destroy decimator object */ \
void FIRDECIM(_destroy)(FIRDECIM() _q); \
\
/* print decimator object propreties to stdout */ \
void FIRDECIM(_print)(FIRDECIM() _q); \
\
/* reset decimator object internal state */ \
void FIRDECIM(_reset)(FIRDECIM() _q); \
\
/* execute decimator on _M input samples */ \
/* _q : decimator object */ \
/* _x : input samples [size: _M x 1] */ \
/* _y : output sample pointer */ \
void FIRDECIM(_execute)(FIRDECIM() _q, \
TI * _x, \
TO * _y); \
\
/* execute decimator on block of _n*_M input samples */ \
/* _q : decimator object */ \
/* _x : input array [size: _n*_M x 1] */ \
/* _n : number of _output_ samples */ \
/* _y : output array [_size: _n x 1] */ \
void FIRDECIM(_execute_block)(FIRDECIM() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
LIQUID_FIRDECIM_DEFINE_API(FIRDECIM_MANGLE_RRRF,
float,
float,
float)
LIQUID_FIRDECIM_DEFINE_API(FIRDECIM_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_FIRDECIM_DEFINE_API(FIRDECIM_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
// iirdecim : infinite impulse response decimator
#define IIRDECIM_MANGLE_RRRF(name) LIQUID_CONCAT(iirdecim_rrrf,name)
#define IIRDECIM_MANGLE_CRCF(name) LIQUID_CONCAT(iirdecim_crcf,name)
#define IIRDECIM_MANGLE_CCCF(name) LIQUID_CONCAT(iirdecim_cccf,name)
#define LIQUID_IIRDECIM_DEFINE_API(IIRDECIM,TO,TC,TI) \
typedef struct IIRDECIM(_s) * IIRDECIM(); \
\
/* create decimator from external coefficients */ \
/* _M : decimation factor */ \
/* _b : feed-back coefficients [size: _nb x 1] */ \
/* _nb : feed-back coefficients length */ \
/* _a : feed-forward coefficients [size: _na x 1] */ \
/* _na : feed-forward coefficients length */ \
IIRDECIM() IIRDECIM(_create)(unsigned int _M, \
TC * _b, \
unsigned int _nb, \
TC * _a, \
unsigned int _na); \
\
/* create decimator with default Butterworth prototype */ \
/* _M : decimation factor */ \
/* _order : filter order */ \
IIRDECIM() IIRDECIM(_create_default)(unsigned int _M, \
unsigned int _order); \
\
/* create decimator from prototype */ \
/* _M : decimation factor */ \
/* _ftype : filter type (e.g. LIQUID_IIRDES_BUTTER) */ \
/* _btype : band type (e.g. LIQUID_IIRDES_BANDPASS) */ \
/* _format : coefficients format (e.g. LIQUID_IIRDES_SOS) */ \
/* _n : filter order */ \
/* _fc : low-pass prototype cut-off frequency */ \
/* _f0 : center frequency (band-pass, band-stop) */ \
/* _Ap : pass-band ripple in dB */ \
/* _As : stop-band ripple in dB */ \
IIRDECIM() IIRDECIM(_create_prototype)( \
unsigned int _M, \
liquid_iirdes_filtertype _ftype, \
liquid_iirdes_bandtype _btype, \
liquid_iirdes_format _format, \
unsigned int _order, \
float _fc, \
float _f0, \
float _Ap, \
float _As); \
\
/* destroy decimator object and free internal memory */ \
void IIRDECIM(_destroy)(IIRDECIM() _q); \
\
/* print decimator object internals */ \
void IIRDECIM(_print)(IIRDECIM() _q); \
\
/* reset decimator object */ \
void IIRDECIM(_reset)(IIRDECIM() _q); \
\
/* execute decimator on _M input samples */ \
/* _q : decimator object */ \
/* _x : input samples [size: _M x 1] */ \
/* _y : output sample pointer */ \
void IIRDECIM(_execute)(IIRDECIM() _q, \
TI * _x, \
TO * _y); \
\
/* execute decimator on block of _n*_M input samples */ \
/* _q : decimator object */ \
/* _x : input array [size: _n*_M x 1] */ \
/* _n : number of _output_ samples */ \
/* _y : output array [_sze: _n x 1] */ \
void IIRDECIM(_execute_block)(IIRDECIM() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
\
/* get system group delay at frequency _fc */ \
float IIRDECIM(_groupdelay)(IIRDECIM() _q, float _fc); \
LIQUID_IIRDECIM_DEFINE_API(IIRDECIM_MANGLE_RRRF,
float,
float,
float)
LIQUID_IIRDECIM_DEFINE_API(IIRDECIM_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_IIRDECIM_DEFINE_API(IIRDECIM_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Half-band resampler
//
#define RESAMP2_MANGLE_RRRF(name) LIQUID_CONCAT(resamp2_rrrf,name)
#define RESAMP2_MANGLE_CRCF(name) LIQUID_CONCAT(resamp2_crcf,name)
#define RESAMP2_MANGLE_CCCF(name) LIQUID_CONCAT(resamp2_cccf,name)
#define LIQUID_RESAMP2_DEFINE_API(RESAMP2,TO,TC,TI) \
typedef struct RESAMP2(_s) * RESAMP2(); \
\
/* create half-band resampler */ \
/* _m : filter semi-length (h_len = 4*m+1) */ \
/* _f0 : filter center frequency */ \
/* _As : stop-band attenuation [dB] */ \
RESAMP2() RESAMP2(_create)(unsigned int _m, \
float _f0, \
float _As); \
\
/* re-create half-band resampler with new properties */ \
/* _q : original half-band resampler object */ \
/* _m : filter semi-length (h_len = 4*m+1) */ \
/* _f0 : filter center frequency */ \
/* _As : stop-band attenuation [dB] */ \
RESAMP2() RESAMP2(_recreate)(RESAMP2() _q, \
unsigned int _m, \
float _f0, \
float _As); \
\
/* destroy half-band resampler */ \
void RESAMP2(_destroy)(RESAMP2() _q); \
\
/* print resamp2 object's internals */ \
void RESAMP2(_print)(RESAMP2() _q); \
\
/* reset internal buffer */ \
void RESAMP2(_reset)(RESAMP2() _q); \
\
/* get resampler filter delay (semi-length m) */ \
unsigned int RESAMP2(_get_delay)(RESAMP2() _q); \
\
/* execute resamp2 as half-band filter */ \
/* _q : resamp2 object */ \
/* _x : input sample */ \
/* _y0 : output sample pointer (low frequency) */ \
/* _y1 : output sample pointer (high frequency) */ \
void RESAMP2(_filter_execute)(RESAMP2() _q, \
TI _x, \
TO * _y0, \
TO * _y1); \
\
/* execute resamp2 as half-band analysis filterbank */ \
/* _q : resamp2 object */ \
/* _x : input array [size: 2 x 1] */ \
/* _y : output array [size: 2 x 1] */ \
void RESAMP2(_analyzer_execute)(RESAMP2() _q, \
TI * _x, \
TO * _y); \
\
/* execute resamp2 as half-band synthesis filterbank */ \
/* _q : resamp2 object */ \
/* _x : input array [size: 2 x 1] */ \
/* _y : output array [size: 2 x 1] */ \
void RESAMP2(_synthesizer_execute)(RESAMP2() _q, \
TI * _x, \
TO * _y); \
\
/* execute resamp2 as half-band decimator */ \
/* _q : resamp2 object */ \
/* _x : input array [size: 2 x 1] */ \
/* _y : output sample pointer */ \
void RESAMP2(_decim_execute)(RESAMP2() _q, \
TI * _x, \
TO * _y); \
\
/* execute resamp2 as half-band interpolator */ \
/* _q : resamp2 object */ \
/* _x : input sample */ \
/* _y : output array [size: 2 x 1] */ \
void RESAMP2(_interp_execute)(RESAMP2() _q, \
TI _x, \
TO * _y); \
LIQUID_RESAMP2_DEFINE_API(RESAMP2_MANGLE_RRRF,
float,
float,
float)
LIQUID_RESAMP2_DEFINE_API(RESAMP2_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_RESAMP2_DEFINE_API(RESAMP2_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Arbitrary resampler
//
#define RESAMP_MANGLE_RRRF(name) LIQUID_CONCAT(resamp_rrrf,name)
#define RESAMP_MANGLE_CRCF(name) LIQUID_CONCAT(resamp_crcf,name)
#define RESAMP_MANGLE_CCCF(name) LIQUID_CONCAT(resamp_cccf,name)
#define LIQUID_RESAMP_DEFINE_API(RESAMP,TO,TC,TI) \
typedef struct RESAMP(_s) * RESAMP(); \
\
/* create arbitrary resampler object */ \
/* _rate : arbitrary resampling rate */ \
/* _m : filter semi-length (delay) */ \
/* _fc : filter cutoff frequency, 0 < _fc < 0.5 */ \
/* _As : filter stop-band attenuation [dB] */ \
/* _npfb : number of filters in the bank */ \
RESAMP() RESAMP(_create)(float _rate, \
unsigned int _m, \
float _fc, \
float _As, \
unsigned int _npfb); \
\
/* create arbitrary resampler object with a specified input */ \
/* resampling rate and default parameters */ \
/* m (filter semi-length) = 7 */ \
/* fc (filter cutoff frequency) = 0.25 */ \
/* As (filter stop-band attenuation) = 60 dB */ \
/* npfb (number of filters in the bank) = 64 */ \
RESAMP() RESAMP(_create_default)(float _rate); \
\
/* destroy arbitrary resampler object */ \
void RESAMP(_destroy)(RESAMP() _q); \
\
/* print resamp object internals to stdout */ \
void RESAMP(_print)(RESAMP() _q); \
\
/* reset resamp object internals */ \
void RESAMP(_reset)(RESAMP() _q); \
\
/* get resampler delay (output samples) */ \
unsigned int RESAMP(_get_delay)(RESAMP() _q); \
\
/* set rate of arbitrary resampler */ \
void RESAMP(_set_rate)(RESAMP() _q, float _rate); \
\
/* adjust rate of arbitrary resampler */ \
void RESAMP(_adjust_rate)(RESAMP() _q, float _delta); \
\
/* execute arbitrary resampler */ \
/* _q : resamp object */ \
/* _x : single input sample */ \
/* _y : output sample array (pointer) */ \
/* _num_written : number of samples written to _y */ \
void RESAMP(_execute)(RESAMP() _q, \
TI _x, \
TO * _y, \
unsigned int * _num_written); \
\
/* execute arbitrary resampler on a block of samples */ \
/* _q : resamp object */ \
/* _x : input buffer [size: _nx x 1] */ \
/* _nx : input buffer */ \
/* _y : output sample array (pointer) */ \
/* _ny : number of samples written to _y */ \
void RESAMP(_execute_block)(RESAMP() _q, \
TI * _x, \
unsigned int _nx, \
TO * _y, \
unsigned int * _ny); \
LIQUID_RESAMP_DEFINE_API(RESAMP_MANGLE_RRRF,
float,
float,
float)
LIQUID_RESAMP_DEFINE_API(RESAMP_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_RESAMP_DEFINE_API(RESAMP_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Multi-stage half-band resampler
//
// resampling type (interpolator/decimator)
typedef enum {
LIQUID_RESAMP_INTERP=0, // interpolator
LIQUID_RESAMP_DECIM, // decimator
} liquid_resamp_type;
#define MSRESAMP2_MANGLE_RRRF(name) LIQUID_CONCAT(msresamp2_rrrf,name)
#define MSRESAMP2_MANGLE_CRCF(name) LIQUID_CONCAT(msresamp2_crcf,name)
#define MSRESAMP2_MANGLE_CCCF(name) LIQUID_CONCAT(msresamp2_cccf,name)
#define LIQUID_MSRESAMP2_DEFINE_API(MSRESAMP2,TO,TC,TI) \
typedef struct MSRESAMP2(_s) * MSRESAMP2(); \
\
/* create multi-stage half-band resampler */ \
/* _type : resampler type (e.g. LIQUID_RESAMP_DECIM) */ \
/* _num_stages : number of resampling stages */ \
/* _fc : filter cut-off frequency 0 < _fc < 0.5 */ \
/* _f0 : filter center frequency */ \
/* _As : stop-band attenuation [dB] */ \
MSRESAMP2() MSRESAMP2(_create)(int _type, \
unsigned int _num_stages, \
float _fc, \
float _f0, \
float _As); \
\
/* destroy multi-stage half-bandresampler */ \
void MSRESAMP2(_destroy)(MSRESAMP2() _q); \
\
/* print msresamp object internals to stdout */ \
void MSRESAMP2(_print)(MSRESAMP2() _q); \
\
/* reset msresamp object internal state */ \
void MSRESAMP2(_reset)(MSRESAMP2() _q); \
\
/* get group delay (number of output samples) */ \
float MSRESAMP2(_get_delay)(MSRESAMP2() _q); \
\
/* execute multi-stage resampler, M = 2^num_stages */ \
/* LIQUID_RESAMP_INTERP: input: 1, output: M */ \
/* LIQUID_RESAMP_DECIM: input: M, output: 1 */ \
/* _q : msresamp object */ \
/* _x : input sample array */ \
/* _y : output sample array */ \
void MSRESAMP2(_execute)(MSRESAMP2() _q, \
TI * _x, \
TO * _y); \
LIQUID_MSRESAMP2_DEFINE_API(MSRESAMP2_MANGLE_RRRF,
float,
float,
float)
LIQUID_MSRESAMP2_DEFINE_API(MSRESAMP2_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_MSRESAMP2_DEFINE_API(MSRESAMP2_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Multi-stage arbitrary resampler
//
#define MSRESAMP_MANGLE_RRRF(name) LIQUID_CONCAT(msresamp_rrrf,name)
#define MSRESAMP_MANGLE_CRCF(name) LIQUID_CONCAT(msresamp_crcf,name)
#define MSRESAMP_MANGLE_CCCF(name) LIQUID_CONCAT(msresamp_cccf,name)
#define LIQUID_MSRESAMP_DEFINE_API(MSRESAMP,TO,TC,TI) \
typedef struct MSRESAMP(_s) * MSRESAMP(); \
\
/* create multi-stage arbitrary resampler */ \
/* _r : resampling rate [output/input] */ \
/* _As : stop-band attenuation [dB] */ \
MSRESAMP() MSRESAMP(_create)(float _r, \
float _As); \
\
/* destroy multi-stage arbitrary resampler */ \
void MSRESAMP(_destroy)(MSRESAMP() _q); \
\
/* print msresamp object internals to stdout */ \
void MSRESAMP(_print)(MSRESAMP() _q); \
\
/* reset msresamp object internal state */ \
void MSRESAMP(_reset)(MSRESAMP() _q); \
\
/* get filter delay (output samples) */ \
float MSRESAMP(_get_delay)(MSRESAMP() _q); \
\
/* execute multi-stage resampler */ \
/* _q : msresamp object */ \
/* _x : input sample array [size: _nx x 1] */ \
/* _nx : input sample array size */ \
/* _y : output sample array [size: variable] */ \
/* _ny : number of samples written to _y */ \
void MSRESAMP(_execute)(MSRESAMP() _q, \
TI * _x, \
unsigned int _nx, \
TO * _y, \
unsigned int * _ny); \
LIQUID_MSRESAMP_DEFINE_API(MSRESAMP_MANGLE_RRRF,
float,
float,
float)
LIQUID_MSRESAMP_DEFINE_API(MSRESAMP_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_MSRESAMP_DEFINE_API(MSRESAMP_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Symbol timing recovery (symbol synchronizer)
//
#define SYMSYNC_MANGLE_RRRF(name) LIQUID_CONCAT(symsync_rrrf,name)
#define SYMSYNC_MANGLE_CRCF(name) LIQUID_CONCAT(symsync_crcf,name)
#define LIQUID_SYMSYNC_DEFINE_API(SYMSYNC,TO,TC,TI) \
\
typedef struct SYMSYNC(_s) * SYMSYNC(); \
\
/* create synchronizer object from external coefficients */ \
/* _k : samples per symbol */ \
/* _M : number of filters in the bank */ \
/* _h : matched filter coefficients [size: */ \
/* _h_len : length of matched filter */ \
SYMSYNC() SYMSYNC(_create)(unsigned int _k, \
unsigned int _M, \
TC * _h, \
unsigned int _h_len); \
\
/* create square-root Nyquist symbol synchronizer */ \
/* _type : filter type (e.g. LIQUID_FIRFILT_RRC) */ \
/* _k : samples/symbol */ \
/* _m : symbol delay */ \
/* _beta : rolloff factor, beta in (0,1] */ \
/* _M : number of filters in the bank */ \
SYMSYNC() SYMSYNC(_create_rnyquist)(int _type, \
unsigned int _k, \
unsigned int _m, \
float _beta, \
unsigned int _M); \
\
/* create symsync using Kaiser filter interpolator; useful */ \
/* when the input signal has matched filter applied already */ \
/* _k : input samples/symbol */ \
/* _m : symbol delay */ \
/* _beta : rolloff factor, beta in (0,1] */ \
/* _M : number of filters in the bank */ \
SYMSYNC() SYMSYNC(_create_kaiser)(unsigned int _k, \
unsigned int _m, \
float _beta, \
unsigned int _M); \
\
/* destroy symsync object, freeing all internal memory */ \
void SYMSYNC(_destroy)(SYMSYNC() _q); \
\
/* print symsync object's parameters */ \
void SYMSYNC(_print)(SYMSYNC() _q); \
\
/* reset symsync internal state */ \
void SYMSYNC(_reset)(SYMSYNC() _q); \
\
/* lock/unlock loop control */ \
void SYMSYNC(_lock)( SYMSYNC() _q); \
void SYMSYNC(_unlock)(SYMSYNC() _q); \
\
/* set synchronizer output rate (samples/symbol) */ \
/* _q : synchronizer object */ \
/* _k_out : output samples/symbol */ \
void SYMSYNC(_set_output_rate)(SYMSYNC() _q, \
unsigned int _k_out); \
\
/* set loop-filter bandwidth */ \
/* _q : synchronizer object */ \
/* _bt : loop bandwidth */ \
void SYMSYNC(_set_lf_bw)(SYMSYNC() _q, \
float _bt); \
\
/* return instantaneous fractional timing offset estimate */ \
float SYMSYNC(_get_tau)(SYMSYNC() _q); \
\
/* execute synchronizer on input data array */ \
/* _q : synchronizer object */ \
/* _x : input data array */ \
/* _nx : number of input samples */ \
/* _y : output data array */ \
/* _ny : number of samples written to output buffer */ \
void SYMSYNC(_execute)(SYMSYNC() _q, \
TI * _x, \
unsigned int _nx, \
TO * _y, \
unsigned int * _ny); \
LIQUID_SYMSYNC_DEFINE_API(SYMSYNC_MANGLE_RRRF,
float,
float,
float)
LIQUID_SYMSYNC_DEFINE_API(SYMSYNC_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
//
// Finite impulse response Farrow filter
//
#define FIRFARROW_MANGLE_RRRF(name) LIQUID_CONCAT(firfarrow_rrrf,name)
#define FIRFARROW_MANGLE_CRCF(name) LIQUID_CONCAT(firfarrow_crcf,name)
//#define FIRFARROW_MANGLE_CCCF(name) LIQUID_CONCAT(firfarrow_cccf,name)
// Macro:
// FIRFARROW : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_FIRFARROW_DEFINE_API(FIRFARROW,TO,TC,TI) \
\
typedef struct FIRFARROW(_s) * FIRFARROW(); \
\
/* create firfarrow object */ \
/* _h_len : filter length */ \
/* _p : polynomial order */ \
/* _fc : filter cutoff frequency */ \
/* _As : stopband attenuation [dB] */ \
FIRFARROW() FIRFARROW(_create)(unsigned int _h_len, \
unsigned int _p, \
float _fc, \
float _As); \
\
/* destroy firfarrow object, freeing all internal memory */ \
void FIRFARROW(_destroy)(FIRFARROW() _q); \
\
/* print firfarrow object's internal properties */ \
void FIRFARROW(_print)(FIRFARROW() _q); \
\
/* reset firfarrow object's internal state */ \
void FIRFARROW(_reset)(FIRFARROW() _q); \
\
/* push sample into firfarrow object */ \
/* _q : firfarrow object */ \
/* _x : input sample */ \
void FIRFARROW(_push)(FIRFARROW() _q, \
TI _x); \
\
/* set fractional delay of firfarrow object */ \
/* _q : firfarrow object */ \
/* _mu : fractional sample delay */ \
void FIRFARROW(_set_delay)(FIRFARROW() _q, \
float _mu); \
\
/* execute firfarrow internal dot product */ \
/* _q : firfarrow object */ \
/* _y : output sample pointer */ \
void FIRFARROW(_execute)(FIRFARROW() _q, \
TO * _y); \
\
/* compute firfarrow filter on block of samples; the input */ \
/* and output arrays may have the same pointer */ \
/* _q : firfarrow object */ \
/* _x : input array [size: _n x 1] */ \
/* _n : input, output array size */ \
/* _y : output array [size: _n x 1] */ \
void FIRFARROW(_execute_block)(FIRFARROW() _q, \
TI * _x, \
unsigned int _n, \
TO * _y); \
\
/* get length of firfarrow object (number of filter taps) */ \
unsigned int FIRFARROW(_get_length)(FIRFARROW() _q); \
\
/* get coefficients of firfarrow object */ \
/* _q : firfarrow object */ \
/* _h : output coefficients pointer */ \
void FIRFARROW(_get_coefficients)(FIRFARROW() _q, \
float * _h); \
\
/* compute complex frequency response */ \
/* _q : filter object */ \
/* _fc : frequency */ \
/* _H : output frequency response */ \
void FIRFARROW(_freqresponse)(FIRFARROW() _q, \
float _fc, \
liquid_float_complex * _H); \
\
/* compute group delay [samples] */ \
/* _q : filter object */ \
/* _fc : frequency */ \
float FIRFARROW(_groupdelay)(FIRFARROW() _q, \
float _fc); \
LIQUID_FIRFARROW_DEFINE_API(FIRFARROW_MANGLE_RRRF,
float,
float,
float)
LIQUID_FIRFARROW_DEFINE_API(FIRFARROW_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
//
// MODULE : framing
//
// framesyncstats : generic frame synchronizer statistic structure
typedef struct {
// signal quality
float evm; // error vector magnitude [dB]
float rssi; // received signal strength indicator [dB]
float cfo; // carrier frequency offset (f/Fs)
// demodulated frame symbols
liquid_float_complex * framesyms; // pointer to array [size: framesyms x 1]
unsigned int num_framesyms; // length of framesyms
// modulation/coding scheme etc.
unsigned int mod_scheme; // modulation scheme
unsigned int mod_bps; // modulation depth (bits/symbol)
unsigned int check; // data validity check (crc, checksum)
unsigned int fec0; // forward error-correction (inner)
unsigned int fec1; // forward error-correction (outer)
} framesyncstats_s;
// external framesyncstats default object
extern framesyncstats_s framesyncstats_default;
// initialize framesyncstats object on default
void framesyncstats_init_default(framesyncstats_s * _stats);
// print framesyncstats object
void framesyncstats_print(framesyncstats_s * _stats);
// framedatastats : gather frame data
typedef struct {
unsigned int num_frames_detected;
unsigned int num_headers_valid;
unsigned int num_payloads_valid;
unsigned long int num_bytes_received;
} framedatastats_s;
// reset framedatastats object
void framedatastats_reset(framedatastats_s * _stats);
// print framedatastats object
void framedatastats_print(framedatastats_s * _stats);
// Generic frame synchronizer callback function type
// _header : header data [size: 8 bytes]
// _header_valid : is header valid? (0:no, 1:yes)
// _payload : payload data [size: _payload_len]
// _payload_len : length of payload (bytes)
// _payload_valid : is payload valid? (0:no, 1:yes)
// _stats : frame statistics object
// _userdata : pointer to userdata
typedef int (*framesync_callback)(unsigned char * _header,
int _header_valid,
unsigned char * _payload,
unsigned int _payload_len,
int _payload_valid,
framesyncstats_s _stats,
void * _userdata);
// framesync csma callback functions invoked when signal levels is high or low
// _userdata : user-defined data pointer
typedef void (*framesync_csma_callback)(void * _userdata);
//
// packet encoder/decoder
//
typedef struct qpacketmodem_s * qpacketmodem;
// create packet encoder
qpacketmodem qpacketmodem_create ();
void qpacketmodem_destroy(qpacketmodem _q);
void qpacketmodem_reset (qpacketmodem _q);
void qpacketmodem_print (qpacketmodem _q);
int qpacketmodem_configure(qpacketmodem _q,
unsigned int _payload_len,
crc_scheme _check,
fec_scheme _fec0,
fec_scheme _fec1,
int _ms);
// get length of encoded frame in symbols
unsigned int qpacketmodem_get_frame_len(qpacketmodem _q);
// get unencoded/decoded payload length (bytes)
unsigned int qpacketmodem_get_payload_len(qpacketmodem _q);
// regular access methods
unsigned int qpacketmodem_get_crc (qpacketmodem _q);
unsigned int qpacketmodem_get_fec0 (qpacketmodem _q);
unsigned int qpacketmodem_get_fec1 (qpacketmodem _q);
unsigned int qpacketmodem_get_modscheme(qpacketmodem _q);
// encode packet into un-modulated frame symbol indices
// _q : qpacketmodem object
// _payload : unencoded payload bytes
// _syms : encoded but un-modulated payload symbol indices
void qpacketmodem_encode_syms(qpacketmodem _q,
const unsigned char * _payload,
unsigned char * _syms);
// decode packet from demodulated frame symbol indices (hard-decision decoding)
// _q : qpacketmodem object
// _syms : received hard-decision symbol indices [size: frame_len x 1]
// _payload : recovered decoded payload bytes
int qpacketmodem_decode_syms(qpacketmodem _q,
unsigned char * _syms,
unsigned char * _payload);
// decode packet from demodulated frame bits (soft-decision decoding)
// _q : qpacketmodem object
// _bits : received soft-decision bits, [size: bps*frame_len x 1]
// _payload : recovered decoded payload bytes
int qpacketmodem_decode_bits(qpacketmodem _q,
unsigned char * _bits,
unsigned char * _payload);
// encode and modulate packet into modulated frame samples
// _q : qpacketmodem object
// _payload : unencoded payload bytes
// _frame : encoded/modulated payload symbols
void qpacketmodem_encode(qpacketmodem _q,
const unsigned char * _payload,
liquid_float_complex * _frame);
// decode packet from modulated frame samples, returning flag if CRC passed
// NOTE: hard-decision decoding
// _q : qpacketmodem object
// _frame : encoded/modulated payload symbols
// _payload : recovered decoded payload bytes
int qpacketmodem_decode(qpacketmodem _q,
liquid_float_complex * _frame,
unsigned char * _payload);
// decode packet from modulated frame samples, returning flag if CRC passed
// NOTE: soft-decision decoding
// _q : qpacketmodem object
// _frame : encoded/modulated payload symbols
// _payload : recovered decoded payload bytes
int qpacketmodem_decode_soft(qpacketmodem _q,
liquid_float_complex * _frame,
unsigned char * _payload);
//
// pilot generator for streaming applications
//
typedef struct qpilotgen_s * qpilotgen;
// create packet encoder
qpilotgen qpilotgen_create(unsigned int _payload_len,
unsigned int _pilot_spacing);
qpilotgen qpilotgen_recreate(qpilotgen _q,
unsigned int _payload_len,
unsigned int _pilot_spacing);
void qpilotgen_destroy(qpilotgen _q);
void qpilotgen_reset( qpilotgen _q);
void qpilotgen_print( qpilotgen _q);
unsigned int qpilotgen_get_frame_len(qpilotgen _q);
// insert pilot symbols
void qpilotgen_execute(qpilotgen _q,
liquid_float_complex * _payload,
liquid_float_complex * _frame);
//
// pilot synchronizer for streaming applications
//
typedef struct qpilotsync_s * qpilotsync;
// create packet encoder
qpilotsync qpilotsync_create(unsigned int _payload_len,
unsigned int _pilot_spacing);
qpilotsync qpilotsync_recreate(qpilotsync _q,
unsigned int _payload_len,
unsigned int _pilot_spacing);
void qpilotsync_destroy(qpilotsync _q);
void qpilotsync_reset( qpilotsync _q);
void qpilotsync_print( qpilotsync _q);
unsigned int qpilotsync_get_frame_len(qpilotsync _q);
// recover frame symbols from received frame
void qpilotsync_execute(qpilotsync _q,
liquid_float_complex * _frame,
liquid_float_complex * _payload);
// get estimates
float qpilotsync_get_dphi(qpilotsync _q);
float qpilotsync_get_phi (qpilotsync _q);
float qpilotsync_get_gain(qpilotsync _q);
//
// Basic frame generator (64 bytes data payload)
//
// frame length in samples
#define LIQUID_FRAME64_LEN (1440)
typedef struct framegen64_s * framegen64;
// create frame generator
framegen64 framegen64_create();
// destroy frame generator
void framegen64_destroy(framegen64 _q);
// print frame generator internal properties
void framegen64_print(framegen64 _q);
// generate frame
// _q : frame generator object
// _header : 8-byte header data
// _payload : 64-byte payload data
// _frame : output frame samples [size: LIQUID_FRAME64_LEN x 1]
void framegen64_execute(framegen64 _q,
unsigned char * _header,
unsigned char * _payload,
liquid_float_complex * _frame);
typedef struct framesync64_s * framesync64;
// create framesync64 object
// _callback : callback function
// _userdata : user data pointer passed to callback function
framesync64 framesync64_create(framesync_callback _callback,
void * _userdata);
// destroy frame synchronizer
void framesync64_destroy(framesync64 _q);
// print frame synchronizer internal properties
void framesync64_print(framesync64 _q);
// reset frame synchronizer internal state
void framesync64_reset(framesync64 _q);
// push samples through frame synchronizer
// _q : frame synchronizer object
// _x : input samples [size: _n x 1]
// _n : number of input samples
void framesync64_execute(framesync64 _q,
liquid_float_complex * _x,
unsigned int _n);
// enable/disable debugging
void framesync64_debug_enable(framesync64 _q);
void framesync64_debug_disable(framesync64 _q);
void framesync64_debug_print(framesync64 _q, const char * _filename);
#if 0
// advanced modes
void framesync64_set_csma_callbacks(framesync64 _q,
framesync_csma_callback _csma_lock,
framesync_csma_callback _csma_unlock,
void * _csma_userdata);
#endif
//
// Flexible frame : adjustable payload, mod scheme, etc., but bring
// your own error correction, redundancy check
//
// frame generator
typedef struct {
unsigned int check; // data validity check
unsigned int fec0; // forward error-correction scheme (inner)
unsigned int fec1; // forward error-correction scheme (outer)
unsigned int mod_scheme; // modulation scheme
} flexframegenprops_s;
void flexframegenprops_init_default(flexframegenprops_s * _fgprops);
typedef struct flexframegen_s * flexframegen;
// create flexframegen object
// _props : frame properties (modulation scheme, etc.)
flexframegen flexframegen_create(flexframegenprops_s * _props);
// destroy flexframegen object
void flexframegen_destroy(flexframegen _q);
// print flexframegen object internals
void flexframegen_print(flexframegen _q);
// reset flexframegen object internals
void flexframegen_reset(flexframegen _q);
// is frame assembled?
int flexframegen_is_assembled(flexframegen _q);
// get frame properties
void flexframegen_getprops(flexframegen _q, flexframegenprops_s * _props);
// set frame properties
int flexframegen_setprops(flexframegen _q, flexframegenprops_s * _props);
// get length of assembled frame (samples)
unsigned int flexframegen_getframelen(flexframegen _q);
// assemble a frame from an array of data
// _q : frame generator object
// _header : frame header
// _payload : payload data [size: _payload_len x 1]
// _payload_len : payload data length
void flexframegen_assemble(flexframegen _q,
const unsigned char * _header,
const unsigned char * _payload,
unsigned int _payload_len);
// write samples of assembled frame, two samples at a time, returning
// '1' when frame is complete, '0' otherwise. Zeros will be written
// to the buffer if the frame is not assembled
// _q : frame generator object
// _buffer : output buffer [size: _buffer_len x 1]
// _buffer_len : output buffer length
int flexframegen_write_samples(flexframegen _q,
liquid_float_complex * _buffer,
unsigned int _buffer_len);
// frame synchronizer
typedef struct flexframesync_s * flexframesync;
// create flexframesync object
// _callback : callback function
// _userdata : user data pointer passed to callback function
flexframesync flexframesync_create(framesync_callback _callback,
void * _userdata);
// destroy frame synchronizer
void flexframesync_destroy(flexframesync _q);
// print frame synchronizer internal properties
void flexframesync_print(flexframesync _q);
// reset frame synchronizer internal state
void flexframesync_reset(flexframesync _q);
// push samples through frame synchronizer
// _q : frame synchronizer object
// _x : input samples [size: _n x 1]
// _n : number of input samples
void flexframesync_execute(flexframesync _q,
liquid_float_complex * _x,
unsigned int _n);
// frame data statistics
void flexframesync_reset_framedatastats(flexframesync _q);
framedatastats_s flexframesync_get_framedatastats (flexframesync _q);
// enable/disable debugging
void flexframesync_debug_enable(flexframesync _q);
void flexframesync_debug_disable(flexframesync _q);
void flexframesync_debug_print(flexframesync _q,
const char * _filename);
//
// bpacket : binary packet suitable for data streaming
//
//
// bpacket generator/encoder
//
typedef struct bpacketgen_s * bpacketgen;
// create bpacketgen object
// _m : p/n sequence length (ignored)
// _dec_msg_len : decoded message length (original uncoded data)
// _crc : data validity check (e.g. cyclic redundancy check)
// _fec0 : inner forward error-correction code scheme
// _fec1 : outer forward error-correction code scheme
bpacketgen bpacketgen_create(unsigned int _m,
unsigned int _dec_msg_len,
int _crc,
int _fec0,
int _fec1);
// re-create bpacketgen object from old object
// _q : old bpacketgen object
// _m : p/n sequence length (ignored)
// _dec_msg_len : decoded message length (original uncoded data)
// _crc : data validity check (e.g. cyclic redundancy check)
// _fec0 : inner forward error-correction code scheme
// _fec1 : outer forward error-correction code scheme
bpacketgen bpacketgen_recreate(bpacketgen _q,
unsigned int _m,
unsigned int _dec_msg_len,
int _crc,
int _fec0,
int _fec1);
// destroy bpacketgen object, freeing all internally-allocated memory
void bpacketgen_destroy(bpacketgen _q);
// print bpacketgen internals
void bpacketgen_print(bpacketgen _q);
// return length of full packet
unsigned int bpacketgen_get_packet_len(bpacketgen _q);
// encode packet
void bpacketgen_encode(bpacketgen _q,
unsigned char * _msg_dec,
unsigned char * _packet);
//
// bpacket synchronizer/decoder
//
typedef struct bpacketsync_s * bpacketsync;
typedef int (*bpacketsync_callback)(unsigned char * _payload,
int _payload_valid,
unsigned int _payload_len,
framesyncstats_s _stats,
void * _userdata);
bpacketsync bpacketsync_create(unsigned int _m,
bpacketsync_callback _callback,
void * _userdata);
void bpacketsync_destroy(bpacketsync _q);
void bpacketsync_print(bpacketsync _q);
void bpacketsync_reset(bpacketsync _q);
// run synchronizer on array of input bytes
// _q : bpacketsync object
// _bytes : input data array [size: _n x 1]
// _n : input array size
void bpacketsync_execute(bpacketsync _q,
unsigned char * _bytes,
unsigned int _n);
// run synchronizer on input byte
// _q : bpacketsync object
// _byte : input byte
void bpacketsync_execute_byte(bpacketsync _q,
unsigned char _byte);
// run synchronizer on input symbol
// _q : bpacketsync object
// _sym : input symbol with _bps significant bits
// _bps : number of bits in input symbol
void bpacketsync_execute_sym(bpacketsync _q,
unsigned char _sym,
unsigned int _bps);
// execute one bit at a time
void bpacketsync_execute_bit(bpacketsync _q,
unsigned char _bit);
//
// GMSK frame generator
//
typedef struct gmskframegen_s * gmskframegen;
// create GMSK frame generator
gmskframegen gmskframegen_create();
void gmskframegen_destroy (gmskframegen _q);
int gmskframegen_is_assembled (gmskframegen _q);
void gmskframegen_print (gmskframegen _q);
void gmskframegen_reset (gmskframegen _q);
void gmskframegen_assemble (gmskframegen _q,
const unsigned char * _header,
const unsigned char * _payload,
unsigned int _payload_len,
crc_scheme _check,
fec_scheme _fec0,
fec_scheme _fec1);
unsigned int gmskframegen_getframelen(gmskframegen _q);
int gmskframegen_write_samples(gmskframegen _q,
liquid_float_complex * _y);
//
// GMSK frame synchronizer
//
typedef struct gmskframesync_s * gmskframesync;
// create GMSK frame synchronizer
// _callback : callback function
// _userdata : user data pointer passed to callback function
gmskframesync gmskframesync_create(framesync_callback _callback,
void * _userdata);
void gmskframesync_destroy(gmskframesync _q);
void gmskframesync_print(gmskframesync _q);
void gmskframesync_reset(gmskframesync _q);
void gmskframesync_execute(gmskframesync _q,
liquid_float_complex * _x,
unsigned int _n);
// debugging
void gmskframesync_debug_enable(gmskframesync _q);
void gmskframesync_debug_disable(gmskframesync _q);
void gmskframesync_debug_print(gmskframesync _q, const char * _filename);
//
// OFDM flexframe generator
//
// ofdm frame generator properties
typedef struct {
unsigned int check; // data validity check
unsigned int fec0; // forward error-correction scheme (inner)
unsigned int fec1; // forward error-correction scheme (outer)
unsigned int mod_scheme; // modulation scheme
//unsigned int block_size; // framing block size
} ofdmflexframegenprops_s;
void ofdmflexframegenprops_init_default(ofdmflexframegenprops_s * _props);
typedef struct ofdmflexframegen_s * ofdmflexframegen;
// create OFDM flexible framing generator object
// _M : number of subcarriers, >10 typical
// _cp_len : cyclic prefix length
// _taper_len : taper length (OFDM symbol overlap)
// _p : subcarrier allocation (null, pilot, data), [size: _M x 1]
// _fgprops : frame properties (modulation scheme, etc.)
ofdmflexframegen ofdmflexframegen_create(unsigned int _M,
unsigned int _cp_len,
unsigned int _taper_len,
unsigned char * _p,
ofdmflexframegenprops_s * _fgprops);
// destroy ofdmflexframegen object
void ofdmflexframegen_destroy(ofdmflexframegen _q);
// print parameters, properties, etc.
void ofdmflexframegen_print(ofdmflexframegen _q);
// reset ofdmflexframegen object internals
void ofdmflexframegen_reset(ofdmflexframegen _q);
// is frame assembled?
int ofdmflexframegen_is_assembled(ofdmflexframegen _q);
// get properties
void ofdmflexframegen_getprops(ofdmflexframegen _q,
ofdmflexframegenprops_s * _props);
// set properties
void ofdmflexframegen_setprops(ofdmflexframegen _q,
ofdmflexframegenprops_s * _props);
// get length of frame (symbols)
// _q : OFDM frame generator object
unsigned int ofdmflexframegen_getframelen(ofdmflexframegen _q);
// assemble a frame from an array of data (NULL pointers will use random data)
// _q : OFDM frame generator object
// _header : frame header [8 bytes]
// _payload : payload data [size: _payload_len x 1]
// _payload_len : payload data length
void ofdmflexframegen_assemble(ofdmflexframegen _q,
const unsigned char * _header,
const unsigned char * _payload,
unsigned int _payload_len);
// write samples of assembled frame
// _q : OFDM frame generator object
// _buf : output buffer [size: _buf_len x 1]
// _buf_len : output buffer length
int ofdmflexframegen_write(ofdmflexframegen _q,
liquid_float_complex * _buf,
unsigned int _buf_len);
//
// OFDM flex frame synchronizer
//
typedef struct ofdmflexframesync_s * ofdmflexframesync;
// create OFDM flexible framing synchronizer object
// _M : number of subcarriers
// _cp_len : cyclic prefix length
// _taper_len : taper length (OFDM symbol overlap)
// _p : subcarrier allocation (null, pilot, data), [size: _M x 1]
// _callback : user-defined callback function
// _userdata : user-defined data pointer
ofdmflexframesync ofdmflexframesync_create(unsigned int _M,
unsigned int _cp_len,
unsigned int _taper_len,
unsigned char * _p,
framesync_callback _callback,
void * _userdata);
void ofdmflexframesync_destroy(ofdmflexframesync _q);
void ofdmflexframesync_print(ofdmflexframesync _q);
void ofdmflexframesync_reset(ofdmflexframesync _q);
void ofdmflexframesync_execute(ofdmflexframesync _q,
liquid_float_complex * _x,
unsigned int _n);
// query the received signal strength indication
float ofdmflexframesync_get_rssi(ofdmflexframesync _q);
// query the received carrier offset estimate
float ofdmflexframesync_get_cfo(ofdmflexframesync _q);
// enable/disable debugging
void ofdmflexframesync_debug_enable(ofdmflexframesync _q);
void ofdmflexframesync_debug_disable(ofdmflexframesync _q);
void ofdmflexframesync_debug_print(ofdmflexframesync _q,
const char * _filename);
//
// Binary P/N synchronizer
//
#define BSYNC_MANGLE_RRRF(name) LIQUID_CONCAT(bsync_rrrf,name)
#define BSYNC_MANGLE_CRCF(name) LIQUID_CONCAT(bsync_crcf,name)
#define BSYNC_MANGLE_CCCF(name) LIQUID_CONCAT(bsync_cccf,name)
// Macro:
// BSYNC : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_BSYNC_DEFINE_API(BSYNC,TO,TC,TI) \
typedef struct BSYNC(_s) * BSYNC(); \
\
BSYNC() BSYNC(_create)(unsigned int _n, TC * _v); \
\
/* create binary synchronizer from m-sequence */ \
/* _g : m-sequence generator polynomial */ \
/* _k : samples/symbol (over-sampling factor) */ \
BSYNC() BSYNC(_create_msequence)(unsigned int _g, \
unsigned int _k); \
void BSYNC(_destroy)(BSYNC() _fs); \
void BSYNC(_print)(BSYNC() _fs); \
void BSYNC(_correlate)(BSYNC() _fs, TI _sym, TO * _y);
LIQUID_BSYNC_DEFINE_API(BSYNC_MANGLE_RRRF,
float,
float,
float)
LIQUID_BSYNC_DEFINE_API(BSYNC_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_BSYNC_DEFINE_API(BSYNC_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Pre-demodulation synchronizers (binary and otherwise)
//
#define PRESYNC_MANGLE_CCCF(name) LIQUID_CONCAT( presync_cccf,name)
#define BPRESYNC_MANGLE_CCCF(name) LIQUID_CONCAT(bpresync_cccf,name)
// Macro:
// PRESYNC : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_PRESYNC_DEFINE_API(PRESYNC,TO,TC,TI) \
typedef struct PRESYNC(_s) * PRESYNC(); \
\
/* create pre-demod synchronizer */ \
/* _v : baseband sequence */ \
/* _n : baseband sequence length */ \
/* _dphi_max : maximum absolute frequency deviation */ \
/* _m : number of correlators */ \
PRESYNC() PRESYNC(_create)(TC * _v, \
unsigned int _n, \
float _dphi_max, \
unsigned int _m); \
\
/* destroy pre-demod synchronizer */ \
void PRESYNC(_destroy)(PRESYNC() _q); \
\
/* print pre-demod synchronizer internal state */ \
void PRESYNC(_print)(PRESYNC() _q); \
\
/* reset pre-demod synchronizer internal state */ \
void PRESYNC(_reset)(PRESYNC() _q); \
\
/* push input sample into pre-demod synchronizer */ \
/* _q : pre-demod synchronizer object */ \
/* _x : input sample */ \
void PRESYNC(_push)(PRESYNC() _q, \
TI _x); \
\
/* correlate input sequence */ \
/* _q : pre-demod synchronizer object */ \
/* _rxy : output cross correlation */ \
/* _dphi_hat : output frequency offset estiamte */ \
void PRESYNC(_correlate)(PRESYNC() _q, \
TO * _rxy, \
float * _dphi_hat); \
// non-binary pre-demodulation synchronizer
LIQUID_PRESYNC_DEFINE_API(PRESYNC_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
// binary pre-demodulation synchronizer
LIQUID_PRESYNC_DEFINE_API(BPRESYNC_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Frame detector
//
typedef struct qdetector_cccf_s * qdetector_cccf;
// create detector with generic sequence
// _s : sample sequence
// _s_len : length of sample sequence
qdetector_cccf qdetector_cccf_create(liquid_float_complex * _s,
unsigned int _s_len);
// create detector from sequence of symbols using internal linear interpolator
// _sequence : symbol sequence
// _sequence_len : length of symbol sequence
// _ftype : filter prototype (e.g. LIQUID_FIRFILT_RRC)
// _k : samples/symbol
// _m : filter delay
// _beta : excess bandwidth factor
qdetector_cccf qdetector_cccf_create_linear(liquid_float_complex * _sequence,
unsigned int _sequence_len,
int _ftype,
unsigned int _k,
unsigned int _m,
float _beta);
// create detector from sequence of GMSK symbols
// _sequence : bit sequence
// _sequence_len : length of bit sequence
// _k : samples/symbol
// _m : filter delay
// _beta : excess bandwidth factor
qdetector_cccf qdetector_cccf_create_gmsk(unsigned char * _sequence,
unsigned int _sequence_len,
unsigned int _k,
unsigned int _m,
float _beta);
void qdetector_cccf_destroy(qdetector_cccf _q);
void qdetector_cccf_print (qdetector_cccf _q);
void qdetector_cccf_reset (qdetector_cccf _q);
// run detector, looking for sequence; return pointer to aligned, buffered samples
void * qdetector_cccf_execute(qdetector_cccf _q,
liquid_float_complex _x);
// set detection threshold (should be between 0 and 1, good starting point is 0.5)
void qdetector_cccf_set_threshold(qdetector_cccf _q,
float _threshold);
// set carrier offset search range
void qdetector_cccf_set_range(qdetector_cccf _q,
float _dphi_max);
// access methods
unsigned int qdetector_cccf_get_seq_len (qdetector_cccf _q); // sequence length
const void * qdetector_cccf_get_sequence(qdetector_cccf _q); // pointer to sequence
unsigned int qdetector_cccf_get_buf_len (qdetector_cccf _q); // buffer length
float qdetector_cccf_get_tau (qdetector_cccf _q); // fractional timing offset estimate
float qdetector_cccf_get_gamma (qdetector_cccf _q); // channel gain
float qdetector_cccf_get_dphi (qdetector_cccf _q); // carrier frequency offset estimate
float qdetector_cccf_get_phi (qdetector_cccf _q); // carrier phase offset estimate
//
// Pre-demodulation detector
//
typedef struct detector_cccf_s * detector_cccf;
// create pre-demod detector
// _s : sequence
// _n : sequence length
// _threshold : detection threshold (default: 0.7)
// _dphi_max : maximum carrier offset
detector_cccf detector_cccf_create(liquid_float_complex * _s,
unsigned int _n,
float _threshold,
float _dphi_max);
// destroy pre-demo detector object
void detector_cccf_destroy(detector_cccf _q);
// print pre-demod detector internal state
void detector_cccf_print(detector_cccf _q);
// reset pre-demod detector internal state
void detector_cccf_reset(detector_cccf _q);
// Run sample through pre-demod detector's correlator.
// Returns '1' if signal was detected, '0' otherwise
// _q : pre-demod detector
// _x : input sample
// _tau_hat : fractional sample offset estimate (set when detected)
// _dphi_hat : carrier frequency offset estimate (set when detected)
// _gamma_hat : channel gain estimate (set when detected)
int detector_cccf_correlate(detector_cccf _q,
liquid_float_complex _x,
float * _tau_hat,
float * _dphi_hat,
float * _gamma_hat);
//
// symbol streaming for testing (no meaningful data, just symbols)
//
#define SYMSTREAM_MANGLE_CFLOAT(name) LIQUID_CONCAT(symstreamcf,name)
#define LIQUID_SYMSTREAM_DEFINE_API(SYMSTREAM,TO) \
\
typedef struct SYMSTREAM(_s) * SYMSTREAM(); \
\
/* create default symstream object */ \
/* (LIQUID_RNYQUIST_ARKAISER, k=2, m=7, beta=0.3, QPSK) */ \
SYMSTREAM() SYMSTREAM(_create)(void); \
\
/* create symstream object with linear modulation */ \
/* _ftype : filter type (e.g. LIQUID_RNYQUIST_RRC) */ \
/* _k : samples per symbol */ \
/* _m : filter delay (symbols) */ \
/* _beta : filter excess bandwidth */ \
/* _ms : modulation scheme (e.g. LIQUID_MODEM_QPSK) */ \
SYMSTREAM() SYMSTREAM(_create_linear)(int _ftype, \
unsigned int _k, \
unsigned int _m, \
float _beta, \
int _ms); \
\
/* destroy symstream object, freeing all internal memory */ \
void SYMSTREAM(_destroy)(SYMSTREAM() _q); \
\
/* print symstream object's parameters */ \
void SYMSTREAM(_print)(SYMSTREAM() _q); \
\
/* reset symstream internal state */ \
void SYMSTREAM(_reset)(SYMSTREAM() _q); \
\
/* write block of samples to output buffer */ \
/* _q : synchronizer object */ \
/* _buf : output buffer [size: _buf_len x 1] */ \
/* _buf_len: output buffer size */ \
void SYMSTREAM(_write_samples)(SYMSTREAM() _q, \
TO * _buf, \
unsigned int _buf_len); \
LIQUID_SYMSTREAM_DEFINE_API(SYMSTREAM_MANGLE_CFLOAT, liquid_float_complex)
//
// multi-signal source for testing (no meaningful data, just signals)
//
#define MSOURCE_MANGLE_CFLOAT(name) LIQUID_CONCAT(msourcecf,name)
#define LIQUID_MSOURCE_DEFINE_API(MSOURCE,TO) \
\
typedef struct MSOURCE(_s) * MSOURCE(); \
\
/* create default msource object */ \
MSOURCE() MSOURCE(_create)(void); \
\
/* destroy msource object */ \
void MSOURCE(_destroy)(MSOURCE() _q); \
\
/* print msrouce object */ \
void MSOURCE(_print)(MSOURCE() _q); \
\
/* reset msrouce object */ \
void MSOURCE(_reset)(MSOURCE() _q); \
\
/* add signal sources */ \
int MSOURCE(_add_tone) (MSOURCE() _q); \
int MSOURCE(_add_noise)(MSOURCE() _q, float _bandwidth); \
int MSOURCE(_add_modem)(MSOURCE() _q, \
int _ms, \
unsigned int _k, \
unsigned int _m, \
float _beta); \
\
/* remove signal */ \
void MSOURCE(_remove)(MSOURCE() _q, int _id); \
\
/* enable/disable signal */ \
void MSOURCE(_enable) (MSOURCE() _q, int _id); \
void MSOURCE(_disable)(MSOURCE() _q, int _id); \
\
/* set signal gain */ \
/* _q : msource object */ \
/* _id : source id */ \
/* _gain : signal gain */ \
void MSOURCE(_set_gain)(MSOURCE() _q, \
int _id, \
float _gain_dB); \
\
/* set carrier offset to signal */ \
/* _q : msource object */ \
/* _id : source id */ \
/* _fc : carrier offset, fc in [-0.5,0.5] */ \
void MSOURCE(_set_frequency)(MSOURCE() _q, \
int _id, \
float _dphi); \
\
/* write block of samples to output buffer */ \
/* _q : synchronizer object */ \
/* _buf : output buffer [size: _buf_len x 1] */ \
/* _buf_len: output buffer size */ \
void MSOURCE(_write_samples)(MSOURCE() _q, \
TO * _buf, \
unsigned int _buf_len); \
LIQUID_MSOURCE_DEFINE_API(MSOURCE_MANGLE_CFLOAT, liquid_float_complex)
//
// Symbol tracking: AGC > symsync > EQ > carrier recovery
//
#define SYMTRACK_MANGLE_RRRF(name) LIQUID_CONCAT(symtrack_rrrf,name)
#define SYMTRACK_MANGLE_CCCF(name) LIQUID_CONCAT(symtrack_cccf,name)
// large macro
// SYMTRACK : name-mangling macro
// T : data type, primitive
// TO : data type, output
// TC : data type, coefficients
// TI : data type, input
#define LIQUID_SYMTRACK_DEFINE_API(SYMTRACK,T,TO,TC,TI) \
\
typedef struct SYMTRACK(_s) * SYMTRACK(); \
\
/* create symtrack object with default parameters */ \
/* _ftype : filter type (e.g. LIQUID_RNYQUIST_RRC) */ \
/* _k : samples per symbol */ \
/* _m : filter delay (symbols) */ \
/* _beta : filter excess bandwidth */ \
/* _ms : modulation scheme (e.g. LIQUID_MODEM_QPSK) */ \
SYMTRACK() SYMTRACK(_create)(int _ftype, \
unsigned int _k, \
unsigned int _m, \
float _beta, \
int _ms); \
\
/* create symtrack object using default parameters */ \
SYMTRACK() SYMTRACK(_create_default)(); \
\
/* destroy symtrack object, freeing all internal memory */ \
void SYMTRACK(_destroy)(SYMTRACK() _q); \
\
/* print symtrack object's parameters */ \
void SYMTRACK(_print)(SYMTRACK() _q); \
\
/* reset symtrack internal state */ \
void SYMTRACK(_reset)(SYMTRACK() _q); \
\
/* set symtrack modulation scheme */ \
void SYMTRACK(_set_modscheme)(SYMTRACK() _q, int _ms); \
\
/* set symtrack internal bandwidth */ \
void SYMTRACK(_set_bandwidth)(SYMTRACK() _q, float _bw); \
\
/* adjust internal nco by requested phase */ \
void SYMTRACK(_adjust_phase)(SYMTRACK() _q, T _dphi); \
\
/* execute synchronizer on single input sample */ \
/* _q : synchronizer object */ \
/* _x : input data sample */ \
/* _y : output data array */ \
/* _ny : number of samples written to output buffer */ \
void SYMTRACK(_execute)(SYMTRACK() _q, \
TI _x, \
TO * _y, \
unsigned int * _ny); \
\
/* execute synchronizer on input data array */ \
/* _q : synchronizer object */ \
/* _x : input data array */ \
/* _nx : number of input samples */ \
/* _y : output data array */ \
/* _ny : number of samples written to output buffer */ \
void SYMTRACK(_execute_block)(SYMTRACK() _q, \
TI * _x, \
unsigned int _nx, \
TO * _y, \
unsigned int * _ny); \
LIQUID_SYMTRACK_DEFINE_API(SYMTRACK_MANGLE_RRRF,
float,
float,
float,
float)
LIQUID_SYMTRACK_DEFINE_API(SYMTRACK_MANGLE_CCCF,
float,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// MODULE : math
//
// ln( Gamma(z) )
float liquid_lngammaf(float _z);
// Gamma(z)
float liquid_gammaf(float _z);
// ln( gamma(z,alpha) ) : lower incomplete gamma function
float liquid_lnlowergammaf(float _z, float _alpha);
// ln( Gamma(z,alpha) ) : upper incomplete gamma function
float liquid_lnuppergammaf(float _z, float _alpha);
// gamma(z,alpha) : lower incomplete gamma function
float liquid_lowergammaf(float _z, float _alpha);
// Gamma(z,alpha) : upper incomplete gamma function
float liquid_uppergammaf(float _z, float _alpha);
// n!
float liquid_factorialf(unsigned int _n);
// ln(I_v(z)) : log Modified Bessel function of the first kind
float liquid_lnbesselif(float _nu, float _z);
// I_v(z) : Modified Bessel function of the first kind
float liquid_besselif(float _nu, float _z);
// I_0(z) : Modified Bessel function of the first kind (order zero)
float liquid_besseli0f(float _z);
// J_v(z) : Bessel function of the first kind
float liquid_besseljf(float _nu, float _z);
// J_0(z) : Bessel function of the first kind (order zero)
float liquid_besselj0f(float _z);
// Q function
float liquid_Qf(float _z);
// Marcum Q-function
float liquid_MarcumQf(int _M,
float _alpha,
float _beta);
// Marcum Q-function (M=1)
float liquid_MarcumQ1f(float _alpha,
float _beta);
// sin(pi x) / (pi x)
float sincf(float _x);
// next power of 2 : y = ceil(log2(_x))
unsigned int liquid_nextpow2(unsigned int _x);
// (n choose k) = n! / ( k! (n-k)! )
float liquid_nchoosek(unsigned int _n, unsigned int _k);
//
// Windowing functions
//
// Modulation schemes available
#define LIQUID_WINDOW_NUM_FUNCTIONS (10)
// prototypes
typedef enum {
LIQUID_WINDOW_UNKNOWN=0, // unknown/unsupported scheme
LIQUID_WINDOW_HAMMING, // Hamming
LIQUID_WINDOW_HANN, // Hann
LIQUID_WINDOW_BLACKMANHARRIS, // Blackman-harris (4-term)
LIQUID_WINDOW_BLACKMANHARRIS7, // Blackman-harris (7-term)
LIQUID_WINDOW_KAISER, // Kaiser (beta factor unspecified)
LIQUID_WINDOW_FLATTOP, // flat top (includes negative values)
LIQUID_WINDOW_TRIANGULAR, // triangular
LIQUID_WINDOW_RCOSTAPER, // raised-cosine taper (taper size unspecified)
LIQUID_WINDOW_KBD, // Kaiser-Bessel derived window (beta factor unspecified)
} liquid_window_type;
// pretty names for window
extern const char * liquid_window_str[LIQUID_WINDOW_NUM_FUNCTIONS][2];
// Print compact list of existing and available windowing functions
void liquid_print_windows();
// returns modulation_scheme based on input string
liquid_window_type liquid_getopt_str2window(const char * _str);
// Kaiser-Bessel derived window (single sample)
// _n : index (0 <= _n < _N)
// _N : length of filter (must be even)
// _beta : Kaiser window parameter (_beta > 0)
float liquid_kbd(unsigned int _n, unsigned int _N, float _beta);
// Kaiser-Bessel derived window (full window)
// _n : length of filter (must be even)
// _beta : Kaiser window parameter (_beta > 0)
// _w : resulting window
void liquid_kbd_window(unsigned int _n, float _beta, float * _w);
// Kaiser window
// _n : window index
// _N : full window length
// _beta : Kaiser-Bessel window shape parameter
// _dt : fractional sample offset
float kaiser(unsigned int _n,
unsigned int _N,
float _beta,
float _dt);
// Hamming window
// _n : window index
// _N : full window length
float hamming(unsigned int _n,
unsigned int _N);
// Hann window
// _n : window index
// _N : full window length
float hann(unsigned int _n,
unsigned int _N);
// Blackman-harris window
// _n : window index
// _N : full window length
float blackmanharris(unsigned int _n,
unsigned int _N);
// 7th order Blackman-harris window
// _n : window index
// _N : full window length
float blackmanharris7(unsigned int _n,
unsigned int _N);
// Flat-top window
// _n : window index
// _N : full window length
float flattop(unsigned int _n,
unsigned int _N);
// Triangular window
// _n : window index
// _N : full window length
// _L : triangle length, _L in {_N, _N+1, _N-1}
float triangular(unsigned int _n,
unsigned int _N,
unsigned int _L);
// raised-cosine tapering window
// _n : window index
// _t : taper length
// _N : full window length
float liquid_rcostaper_windowf(unsigned int _n,
unsigned int _t,
unsigned int _N);
// polynomials
#define POLY_MANGLE_DOUBLE(name) LIQUID_CONCAT(poly, name)
#define POLY_MANGLE_FLOAT(name) LIQUID_CONCAT(polyf, name)
#define POLY_MANGLE_CDOUBLE(name) LIQUID_CONCAT(polyc, name)
#define POLY_MANGLE_CFLOAT(name) LIQUID_CONCAT(polycf, name)
// large macro
// POLY : name-mangling macro
// T : data type
// TC : data type (complex)
#define LIQUID_POLY_DEFINE_API(POLY,T,TC) \
/* evaluate polynomial _p (order _k-1) at value _x */ \
T POLY(_val)(T * _p, unsigned int _k, T _x); \
\
/* least-squares polynomial fit (order _k-1) */ \
void POLY(_fit)(T * _x, \
T * _y, \
unsigned int _n, \
T * _p, \
unsigned int _k); \
\
/* Lagrange polynomial exact fit (order _n-1) */ \
void POLY(_fit_lagrange)(T * _x, \
T * _y, \
unsigned int _n, \
T * _p); \
\
/* Lagrange polynomial interpolation */ \
T POLY(_interp_lagrange)(T * _x, \
T * _y, \
unsigned int _n, \
T _x0); \
\
/* Lagrange polynomial fit (barycentric form) */ \
void POLY(_fit_lagrange_barycentric)(T * _x, \
unsigned int _n, \
T * _w); \
\
/* Lagrange polynomial interpolation (barycentric form) */ \
T POLY(_val_lagrange_barycentric)(T * _x, \
T * _y, \
T * _w, \
T _x0, \
unsigned int _n); \
\
/* expands the polynomial: \
* P_n(x) = (1+x)^n \
* as \
* P_n(x) = p[0] + p[1]*x + p[2]*x^2 + ... + p[n]x^n \
* NOTE: _p has order n=m+k (array is length n+1) \
*/ \
void POLY(_expandbinomial)(unsigned int _n, \
T * _p); \
\
/* expands the polynomial: \
* P_n(x) = (1+x)^m * (1-x)^k \
* as \
* P_n(x) = p[0] + p[1]*x + p[2]*x^2 + ... + p[n]x^n \
* NOTE: _p has order n=m+k (array is length n+1) \
*/ \
void POLY(_expandbinomial_pm)(unsigned int _m, \
unsigned int _k, \
T * _p); \
\
/* expands the polynomial: \
* P_n(x) = (x-r[0]) * (x-r[1]) * ... * (x-r[n-1]) \
* as \
* P_n(x) = c[0] + c[1]*x + ... + c[n]*x^n \
* where r[0],r[1],...,r[n-1] are the roots of P_n(x) \
* NOTE: _c has order _n (array is length _n+1) \
*/ \
void POLY(_expandroots)(T * _a, \
unsigned int _n, \
T * _c); \
\
/* expands the polynomial: \
* P_n(x) = \
* (x*b[0]-a[0]) * (x*b[1]-a[1]) * ... * (x*b[n-1]-a[n-1]) \
* as \
* P_n(x) = c[0] + c[1]*x + ... + c[n]*x^n \
* NOTE: _c has order _n (array is length _n+1) \
*/ \
void POLY(_expandroots2)(T * _a, \
T * _b, \
unsigned int _n, \
T * _c); \
\
/* find roots of the polynomial (complex) */ \
/* _poly : poly array, ascending powers [size: _k x 1] */ \
/* _k : poly length (poly order = _k - 1) */ \
/* _roots : resulting complex roots [size: _k-1 x 1] */ \
void POLY(_findroots)(T * _poly, \
unsigned int _n, \
TC * _roots); \
\
/* find the complex roots of the polynomial using the */ \
/* Durand-Kerner method */ \
void POLY(_findroots_durandkerner)(T * _poly, \
unsigned int _k, \
TC * _roots); \
\
/* find the complex roots of the polynomial using */ \
/* Bairstow's method */ \
void POLY(_findroots_bairstow)(T * _poly, \
unsigned int _k, \
TC * _roots); \
\
/* expands the multiplication of two polynomials */ \
void POLY(_mul)(T * _a, \
unsigned int _order_a, \
T * _b, \
unsigned int _order_b, \
T * _c); \
LIQUID_POLY_DEFINE_API(POLY_MANGLE_DOUBLE,
double,
liquid_double_complex)
LIQUID_POLY_DEFINE_API(POLY_MANGLE_FLOAT,
float,
liquid_float_complex)
LIQUID_POLY_DEFINE_API(POLY_MANGLE_CDOUBLE,
liquid_double_complex,
liquid_double_complex)
LIQUID_POLY_DEFINE_API(POLY_MANGLE_CFLOAT,
liquid_float_complex,
liquid_float_complex)
#if 0
// expands the polynomial: (1+x)^n
void poly_binomial_expand(unsigned int _n, int * _c);
// expands the polynomial: (1+x)^k * (1-x)^(n-k)
void poly_binomial_expand_pm(unsigned int _n,
unsigned int _k,
int * _c);
#endif
//
// modular arithmetic, etc.
//
// maximum number of factors
#define LIQUID_MAX_FACTORS (40)
// is number prime?
int liquid_is_prime(unsigned int _n);
// compute number's prime factors
// _n : number to factor
// _factors : pre-allocated array of factors [size: LIQUID_MAX_FACTORS x 1]
// _num_factors: number of factors found, sorted ascending
void liquid_factor(unsigned int _n,
unsigned int * _factors,
unsigned int * _num_factors);
// compute number's unique prime factors
// _n : number to factor
// _factors : pre-allocated array of factors [size: LIQUID_MAX_FACTORS x 1]
// _num_factors: number of unique factors found, sorted ascending
void liquid_unique_factor(unsigned int _n,
unsigned int * _factors,
unsigned int * _num_factors);
// compute c = base^exp (mod n)
unsigned int liquid_modpow(unsigned int _base,
unsigned int _exp,
unsigned int _n);
// find smallest primitive root of _n
unsigned int liquid_primitive_root(unsigned int _n);
// find smallest primitive root of _n, assuming _n is prime
unsigned int liquid_primitive_root_prime(unsigned int _n);
// Euler's totient function
unsigned int liquid_totient(unsigned int _n);
//
// MODULE : matrix
//
#define MATRIX_MANGLE_DOUBLE(name) LIQUID_CONCAT(matrix, name)
#define MATRIX_MANGLE_FLOAT(name) LIQUID_CONCAT(matrixf, name)
#define MATRIX_MANGLE_CDOUBLE(name) LIQUID_CONCAT(matrixc, name)
#define MATRIX_MANGLE_CFLOAT(name) LIQUID_CONCAT(matrixcf, name)
// large macro
// MATRIX : name-mangling macro
// T : data type
#define LIQUID_MATRIX_DEFINE_API(MATRIX,T) \
\
/* print array as matrix */ \
/* _x : input matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_print)(T * _x, \
unsigned int _r, \
unsigned int _c); \
\
/* add two matrices _x and _y saving the result in _z */ \
/* _x : input matrix [size: _r x _c] */ \
/* _y : input matrix [size: _r x _c] */ \
/* _z : output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_add)(T * _x, \
T * _y, \
T * _z, \
unsigned int _r, \
unsigned int _c); \
\
/* subtract two matrices _x and _y saving the result in _z */ \
/* _x : input matrix [size: _r x _c] */ \
/* _y : input matrix [size: _r x _c] */ \
/* _z : output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_sub)(T * _x, \
T * _y, \
T * _z, \
unsigned int _r, \
unsigned int _c); \
\
/* perform point-wise multiplication of two matrices _x */ \
/* and _y saving the result in _z */ \
/* _x : input matrix [size: _r x _c] */ \
/* _y : input matrix [size: _r x _c] */ \
/* _z : output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_pmul)(T * _x, \
T * _y, \
T * _z, \
unsigned int _r, \
unsigned int _c); \
\
/* perform point-wise division of two matrices _x and _y */ \
/* saving the result in _z */ \
/* _x : input matrix [size: _r x _c] */ \
/* _y : input matrix [size: _r x _c] */ \
/* _z : output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_pdiv)(T * _x, \
T * _y, \
T * _z, \
unsigned int _r, \
unsigned int _c); \
\
/* multiply two matrices _x and _y storing the result in _z */ \
/* NOTE: _rz = _rx, _cz = _cy, and _cx = _ry */ \
/* _x : input matrix [size: _rx x _cx] */ \
/* _y : input matrix [size: _ry x _cy] */ \
/* _z : output matrix [size: _rz x _cz] */ \
void MATRIX(_mul)(T * _x, unsigned int _rx, unsigned int _cx, \
T * _y, unsigned int _ry, unsigned int _cy, \
T * _z, unsigned int _rz, unsigned int _cz); \
\
/* solve _x = _y*_z for _z for square matrices of size _n */ \
/* _x : input matrix [size: _n x _n] */ \
/* _y : input matrix [size: _n x _n] */ \
/* _z : output matrix [size: _n x _n] */ \
void MATRIX(_div)(T * _x, \
T * _y, \
T * _z, \
unsigned int _n); \
\
/* compute the determinant of a square matrix _x */ \
/* _x : input matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
T MATRIX(_det)(T * _x, \
unsigned int _r, \
unsigned int _c); \
\
/* compute the in-place transpose of the matrix _x */ \
/* _x : input matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_trans)(T * _x, \
unsigned int _r, \
unsigned int _c); \
\
/* compute the in-place Hermitian transpose of _x */ \
/* _x : input matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_hermitian)(T * _x, \
unsigned int _r, \
unsigned int _c); \
\
/* compute x*x' on [m x n] matrix, result: [m x m] */ \
/* _x : input matrix [size: _m x _n] */ \
/* _m : input rows */ \
/* _n : input columns */ \
/* _xxT : output matrix [size: _m x _m] */ \
void MATRIX(_mul_transpose)(T * _x, \
unsigned int _m, \
unsigned int _n, \
T * _xxT); \
\
/* compute x'*x on [m x n] matrix, result: [n x n] */ \
/* _x : input matrix [size: _m x _n] */ \
/* _m : input rows */ \
/* _n : input columns */ \
/* _xTx : output matrix [size: _n x _n] */ \
void MATRIX(_transpose_mul)(T * _x, \
unsigned int _m, \
unsigned int _n, \
T * _xTx); \
\
/* compute x*x.' on [m x n] matrix, result: [m x m] */ \
/* _x : input matrix [size: _m x _n] */ \
/* _m : input rows */ \
/* _n : input columns */ \
/* _xxH : output matrix [size: _m x _m] */ \
void MATRIX(_mul_hermitian)(T * _x, \
unsigned int _m, \
unsigned int _n, \
T * _xxH); \
\
/* compute x.'*x on [m x n] matrix, result: [n x n] */ \
/* _x : input matrix [size: _m x _n] */ \
/* _m : input rows */ \
/* _n : input columns */ \
/* _xHx : output matrix [size: _n x _n] */ \
void MATRIX(_hermitian_mul)(T * _x, \
unsigned int _m, \
unsigned int _n, \
T * _xHx); \
\
\
/* augment two matrices _x and _y storing the result in _z */ \
/* NOTE: _rz = _rx = _ry, _rx = _ry, and _cz = _cx + _cy */ \
/* _x : input matrix [size: _rx x _cx] */ \
/* _y : input matrix [size: _ry x _cy] */ \
/* _z : output matrix [size: _rz x _cz] */ \
void MATRIX(_aug)(T * _x, unsigned int _rx, unsigned int _cx, \
T * _y, unsigned int _ry, unsigned int _cy, \
T * _z, unsigned int _rz, unsigned int _cz); \
\
/* compute the inverse of a square matrix _x */ \
/* _x : input/output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_inv)(T * _x, \
unsigned int _r, \
unsigned int _c); \
\
/* generate the identity square matrix of size _n */ \
/* _x : output matrix [size: _n x _n] */ \
/* _n : dimensions of _x */ \
void MATRIX(_eye)(T * _x, \
unsigned int _n); \
\
/* generate the all-ones matrix of size _n */ \
/* _x : output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_ones)(T * _n, \
unsigned int _r, \
unsigned int _c); \
\
/* generate the all-zeros matrix of size _n */ \
/* _x : output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_zeros)(T * _x, \
unsigned int _r, \
unsigned int _c); \
\
/* perform Gauss-Jordan elimination on matrix _x */ \
/* _x : input/output matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
void MATRIX(_gjelim)(T * _x, \
unsigned int _r, \
unsigned int _c); \
\
/* pivot on element _x[_r,_c] */ \
/* _x : output matrix [size: _rx x _cx] */ \
/* _rx : rows of _x */ \
/* _cx : columns of _x */ \
/* _r : pivot row */ \
/* _c : pivot column */ \
void MATRIX(_pivot)(T * _x, \
unsigned int _rx, \
unsigned int _cx, \
unsigned int _r, \
unsigned int _c); \
\
/* swap rows _r1 and _r2 of matrix _x */ \
/* _x : input/output matrix [size: _rx x _cx] */ \
/* _rx : rows of _x */ \
/* _cx : columns of _x */ \
/* _r1 : first row to swap */ \
/* _r2 : second row to swap */ \
void MATRIX(_swaprows)(T * _x, \
unsigned int _rx, \
unsigned int _cx, \
unsigned int _r1, \
unsigned int _r2); \
\
/* solve linear system of _n equations: _A*_x = _b */ \
/* _A : system matrix [size: _n x _n] */ \
/* _n : system size */ \
/* _b : equality vector [size: _n x 1] */ \
/* _x : solution vector [size: _n x 1] */ \
/* _opts : options (ignored for now) */ \
void MATRIX(_linsolve)(T * _A, \
unsigned int _n, \
T * _b, \
T * _x, \
void * _opts); \
\
/* solve linear system of equations using conjugate */ \
/* gradient method */ \
/* _A : symmetric positive definite square matrix */ \
/* _n : system dimension */ \
/* _b : equality [size: _n x 1] */ \
/* _x : solution estimate [size: _n x 1] */ \
/* _opts : options (ignored for now) */ \
void MATRIX(_cgsolve)(T * _A, \
unsigned int _n, \
T * _b, \
T * _x, \
void * _opts); \
\
/* L/U/P decomposition, Crout's method */ \
/* _x : input/output matrix [size: _rx x _cx] */ \
/* _rx : rows of _x */ \
/* _cx : columns of _x */ \
/* _L : first row to swap */ \
/* _U : first row to swap */ \
/* _P : first row to swap */ \
void MATRIX(_ludecomp_crout)(T * _x, \
unsigned int _rx, \
unsigned int _cx, \
T * _L, \
T * _U, \
T * _P); \
\
/* L/U/P decomposition, Doolittle's method */ \
/* _x : input/output matrix [size: _rx x _cx] */ \
/* _rx : rows of _x */ \
/* _cx : columns of _x */ \
/* _L : first row to swap */ \
/* _U : first row to swap */ \
/* _P : first row to swap */ \
void MATRIX(_ludecomp_doolittle)(T * _x, \
unsigned int _rx, \
unsigned int _cx, \
T * _L, \
T * _U, \
T * _P); \
\
/* Orthnormalization using the Gram-Schmidt algorithm */ \
/* _A : input matrix [size: _r x _c] */ \
/* _r : rows */ \
/* _c : columns */ \
/* _v : output matrix */ \
void MATRIX(_gramschmidt)(T * _A, \
unsigned int _r, \
unsigned int _c, \
T * _v); \
\
/* Q/R decomposition using the Gram-Schmidt algorithm such */ \
/* that _A = _Q*_R and _Q^T * _Q = _In and _R is a diagonal */ \
/* matrix */ \
/* NOTE: all matrices are square */ \
/* _A : input matrix [size: _m x _m] */ \
/* _m : rows */ \
/* _n : columns (same as cols) */ \
/* _Q : output matrix [size: _m x _m] */ \
/* _R : output matrix [size: _m x _m] */ \
void MATRIX(_qrdecomp_gramschmidt)(T * _A, \
unsigned int _m, \
unsigned int _n, \
T * _Q, \
T * _R); \
\
/* Compute Cholesky decomposition of a symmetric/Hermitian */ \
/* positive-definite matrix as A = L * L^T */ \
/* _A : input square matrix [size: _n x _n] */ \
/* _n : input matrix dimension */ \
/* _L : output lower-triangular matrix */ \
void MATRIX(_chol)(T * _A, \
unsigned int _n, \
T * _L); \
#define matrix_access(X,R,C,r,c) ((X)[(r)*(C)+(c)])
#define matrixc_access(X,R,C,r,c) matrix_access(X,R,C,r,c)
#define matrixf_access(X,R,C,r,c) matrix_access(X,R,C,r,c)
#define matrixcf_access(X,R,C,r,c) matrix_access(X,R,C,r,c)
LIQUID_MATRIX_DEFINE_API(MATRIX_MANGLE_FLOAT, float)
LIQUID_MATRIX_DEFINE_API(MATRIX_MANGLE_DOUBLE, double)
LIQUID_MATRIX_DEFINE_API(MATRIX_MANGLE_CFLOAT, liquid_float_complex)
LIQUID_MATRIX_DEFINE_API(MATRIX_MANGLE_CDOUBLE, liquid_double_complex)
#define SMATRIX_MANGLE_BOOL(name) LIQUID_CONCAT(smatrixb, name)
#define SMATRIX_MANGLE_FLOAT(name) LIQUID_CONCAT(smatrixf, name)
#define SMATRIX_MANGLE_INT(name) LIQUID_CONCAT(smatrixi, name)
// sparse 'alist' matrix type (similar to MacKay, Davey Lafferty convention)
// large macro
// SMATRIX : name-mangling macro
// T : primitive data type
#define LIQUID_SMATRIX_DEFINE_API(SMATRIX,T) \
typedef struct SMATRIX(_s) * SMATRIX(); \
\
/* create _M x _N matrix, initialized with zeros */ \
SMATRIX() SMATRIX(_create)(unsigned int _M, \
unsigned int _N); \
\
/* create _M x _N matrix, initialized on array */ \
SMATRIX() SMATRIX(_create_array)(T * _x, \
unsigned int _m, \
unsigned int _n); \
\
/* destroy object */ \
void SMATRIX(_destroy)(SMATRIX() _q); \
\
/* print compact form */ \
void SMATRIX(_print)(SMATRIX() _q); \
\
/* print expanded form */ \
void SMATRIX(_print_expanded)(SMATRIX() _q); \
\
/* query properties methods */ \
void SMATRIX(_size)(SMATRIX() _q, \
unsigned int * _m, \
unsigned int * _n); \
\
/* zero all elements */ \
void SMATRIX(_clear)(SMATRIX() _q); /* zero and keep memory */ \
void SMATRIX(_reset)(SMATRIX() _q); /* zero and clear memory */ \
\
/* determine if value has been set (allocated memory) */ \
int SMATRIX(_isset)(SMATRIX() _q, \
unsigned int _m, \
unsigned int _n); \
\
/* inserts/deletes element at index (memory allocation) */ \
void SMATRIX(_insert)(SMATRIX() _q, \
unsigned int _m, \
unsigned int _n, \
T _v); \
void SMATRIX(_delete)(SMATRIX() _q, \
unsigned int _m, \
unsigned int _n); \
\
/* sets/gets the value (with memory allocation if needed) */ \
void SMATRIX(_set)(SMATRIX() _q, \
unsigned int _m, \
unsigned int _n, \
T _v); \
T SMATRIX(_get)(SMATRIX() _q, \
unsigned int _m, \
unsigned int _n); \
\
/* initialize to identity matrix */ \
void SMATRIX(_eye)(SMATRIX() _q); \
\
/* multiply two sparse binary matrices */ \
void SMATRIX(_mul)(SMATRIX() _x, \
SMATRIX() _y, \
SMATRIX() _z); \
\
/* multiply sparse matrix by vector */ \
/* _q : sparse matrix */ \
/* _x : input vector [size: _N x 1] */ \
/* _y : output vector [size: _M x 1] */ \
void SMATRIX(_vmul)(SMATRIX() _q, \
T * _x, \
T * _y); \
LIQUID_SMATRIX_DEFINE_API(SMATRIX_MANGLE_BOOL, unsigned char)
LIQUID_SMATRIX_DEFINE_API(SMATRIX_MANGLE_FLOAT, float)
LIQUID_SMATRIX_DEFINE_API(SMATRIX_MANGLE_INT, short int)
//
// smatrix cross methods
//
// multiply sparse binary matrix by floating-point matrix
// _q : sparse matrix [size: A->M x A->N]
// _x : input vector [size: mx x nx ]
// _y : output vector [size: my x ny ]
void smatrixb_mulf(smatrixb _A,
float * _x,
unsigned int _mx,
unsigned int _nx,
float * _y,
unsigned int _my,
unsigned int _ny);
// multiply sparse binary matrix by floating-point vector
// _q : sparse matrix
// _x : input vector [size: _N x 1]
// _y : output vector [size: _M x 1]
void smatrixb_vmulf(smatrixb _q,
float * _x,
float * _y);
//
// MODULE : modem (modulator/demodulator)
//
// Maximum number of allowed bits per symbol
#define MAX_MOD_BITS_PER_SYMBOL 8
// Modulation schemes available
#define LIQUID_MODEM_NUM_SCHEMES (52)
typedef enum {
LIQUID_MODEM_UNKNOWN=0, // Unknown modulation scheme
// Phase-shift keying (PSK)
LIQUID_MODEM_PSK2, LIQUID_MODEM_PSK4,
LIQUID_MODEM_PSK8, LIQUID_MODEM_PSK16,
LIQUID_MODEM_PSK32, LIQUID_MODEM_PSK64,
LIQUID_MODEM_PSK128, LIQUID_MODEM_PSK256,
// Differential phase-shift keying (DPSK)
LIQUID_MODEM_DPSK2, LIQUID_MODEM_DPSK4,
LIQUID_MODEM_DPSK8, LIQUID_MODEM_DPSK16,
LIQUID_MODEM_DPSK32, LIQUID_MODEM_DPSK64,
LIQUID_MODEM_DPSK128, LIQUID_MODEM_DPSK256,
// amplitude-shift keying
LIQUID_MODEM_ASK2, LIQUID_MODEM_ASK4,
LIQUID_MODEM_ASK8, LIQUID_MODEM_ASK16,
LIQUID_MODEM_ASK32, LIQUID_MODEM_ASK64,
LIQUID_MODEM_ASK128, LIQUID_MODEM_ASK256,
// rectangular quadrature amplitude-shift keying (QAM)
LIQUID_MODEM_QAM4,
LIQUID_MODEM_QAM8, LIQUID_MODEM_QAM16,
LIQUID_MODEM_QAM32, LIQUID_MODEM_QAM64,
LIQUID_MODEM_QAM128, LIQUID_MODEM_QAM256,
// amplitude phase-shift keying (APSK)
LIQUID_MODEM_APSK4,
LIQUID_MODEM_APSK8, LIQUID_MODEM_APSK16,
LIQUID_MODEM_APSK32, LIQUID_MODEM_APSK64,
LIQUID_MODEM_APSK128, LIQUID_MODEM_APSK256,
// specific modem types
LIQUID_MODEM_BPSK, // Specific: binary PSK
LIQUID_MODEM_QPSK, // specific: quaternary PSK
LIQUID_MODEM_OOK, // Specific: on/off keying
LIQUID_MODEM_SQAM32, // 'square' 32-QAM
LIQUID_MODEM_SQAM128, // 'square' 128-QAM
LIQUID_MODEM_V29, // V.29 star constellation
LIQUID_MODEM_ARB16OPT, // optimal 16-QAM
LIQUID_MODEM_ARB32OPT, // optimal 32-QAM
LIQUID_MODEM_ARB64OPT, // optimal 64-QAM
LIQUID_MODEM_ARB128OPT, // optimal 128-QAM
LIQUID_MODEM_ARB256OPT, // optimal 256-QAM
LIQUID_MODEM_ARB64VT, // Virginia Tech logo
// arbitrary modem type
LIQUID_MODEM_ARB // arbitrary QAM
} modulation_scheme;
// structure for holding full modulation type descriptor
struct modulation_type_s {
const char * name; // short name (e.g. 'bpsk')
const char * fullname; // full name (e.g. 'binary phase-shift keying')
modulation_scheme scheme; // modulation scheme (e.g. LIQUID_MODEM_BPSK)
unsigned int bps; // modulation depth (e.g. 1)
};
// full modulation type descriptor
extern const struct modulation_type_s modulation_types[LIQUID_MODEM_NUM_SCHEMES];
// Print compact list of existing and available modulation schemes
void liquid_print_modulation_schemes();
// returns modulation_scheme based on input string
modulation_scheme liquid_getopt_str2mod(const char * _str);
// query basic modulation types
int liquid_modem_is_psk(modulation_scheme _ms);
int liquid_modem_is_dpsk(modulation_scheme _ms);
int liquid_modem_is_ask(modulation_scheme _ms);
int liquid_modem_is_qam(modulation_scheme _ms);
int liquid_modem_is_apsk(modulation_scheme _ms);
// useful functions
// counts the number of different bits between two symbols
unsigned int count_bit_errors(unsigned int _s1, unsigned int _s2);
// counts the number of different bits between two arrays of symbols
// _msg0 : original message [size: _n x 1]
// _msg1 : copy of original message [size: _n x 1]
// _n : message size
unsigned int count_bit_errors_array(unsigned char * _msg0,
unsigned char * _msg1,
unsigned int _n);
// converts binary-coded decimal (BCD) to gray, ensuring successive values
// differ by exactly one bit
unsigned int gray_encode(unsigned int symbol_in);
// converts a gray-encoded symbol to binary-coded decimal (BCD)
unsigned int gray_decode(unsigned int symbol_in);
// pack soft bits into symbol
// _soft_bits : soft input bits [size: _bps x 1]
// _bps : bits per symbol
// _sym_out : output symbol, value in [0,2^_bps)
void liquid_pack_soft_bits(unsigned char * _soft_bits,
unsigned int _bps,
unsigned int * _sym_out);
// unpack soft bits into symbol
// _sym_in : input symbol, value in [0,2^_bps)
// _bps : bits per symbol
// _soft_bits : soft output bits [size: _bps x 1]
void liquid_unpack_soft_bits(unsigned int _sym_in,
unsigned int _bps,
unsigned char * _soft_bits);
//
// Linear modem
//
#define LIQUID_MODEM_MANGLE_FLOAT(name) LIQUID_CONCAT(modem,name)
// Macro : MODEM
// MODEM : name-mangling macro
// T : primitive data type
// TC : primitive data type (complex)
#define LIQUID_MODEM_DEFINE_API(MODEM,T,TC) \
\
/* define struct pointer */ \
typedef struct MODEM(_s) * MODEM(); \
\
/* create digital modem object */ \
MODEM() MODEM(_create)(modulation_scheme _scheme); \
\
/* create arbitrary digital modem object */ \
/* _table : array of complex constellation points */ \
/* _M : modulation order and table size */ \
MODEM() MODEM(_create_arbitrary)(TC * _table, \
unsigned int _M); \
\
/* recreate modulation scheme, re-allocating memory as */ \
/* necessary */ \
MODEM() MODEM(_recreate)(MODEM() _q, \
modulation_scheme _scheme); \
\
void MODEM(_destroy)(MODEM() _q); \
void MODEM(_print)( MODEM() _q); \
void MODEM(_reset)( MODEM() _q); \
\
/* generate random symbol */ \
unsigned int MODEM(_gen_rand_sym)(MODEM() _q); \
\
/* Accessor functions */ \
unsigned int MODEM(_get_bps) (MODEM() _q); \
modulation_scheme MODEM(_get_scheme)(MODEM() _q); \
\
/* generic modulate function; simply queries modem scheme */ \
/* and calls appropriate subroutine */ \
/* _q : modem object */ \
/* _s : input symbol */ \
/* _x : output sample */ \
void MODEM(_modulate)(MODEM() _q, \
unsigned int _s, \
TC * _y); \
\
/* generic hard-decision demodulation function */ \
/* _q : modem object */ \
/* _x : input sample */ \
/* _s : output symbol */ \
void MODEM(_demodulate)(MODEM() _q, \
TC _x, \
unsigned int * _s); \
\
/* generic soft-decision demodulation function */ \
/* _q : modem object */ \
/* _x : input sample */ \
/* _s : output hard symbol */ \
/* _soft_bits : output soft bits */ \
void MODEM(_demodulate_soft)(MODEM() _q, \
TC _x, \
unsigned int * _s, \
unsigned char * _soft_bits); \
\
/* get demodulator's estimated transmit sample */ \
void MODEM(_get_demodulator_sample)(MODEM() _q, \
TC * _x_hat); \
\
/* get demodulator phase error */ \
float MODEM(_get_demodulator_phase_error)(MODEM() _q); \
\
/* get demodulator error vector magnitude */ \
float MODEM(_get_demodulator_evm)(MODEM() _q); \
// define modem APIs
LIQUID_MODEM_DEFINE_API(LIQUID_MODEM_MANGLE_FLOAT,float,liquid_float_complex)
//
// continuous-phase modulation
//
// gmskmod : GMSK modulator
typedef struct gmskmod_s * gmskmod;
// create gmskmod object
// _k : samples/symbol
// _m : filter delay (symbols)
// _BT : excess bandwidth factor
gmskmod gmskmod_create(unsigned int _k,
unsigned int _m,
float _BT);
void gmskmod_destroy(gmskmod _q);
void gmskmod_print(gmskmod _q);
void gmskmod_reset(gmskmod _q);
void gmskmod_modulate(gmskmod _q,
unsigned int _sym,
liquid_float_complex * _y);
// gmskdem : GMSK demodulator
typedef struct gmskdem_s * gmskdem;
// create gmskdem object
// _k : samples/symbol
// _m : filter delay (symbols)
// _BT : excess bandwidth factor
gmskdem gmskdem_create(unsigned int _k,
unsigned int _m,
float _BT);
void gmskdem_destroy(gmskdem _q);
void gmskdem_print(gmskdem _q);
void gmskdem_reset(gmskdem _q);
void gmskdem_set_eq_bw(gmskdem _q, float _bw);
void gmskdem_demodulate(gmskdem _q,
liquid_float_complex * _y,
unsigned int * _sym);
//
// continuous phase frequency-shift keying (CP-FSK) modems
//
// CP-FSK filter prototypes
typedef enum {
LIQUID_CPFSK_SQUARE=0, // square pulse
LIQUID_CPFSK_RCOS_FULL, // raised-cosine (full response)
LIQUID_CPFSK_RCOS_PARTIAL, // raised-cosine (partial response)
LIQUID_CPFSK_GMSK, // Gauss minimum-shift keying pulse
} liquid_cpfsk_filter;
// CP-FSK modulator
typedef struct cpfskmod_s * cpfskmod;
// create cpfskmod object (frequency modulator)
// _bps : bits per symbol, _bps > 0
// _h : modulation index, _h > 0
// _k : samples/symbol, _k > 1, _k even
// _m : filter delay (symbols), _m > 0
// _beta : filter bandwidth parameter, _beta > 0
// _type : filter type (e.g. LIQUID_CPFSK_SQUARE)
cpfskmod cpfskmod_create(unsigned int _bps,
float _h,
unsigned int _k,
unsigned int _m,
float _beta,
int _type);
//cpfskmod cpfskmod_create_msk(unsigned int _k);
//cpfskmod cpfskmod_create_gmsk(unsigned int _k, float _BT);
// destroy cpfskmod object
void cpfskmod_destroy(cpfskmod _q);
// print cpfskmod object internals
void cpfskmod_print(cpfskmod _q);
// reset state
void cpfskmod_reset(cpfskmod _q);
// get transmit delay [symbols]
unsigned int cpfskmod_get_delay(cpfskmod _q);
// modulate sample
// _q : frequency modulator object
// _s : input symbol
// _y : output sample array [size: _k x 1]
void cpfskmod_modulate(cpfskmod _q,
unsigned int _s,
liquid_float_complex * _y);
// CP-FSK demodulator
typedef struct cpfskdem_s * cpfskdem;
// create cpfskdem object (frequency modulator)
// _bps : bits per symbol, _bps > 0
// _h : modulation index, _h > 0
// _k : samples/symbol, _k > 1, _k even
// _m : filter delay (symbols), _m > 0
// _beta : filter bandwidth parameter, _beta > 0
// _type : filter type (e.g. LIQUID_CPFSK_SQUARE)
cpfskdem cpfskdem_create(unsigned int _bps,
float _h,
unsigned int _k,
unsigned int _m,
float _beta,
int _type);
//cpfskdem cpfskdem_create_msk(unsigned int _k);
//cpfskdem cpfskdem_create_gmsk(unsigned int _k, float _BT);
// destroy cpfskdem object
void cpfskdem_destroy(cpfskdem _q);
// print cpfskdem object internals
void cpfskdem_print(cpfskdem _q);
// reset state
void cpfskdem_reset(cpfskdem _q);
// get receive delay [symbols]
unsigned int cpfskdem_get_delay(cpfskdem _q);
#if 0
// demodulate array of samples
// _q : continuous-phase frequency demodulator object
// _y : input sample array [size: _n x 1]
// _n : input sample array length
// _s : output symbol array
// _nw : number of output symbols written
void cpfskdem_demodulate(cpfskdem _q,
liquid_float_complex * _y,
unsigned int _n,
unsigned int * _s,
unsigned int * _nw);
#else
// demodulate array of samples, assuming perfect timing
// _q : continuous-phase frequency demodulator object
// _y : input sample array [size: _k x 1]
unsigned int cpfskdem_demodulate(cpfskdem _q,
liquid_float_complex * _y);
#endif
//
// M-ary frequency-shift keying (MFSK) modems
//
// FSK modulator
typedef struct fskmod_s * fskmod;
// create fskmod object (frequency modulator)
// _m : bits per symbol, _bps > 0
// _k : samples/symbol, _k >= 2^_m
// _bandwidth : total signal bandwidth, (0,0.5)
fskmod fskmod_create(unsigned int _m,
unsigned int _k,
float _bandwidth);
// destroy fskmod object
void fskmod_destroy(fskmod _q);
// print fskmod object internals
void fskmod_print(fskmod _q);
// reset state
void fskmod_reset(fskmod _q);
// modulate sample
// _q : frequency modulator object
// _s : input symbol
// _y : output sample array [size: _k x 1]
void fskmod_modulate(fskmod _q,
unsigned int _s,
liquid_float_complex * _y);
// FSK demodulator
typedef struct fskdem_s * fskdem;
// create fskdem object (frequency demodulator)
// _m : bits per symbol, _bps > 0
// _k : samples/symbol, _k >= 2^_m
// _bandwidth : total signal bandwidth, (0,0.5)
fskdem fskdem_create(unsigned int _m,
unsigned int _k,
float _bandwidth);
// destroy fskdem object
void fskdem_destroy(fskdem _q);
// print fskdem object internals
void fskdem_print(fskdem _q);
// reset state
void fskdem_reset(fskdem _q);
// demodulate symbol, assuming perfect symbol timing
// _q : fskdem object
// _y : input sample array [size: _k x 1]
unsigned int fskdem_demodulate(fskdem _q,
liquid_float_complex * _y);
// get demodulator frequency error
float fskdem_get_frequency_error(fskdem _q);
//
// Analog frequency modulator
//
#define LIQUID_FREQMOD_MANGLE_FLOAT(name) LIQUID_CONCAT(freqmod,name)
// Macro : FREQMOD (analog frequency modulator)
// FREQMOD : name-mangling macro
// T : primitive data type
// TC : primitive data type (complex)
#define LIQUID_FREQMOD_DEFINE_API(FREQMOD,T,TC) \
\
/* define struct pointer */ \
typedef struct FREQMOD(_s) * FREQMOD(); \
\
/* create freqmod object (frequency modulator) */ \
/* _kf : modulation factor */ \
FREQMOD() FREQMOD(_create)(float _kf); \
\
/* destroy freqmod object */ \
void FREQMOD(_destroy)(FREQMOD() _q); \
\
/* print freqmod object internals */ \
void FREQMOD(_print)(FREQMOD() _q); \
\
/* reset state */ \
void FREQMOD(_reset)(FREQMOD() _q); \
\
/* modulate single sample */ \
/* _q : frequency modulator object */ \
/* _m : message signal m(t) */ \
/* _s : complex baseband signal s(t) */ \
void FREQMOD(_modulate)(FREQMOD() _q, \
T _m, \
TC * _s); \
\
/* modulate block of samples */ \
/* _q : frequency modulator object */ \
/* _m : message signal m(t), [size: _n x 1] */ \
/* _n : number of input, output samples */ \
/* _s : complex baseband signal s(t) [size: _n x 1] */ \
void FREQMOD(_modulate_block)(FREQMOD() _q, \
T * _m, \
unsigned int _n, \
TC * _s); \
// define freqmod APIs
LIQUID_FREQMOD_DEFINE_API(LIQUID_FREQMOD_MANGLE_FLOAT,float,liquid_float_complex)
//
// Analog frequency demodulator
//
#define LIQUID_FREQDEM_MANGLE_FLOAT(name) LIQUID_CONCAT(freqdem,name)
// Macro : FREQDEM (analog frequency modulator)
// FREQDEM : name-mangling macro
// T : primitive data type
// TC : primitive data type (complex)
#define LIQUID_FREQDEM_DEFINE_API(FREQDEM,T,TC) \
\
/* define struct pointer */ \
typedef struct FREQDEM(_s) * FREQDEM(); \
\
/* create freqdem object (frequency modulator) */ \
/* _kf : modulation factor */ \
FREQDEM() FREQDEM(_create)(float _kf); \
\
/* destroy freqdem object */ \
void FREQDEM(_destroy)(FREQDEM() _q); \
\
/* print freqdem object internals */ \
void FREQDEM(_print)(FREQDEM() _q); \
\
/* reset state */ \
void FREQDEM(_reset)(FREQDEM() _q); \
\
/* demodulate sample */ \
/* _q : frequency modulator object */ \
/* _r : received signal r(t) */ \
/* _m : output message signal m(t) */ \
void FREQDEM(_demodulate)(FREQDEM() _q, \
TC _r, \
T * _m); \
\
/* demodulate block of samples */ \
/* _q : frequency demodulator object */ \
/* _r : received signal r(t) [size: _n x 1] */ \
/* _n : number of input, output samples */ \
/* _m : message signal m(t), [size: _n x 1] */ \
void FREQDEM(_demodulate_block)(FREQDEM() _q, \
TC * _r, \
unsigned int _n, \
T * _m); \
// define freqdem APIs
LIQUID_FREQDEM_DEFINE_API(LIQUID_FREQDEM_MANGLE_FLOAT,float,liquid_float_complex)
// amplitude modulation types
typedef enum {
LIQUID_AMPMODEM_DSB=0, // double side-band
LIQUID_AMPMODEM_USB, // single side-band (upper)
LIQUID_AMPMODEM_LSB // single side-band (lower)
} liquid_ampmodem_type;
typedef struct ampmodem_s * ampmodem;
// create ampmodem object
// _m : modulation index
// _fc : carrier frequency, range: [-0.5,0.5]
// _type : AM type (e.g. LIQUID_AMPMODEM_DSB)
// _suppressed_carrier : carrier suppression flag
ampmodem ampmodem_create(float _m,
float _fc,
liquid_ampmodem_type _type,
int _suppressed_carrier);
// destroy ampmodem object
void ampmodem_destroy(ampmodem _fm);
// print ampmodem object internals
void ampmodem_print(ampmodem _fm);
// reset ampmodem object state
void ampmodem_reset(ampmodem _fm);
// modulate sample
void ampmodem_modulate(ampmodem _fm,
float _x,
liquid_float_complex *_y);
void ampmodem_modulate_block(ampmodem _q,
float * _m,
unsigned int _n,
liquid_float_complex *_s);
// demodulate sample
void ampmodem_demodulate(ampmodem _fm,
liquid_float_complex _y,
float *_x);
void ampmodem_demodulate_block(ampmodem _q,
liquid_float_complex * _r,
unsigned int _n,
float * _m);
//
// MODULE : multichannel
//
#define FIRPFBCH_NYQUIST 0
#define FIRPFBCH_ROOTNYQUIST 1
#define LIQUID_ANALYZER 0
#define LIQUID_SYNTHESIZER 1
//
// Finite impulse response polyphase filterbank channelizer
//
#define FIRPFBCH_MANGLE_CRCF(name) LIQUID_CONCAT(firpfbch_crcf,name)
#define FIRPFBCH_MANGLE_CCCF(name) LIQUID_CONCAT(firpfbch_cccf,name)
// Macro:
// FIRPFBCH : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_FIRPFBCH_DEFINE_API(FIRPFBCH,TO,TC,TI) \
typedef struct FIRPFBCH(_s) * FIRPFBCH(); \
\
/* create finite impulse response polyphase filter-bank */ \
/* channelizer object from external coefficients */ \
/* _type : channelizer type, e.g. LIQUID_ANALYZER */ \
/* _M : number of channels */ \
/* _p : number of coefficients for each channel */ \
/* _h : coefficients [size: _M*_p x 1] */ \
FIRPFBCH() FIRPFBCH(_create)(int _type, \
unsigned int _M, \
unsigned int _p, \
TC * _h); \
\
/* create FIR polyphase filterbank channelizer object with */ \
/* prototype filter based on windowed Kaiser design */ \
/* _type : type (LIQUID_ANALYZER | LIQUID_SYNTHESIZER) */ \
/* _M : number of channels */ \
/* _m : filter delay (symbols) */ \
/* _As : stop-band attentuation [dB] */ \
FIRPFBCH() FIRPFBCH(_create_kaiser)(int _type, \
unsigned int _M, \
unsigned int _m, \
float _As); \
\
/* create FIR polyphase filterbank channelizer object with */ \
/* prototype root-Nyquist filter */ \
/* _type : type (LIQUID_ANALYZER | LIQUID_SYNTHESIZER) */ \
/* _M : number of channels */ \
/* _m : filter delay (symbols) */ \
/* _beta : filter excess bandwidth factor, in [0,1] */ \
/* _ftype : filter prototype (rrcos, rkaiser, etc.) */ \
FIRPFBCH() FIRPFBCH(_create_rnyquist)(int _type, \
unsigned int _M, \
unsigned int _m, \
float _beta, \
int _ftype); \
\
/* destroy firpfbch object */ \
void FIRPFBCH(_destroy)(FIRPFBCH() _q); \
\
/* clear/reset firpfbch internal state */ \
void FIRPFBCH(_reset)(FIRPFBCH() _q); \
\
/* print firpfbch internal parameters to stdout */ \
void FIRPFBCH(_print)(FIRPFBCH() _q); \
\
/* execute filterbank as synthesizer on block of samples */ \
/* _q : filterbank channelizer object */ \
/* _x : channelized input, [size: num_channels x 1] */ \
/* _y : output time series, [size: num_channels x 1] */ \
void FIRPFBCH(_synthesizer_execute)(FIRPFBCH() _q, \
TI * _x, \
TO * _y); \
\
/* execute filterbank as analyzer on block of samples */ \
/* _q : filterbank channelizer object */ \
/* _x : input time series, [size: num_channels x 1] */ \
/* _y : channelized output, [size: num_channels x 1] */ \
void FIRPFBCH(_analyzer_execute)(FIRPFBCH() _q, \
TI * _x, \
TO * _y); \
LIQUID_FIRPFBCH_DEFINE_API(FIRPFBCH_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
LIQUID_FIRPFBCH_DEFINE_API(FIRPFBCH_MANGLE_CCCF,
liquid_float_complex,
liquid_float_complex,
liquid_float_complex)
//
// Finite impulse response polyphase filterbank channelizer
// with output rate 2 Fs / M
//
#define FIRPFBCH2_MANGLE_CRCF(name) LIQUID_CONCAT(firpfbch2_crcf,name)
// Macro:
// FIRPFBCH2 : name-mangling macro
// TO : output data type
// TC : coefficients data type
// TI : input data type
#define LIQUID_FIRPFBCH2_DEFINE_API(FIRPFBCH2,TO,TC,TI) \
typedef struct FIRPFBCH2(_s) * FIRPFBCH2(); \
\
/* create firpfbch2 object */ \
/* _type : channelizer type (e.g. LIQUID_ANALYZER) */ \
/* _M : number of channels (must be even) */ \
/* _m : prototype filter semi-lenth, length=2*M*m */ \
/* _h : prototype filter coefficient array */ \
FIRPFBCH2() FIRPFBCH2(_create)(int _type, \
unsigned int _M, \
unsigned int _m, \
TC * _h); \
\
/* create firpfbch2 object using Kaiser window prototype */ \
/* _type : channelizer type (e.g. LIQUID_ANALYZER) */ \
/* _M : number of channels (must be even) */ \
/* _m : prototype filter semi-lenth, length=2*M*m+1 */ \
/* _As : filter stop-band attenuation [dB] */ \
FIRPFBCH2() FIRPFBCH2(_create_kaiser)(int _type, \
unsigned int _M, \
unsigned int _m, \
float _As); \
\
/* destroy firpfbch2 object, freeing internal memory */ \
void FIRPFBCH2(_destroy)(FIRPFBCH2() _q); \
\
/* reset firpfbch2 object internals */ \
void FIRPFBCH2(_reset)(FIRPFBCH2() _q); \
\
/* print firpfbch2 object internals */ \
void FIRPFBCH2(_print)(FIRPFBCH2() _q); \
\
/* execute filterbank channelizer */ \
/* LIQUID_ANALYZER: input: M/2, output: M */ \
/* LIQUID_SYNTHESIZER: input: M, output: M/2 */ \
/* _x : channelizer input */ \
/* _y : channelizer output */ \
void FIRPFBCH2(_execute)(FIRPFBCH2() _q, \
TI * _x, \
TO * _y); \
LIQUID_FIRPFBCH2_DEFINE_API(FIRPFBCH2_MANGLE_CRCF,
liquid_float_complex,
float,
liquid_float_complex)
#define OFDMFRAME_SCTYPE_NULL 0
#define OFDMFRAME_SCTYPE_PILOT 1
#define OFDMFRAME_SCTYPE_DATA 2
// initialize default subcarrier allocation
// _M : number of subcarriers
// _p : output subcarrier allocation array, [size: _M x 1]
void ofdmframe_init_default_sctype(unsigned int _M,
unsigned char * _p);
// initialize default subcarrier allocation
// _M : number of subcarriers
// _f0 : lower frequency band, _f0 in [-0.5,0.5]
// _f1 : upper frequency band, _f1 in [-0.5,0.5]
// _p : output subcarrier allocation array, [size: _M x 1]
void ofdmframe_init_sctype_range(unsigned int _M,
float _f0,
float _f1,
unsigned char * _p);
// validate subcarrier type (count number of null, pilot, and data
// subcarriers in the allocation)
// _p : subcarrier allocation array, [size: _M x 1]
// _M : number of subcarriers
// _M_null : output number of null subcarriers
// _M_pilot : output number of pilot subcarriers
// _M_data : output number of data subcarriers
void ofdmframe_validate_sctype(unsigned char * _p,
unsigned int _M,
unsigned int * _M_null,
unsigned int * _M_pilot,
unsigned int * _M_data);
// print subcarrier allocation to screen
// _p : output subcarrier allocation array, [size: _M x 1]
// _M : number of subcarriers
void ofdmframe_print_sctype(unsigned char * _p,
unsigned int _M);
//
// OFDM frame (symbol) generator
//
typedef struct ofdmframegen_s * ofdmframegen;
// create OFDM framing generator object
// _M : number of subcarriers, >10 typical
// _cp_len : cyclic prefix length
// _taper_len : taper length (OFDM symbol overlap)
// _p : subcarrier allocation (null, pilot, data), [size: _M x 1]
ofdmframegen ofdmframegen_create(unsigned int _M,
unsigned int _cp_len,
unsigned int _taper_len,
unsigned char * _p);
void ofdmframegen_destroy(ofdmframegen _q);
void ofdmframegen_print(ofdmframegen _q);
void ofdmframegen_reset(ofdmframegen _q);
// write first S0 symbol
void ofdmframegen_write_S0a(ofdmframegen _q,
liquid_float_complex *_y);
// write second S0 symbol
void ofdmframegen_write_S0b(ofdmframegen _q,
liquid_float_complex *_y);
// write S1 symbol
void ofdmframegen_write_S1(ofdmframegen _q,
liquid_float_complex *_y);
// write data symbol
void ofdmframegen_writesymbol(ofdmframegen _q,
liquid_float_complex * _x,
liquid_float_complex *_y);
// write tail
void ofdmframegen_writetail(ofdmframegen _q,
liquid_float_complex * _x);
//
// OFDM frame (symbol) synchronizer
//
typedef int (*ofdmframesync_callback)(liquid_float_complex * _y,
unsigned char * _p,
unsigned int _M,
void * _userdata);
typedef struct ofdmframesync_s * ofdmframesync;
// create OFDM framing synchronizer object
// _M : number of subcarriers, >10 typical
// _cp_len : cyclic prefix length
// _taper_len : taper length (OFDM symbol overlap)
// _p : subcarrier allocation (null, pilot, data), [size: _M x 1]
// _callback : user-defined callback function
// _userdata : user-defined data pointer
ofdmframesync ofdmframesync_create(unsigned int _M,
unsigned int _cp_len,
unsigned int _taper_len,
unsigned char * _p,
ofdmframesync_callback _callback,
void * _userdata);
void ofdmframesync_destroy(ofdmframesync _q);
void ofdmframesync_print(ofdmframesync _q);
void ofdmframesync_reset(ofdmframesync _q);
void ofdmframesync_execute(ofdmframesync _q,
liquid_float_complex * _x,
unsigned int _n);
// query methods
float ofdmframesync_get_rssi(ofdmframesync _q); // received signal strength indication
float ofdmframesync_get_cfo(ofdmframesync _q); // carrier offset estimate
// debugging
void ofdmframesync_debug_enable(ofdmframesync _q);
void ofdmframesync_debug_disable(ofdmframesync _q);
void ofdmframesync_debug_print(ofdmframesync _q, const char * _filename);
//
// MODULE : nco (numerically-controlled oscillator)
//
// oscillator type
// LIQUID_NCO : numerically-controlled oscillator (fast)
// LIQUID_VCO : "voltage"-controlled oscillator (precise)
typedef enum {
LIQUID_NCO=0,
LIQUID_VCO
} liquid_ncotype;
#define NCO_MANGLE_FLOAT(name) LIQUID_CONCAT(nco_crcf, name)
// large macro
// NCO : name-mangling macro
// T : primitive data type
// TC : input/output data type
#define LIQUID_NCO_DEFINE_API(NCO,T,TC) \
typedef struct NCO(_s) * NCO(); \
\
NCO() NCO(_create)(liquid_ncotype _type); \
void NCO(_destroy)(NCO() _q); \
void NCO(_print)(NCO() _q); \
\
/* set phase/frequency to zero, reset pll filter */ \
void NCO(_reset)(NCO() _q); \
\
/* get/set/adjust internal frequency/phase */ \
T NCO(_get_frequency)( NCO() _q); \
void NCO(_set_frequency)( NCO() _q, T _f); \
void NCO(_adjust_frequency)(NCO() _q, T _df); \
T NCO(_get_phase)( NCO() _q); \
void NCO(_set_phase)( NCO() _q, T _phi); \
void NCO(_adjust_phase)( NCO() _q, T _dphi); \
\
/* increment phase by internal phase step (frequency) */ \
void NCO(_step)(NCO() _q); \
\
/* compute trigonometric functions */ \
T NCO(_sin)(NCO() _q); \
T NCO(_cos)(NCO() _q); \
void NCO(_sincos)(NCO() _q, T* _s, T* _c); \
void NCO(_cexpf)(NCO() _q, TC * _y); \
\
/* pll : phase-locked loop */ \
void NCO(_pll_set_bandwidth)(NCO() _q, T _bandwidth); \
void NCO(_pll_step)(NCO() _q, T _dphi); \
\
/* Rotate input sample up by NCO angle (no stepping) */ \
void NCO(_mix_up)(NCO() _q, TC _x, TC *_y); \
\
/* Rotate input sample down by NCO angle (no stepping) */ \
void NCO(_mix_down)(NCO() _q, TC _x, TC *_y); \
\
/* Rotate input vector up by NCO angle (stepping) */ \
/* _q : nco object */ \
/* _x : input vector [size: _N x 1] */ \
/* _y : output vector [size: _N x 1] */ \
/* _N : vector size */ \
void NCO(_mix_block_up)(NCO() _q, \
TC *_x, \
TC *_y, \
unsigned int _N); \
\
/* Rotate input vector down by NCO angle (stepping) */ \
/* _q : nco object */ \
/* _x : input vector [size: _N x 1] */ \
/* _y : output vector [size: _N x 1] */ \
/* _N : vector size */ \
void NCO(_mix_block_down)(NCO() _q, \
TC *_x, \
TC *_y, \
unsigned int _N); \
// Define nco APIs
LIQUID_NCO_DEFINE_API(NCO_MANGLE_FLOAT, float, liquid_float_complex)
// nco utilities
// unwrap phase of array (basic)
void liquid_unwrap_phase(float * _theta, unsigned int _n);
// unwrap phase of array (advanced)
void liquid_unwrap_phase2(float * _theta, unsigned int _n);
//
// MODULE : optimization
//
// utility function pointer definition
typedef float (*utility_function)(void * _userdata,
float * _v,
unsigned int _n);
// n-dimensional Rosenbrock utility function (minimum at _v = {1,1,1...}
// _userdata : user-defined data structure (convenience)
// _v : input vector [size: _n x 1]
// _n : input vector size
float liquid_rosenbrock(void * _userdata,
float * _v,
unsigned int _n);
// n-dimensional inverse Gauss utility function (minimum at _v = {0,0,0...}
// _userdata : user-defined data structure (convenience)
// _v : input vector [size: _n x 1]
// _n : input vector size
float liquid_invgauss(void * _userdata,
float * _v,
unsigned int _n);
// n-dimensional multimodal utility function (minimum at _v = {0,0,0...}
// _userdata : user-defined data structure (convenience)
// _v : input vector [size: _n x 1]
// _n : input vector size
float liquid_multimodal(void * _userdata,
float * _v,
unsigned int _n);
// n-dimensional spiral utility function (minimum at _v = {0,0,0...}
// _userdata : user-defined data structure (convenience)
// _v : input vector [size: _n x 1]
// _n : input vector size
float liquid_spiral(void * _userdata,
float * _v,
unsigned int _n);
//
// Gradient search
//
#define LIQUID_OPTIM_MINIMIZE (0)
#define LIQUID_OPTIM_MAXIMIZE (1)
typedef struct gradsearch_s * gradsearch;
// Create a gradient search object
// _userdata : user data object pointer
// _v : array of parameters to optimize
// _num_parameters : array length (number of parameters to optimize)
// _u : utility function pointer
// _direction : search direction (e.g. LIQUID_OPTIM_MAXIMIZE)
gradsearch gradsearch_create(void * _userdata,
float * _v,
unsigned int _num_parameters,
utility_function _utility,
int _direction);
// Destroy a gradsearch object
void gradsearch_destroy(gradsearch _q);
// Prints current status of search
void gradsearch_print(gradsearch _q);
// Iterate once
float gradsearch_step(gradsearch _q);
// Execute the search
float gradsearch_execute(gradsearch _q,
unsigned int _max_iterations,
float _target_utility);
// quasi-Newton search
typedef struct qnsearch_s * qnsearch;
// Create a simple qnsearch object; parameters are specified internally
// _userdata : userdata
// _v : array of parameters to optimize
// _num_parameters : array length
// _get_utility : utility function pointer
// _direction : search direction (e.g. LIQUID_OPTIM_MAXIMIZE)
qnsearch qnsearch_create(void * _userdata,
float * _v,
unsigned int _num_parameters,
utility_function _u,
int _direction);
// Destroy a qnsearch object
void qnsearch_destroy(qnsearch _g);
// Prints current status of search
void qnsearch_print(qnsearch _g);
// Resets internal state
void qnsearch_reset(qnsearch _g);
// Iterate once
void qnsearch_step(qnsearch _g);
// Execute the search
float qnsearch_execute(qnsearch _g,
unsigned int _max_iterations,
float _target_utility);
//
// chromosome (for genetic algorithm search)
//
typedef struct chromosome_s * chromosome;
// create a chromosome object, variable bits/trait
chromosome chromosome_create(unsigned int * _bits_per_trait,
unsigned int _num_traits);
// create a chromosome object, all traits same resolution
chromosome chromosome_create_basic(unsigned int _num_traits,
unsigned int _bits_per_trait);
// create a chromosome object, cloning a parent
chromosome chromosome_create_clone(chromosome _parent);
// copy existing chromosomes' internal traits (all other internal
// parameters must be equal)
void chromosome_copy(chromosome _parent, chromosome _child);
// Destroy a chromosome object
void chromosome_destroy(chromosome _c);
// get number of traits in chromosome
unsigned int chromosome_get_num_traits(chromosome _c);
// Print chromosome values to screen (binary representation)
void chromosome_print(chromosome _c);
// Print chromosome values to screen (floating-point representation)
void chromosome_printf(chromosome _c);
// clear chromosome (set traits to zero)
void chromosome_reset(chromosome _c);
// initialize chromosome on integer values
void chromosome_init(chromosome _c,
unsigned int * _v);
// initialize chromosome on floating-point values
void chromosome_initf(chromosome _c,
float * _v);
// Mutates chromosome _c at _index
void chromosome_mutate(chromosome _c, unsigned int _index);
// Resulting chromosome _c is a crossover of parents _p1 and _p2 at _threshold
void chromosome_crossover(chromosome _p1,
chromosome _p2,
chromosome _c,
unsigned int _threshold);
// Initializes chromosome to random value
void chromosome_init_random(chromosome _c);
// Returns integer representation of chromosome
unsigned int chromosome_value(chromosome _c,
unsigned int _index);
// Returns floating-point representation of chromosome
float chromosome_valuef(chromosome _c,
unsigned int _index);
//
// genetic algorithm search
//
typedef struct gasearch_s * gasearch;
typedef float (*gasearch_utility)(void * _userdata, chromosome _c);
// Create a simple gasearch object; parameters are specified internally
// _utility : chromosome fitness utility function
// _userdata : user data, void pointer passed to _get_utility() callback
// _parent : initial population parent chromosome, governs precision, etc.
// _minmax : search direction
gasearch gasearch_create(gasearch_utility _u,
void * _userdata,
chromosome _parent,
int _minmax);
// Create a gasearch object, specifying search parameters
// _utility : chromosome fitness utility function
// _userdata : user data, void pointer passed to _get_utility() callback
// _parent : initial population parent chromosome, governs precision, etc.
// _minmax : search direction
// _population_size : number of chromosomes in population
// _mutation_rate : probability of mutating chromosomes
gasearch gasearch_create_advanced(gasearch_utility _utility,
void * _userdata,
chromosome _parent,
int _minmax,
unsigned int _population_size,
float _mutation_rate);
// Destroy a gasearch object
void gasearch_destroy(gasearch _q);
// print search parameter internals
void gasearch_print(gasearch _q);
// set mutation rate
void gasearch_set_mutation_rate(gasearch _q,
float _mutation_rate);
// set population/selection size
// _q : ga search object
// _population_size : new population size (number of chromosomes)
// _selection_size : selection size (number of parents for new generation)
void gasearch_set_population_size(gasearch _q,
unsigned int _population_size,
unsigned int _selection_size);
// Execute the search
// _q : ga search object
// _max_iterations : maximum number of iterations to run before bailing
// _target_utility : target utility
float gasearch_run(gasearch _q,
unsigned int _max_iterations,
float _target_utility);
// iterate over one evolution of the search algorithm
void gasearch_evolve(gasearch _q);
// get optimal chromosome
// _q : ga search object
// _c : output optimal chromosome
// _utility_opt : fitness of _c
void gasearch_getopt(gasearch _q,
chromosome _c,
float * _utility_opt);
//
// MODULE : quantization
//
float compress_mulaw(float _x, float _mu);
float expand_mulaw(float _x, float _mu);
void compress_cf_mulaw(liquid_float_complex _x, float _mu, liquid_float_complex * _y);
void expand_cf_mulaw(liquid_float_complex _y, float _mu, liquid_float_complex * _x);
//float compress_alaw(float _x, float _a);
//float expand_alaw(float _x, float _a);
// inline quantizer: 'analog' signal in [-1, 1]
unsigned int quantize_adc(float _x, unsigned int _num_bits);
float quantize_dac(unsigned int _s, unsigned int _num_bits);
// structured quantizer
typedef enum {
LIQUID_COMPANDER_NONE=0,
LIQUID_COMPANDER_LINEAR,
LIQUID_COMPANDER_MULAW,
LIQUID_COMPANDER_ALAW
} liquid_compander_type;
#define QUANTIZER_MANGLE_FLOAT(name) LIQUID_CONCAT(quantizerf, name)
#define QUANTIZER_MANGLE_CFLOAT(name) LIQUID_CONCAT(quantizercf, name)
// large macro
// QUANTIZER : name-mangling macro
// T : data type
#define LIQUID_QUANTIZER_DEFINE_API(QUANTIZER,T) \
typedef struct QUANTIZER(_s) * QUANTIZER(); \
QUANTIZER() QUANTIZER(_create)(liquid_compander_type _ctype, \
float _range, \
unsigned int _num_bits); \
void QUANTIZER(_destroy)(QUANTIZER() _q); \
void QUANTIZER(_print)(QUANTIZER() _q); \
void QUANTIZER(_execute_adc)(QUANTIZER() _q, \
T _x, \
unsigned int * _sample); \
void QUANTIZER(_execute_dac)(QUANTIZER() _q, \
unsigned int _sample, \
T * _x);
LIQUID_QUANTIZER_DEFINE_API(QUANTIZER_MANGLE_FLOAT, float)
LIQUID_QUANTIZER_DEFINE_API(QUANTIZER_MANGLE_CFLOAT, liquid_float_complex)
//
// MODULE : random (number generators)
//
// Uniform random number generator, (0,1]
float randf();
float randf_pdf(float _x);
float randf_cdf(float _x);
// Gauss random number generator, N(0,1)
// f(x) = 1/sqrt(2*pi*sigma^2) * exp{-(x-eta)^2/(2*sigma^2)}
//
// where
// eta = mean
// sigma = standard deviation
//
float randnf();
void awgn(float *_x, float _nstd);
void crandnf(liquid_float_complex *_y);
void cawgn(liquid_float_complex *_x, float _nstd);
float randnf_pdf(float _x, float _eta, float _sig);
float randnf_cdf(float _x, float _eta, float _sig);
// Exponential
// f(x) = lambda exp{ -lambda x }
// where
// lambda = spread parameter, lambda > 0
// x >= 0
float randexpf(float _lambda);
float randexpf_pdf(float _x, float _lambda);
float randexpf_cdf(float _x, float _lambda);
// Weibull
// f(x) = (a/b) (x/b)^(a-1) exp{ -(x/b)^a }
// where
// a = alpha : shape parameter
// b = beta : scaling parameter
// g = gamma : location (threshold) parameter
//
float randweibf(float _alpha, float _beta, float _gamma);
float randweibf_pdf(float _x, float _a, float _b, float _g);
float randweibf_cdf(float _x, float _a, float _b, float _g);
// Gamma
// x^(a-1) exp(-x/b)
// f(x) = -------------------
// Gamma(a) b^a
// where
// a = alpha : shape parameter, a > 0
// b = beta : scale parameter, b > 0
// Gamma(z) = regular gamma function
// x >= 0
float randgammaf(float _alpha, float _beta);
float randgammaf_pdf(float _x, float _alpha, float _beta);
float randgammaf_cdf(float _x, float _alpha, float _beta);
// Nakagami-m
// f(x) = (2/Gamma(m)) (m/omega)^m x^(2m-1) exp{-(m/omega)x^2}
// where
// m : shape parameter, m >= 0.5
// omega : spread parameter, omega > 0
// Gamma(z): regular complete gamma function
// x >= 0
float randnakmf(float _m, float _omega);
float randnakmf_pdf(float _x, float _m, float _omega);
float randnakmf_cdf(float _x, float _m, float _omega);
// Rice-K
// f(x) = (x/sigma^2) exp{ -(x^2+s^2)/(2sigma^2) } I0( x s / sigma^2 )
// where
// s = sqrt( omega*K/(K+1) )
// sigma = sqrt(0.5 omega/(K+1))
// and
// K = shape parameter
// omega = spread parameter
// I0 = modified Bessel function of the first kind
// x >= 0
float randricekf(float _K, float _omega);
float randricekf_cdf(float _x, float _K, float _omega);
float randricekf_pdf(float _x, float _K, float _omega);
// Data scrambler : whiten data sequence
void scramble_data(unsigned char * _x, unsigned int _len);
void unscramble_data(unsigned char * _x, unsigned int _len);
void unscramble_data_soft(unsigned char * _x, unsigned int _len);
//
// MODULE : sequence
//
// Binary sequence (generic)
typedef struct bsequence_s * bsequence;
// Create a binary sequence of a specific length (number of bits)
bsequence bsequence_create(unsigned int num_bits);
// Free memory in a binary sequence
void bsequence_destroy(bsequence _bs);
// Clear binary sequence (set to 0's)
void bsequence_reset(bsequence _bs);
// initialize sequence on external array
void bsequence_init(bsequence _bs,
unsigned char * _v);
// Print sequence to the screen
void bsequence_print(bsequence _bs);
// Push bit into to back of a binary sequence
void bsequence_push(bsequence _bs,
unsigned int _bit);
// circular shift (left)
void bsequence_circshift(bsequence _bs);
// Correlate two binary sequences together
int bsequence_correlate(bsequence _bs1, bsequence _bs2);
// compute the binary addition of two bit sequences
void bsequence_add(bsequence _bs1, bsequence _bs2, bsequence _bs3);
// compute the binary multiplication of two bit sequences
void bsequence_mul(bsequence _bs1, bsequence _bs2, bsequence _bs3);
// accumulate the 1's in a binary sequence
unsigned int bsequence_accumulate(bsequence _bs);
// accessor functions
unsigned int bsequence_get_length(bsequence _bs);
unsigned int bsequence_index(bsequence _bs, unsigned int _i);
// Complementary codes
// intialize two sequences to complementary codes. sequences must
// be of length at least 8 and a power of 2 (e.g. 8, 16, 32, 64,...)
// _a : sequence 'a' (bsequence object)
// _b : sequence 'b' (bsequence object)
void bsequence_create_ccodes(bsequence _a,
bsequence _b);
// M-Sequence
#define LIQUID_MAX_MSEQUENCE_LENGTH 32767
// default m-sequence generators: g (hex) m n g (oct) g (binary)
#define LIQUID_MSEQUENCE_GENPOLY_M2 0x0007 // 2 3 7 111
#define LIQUID_MSEQUENCE_GENPOLY_M3 0x000B // 3 7 13 1011
#define LIQUID_MSEQUENCE_GENPOLY_M4 0x0013 // 4 15 23 10011
#define LIQUID_MSEQUENCE_GENPOLY_M5 0x0025 // 5 31 45 100101
#define LIQUID_MSEQUENCE_GENPOLY_M6 0x0043 // 6 63 103 1000011
#define LIQUID_MSEQUENCE_GENPOLY_M7 0x0089 // 7 127 211 10001001
#define LIQUID_MSEQUENCE_GENPOLY_M8 0x011D // 8 255 435 100101101
#define LIQUID_MSEQUENCE_GENPOLY_M9 0x0211 // 9 511 1021 1000010001
#define LIQUID_MSEQUENCE_GENPOLY_M10 0x0409 // 10 1023 2011 10000001001
#define LIQUID_MSEQUENCE_GENPOLY_M11 0x0805 // 11 2047 4005 100000000101
#define LIQUID_MSEQUENCE_GENPOLY_M12 0x1053 // 12 4095 10123 1000001010011
#define LIQUID_MSEQUENCE_GENPOLY_M13 0x201b // 13 8191 20033 10000000011011
#define LIQUID_MSEQUENCE_GENPOLY_M14 0x402b // 14 16383 40053 100000000101011
#define LIQUID_MSEQUENCE_GENPOLY_M15 0x8003 // 15 32767 100003 1000000000000011
typedef struct msequence_s * msequence;
// create a maximal-length sequence (m-sequence) object with
// an internal shift register length of _m bits.
// _m : generator polynomial length, sequence length is (2^m)-1
// _g : generator polynomial, starting with most-significant bit
// _a : initial shift register state, default: 000...001
msequence msequence_create(unsigned int _m,
unsigned int _g,
unsigned int _a);
// create a maximal-length sequence (m-sequence) object from a generator polynomial
msequence msequence_create_genpoly(unsigned int _g);
// creates a default maximal-length sequence
msequence msequence_create_default(unsigned int _m);
// destroy an msequence object, freeing all internal memory
void msequence_destroy(msequence _m);
// prints the sequence's internal state to the screen
void msequence_print(msequence _m);
// advance msequence on shift register, returning output bit
unsigned int msequence_advance(msequence _ms);
// generate pseudo-random symbol from shift register by
// advancing _bps bits and returning compacted symbol
// _ms : m-sequence object
// _bps : bits per symbol of output
unsigned int msequence_generate_symbol(msequence _ms,
unsigned int _bps);
// reset msequence shift register to original state, typically '1'
void msequence_reset(msequence _ms);
// initialize a bsequence object on an msequence object
// _bs : bsequence object
// _ms : msequence object
void bsequence_init_msequence(bsequence _bs,
msequence _ms);
// get the length of the sequence
unsigned int msequence_get_length(msequence _ms);
// get the internal state of the sequence
unsigned int msequence_get_state(msequence _ms);
// set the internal state of the sequence
void msequence_set_state(msequence _ms,
unsigned int _a);
//
// MODULE : utility
//
// pack binary array with symbol(s)
// _src : source array [size: _n x 1]
// _n : input source array length
// _k : bit index to write in _src
// _b : number of bits in input symbol
// _sym_in : input symbol
void liquid_pack_array(unsigned char * _src,
unsigned int _n,
unsigned int _k,
unsigned int _b,
unsigned char _sym_in);
// unpack symbols from binary array
// _src : source array [size: _n x 1]
// _n : input source array length
// _k : bit index to write in _src
// _b : number of bits in output symbol
// _sym_out : output symbol
void liquid_unpack_array(unsigned char * _src,
unsigned int _n,
unsigned int _k,
unsigned int _b,
unsigned char * _sym_out);
// pack one-bit symbols into bytes (8-bit symbols)
// _sym_in : input symbols array [size: _sym_in_len x 1]
// _sym_in_len : number of input symbols
// _sym_out : output symbols
// _sym_out_len : number of bytes allocated to output symbols array
// _num_written : number of output symbols actually written
void liquid_pack_bytes(unsigned char * _sym_in,
unsigned int _sym_in_len,
unsigned char * _sym_out,
unsigned int _sym_out_len,
unsigned int * _num_written);
// unpack 8-bit symbols (full bytes) into one-bit symbols
// _sym_in : input symbols array [size: _sym_in_len x 1]
// _sym_in_len : number of input symbols
// _sym_out : output symbols array
// _sym_out_len : number of bytes allocated to output symbols array
// _num_written : number of output symbols actually written
void liquid_unpack_bytes(unsigned char * _sym_in,
unsigned int _sym_in_len,
unsigned char * _sym_out,
unsigned int _sym_out_len,
unsigned int * _num_written);
// repack bytes with arbitrary symbol sizes
// _sym_in : input symbols array [size: _sym_in_len x 1]
// _sym_in_bps : number of bits per input symbol
// _sym_in_len : number of input symbols
// _sym_out : output symbols array
// _sym_out_bps : number of bits per output symbol
// _sym_out_len : number of bytes allocated to output symbols array
// _num_written : number of output symbols actually written
void liquid_repack_bytes(unsigned char * _sym_in,
unsigned int _sym_in_bps,
unsigned int _sym_in_len,
unsigned char * _sym_out,
unsigned int _sym_out_bps,
unsigned int _sym_out_len,
unsigned int * _num_written);
// shift array to the left _b bits, filling in zeros
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bits to shift
void liquid_lbshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// shift array to the right _b bits, filling in zeros
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bits to shift
void liquid_rbshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// circularly shift array to the left _b bits
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bits to shift
void liquid_lbcircshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// circularly shift array to the right _b bits
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bits to shift
void liquid_rbcircshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// shift array to the left _b bytes, filling in zeros
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bytes to shift
void liquid_lshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// shift array to the right _b bytes, filling in zeros
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bytes to shift
void liquid_rshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// circular shift array to the left _b bytes
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bytes to shift
void liquid_lcircshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// circular shift array to the right _b bytes
// _src : source address [size: _n x 1]
// _n : input data array size
// _b : number of bytes to shift
void liquid_rcircshift(unsigned char * _src,
unsigned int _n,
unsigned int _b);
// Count the number of ones in an integer
unsigned int liquid_count_ones(unsigned int _x);
// count number of ones in an integer, modulo 2
unsigned int liquid_count_ones_mod2(unsigned int _x);
// compute bindary dot-product between two integers
unsigned int liquid_bdotprod(unsigned int _x,
unsigned int _y);
// Count leading zeros in an integer
unsigned int liquid_count_leading_zeros(unsigned int _x);
// Most-significant bit index
unsigned int liquid_msb_index(unsigned int _x);
// Print string of bits to stdout
void liquid_print_bitstring(unsigned int _x,
unsigned int _n);
// reverse byte, word, etc.
unsigned char liquid_reverse_byte( unsigned char _x);
unsigned int liquid_reverse_uint16(unsigned int _x);
unsigned int liquid_reverse_uint24(unsigned int _x);
unsigned int liquid_reverse_uint32(unsigned int _x);
//
// MODULE : vector
//
#define VECTOR_MANGLE_RF(name) LIQUID_CONCAT(liquid_vectorf, name)
#define VECTOR_MANGLE_CF(name) LIQUID_CONCAT(liquid_vectorcf,name)
// large macro
// VECTOR : name-mangling macro
// T : data type
// TP : data type (primitive)
#define LIQUID_VECTOR_DEFINE_API(VECTOR,T,TP) \
\
/* initialize vector with scalar: x[i] = c (scalar) */ \
void VECTOR(_init)(T _c, \
T * _x, \
unsigned int _n); \
\
/* add each element: z[i] = x[i] + y[i] */ \
void VECTOR(_add)(T * _x, \
T * _y, \
unsigned int _n, \
T * _z); \
/* add scalar to each element: y[i] = x[i] + c */ \
void VECTOR(_addscalar)(T * _x, \
unsigned int _n, \
T _c, \
T * _y); \
\
/* multiply each element: z[i] = x[i] * y[i] */ \
void VECTOR(_mul)(T * _x, \
T * _y, \
unsigned int _n, \
T * _z); \
/* multiply each element with scalar: y[i] = x[i] * c */ \
void VECTOR(_mulscalar)(T * _x, \
unsigned int _n, \
T _c, \
T * _y); \
\
/* compute complex phase rotation: x[i] = exp{j theta[i]} */ \
void VECTOR(_cexpj)(TP * _theta, \
unsigned int _n, \
T * _x); \
/* compute angle of each element: theta[i] = arg{ x[i] } */ \
void VECTOR(_carg)(T * _x, \
unsigned int _n, \
TP * _theta); \
/* compute absolute value of each element: y[i] = |x[i]| */ \
void VECTOR(_abs)(T * _x, \
unsigned int _n, \
TP * _y); \
\
/* compute sum of squares: sum{ |x|^2 } */ \
TP VECTOR(_sumsq)(T * _x, \
unsigned int _n); \
\
/* compute l-2 norm: sqrt{ sum{ |x|^2 } } */ \
TP VECTOR(_norm)(T * _x, \
unsigned int _n); \
\
/* compute l-p norm: { sum{ |x|^p } }^(1/p) */ \
TP VECTOR(_pnorm)(T * _x, \
unsigned int _n, \
TP _p); \
\
/* scale vector elements by l-2 norm: y[i] = x[i]/norm(x) */ \
void VECTOR(_normalize)(T * _x, \
unsigned int _n, \
T * _y); \
LIQUID_VECTOR_DEFINE_API(VECTOR_MANGLE_RF, float, float);
LIQUID_VECTOR_DEFINE_API(VECTOR_MANGLE_CF, liquid_float_complex, float);
//
// mixed types
//
#if 0
void liquid_vectorf_add(float * _a,
float * _b,
unsigned int _n,
float * _c);
#endif
#ifdef __cplusplus
} //extern "C"
#endif // __cplusplus
#ifdef _MSC_VER
#pragma warning( pop )
#endif
#endif // __LIQUID_H__
Reference in New Issue View Git Blame Copy Permalink