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

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

3G开发

开发平台：

Visual C++

gausrand.cpp：源码内容
  1. //
  2. //  File = gausrand.cpp
  3. //
  5. #include <math.h>
  6. #include "uni_rand.h"
  7. #include "gausrand.h"
  8.  
  9. void GaussRandom(long *seed, float *result)
  10. {
  11.  double x1, x2, g1; 
  12.  double XA, XB, big_S;
  13.  double radical;
  15.  big_S=1.0;
  16.  while(big_S >= 1.0){
  17.    XA = DoubleUniformRandom(seed);
  18.    XB = DoubleUniformRandom(seed);
  19.    x1 = 1.0 - 2.0 * XA;
  20.    x2 = 1.0 - 2.0 * XB;
  21.  big_S = x1*x1 + x2*x2;
  22.  }
  23.  radical = sqrt(-2.0*log(big_S)/big_S);
  24.  g1 = x1 * radical;
  25.  *result = float(g1);
  26. }
  27. //------------------------------------------------
  28. void GaussRandom(long *seed, double *result)
  29. {
  30.  double x1, x2, g1; 
  31.  double XA, XB, big_S;
  32.  double radical;
  34.  big_S=1.0;
  35.  while(big_S >= 1.0){
  36.    XA = DoubleUniformRandom(seed);
  37.    XB = DoubleUniformRandom(seed);
  38.    x1 = 1.0 - 2.0 * XA;
  39.    x2 = 1.0 - 2.0 * XB;
  40.    big_S = x1*x1 + x2*x2;
  41.  }
  42.  radical = sqrt(-2.0*log(big_S)/big_S);
  43.  g1 = x1 * radical;
  44.  *result = g1;
  45. }
  46. //------------------------------------------------
  47. void GaussRandom(long *seed, std::complex<double> *result)
  48. {
  49.  double x1, x2, g1, g2; 
  50.  double XA, XB, big_S;
  51.  double radical;
  53.  big_S=1.0;
  54.  while(big_S >= 1.0){
  55.    XA = DoubleUniformRandom(seed);
  56.    XB = DoubleUniformRandom(seed);
  57.    x1 = 1.0 - 2.0 * XA;
  58.    x2 = 1.0 - 2.0 * XB;
  59.    big_S = x1*x1 + x2*x2;
  60.  }
  61.  radical = sqrt(-2.0*log(big_S)/big_S);
  62.  g1 = x1 * radical;
  63.  g2 = x2 * radical;
  64.  *result = std::complex<double>(g1, g2);
  65. }
  66. //------------------------------------------------
  67. void GaussRandom(long *seed, std::complex<float> *result)
  68. {
  69.  double x1, x2;
  70.  float g1, g2; 
  71.  double XA, XB, big_S;
  72.  double radical;
  74.  big_S=1.0;
  75.  while(big_S >= 1.0){
  76.    XA = DoubleUniformRandom(seed);
  77.    XB = DoubleUniformRandom(seed);
  78.    x1 = 1.0 - 2.0 * XA;
  79.    x2 = 1.0 - 2.0 * XB;
  80.    big_S = x1*x1 + x2*x2;
  81.  }
  82.  radical = sqrt(-2.0*log(big_S)/big_S);
  83.  g1 = float(x1 * radical);
  84.  g2 = float(x2 * radical);
  85.  *result = std::complex<float>(g1, g2);
  86. }