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

3G开发

开发平台:

Visual C++

  1. //
  2. //  File = qamoptimdem.cpp
  3. //
  4. #include <stdlib.h>
  5. //#include <fstream>
  6. #include "parmfile.h"
  7. #include "qamoptimdem.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. QamOptimalDemod::QamOptimalDemod( char* instance_name,
  17.                                   PracSimModel* outer_model,
  18.                                   Signal< std::complex< float > >* in_sig,
  19.                                   Signal< bit_t >* symb_clock_in,
  20.                                   Signal< byte_t >* out_sig )
  21.                 :PracSimModel(instance_name,
  22.                               outer_model)
  23. {
  24.   MODEL_NAME(QamOptimalDemod);
  25.   ENABLE_MULTIRATE;
  26.   //-----------------------------------------
  27.   //  Read model config parms
  28.   OPEN_PARM_BLOCK;
  29.   GET_INT_PARM(Bits_Per_Symb);
  30.   GET_INT_PARM(Samps_Per_Symb);
  31.   //--------------------------------------
  32.   //  Connect input and output signals
  33.   Out_Sig = out_sig;
  34.   Symb_Clock_In = symb_clock_in;
  35.   In_Sig = in_sig;
  36.   MAKE_OUTPUT( Out_Sig );
  37.   MAKE_INPUT( Symb_Clock_In );
  38.   MAKE_INPUT( In_Sig );
  39.   double resamp_rate = 1.0/double(Samps_Per_Symb);
  40.   CHANGE_RATE( In_Sig, Out_Sig, resamp_rate );
  41.   CHANGE_RATE( Symb_Clock_In, Out_Sig, resamp_rate );
  42.   //------------------
  43.   //  compute decision boundaries
  44.   Num_Symb_Rows = 1;
  45.   for(int i=1; i<=Bits_Per_Symb/2; i++)
  46.     Num_Symb_Rows *=2;
  47.   I_Boundary = new double[Num_Symb_Rows-1];
  48.   Q_Boundary = new double[Num_Symb_Rows-1];
  49.   for(int isymb=0; isymb<Num_Symb_Rows; isymb++)
  50.     {
  51.     I_Boundary[isymb] = double(Samps_Per_Symb*(2 - Num_Symb_Rows + 2 * isymb));
  52.     Q_Boundary[isymb] = double(Samps_Per_Symb*(2 - Num_Symb_Rows + 2 * isymb));
  53.     }
  54. }
  55. //==============================================
  56. QamOptimalDemod::~QamOptimalDemod( void ){ };
  57. //==============================================
  58. void QamOptimalDemod::Initialize(void)
  59. {
  60.   Block_Size = In_Sig->GetBlockSize();
  61.   Out_Samp_Intvl = Out_Sig->GetSampIntvl();
  62.   I_Integ_Val = 0.0;
  63.   Q_Integ_Val = 0.0;
  64. }
  65. //============================================
  66. int QamOptimalDemod::Execute()
  67. {
  68.   byte_t *out_sig_ptr;
  69.   byte_t i_decis, q_decis;
  70.   std::complex<float> *in_sig_ptr;
  71.   bit_t *symb_clock_in_ptr;
  72.   std::complex<float> in_val;
  73.   double i_integ_val, q_integ_val;
  74.   double max_val=0.0;
  75.   int is;
  76.   byte_t isymb, symb_decis;
  77.   #ifdef _DEBUG
  78.     *DebugFile << "In QamOptimalDemod::Execute" << endl;
  79.   #endif
  80.   in_sig_ptr = GET_INPUT_PTR( In_Sig );
  81.   symb_clock_in_ptr = GET_INPUT_PTR( Symb_Clock_In );
  82.   out_sig_ptr = GET_OUTPUT_PTR( Out_Sig );
  83.   i_integ_val = I_Integ_Val;
  84.   q_integ_val = Q_Integ_Val;
  85.   for (is=0; is<Block_Size; is++)
  86.     {
  87.     in_val = *in_sig_ptr++;
  88.     I_Integ_Val += std::real(in_val);
  89.     Q_Integ_Val += std::imag(in_val);
  90.     if(*symb_clock_in_ptr != 0)
  91.       {
  92.       // time to make a decision
  93.       i_decis = 0;
  94.       for(isymb=0; isymb<Num_Symb_Rows-1; isymb++)
  95.         {
  96.         if(I_Integ_Val < I_Boundary[isymb]) break;
  97.         i_decis++;
  98.         }
  99.       q_decis = 0;
  100.       for(isymb=0; isymb<Num_Symb_Rows-1; isymb++)
  101.         {
  102.         if(Q_Integ_Val < Q_Boundary[isymb]) break;
  103.         q_decis++;
  104.         }
  105.       symb_decis = q_decis + Num_Symb_Rows * i_decis;
  106.       *out_sig_ptr = symb_decis;
  107.       out_sig_ptr++;
  108.       I_Integ_Val = 0.0;
  109.       Q_Integ_Val = 0.0;
  110.       }
  111.     symb_clock_in_ptr++;
  112.     }
  113.   //Integ_Val = integ_val;
  114.   return(_MES_AOK);
  115. }