rayl_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 "rayl_gen.h"
  9. #include "model_graph.h"
  10. #include "syst_graph.h"
  11. #include "misdefs.h"
  12. #include "raylrand.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. template <class T>
  22. RayleighNoiseGenerator<T>::RayleighNoiseGenerator( char* instance_name,
  23.                                               PracSimModel* outer_model,
  24.                                               Signal<T>* noise_sig)
  25.                       :PracSimModel(instance_name,
  26.                                     outer_model)
  27. {
  28.   MODEL_NAME(RayleighNoiseGenerator);
  29.   Noise_Sig = noise_sig;
  30.   OPEN_PARM_BLOCK;
  31.   GET_INT_PARM(Seed);
  32.   GET_FLOAT_PARM(Noise_Sigma);
  33.   //----------------------------------------------
  34.   MAKE_OUTPUT(Noise_Sig);
  35. }
  36. //=============================================
  37. template <class T>
  38. RayleighNoiseGenerator<T>::~RayleighNoiseGenerator( void ){ };
  39. //===========================================
  40. template <class T>
  41. void RayleighNoiseGenerator<T>::Initialize(void)
  42. {
  43.   *DebugFile << "Now in RayleighNoiseGenerator::Initialize()" << endl;
  44.   Proc_Block_Size = Noise_Sig->GetBlockSize();
  45. }
  46. //=============================================
  47. template <class T>
  48. int RayleighNoiseGenerator<T>::Execute(void)
  49. {
  50.   int is;
  51.   float *noise_sig_ptr;
  52.   float noise_sigma;
  53.   float rand_var;
  54.   float noise_sig_val;
  55.   //--------------------------------------------
  56.   noise_sigma = Noise_Sigma;
  57.   long seed = Seed;
  58.   //---------------------------------------------
  59.   noise_sig_ptr = GET_OUTPUT_PTR(Noise_Sig);
  60.   //-------------------------------------------------
  61.   //  main loop
  62.   for(is=0; is<Proc_Block_Size; is++)
  63.     {
  64.     // generate gaussian RV
  65.     RayleighRandom(&seed, &rand_var);
  66.     noise_sig_val = noise_sigma * rand_var;
  67.     *noise_sig_ptr++ = noise_sig_val;
  68.     }// end of main loop
  69.   // put back variables that have changed
  70.   Seed = seed;
  71.   Noise_Sig->SetValidBlockSize(Proc_Block_Size);
  72.   //----------------------------------------------
  73.   return(_MES_AOK);
  74. };
  75. template RayleighNoiseGenerator<float>;