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

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

3G开发

开发平台：

Visual C++

nrz_theory.cpp：源码内容
  1. //
  2. //  File = nrz_theory.cpp
  3. //
  5. #include <stdlib.h>
  6. #include <math.h>
  7. #include "nrz_theory.h"
  8. #include "sinc.h"
  9. #include <fstream>
  10. #ifndef PI
  11.   #define PI 3.1415926535897932
  12. #endif
  13. using namespace std;
  15. void NrzPsd(  double bit_energy,
  16.               double bit_rate,
  17.               double beg_freq,
  18.               double end_freq,
  19.               int num_pts,
  20.               bool plot_in_db,
  21.               char* out_filename)
  22. {
  23.   double freq, val;
  24.   ofstream *out_file;
  25.   out_file = new ofstream(out_filename, ios::out);
  27.   double delta_freq = (end_freq - beg_freq)/double(num_pts-1);
  28.   for(int n=0; n<num_pts; n++)
  29.     {
  30.     freq = beg_freq + n * delta_freq;
  31.     val = bit_energy * sinc_sqrd(freq/bit_rate);
  32.     if( plot_in_db )
  33.       {
  34.       if(val>0.0)
  35.         {
  36.         val = 10.0 * log10(val);
  37.         }
  38.       else
  39.         {
  40.         val = -200.0;
  41.         }
  42.       }
  43.     (*out_file) << freq << ", " << val << endl;
  44.     }
  45.   delete out_file;
  46. }
  47. //================================================================
  48. void BiphasePsd(  double bit_energy,
  49.               double bit_rate,
  50.               double beg_freq,
  51.               double end_freq,
  52.               int num_pts,
  53.               bool plot_in_db,
  54.               char* out_filename)
  55. {
  56.   double freq, val, argument;
  57.   ofstream *out_file;
  58.   out_file = new ofstream(out_filename, ios::out);
  60.   double delta_freq = (end_freq - beg_freq)/double(num_pts-1);
  61.   for(int n=0; n<num_pts; n++)
  62.     {
  63.     freq = beg_freq + n * delta_freq;
  64.     argument = freq/bit_rate/2.0;
  65.     val = sin(PI*argument);
  66.     val *= val;
  67.     val *= bit_energy * sinc_sqrd(argument);
  68.     if( plot_in_db )
  69.       {
  70.       if(val>0.0)
  71.         {
  72.         val = 10.0 * log10(val);
  73.         }
  74.       else
  75.         {
  76.         val = -200.0;
  77.         }
  78.       }
  79.     (*out_file) << freq << ", " << val << endl;
  80.     }
  81.   delete out_file;
  82. }
  84. //================================================================
  85. void MillerPsd(  double bit_energy,
  86.               double bit_rate,
  87.               double beg_freq,
  88.               double end_freq,
  89.               int num_pts,
  90.               bool plot_in_db,
  91.               char* out_filename)
  92. {
  93.   double freq, val, phi;
  94.   ofstream *out_file;
  95.   out_file = new ofstream(out_filename, ios::out);
  97.   double delta_freq = (end_freq - beg_freq)/double(num_pts-1);
  98.   for(int n=0; n<num_pts; n++)
  99.     {
  100.     freq = beg_freq + n * delta_freq;
  101.     phi = PI*freq/bit_rate;
  102.     val = 23 - 2.0*cos(phi) - 22.0*cos(2.0*phi) - 12.0*cos(3.0*phi) + 5.0*cos(4.0*phi);
  103.     val += 12.0*cos(5.0*phi) + 2.0*cos(6.0*phi) - 8.0*cos(7.0*phi) + 2.0*cos(8.0*phi);
  104.     val /= phi * phi * (34.0 + 16.0*cos(8.0*phi));
  105.     val *= bit_energy;
  106.     if( plot_in_db )
  107.       {
  108.       if(val>0.0)
  109.         {
  110.         val = 10.0 * log10(val);
  111.         }
  112.       else
  113.         {
  114.         val = -200.0;
  115.         }
  116.       }
  117.     (*out_file) << freq << ", " << val << endl;
  118.     }
  119.   delete out_file;
  120. }
  122. //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++