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

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

3G开发

开发平台：

Visual C++

disc_auto_cov.cpp：源码内容
  1. //
  2. //  File = disc_auto_corr.cpp
  3. //
  5. #include <stdlib.h>
  6. #include "parmfile.h"
  7. #include "disc_auto_cov.h"
  8. #include "misdefs.h"
  9. #include "model_graph.h"
  10. #include "sigplot.h"
  11. extern ParmFile* ParmInput;
  12. extern int PassNumber;
  13. #ifdef _DEBUG
  14.   extern ofstream *DebugFile;
  15. #endif
  16. ofstream CovarFile("covar_res.txt", ios::out);
  18. //======================================================
  20. DiscreteAutoCovar::DiscreteAutoCovar( char* instance_name,
  21.                                 PracSimModel* outer_model,
  22.                                 Signal<byte_t>* in_sig)
  23.             :PracSimModel(instance_name,
  24.                           outer_model)
  25. {
  26.   In_Sig = in_sig;
  28.   OPEN_PARM_BLOCK;
  30.   GET_INT_PARM(Num_Corr_Passes);
  31.   GET_DOUBLE_PARM(Mean);
  32.   Mean_Sqrd = Mean * Mean;
  34.   MAKE_INPUT(In_Sig);
  36. }
  37. DiscreteAutoCovar::~DiscreteAutoCovar( void ){ };
  39. void DiscreteAutoCovar::Initialize(void)
  40. {
  41.    Proc_Block_Size = In_Sig->GetBlockSize();
  42.    Samp_Intvl = In_Sig->GetSampIntvl();
  43.    Corr_Buf = new double[Proc_Block_Size/2];
  44.    for(int i=0; i<Proc_Block_Size/2; i++) Corr_Buf[i]=0.0;
  46. //   Max_Corr = 0.0;
  47. //   Max_Corr_Time = 0.0;
  49.    //INITIALIZATION_REPORT(BasicResults);
  50.    return;
  51. }
  53. int DiscreteAutoCovar::Execute()
  54. {
  55.    byte_t *in_sig_beg;
  56.    byte_t *in_sig_ptr;
  57.    byte_t *slide_sig_beg;
  58.    byte_t *slide_sig_ptr;
  59.    int i, j;
  60.    int sum;
  61.    double covar_val;
  63.    if(PassNumber > Num_Corr_Passes) return(_MES_AOK);
  65.    //-------------------------------------------------------
  66.    //  Copy frequently accessed member vars into local vars
  67.    int proc_block_size = Proc_Block_Size;
  68.    int slide_block_size = proc_block_size/2;
  69.    double samp_intvl = Samp_Intvl;
  71.    //----------------------------------------
  72.    // Get pointers for input and output
  73.    in_sig_beg = GET_INPUT_PTR(In_Sig);
  74.    slide_sig_beg = in_sig_beg + proc_block_size/4;
  76.    //----------------------------------------
  78.    for( j=0; j<slide_block_size; j++)
  79.    {
  80.       sum = 0;
  81.       in_sig_ptr = in_sig_beg;
  82.       slide_sig_ptr = slide_sig_beg;
  83.       for( i=0; i<slide_block_size; i++)
  84.       {
  85.          sum += (*in_sig_ptr )*(*slide_sig_ptr );
  86.          slide_sig_ptr++;
  87.          in_sig_ptr++;
  88.       }
  89.       Corr_Buf[j] += ((float(sum)/slide_block_size)- Mean_Sqrd);
  90.       in_sig_beg++;
  91.    }
  92.    if(PassNumber == Num_Corr_Passes)
  93.    {
  94.       for( j=0; j<slide_block_size; j++)
  95.       {
  96.          covar_val = fabs((Corr_Buf[j]/Num_Corr_Passes));
  97.          if(covar_val <1e-20) covar_val = 1e-20;
  98.          CovarFile << j-slide_block_size/2 << ", " << 10.0*log10(covar_val) << endl;
  99.          //CovarFile << j-slide_block_size/2 << ", " << Corr_Buf[j] << endl;
  100.       }
  101.       CovarFile.close();
  102.       exit(22);
  103.    }
  104.    //---------------------
  105.    return(_MES_AOK);
  106. }