fsk_genie.cpp

this is the simulation model for wireless communication

//
//  File = fsk_genie.cpp
//

#include <stdlib.h>
//#include <fstream>
//#include <strstream>
#include "parmfile.h"
#include "misdefs.h"
#include "fsk_genie.h"
#include "model_graph.h"
extern ParmFile *ParmInput;
extern int PassNumber;

//======================================================

FskCarrierGenie::FskCarrierGenie( char* instance_name,
                                PracSimModel* outer_model,
                                Signal< std::complex<float> > *lo_ref_sig,
                                Signal< std::complex<float> > *hi_ref_sig )
              :PracSimModel(instance_name,
                            outer_model)
{
  MODEL_NAME(FskCarrierGenie);

  //-----------------------------------
  // Read configuration parameters
  OPEN_PARM_BLOCK;

  GET_DOUBLE_PARM(Carrier_Freq_Hz);
  GET_DOUBLE_PARM(Freq_Offset_Hz);
  GET_INT_PARM( Lo_Ref_Samp_Shift );
  GET_INT_PARM( Hi_Ref_Samp_Shift );

  //-----------------------------------
  //  Signals

  Lo_Ref_Sig = lo_ref_sig;
  Hi_Ref_Sig = hi_ref_sig;

  MAKE_OUTPUT( Lo_Ref_Sig );
  MAKE_OUTPUT( Hi_Ref_Sig );
  //----------------------------------
  //  Compute derived parameters

  Freq_Hi_Hz = Carrier_Freq_Hz + Freq_Offset_Hz;
  Freq_Lo_Hz = Carrier_Freq_Hz - Freq_Offset_Hz;


  return;
  }
//================================================
FskCarrierGenie::~FskCarrierGenie( void ){ };

//=======================================================

void FskCarrierGenie::Initialize()
{
  //---------------------------------------
  //  Initialize derived parameters
  Block_Size = Lo_Ref_Sig->GetBlockSize();
  Samp_Intvl = Lo_Ref_Sig->GetSampIntvl();
  Cumul_Phase_Lo = Lo_Ref_Samp_Shift * Samp_Intvl * Freq_Lo_Hz;
  Cumul_Phase_Hi = Hi_Ref_Samp_Shift * Samp_Intvl * Freq_Hi_Hz;
}

//=======================================
int FskCarrierGenie::Execute()
{
  std::complex<float> *lo_ref_sig_ptr;
  std::complex<float> *hi_ref_sig_ptr;
  double cumul_phase_lo, cumul_phase_hi;
  double freq_lo_hz, freq_hi_hz;
  double samp_intvl;
  int is;

  lo_ref_sig_ptr = GET_OUTPUT_PTR( Lo_Ref_Sig );
  hi_ref_sig_ptr = GET_OUTPUT_PTR( Hi_Ref_Sig );

  cumul_phase_lo = Cumul_Phase_Lo;
  cumul_phase_hi = Cumul_Phase_Hi;
  freq_lo_hz = Freq_Lo_Hz;
  freq_hi_hz = Freq_Hi_Hz;

  samp_intvl = Samp_Intvl;

  for(is=0; is<Block_Size; is++)
    {
    cumul_phase_lo = fmod(cumul_phase_lo + samp_intvl * freq_lo_hz, 1.0);
    cumul_phase_hi = fmod(cumul_phase_hi + samp_intvl * freq_hi_hz, 1.0);
    *lo_ref_sig_ptr++ = std::complex<float>( float(cos(TWO_PI * cumul_phase_lo)),
                                             float(sin(TWO_PI * cumul_phase_lo)));
    *hi_ref_sig_ptr++ = std::complex<float>( float(cos(TWO_PI * cumul_phase_hi)),
                                             float(sin(TWO_PI * cumul_phase_hi)));
    }
  Cumul_Phase_Lo = cumul_phase_lo;
  Cumul_Phase_Hi = cumul_phase_hi;
  return(_MES_AOK);

}