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

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

3G开发

开发平台：

Visual C++

dump_spect.cpp：源码内容
  1. //
  2. //  File = dump_spect.cpp
  3. //
  5. #include <stdlib.h>
  6. #include "dump_spect.h"
  7. #include "math.h"
  9. #ifdef _DEBUG
  10.   extern ofstream *DebugFile;
  11. #endif
  13. void DumpSpectrum(  double *psd_est,
  14.                int nsamps_in_psd,
  15.                double delta_f,
  16.                double freq_norm_factor,
  17.                bool output_in_decibels,
  18.                bool plot_two_sided,
  19.                ofstream *psd_file)
  20. {
  21.    int is;
  22.    //------------------------------------
  23.    // find peak
  25.    double peak_val=0.0;
  26.    int peak_idx=0;
  27.    for(is=0; is<nsamps_in_psd/2; is++)
  28.    {
  29.       if(psd_est[is] <= peak_val) continue;
  30.       //else
  31.       peak_val = psd_est[is];
  32.       peak_idx = is;
  33.    }
  35.    if(output_in_decibels)
  36.    {
  37.       //double offset = 10.0*log10(Psd_Est[peak_idx]);
  38.       double offset = 0.0;
  40.       if(plot_two_sided)
  41.       {
  42.          for(is=-(nsamps_in_psd/2-1); is<0; is++)
  43.          {
  44.             if( psd_est[-is] >0 )
  45.             {
  46.                (*psd_file) << is * delta_f * freq_norm_factor << ", " 
  47.                            << (10.0*log10(psd_est[-is])) << endl;
  48.                           // << (10.0*log10(Psd_Est[-is])-offset) << endl;
  49.             }
  50.             else
  51.             {
  52.                (*psd_file) << is * delta_f * freq_norm_factor << ", -200.0" << endl;
  53.             }
  54.          }
  55.       }
  57.       for(is=0; is<nsamps_in_psd/2; is++)
  58.       {
  59.          if( psd_est[is] >0 )
  60.          {
  61.             (*psd_file) << is * delta_f * freq_norm_factor << ", " 
  62.                         << (10.0*log10(psd_est[is])) << endl;
  63.                        // << (10.0*log10(Psd_Est[is])-offset) << endl;
  64.             psd_est[is] = 0.0;
  65.          }
  66.          else
  67.          {
  68.             (*psd_file) << is * delta_f * freq_norm_factor << ", -200.0" << endl;
  69.          }
  70.       }
  71.    }
  72.    else
  73.    {
  74.       // plot as linear ordinate
  75.       if(plot_two_sided)
  76.       {
  77.          for(is=-(nsamps_in_psd/2-1); is<0; is++)
  78.          {
  79.             (*psd_file) << is * delta_f * freq_norm_factor << ", " 
  80.                         << psd_est[-is] << endl;
  81.          }
  82.       }
  84.       for(is=0; is<nsamps_in_psd/2; is++)
  85.       {
  86.          (*psd_file) << is * delta_f * freq_norm_factor << ", " 
  87.                      << psd_est[is] << endl;
  88.          psd_est[is] = 0.0;
  89.       }
  90.    }
  91. }