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

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

3G开发

开发平台：

Visual C++

mskmod.cpp：源码内容
  1. //
  2. //  File = mskmod.cpp
  3. //
  5. #include <stdlib.h>
  6. #include "parmfile.h"
  7. #include "mskmod.h"
  8. #include "misdefs.h"
  9. #include "model_graph.h"
  10. extern ParmFile *ParmInput;
  11. using std::complex;
  13. //======================================================
  15. MskModulator::MskModulator( 
  16.                char* instance_name,
  17.                PracSimModel* outer_model,
  18.                Signal< float >* i_in_sig,
  19.                Signal< float >* q_in_sig,
  20.                Signal< complex<float> >* cmpx_out_sig,
  21.                Signal< float >* mag_out_sig,
  22.                Signal< float >* phase_out_sig )
  23.           :PracSimModel(instance_name,
  24.                         outer_model)
  25. {
  26.    MODEL_NAME(MskModulator);
  28.    //  Read model config parms
  29.    OPEN_PARM_BLOCK;
  30.    GET_DOUBLE_PARM(Bit_Durat);
  31.    GET_DOUBLE_PARM(Data_Skew);
  32.    GET_DOUBLE_PARM(Subcar_Misalign);
  33.    GET_DOUBLE_PARM(Phase_Unbal);
  34.    GET_DOUBLE_PARM(Amp_Unbal);
  35.    GET_BOOL_PARM(Shaping_Is_Bipolar);
  37.    //  Connect input and output signals
  38.    I_In_Sig = i_in_sig;
  39.    Q_In_Sig = q_in_sig;
  40.    Cmpx_Out_Sig = cmpx_out_sig;
  41.    Mag_Out_Sig = mag_out_sig;
  42.    Phase_Out_Sig = phase_out_sig;
  44.    MAKE_OUTPUT( Cmpx_Out_Sig );
  45.    MAKE_OUTPUT( Mag_Out_Sig );
  46.    MAKE_OUTPUT( Phase_Out_Sig );
  47.    MAKE_INPUT( I_In_Sig );
  48.    MAKE_INPUT( Q_In_Sig );
  50.    // Set up derived parms
  51.    double phase_unbal_rad = PI * Phase_Unbal / 180.0;
  52.    Pi_Over_Bit_Dur = PI/Bit_Durat;
  53.    Phase_Shift = complex<float>( -sin(phase_unbal_rad),
  54.                                  cos(phase_unbal_rad));
  55. }
  56. //======================================================
  57. MskModulator::~MskModulator( void ){ };
  58. //======================================================
  59. void MskModulator::Initialize(void)
  60. {
  61.    Samps_Out_Cnt = 0;
  62.    Block_Size = I_In_Sig->GetBlockSize();
  63.    Samp_Intvl = I_In_Sig->GetSampIntvl();
  64. }
  65. //======================================================
  66. int MskModulator::Execute(void)
  67. {
  68.    float *i_in_sig_ptr, *q_in_sig_ptr;
  69.    float *phase_out_sig_ptr, *mag_out_sig_ptr;
  70.    float subcar_misalign, amp_unbal, data_skew;
  71.    float work, work1;
  72.    std::complex<float> work2;
  73.    std::complex<float> *cmpx_out_sig_ptr;
  74.    int samps_out_cnt;
  75.    double samp_intvl;
  76.    double pi_over_bit_dur, argument;
  77.    std::complex<float> phase_shift;
  78.    long int_mult;
  79.    int shaping_is_bipolar;
  80.    int block_size;
  81.    int is;
  83.    cmpx_out_sig_ptr = GET_OUTPUT_PTR( Cmpx_Out_Sig );
  84.    phase_out_sig_ptr = GET_OUTPUT_PTR( Phase_Out_Sig );
  85.    mag_out_sig_ptr = GET_OUTPUT_PTR( Mag_Out_Sig );
  86.    i_in_sig_ptr = GET_INPUT_PTR( I_In_Sig );
  87.    q_in_sig_ptr = GET_INPUT_PTR( Q_In_Sig );
  89.    samps_out_cnt = Samps_Out_Cnt;
  90.    samp_intvl = Samp_Intvl;
  91.    subcar_misalign = Subcar_Misalign;
  92.    amp_unbal = Amp_Unbal;
  93.    data_skew = Data_Skew;
  94.    pi_over_bit_dur = Pi_Over_Bit_Dur;
  95.    phase_shift = Phase_Shift;
  96.    shaping_is_bipolar = Shaping_Is_Bipolar;
  98.    block_size = I_In_Sig->GetValidBlockSize();
  99.    Cmpx_Out_Sig->SetValidBlockSize(block_size);
  100.    Mag_Out_Sig->SetValidBlockSize(block_size);
  101.    Phase_Out_Sig->SetValidBlockSize(block_size);
  103.    for (is=0; is<block_size; is++){
  104.       argument = pi_over_bit_dur * samps_out_cnt *
  105.                   samp_intvl;
  106.       int_mult = long(argument/TWO_PI);
  107.       argument -= int_mult*TWO_PI;
  109.       if(shaping_is_bipolar){
  110.          work1 = float((*i_in_sig_ptr) 
  111.             * sin(argument-pi_over_bit_dur*
  112.             subcar_misalign));
  113.          work = float(amp_unbal * (*q_in_sig_ptr) 
  114.             * cos(argument-pi_over_bit_dur
  115.             *(subcar_misalign+data_skew)));
  116.       }
  117.       else{
  118.          work1 = float((*i_in_sig_ptr) 
  119.             * fabs(sin(argument-pi_over_bit_dur*
  120.             subcar_misalign)));
  121.          work = float(amp_unbal * (*q_in_sig_ptr) 
  122.             * fabs(cos(argument-pi_over_bit_dur
  123.             *(subcar_misalign+data_skew))));
  124.       }
  125.       work2 = work * phase_shift;
  126.       *cmpx_out_sig_ptr = work1 + work2;
  127.       *phase_out_sig_ptr = float(180.0*
  128.                std::arg<float>(*cmpx_out_sig_ptr)/PI);
  129.       *mag_out_sig_ptr = 
  130.                std::abs<float>(*cmpx_out_sig_ptr);
  131.       samps_out_cnt++;
  132.       i_in_sig_ptr++;
  133.       q_in_sig_ptr++;
  134.       cmpx_out_sig_ptr++;
  135.       phase_out_sig_ptr++;
  136.       mag_out_sig_ptr++;
  137.    }
  138.    Samps_Out_Cnt = samps_out_cnt;
  139.    return(_MES_AOK);
  140. }