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

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

3G开发

开发平台：

Visual C++

nonlinear_amp.cpp：源码内容
  1. //
  2. //  File = nonlinear_amp.cpp
  3. //
  5. #include <stdlib.h>
  6. #include <fstream>
  7. #include <strstream>
  8. #include "parmfile.h"
  9. #include "nonlinear_amp.h"
  10. #include "model_graph.h"
  11. extern ParmFile *ParmInput;
  13. //======================================================
  15. NonlinearAmplifier::NonlinearAmplifier( 
  16.                      char* instance_name,
  17.                      PracSimModel* outer_model,
  18.                      Signal< complex<float> >* in_sig,
  19.                      Signal< complex<float> >* out_sig )
  20.              :PracSimModel(instance_name,
  21.                            outer_model)
  22. {
  23.   MODEL_NAME(NonlinearAmplifier);
  25.   //  Read model config parms
  27.    OPEN_PARM_BLOCK;
  28.    GET_DOUBLE_PARM(Output_Power_Scale_Factor);
  29.    GET_DOUBLE_PARM(Phase_Scale_Factor);
  30.    GET_DOUBLE_PARM(Anticipated_Input_Power);
  31.    GET_DOUBLE_PARM(Operating_Point);
  32.    GET_DOUBLE_PARM(Agc_Time_Constant);
  34.    Input_Power_Scale_Factor = 
  35.          float(Operating_Point/Anticipated_Input_Power);
  37.    Am_Am_Fname = new char[64];
  38.    strcpy(Am_Am_Fname, "");
  39.    GET_STRING_PARM(Am_Am_Fname);
  41.    Am_Pm_Fname = new char[64];
  42.    strcpy(Am_Pm_Fname, "");
  43.    GET_STRING_PARM(Am_Pm_Fname);
  45.   //  Connect input and output signals
  47.   In_Sig = in_sig;
  48.   Out_Sig = out_sig;
  49.   MAKE_OUTPUT( Out_Sig );
  50.   MAKE_INPUT( In_Sig );
  52.   Am_Am_Curve = new SampledCurve(Am_Am_Fname);
  53.   Am_Pm_Curve = new SampledCurve(Am_Pm_Fname);
  54.   }
  56. //======================================================
  57. void NonlinearAmplifier::Initialize(void)
  58. {
  59.   double samp_intvl = Out_Sig->GetSampIntvl();
  60.   Out_Avg_Block_Size = Out_Sig->GetBlockSize();
  61. }
  62. //======================================================
  63. NonlinearAmplifier::~NonlinearAmplifier( void ){ };
  65. //======================================================
  66. int NonlinearAmplifier::Execute()
  67. {
  68.    complex<float>  *out_sig_ptr, out_sig;
  69.    complex<float>  *in_sig_ptr, in_sig;
  70.    complex<float>  agc_in_sig;
  71.    float power, power_out;
  72.    float input_phase;
  73.    double phase_shift;
  74.    double amplitude;
  75.    double phase_out;
  76.    double sum_in, sum_out;
  77.    double amp_sqrd;
  78.    double avg_power_in, avg_power_out;
  79.    int block_size, is;
  81.    block_size = In_Sig->GetValidBlockSize();
  82.    Out_Sig->SetValidBlockSize(block_size);
  83.    //-------------------------------------------------
  85.    out_sig_ptr = GET_OUTPUT_PTR( Out_Sig );
  86.    in_sig_ptr = GET_INPUT_PTR( In_Sig );      
  87.    block_size = In_Sig->GetValidBlockSize();
  88.    Out_Sig->SetValidBlockSize(block_size);
  89.    sum_in = 0.0;
  90.    sum_out = 0.0;
  92.    for(is=0; is<block_size; is++){
  93.       in_sig = *in_sig_ptr++;
  94.       agc_in_sig = in_sig;
  95.       amp_sqrd = std::norm(agc_in_sig);
  96.       sum_in += amp_sqrd;
  97.       power = float(0.5 * Input_Power_Scale_Factor *
  98.                      amp_sqrd);
  99.       if(power == 0.0){
  100.          input_phase = 0.0;
  101.       }
  102.       else{
  103.          input_phase = std::arg(agc_in_sig);
  104.       }
  105.       power_out = float(Output_Power_Scale_Factor *
  106.                      Am_Am_Curve->GetValue(power));
  107.       sum_out += power_out;
  108.       phase_shift = Am_Pm_Curve->GetValue(power);
  109.       amplitude = sqrt(2.0*power_out);
  111.       phase_out = input_phase + 
  112.          Phase_Scale_Factor*phase_shift;
  114.       out_sig = complex<float>(  
  115.             float(amplitude*cos(phase_out)),
  116.             float(amplitude*sin(phase_out)));
  117.       *out_sig_ptr++ = out_sig;
  118.    }
  119.    avg_power_out = sum_out/block_size;
  120.    avg_power_in = sum_in/block_size/2;
  121.    BasicResults 
  122.       << "In NonlinearAmplifier avg input power = " 
  123.       << avg_power_in << endl;
  124.    BasicResults 
  125.       << "                      avg output power = " 
  126.       << avg_power_out << endl;
  128.    return(_MES_AOK);
  129. }