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

3G开发

开发平台:

Visual C++

  1. //
  2. //  File = noise_gen.cpp
  3. //
  4. #include <stdlib.h>
  5. #include <fstream>
  6. #include <math.h>
  7. #include "parmfile.h"
  8. #include "noise_gen.h"
  9. #include "model_graph.h"
  10. #include "syst_graph.h"
  11. #include "misdefs.h"
  12. #include "gausrand.h"
  13. extern int PassNumber;
  14. extern ParmFile *ParmInput;
  15. extern SystemGraph CommSystemGraph;
  16. #ifdef _DEBUG
  17.   extern ofstream *DebugFile;
  18. #endif
  19. //======================================================
  20. // normal constructor
  21. GaussianNoiseGenerator::GaussianNoiseGenerator( char* instance_name,
  22.                                               PracSimModel* outer_model,
  23.                                               Signal<float>* noise_sig)
  24.                       :PracSimModel(instance_name,
  25.                                     outer_model)
  26. {
  27.   MODEL_NAME(GaussianNoiseGenerator);
  28.   Noise_Sig = noise_sig;
  29.   OPEN_PARM_BLOCK;
  30.   GET_INT_PARM(Seed);
  31.   GET_DOUBLE_PARM(Noise_Sigma);
  32.   //----------------------------------------------
  33.   MAKE_OUTPUT(Noise_Sig);
  34. }
  35. //=============================================
  36. GaussianNoiseGenerator::~GaussianNoiseGenerator( void ){ };
  37. //===========================================
  38. void GaussianNoiseGenerator::Initialize(void)
  39. {
  40.   *DebugFile << "Now in GaussianNoiseGenerator::Initialize()" << endl;
  41.   Proc_Block_Size = Noise_Sig->GetBlockSize();
  42. }
  43. //=============================================
  44. int GaussianNoiseGenerator::Execute(void)
  45. {
  46.   int is;
  47.   float *noise_sig_ptr;
  48.   float noise_sigma;
  49.   float rand_var;
  50.   float noise_sig_val;
  51.   //--------------------------------------------
  52.   noise_sigma = Noise_Sigma;
  53.   long seed = Seed;
  54.   //---------------------------------------------
  55.   noise_sig_ptr = GET_OUTPUT_PTR(Noise_Sig);
  56.   //-------------------------------------------------
  57.   //  main loop
  58.   for(is=0; is<Proc_Block_Size; is++)
  59.     {
  60.     // generate gaussian RV
  61.     GaussRandom(&seed, &rand_var);
  62.     noise_sig_val = noise_sigma * rand_var;
  63.     *noise_sig_ptr++ = noise_sig_val;
  64.     }// end of main loop
  65.   // put back variables that have changed
  66.   Seed = seed;
  67.   //----------------------------------------------
  68.   return(_MES_AOK);
  69. };