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

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

3G开发

开发平台：

Visual C++

sines_in_awgn.cpp：源码内容
  1. //
  2. //  File = sines_in_awgn.cpp
  3. //
  5. #include <stdlib.h>
  6. //#include <fstream>
  7. #include "sines_in_awgn.h"
  8. #include "misdefs.h"
  9. #include "parmfile.h"
  10. #include "model_graph.h"
  11. #include "sigstuff.h"
  12. #include "syst_graph.h"
  13. #include "gausrand.h"
  15. extern ParmFile *ParmInput;
  16. //extern ActiveSystemGraph ActSystGraph;
  17. #ifdef _DEBUG
  18.   extern ofstream *DebugFile;
  19. #endif
  22. //======================================================
  24. SinesInAwgn::SinesInAwgn( char* instance_name,
  25.                     PracSimModel* outer_model,
  26.                     Signal<float>* out_sig )
  27.         :PracSimModel( instance_name,
  28.                         outer_model )
  29. {
  30.   MODEL_NAME(SinesInAwgn);
  31.   Out_Sig = out_sig;
  32.   OPEN_PARM_BLOCK;
  34.    GET_DOUBLE_PARM(Awgn_Sigma);
  35.    GET_LONG_PARM(Noise_Seed);
  36.   GET_INT_PARM(Num_Sines);
  38.   Phase_In_Deg = new double[Num_Sines];
  39.   Phase_In_Rad = new double[Num_Sines];
  40.   Arg_Increm = new double[Num_Sines];
  41.   Base_Arg = new double[Num_Sines];
  42.   Freq_In_Hz = new double[Num_Sines];
  43.   Tone_Gain = new double[Num_Sines];
  45.   GET_DOUBLE_PARM_ARRAY(Phase_In_Deg,Num_Sines);
  46.   for(int idx=0; idx<Num_Sines; idx++)
  47.     {
  48.     Phase_In_Rad[idx] = PI * Phase_In_Deg[idx] /180.0;
  49.     cout << "Phase_In_Rad[" << idx << "] = " << Phase_In_Rad[idx] << endl;
  50.     }
  51.   GET_DOUBLE_PARM_ARRAY(Freq_In_Hz, Num_Sines);
  52.   GET_DOUBLE_PARM_ARRAY(Tone_Gain, Num_Sines);
  53.   MAKE_OUTPUT(Out_Sig);
  54. }
  55. //====================================================
  56. SinesInAwgn::~SinesInAwgn( void){};
  57. //====================================================
  58. void SinesInAwgn::Initialize(void)
  59. {
  60.   *DebugFile << "Now in SinesInAwgn::Initialize()" << endl;
  62.   Block_Size = Out_Sig->GetBlockSize();
  63.   Samp_Intvl = Out_Sig->GetSampIntvl();
  65.   for(int isin=0; isin<Num_Sines; isin++)
  66.     {
  67.     Arg_Increm[isin] = TWO_PI * Freq_In_Hz[isin] * Samp_Intvl;
  68.     Base_Arg[isin] = Phase_In_Rad[isin];
  69.     }
  70. }
  71. //====================================================
  72. int SinesInAwgn::Execute()
  73. {
  74.   int is, isin;
  75.   float *out_sig_ptr;
  76.   double out_samp;
  77.   double *base_arg, *arg_increm;
  78.   double *tone_gain;
  79.   long noise_seed;
  80.   double awgn_sigma;
  82.   //--------------------------------
  83.   //  Get pointer for output buffer
  85.   out_sig_ptr = GET_OUTPUT_PTR(Out_Sig);
  86.   base_arg = Base_Arg;
  87.   arg_increm = Arg_Increm;
  88.   tone_gain = Tone_Gain;
  89.   noise_seed = Noise_Seed;
  90.   awgn_sigma = Awgn_Sigma;
  92.   for (is=0; is<Block_Size; is++)
  93.     {
  94.     GaussRandom(&noise_seed, &out_samp);
  95.     out_samp *= awgn_sigma;
  96.     for(isin=0; isin<Num_Sines; isin++)
  97.       {
  98.       out_samp += tone_gain[isin] * sin(base_arg[isin] + (is*arg_increm[isin]));
  99.       }
  100.     if(is < Block_Size/2)
  101.       {
  102.       //*out_sig_ptr++ = 2*is*out_samp/Block_Size;
  103.       *out_sig_ptr++ = out_samp;
  104.       }
  105.     else
  106.       {
  107.       //*out_sig_ptr++ = out_samp * (1.0- 2*(is-Block_Size/2.0)/double(Block_Size));
  108.       *out_sig_ptr++ = out_samp;
  109.       }
  110.     }
  111.   for(isin=0; isin<Num_Sines; isin++)
  112.     {
  113.     Base_Arg[isin] = fmod(Base_Arg[isin] + Block_Size*Arg_Increm[isin],TWO_PI);
  114.     }
  115.   Noise_Seed = noise_seed;
  116.   return(_MES_AOK);
  117. }