gauss.cpp

DSP中输入信号的生成过程。 conio.cpp实现X(n)信号

/*
以下程序是sinwn.cpp程序中调用的一个子程序。
用来产生高斯白噪声。
*/
#include "uniform.c"
double gauss(mean,sigma,s)
double mean,sigma;
long int *s;
{
  int i;
  double x,y;
  double uniform();
  for(x=0,i=0;i<12;i++)
    x+=uniform(0.0,1.0,s);
  x=x-6.0;
  y=mean+x*sigma;
  return(y);
}

double uniform(a,b,seed)
double a,b;
long int *seed;
{ double t;
  *seed=2045*(*seed)+1;
  *seed=*seed-(*seed/1048576)*1048576;
  t=(*seed)/1048576.0;
  t=a+(b-a)*t;
  return(t);
}