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

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

3G开发

开发平台：

Visual C++

ar_proc_gen.cpp：源码内容
  1. //
  2. //  File = ar_proc_gen.cpp
  3. //
  5. #include <stdlib.h>
  6. #include <fstream>
  7. #include <math.h>
  8. #include "parmfile.h"
  9. #include "ar_proc_gen.h"
  10. #include "model_graph.h"
  11. #include "syst_graph.h"
  12. #include "misdefs.h"
  13. #include "gausrand.h"
  14. extern int PassNumber;
  15. extern ParmFile *ParmInput;
  16. extern SystemGraph CommSystemGraph;
  17. #ifdef _DEBUG
  18.   extern ofstream *DebugFile;
  19. #endif
  21. //======================================================
  22. // normal constructor
  24. template <class T>
  25. ArProcessGenerator<T>::ArProcessGenerator( char* instance_name,
  26.                                               PracSimModel* outer_model,
  27.                                               Signal<T>* noise_sig)
  28.                       :PracSimModel(instance_name,
  29.                                     outer_model)
  30. {
  31.   MODEL_NAME(ArProcessGenerator);
  33.   Noise_Sig = noise_sig;
  35.   OPEN_PARM_BLOCK;
  37.   GET_INT_PARM(Noise_Seed);
  38.   GET_DOUBLE_PARM(Driving_Variance);
  39.   Driving_Sigma = float(sqrt(Driving_Variance));
  40.   GET_INT_PARM(Ar_Order);
  41.   int ar_size = Ar_Order+1;
  42.   A_Coeffs = new double[ar_size];
  43.   GET_DOUBLE_PARM_ARRAY(A_Coeffs, ar_size);
  45.   //----------------------------------------------
  46.   MAKE_OUTPUT(Noise_Sig);
  47.   Old_Output = new T[Ar_Order];
  48.   for(int i=0; i<Ar_Order; i++) Old_Output[i] = 0.0;
  50. }
  51. //=============================================
  52. template <class T>
  53. ArProcessGenerator<T>::~ArProcessGenerator( void ){ };
  54. //===========================================
  55. template <class T>
  56. void ArProcessGenerator<T>::Initialize(void)
  57. {
  58.   *DebugFile << "Now in ArProcessGenerator::Initialize()" << endl;
  60.   Proc_Block_Size = Noise_Sig->GetBlockSize();
  61. }
  62. //=============================================
  63. template <class T>
  64. int ArProcessGenerator<T>::Execute(void)
  65. {
  66.    int is, idx;
  67.    float *noise_sig_ptr;
  68.    float driving_sigma;
  69.    float rand_var;
  70.    float noise_sig_val;
  72.    //--------------------------------------------
  73.    driving_sigma = Driving_Sigma;
  74.    long noise_seed = Noise_Seed;
  76.    //---------------------------------------------
  78.    noise_sig_ptr = GET_OUTPUT_PTR(Noise_Sig);
  80.    //-------------------------------------------------
  81.    //  main loop
  83.    for(is=0; is<Proc_Block_Size; is++)
  84.    {
  85.       // generate gaussian RV
  86.       GaussRandom(&noise_seed, &rand_var);
  87.       noise_sig_val = driving_sigma * rand_var;
  89.       for(idx=0; idx<Ar_Order; idx++)
  90.       {
  91.          noise_sig_val -= A_Coeffs[idx+1]*Old_Output[idx];
  92.       }
  93.       for(idx=Ar_Order-1; idx>0; idx--)
  94.       {
  95.          Old_Output[idx] = Old_Output[idx-1];
  96.       }
  98.       Old_Output[0] = noise_sig_val;
  99.       *noise_sig_ptr++ = noise_sig_val;
  100.    }// end of main loop
  102.    // put back variables that have changed
  103.    Noise_Seed = noise_seed;
  105.    //----------------------------------------------
  106.    Noise_Sig->SetValidBlockSize(Proc_Block_Size);
  107.    return(_MES_AOK);
  108. };
  110. template ArProcessGenerator<float>;