fsk_demod_coh_bp.cpp
上传用户:jtjnyq9001
上传日期:2014-11-21
资源大小:3974k
文件大小:4k
源码类别:

3G开发

开发平台:

Visual C++

  1. //
  2. //  File = fsk_demod_coh_bp.cpp
  3. //
  4. #include <stdlib.h>
  5. //#include <fstream>
  6. #include "parmfile.h"
  7. #include "fsk_demod_coh_bp.h"
  8. #include "misdefs.h"
  9. #include "model_graph.h"
  10. extern ParmFile *ParmInput;
  11. extern int PassNumber;
  12. #ifdef _DEBUG
  13.   extern ofstream *DebugFile;
  14. #endif
  15. //======================================================
  16. FskCoherentBandpassDemod::FskCoherentBandpassDemod( char* instance_name,
  17.                                   PracSimModel* outer_model,
  18.                                   Signal< float >* in_sig,
  19.                                   Signal< float >* lo_ref_sig,
  20.                                   Signal< float >* hi_ref_sig,
  21.                                   Signal< bit_t >* symb_clock_in,
  22.                                   Signal< bit_t >* decis_out,
  23.                                   Signal< float >* lo_integ_sig,
  24.                                   Signal< float >* hi_integ_sig )
  25.                 :PracSimModel(instance_name,
  26.                               outer_model)
  27. {
  28.   MODEL_NAME(FskCoherentBandpassDemod);
  29.   ENABLE_MULTIRATE;
  30.   //-----------------------------------------
  31.   //  Read model config parms
  32.   OPEN_PARM_BLOCK;
  33.   GET_INT_PARM(Samps_Per_Symb);
  34.   GET_DOUBLE_PARM(Dly_To_Start);
  35.   //--------------------------------------
  36.   //  Connect input and output signals
  37.   Decis_Out = decis_out;
  38.   Symb_Clock_In = symb_clock_in;
  39.   In_Sig = in_sig;
  40.   Lo_Ref_Sig = lo_ref_sig;
  41.   Hi_Ref_Sig = hi_ref_sig;
  42.   Lo_Integ_Sig = lo_integ_sig;
  43.   Hi_Integ_Sig = hi_integ_sig;
  44.   MAKE_OUTPUT( Decis_Out );
  45.   MAKE_OUTPUT( Lo_Integ_Sig );
  46.   MAKE_OUTPUT( Hi_Integ_Sig );
  47.   MAKE_INPUT( Symb_Clock_In );
  48.   MAKE_INPUT( In_Sig );
  49.   MAKE_INPUT( Lo_Ref_Sig );
  50.   MAKE_INPUT( Hi_Ref_Sig );
  51.   double resamp_rate = 1.0/double(Samps_Per_Symb);
  52.   CHANGE_RATE( In_Sig, Decis_Out, resamp_rate );
  53.   CHANGE_RATE( Lo_Ref_Sig, Decis_Out, resamp_rate );
  54.   CHANGE_RATE( Hi_Ref_Sig, Decis_Out, resamp_rate );
  55.   CHANGE_RATE( Symb_Clock_In, Decis_Out, resamp_rate );
  56.   SAME_RATE( Lo_Ref_Sig, Lo_Integ_Sig );
  57.   SAME_RATE( Hi_Ref_Sig, Hi_Integ_Sig );
  58. }
  59. //==============================================
  60. FskCoherentBandpassDemod::~FskCoherentBandpassDemod( void ){ };
  61. //==============================================
  62. void FskCoherentBandpassDemod::Initialize(void)
  63. {
  64.   Block_Size = In_Sig->GetBlockSize();
  65.   Out_Samp_Intvl = Decis_Out->GetSampIntvl();
  66.   Lo_Integ_Val = 0.0;
  67.   Hi_Integ_Val = 0.0;
  68. }
  69. //============================================
  70. int FskCoherentBandpassDemod::Execute()
  71. {
  72.   bit_t *decis_out_ptr;
  73.   float *in_sig_ptr;
  74.   float *lo_ref_sig_ptr;
  75.   float *hi_ref_sig_ptr;
  76.   float *lo_integ_sig_ptr;
  77.   float *hi_integ_sig_ptr;
  78.   bit_t *symb_clock_in_ptr;
  79.   float in_val;
  80.   double lo_integ_val, hi_integ_val;
  81.   double max_val=0.0;
  82.   int is;
  83.   #ifdef _DEBUG
  84.     *DebugFile << "In FskCoherentBandpassDemod::Execute" << endl;
  85.   #endif
  86.   in_sig_ptr = GET_INPUT_PTR( In_Sig );
  87.   lo_ref_sig_ptr = GET_INPUT_PTR( Lo_Ref_Sig );
  88.   hi_ref_sig_ptr = GET_INPUT_PTR( Hi_Ref_Sig );
  89.   symb_clock_in_ptr = GET_INPUT_PTR( Symb_Clock_In );
  90.   decis_out_ptr = GET_OUTPUT_PTR( Decis_Out );
  91.   lo_integ_sig_ptr = GET_OUTPUT_PTR( Lo_Integ_Sig );
  92.   hi_integ_sig_ptr = GET_OUTPUT_PTR( Hi_Integ_Sig );
  93.   lo_integ_val = Lo_Integ_Val;
  94.   hi_integ_val = Hi_Integ_Val;
  95.   for (is=0; is<Block_Size; is++)
  96.     {
  97.     in_val = *in_sig_ptr++;
  98.     // correlate input signal against recovered carrier
  99.     lo_integ_val += in_val * (*lo_ref_sig_ptr++);
  100.     hi_integ_val += in_val * (*hi_ref_sig_ptr++);
  101.     *lo_integ_sig_ptr++ = float(lo_integ_val/float(Samps_Per_Symb));
  102.     *hi_integ_sig_ptr++ = float(hi_integ_val/float(Samps_Per_Symb));
  103.     if(*symb_clock_in_ptr != 0)
  104.       {
  105.       // time to make a decision
  106.       if(hi_integ_val > lo_integ_val )
  107.         {
  108.         *decis_out_ptr++ = 1;
  109.         }
  110.       else
  111.         {
  112.         *decis_out_ptr++ = 0;
  113.         }
  114.       lo_integ_val = 0.0;
  115.       hi_integ_val = 0.0;
  116.       }
  117.     symb_clock_in_ptr++;
  118.     }
  119.   Lo_Integ_Val = lo_integ_val;
  120.   Hi_Integ_Val = hi_integ_val;
  121.   return(_MES_AOK);
  122. }