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

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

3G开发

开发平台：

Visual C++

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