开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 3G开发 > 查看源码
bpsk_correl_demod.cpp
资源名称：simulation_wireless_c++.rar　[点击查看]
上传用户：jtjnyq9001
上传日期：2014-11-21
资源大小：3974k
文件大小：4k
源码类别：

3G开发

开发平台：

Visual C++

bpsk_correl_demod.cpp：源码内容
  1. //
  2. //  File = bpsk_correl_demod.cpp
  3. //
  5. #include <stdlib.h>
  6. #include <fstream>
  7. #include "parmfile.h"
  8. #include "bpsk_correl_demod.h"
  9. #include "misdefs.h"
  10. #include "model_graph.h"
  11. extern ParmFile *ParmInput;
  12. extern int PassNumber;
  13. #ifdef _DEBUG
  14.   extern ofstream *DebugFile;
  15. #endif
  17. //======================================================
  19. BpskCorrelationDemod::BpskCorrelationDemod( char* instance_name,
  20.                                   PracSimModel* outer_model,
  21.                                   Signal< std::complex< float > >* in_sig,
  22.                                   Signal< std::complex<float> >* phase_ref_sig,
  23.                                   Signal< bit_t >* symb_clock_in,
  24.                                   Signal< bit_t >* decis_out )
  25.                 :PracSimModel(instance_name,
  26.                               outer_model)
  27. {
  28.   MODEL_NAME(BpskCorrelationDemod);
  30.   ENABLE_MULTIRATE;
  32.   //-----------------------------------------
  33.   //  Read model config parms
  35.   OPEN_PARM_BLOCK;
  37.   GET_INT_PARM(Samps_Per_Symb);
  38.   GET_DOUBLE_PARM(Dly_To_Start);
  40.   //--------------------------------------
  41.   //  Connect input and output signals
  43.   Decis_Out = decis_out;
  44.   Symb_Clock_In = symb_clock_in;
  45.   In_Sig = in_sig;
  46.   Phase_Ref_Sig = phase_ref_sig;
  48.   MAKE_OUTPUT( Decis_Out );
  49.   MAKE_INPUT( Symb_Clock_In );
  50.   MAKE_INPUT( Phase_Ref_Sig );
  51.   MAKE_INPUT( In_Sig );
  53.   double resamp_rate = 1.0/double(Samps_Per_Symb);
  54.   CHANGE_RATE( In_Sig, Decis_Out, resamp_rate );
  55.   CHANGE_RATE( Phase_Ref_Sig, Decis_Out, resamp_rate );
  56.   CHANGE_RATE( Symb_Clock_In, Decis_Out, resamp_rate );
  58. }
  59. //==============================================
  60. BpskCorrelationDemod::~BpskCorrelationDemod( void ){ };
  61. //==============================================
  63. void BpskCorrelationDemod::Initialize(void)
  64. {
  65.   Block_Size = In_Sig->GetBlockSize();
  66.   Out_Samp_Intvl = Decis_Out->GetSampIntvl();
  68.   Integ_Val = 0.0;
  70. }
  72. //============================================
  73. int BpskCorrelationDemod::Execute()
  74. {
  75.   bit_t *decis_out_ptr;
  76.   std::complex<float> *in_sig_ptr;
  77.   bit_t *symb_clock_in_ptr;
  78.   std::complex<float> in_val;
  79.   std::complex<float> *phase_ref_sig_ptr;
  80.   std::complex<float> phase_ref_val;
  81.   double integ_val;
  82.   double max_val=0.0;
  83.   int is;
  84.   #ifdef _DEBUG
  85.     *DebugFile << "In BpskCorrelationDemod::Execute" << endl;
  86.   #endif
  88.   int valid_block_size = In_Sig->GetValidBlockSize();
  89.   int out_block_size = valid_block_size/Samps_Per_Symb;
  90.   if(valid_block_size != Samps_Per_Symb*out_block_size) {
  91.       //error
  92.       }
  93.       else {
  94.          Decis_Out->SetValidBlockSize(out_block_size);
  95.       }
  97.   in_sig_ptr = GET_INPUT_PTR( In_Sig );
  98.   symb_clock_in_ptr = GET_INPUT_PTR( Symb_Clock_In );
  99.   phase_ref_sig_ptr = GET_INPUT_PTR( Phase_Ref_Sig );
  100.   decis_out_ptr = GET_OUTPUT_PTR( Decis_Out );
  102.   integ_val = Integ_Val;
  104.   for (is=0; is<valid_block_size; is++)
  105.     {
  106.     in_val = *in_sig_ptr++;
  108.     // correlate input signal against recovered carrier
  109.     phase_ref_val = *phase_ref_sig_ptr++;
  110. //    Integ_Val += std::real(in_val * 
  111. //                  std::complex<float>( cos(phase_ref_val), -sin(phase_ref_val) ));
  112.     Integ_Val += std::real(in_val * phase_ref_val);
  114.     if(*symb_clock_in_ptr != 0)
  115.       {
  116.       // time to make a decision
  118.       if(Integ_Val > 0.0 )
  119.         {
  120.         *decis_out_ptr++ = 1;
  121.         }
  122.       else
  123.         {
  124.         *decis_out_ptr++ = 0;
  125.         }
  127.       Integ_Val = 0.0;
  128.       }
  129.     symb_clock_in_ptr++;
  130.     }
  131.   //Integ_Val = integ_val;
  132.   return(_MES_AOK);
  133. }