// execute synchronizer object on buffer of samples
void ofdmflexframesync_execute(ofdmflexframesync _q,
float complex * _x,
unsigned int _n)
{
// push samples through ofdmframesync object
ofdmframesync_execute(_q->fs, _x, _n);
}
// Helper function to keep code base small
// _num_subcarriers : number of subcarriers
// _cp_len : cyclic prefix lenght
// _taper_len : taper length
void ofdmframesync_acquire_test(unsigned int _num_subcarriers,
unsigned int _cp_len,
unsigned int _taper_len)
{
// options
unsigned int M = _num_subcarriers; // number of subcarriers
unsigned int cp_len = _cp_len; // cyclic prefix lenght
unsigned int taper_len = _taper_len; // taper length
float tol = 1e-2f; // error tolerance
//
float dphi = 1.0f / (float)M; // carrier frequency offset
// subcarrier allocation (initialize to default)
unsigned char p[M];
ofdmframe_init_default_sctype(M, p);
// derived values
unsigned int num_samples = (3 + 1)*(M + cp_len);
// create synthesizer/analyzer objects
ofdmframegen fg = ofdmframegen_create(M, cp_len, taper_len, p);
//ofdmframegen_print(fg);
float complex X[M]; // original data sequence
float complex X_test[M]; // recovered data sequence
ofdmframesync fs = ofdmframesync_create(M,cp_len,taper_len,p,ofdmframesync_autotest_callback,(void*)X_test);
unsigned int i;
float complex s0[M]; // short PLCP sequence
float complex s1[M]; // long PLCP sequence
float complex y[num_samples]; // frame samples
// assemble full frame
unsigned int n=0;
// write first S0 symbol
ofdmframegen_write_S0a(fg, &y[n]);
n += M + cp_len;
// write second S0 symbol
ofdmframegen_write_S0b(fg, &y[n]);
n += M + cp_len;
// write S1 symbol
ofdmframegen_write_S1( fg, &y[n]);
n += M + cp_len;
// generate data symbol (random)
for (i=0; i<M; i++) {
X[i] = cexpf(_Complex_I*2*M_PI*randf());
X_test[i] = 0.0f;
}
// write data symbol
ofdmframegen_writesymbol(fg, X, &y[n]);
n += M + cp_len;
// validate frame length
assert(n == num_samples);
// add carrier offset
for (i=0; i<num_samples; i++)
y[i] *= cexpf(_Complex_I*dphi*i);
// run receiver
ofdmframesync_execute(fs,y,num_samples);
// check output
for (i=0; i<M; i++) {
if (p[i] == OFDMFRAME_SCTYPE_DATA) {
float e = crealf( (X[i] - X_test[i])*conjf(X[i] - X_test[i]) );
CONTEND_DELTA( cabsf(e), 0.0f, tol );
}
}
// destroy objects
ofdmframegen_destroy(fg);
ofdmframesync_destroy(fs);
}
// Helper function to keep code base small
void ofdmframesync_rxsymbol_bench(struct rusage *_start,
struct rusage *_finish,
unsigned long int *_num_iterations,
unsigned int _num_subcarriers,
unsigned int _cp_len)
{
// options
modulation_scheme ms = LIQUID_MODEM_QPSK;
unsigned int M = _num_subcarriers;
unsigned int cp_len = _cp_len;
unsigned int taper_len = 0;
// create synthesizer/analyzer objects
ofdmframegen fg = ofdmframegen_create(M, cp_len, taper_len, NULL);
//ofdmframegen_print(fg);
modem mod = modem_create(ms);
ofdmframesync fs = ofdmframesync_create(M,cp_len,taper_len,NULL,NULL,NULL);
unsigned int i;
float complex X[M]; // channelized symbol
float complex x[M+cp_len]; // time-domain symbol
// synchronize short sequence (first)
ofdmframegen_write_S0a(fg, x);
ofdmframesync_execute(fs, x, M+cp_len);
// synchronize short sequence (second)
ofdmframegen_write_S0b(fg, x);
ofdmframesync_execute(fs, x, M+cp_len);
// synchronize long sequence
ofdmframegen_write_S1(fg, x);
ofdmframesync_execute(fs, x, M+cp_len);
// modulate data symbols (use same symbol, ignore pilot phase)
unsigned int s;
for (i=0; i<M; i++) {
s = modem_gen_rand_sym(mod);
modem_modulate(mod,s,&X[i]);
}
ofdmframegen_writesymbol(fg, X, x);
// add noise
for (i=0; i<M+cp_len; i++)
x[i] += 0.02f*randnf()*cexpf(_Complex_I*2*M_PI*randf());
// normalize number of iterations
*_num_iterations /= M;
// start trials
getrusage(RUSAGE_SELF, _start);
for (i=0; i<(*_num_iterations); i++) {
// receive data symbols (ignoring pilots)
ofdmframesync_execute(fs, x, M+cp_len);
ofdmframesync_execute(fs, x, M+cp_len);
ofdmframesync_execute(fs, x, M+cp_len);
ofdmframesync_execute(fs, x, M+cp_len);
}
getrusage(RUSAGE_SELF, _finish);
*_num_iterations *= 4;
// destroy objects
ofdmframegen_destroy(fg);
ofdmframesync_destroy(fs);
modem_destroy(mod);
}
