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

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

3G开发

开发平台：

Visual C++

qamoptdem_bp.cpp：源码内容
  1. //
  2. //  File = qamoptdem_bp.cpp
  3. //
  5. #include <stdlib.h>
  6. #include <fstream.h>
  7. #include "parmfile.h"
  8. #include "qamoptdem_bp.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. QamOptimalBandpassDemod::QamOptimalBandpassDemod( char* instance_name,
  20.                                   PracSimModel* outer_model,
  21.                                   Signal< float >* in_sig,
  22.                                   Signal< bit_t >* symb_clock_in,
  23.                                   Signal< byte_t >* out_sig )
  24.                 :PracSimModel(instance_name,
  25.                               outer_model)
  26. {
  27.   MODEL_NAME(QamOptimalBandpassDemod);
  29.   ENABLE_MULTIRATE;
  31.   //-----------------------------------------
  32.   //  Read model config parms
  34.   OPEN_PARM_BLOCK;
  36.   GET_DOUBLE_PARM(Phase_Unbal);
  37.   GET_DOUBLE_PARM(Carrier_Freq);
  38.   GET_INT_PARM(Bits_Per_Symb);
  39.   GET_INT_PARM(Samps_Per_Symb);
  40.   GET_DOUBLE_PARM(Dly_To_Start);
  42.   Carrier_Freq_Rad = TWO_PI * Carrier_Freq;
  44.   //int block_size = ParmInput.GetIntParm("block_size");
  45.   //double samp_intvl = ParmInput.GetDoubleParm("samp_intvl");
  47.   //--------------------------------------
  48.   //  Connect input and output signals
  50.   Out_Sig = out_sig;
  51.   Symb_Clock_In = symb_clock_in;
  52.   In_Sig = in_sig;
  54.   MAKE_OUTPUT( Out_Sig );
  55.   MAKE_INPUT( Symb_Clock_In );
  56.   MAKE_INPUT( In_Sig );
  58.   double resamp_rate = 1.0/double(Samps_Per_Symb);
  59.   CHANGE_RATE( In_Sig, Out_Sig, resamp_rate );
  60.   CHANGE_RATE( Symb_Clock_In, Out_Sig, resamp_rate );
  62.   //------------------
  63.   //  compute decision boundaries
  65.   Num_Diff_Symbs = 1;
  66.   for(int i=1; i<=Bits_Per_Symb; i++)
  67.     Num_Diff_Symbs *=2;
  69.   double offset = TWO_PI/(2.0*Num_Diff_Symbs);
  70.   Decis_Bound = new float[Num_Diff_Symbs];
  71.   for(int isymb=0; isymb<Num_Diff_Symbs; isymb++)
  72.     {
  73.     Decis_Bound[isymb] = float(offset + isymb*TWO_PI/double(Num_Diff_Symbs));
  74.     }
  75. }
  76. //==============================================
  77. QamOptimalBandpassDemod::~QamOptimalBandpassDemod( void ){ };
  78. //==============================================
  80. void QamOptimalBandpassDemod::Initialize(void)
  81. {
  82.   Block_Size = In_Sig->GetBlockSize();
  83. //  Out_Samp_Intvl = Out_Sig->GetSampIntvl();
  84.   //
  85.   //  set up table of phase references
  86.   Ref_Angle = new float[Num_Diff_Symbs];
  87.   Integ_Val = new double[Num_Diff_Symbs];
  89.   for( int isymb=0; isymb<Num_Diff_Symbs; isymb++)
  90.     {
  91.     Integ_Val[isymb] = 0.0;
  92.     Ref_Angle[isymb] = TWO_PI * isymb/double(Num_Diff_Symbs);
  93.     }
  94.   Time = 0.0;
  95.   Samp_Intvl = In_Sig->GetSampIntvl();
  96. }
  98. //============================================
  99. int QamOptimalBandpassDemod::Execute()
  100. {
  101.   byte_t *out_sig_ptr;
  102.   float *in_sig_ptr;
  103.   bit_t *symb_clock_in_ptr;
  104.   float in_val;
  105.   double *integ_val;
  106.   double max_val=0.0;
  107.   double omega, time, time_base;
  108.   double samp_intvl;
  109.   int is;
  110.   byte_t isymb, symb_decis;
  111.   #ifdef _DEBUG
  112.     *DebugFile << "In QamOptimalBandpassDemod::Execute" << endl;
  113.   #endif
  115.   in_sig_ptr = GET_INPUT_PTR( In_Sig );
  116.   symb_clock_in_ptr = GET_INPUT_PTR( Symb_Clock_In );
  117.   out_sig_ptr = GET_OUTPUT_PTR( Out_Sig );
  119.   omega = Carrier_Freq_Rad;
  120.   time_base = Time;
  121.   samp_intvl = Samp_Intvl;
  122.   integ_val = Integ_Val;
  124.   for (is=0; is<Block_Size; is++)
  125.     {
  126.     in_val = *in_sig_ptr++;
  128.     time = time_base + is * samp_intvl;
  129.     // correlate input signal against all possible reference phases
  130.     for( isymb=0; isymb<Num_Diff_Symbs; isymb++)
  131.       {
  132.       Integ_Val[isymb] += in_val * cos(omega*time + Ref_Angle[isymb]);
  133.       }
  135.     if(*symb_clock_in_ptr != 0)
  136.       {
  137.       // time to make a decision
  139.       max_val = Integ_Val[0];
  140.       symb_decis = 0;
  141.       for(isymb=1; isymb<Num_Diff_Symbs; isymb++)
  142.         {
  143.         if(Integ_Val[isymb] > max_val)
  144.           {
  145.           max_val = Integ_Val[isymb];
  146.           symb_decis = isymb;
  147.           }
  148.         }
  149.       *out_sig_ptr = symb_decis;
  150.       out_sig_ptr++;
  152.       for(isymb=0; isymb<Num_Diff_Symbs; isymb++)
  153.         {
  154.         Integ_Val[isymb] = 0.0;
  155.         }
  156.       }
  157.     symb_clock_in_ptr++;
  158.     }
  159.   Time = time_base + samp_intvl * Block_Size;
  160.   //Integ_Val = integ_val;
  161.   return(_MES_AOK);
  162. }