sinwn.cpp

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

/*
说明：以下程序是conio.cpp程序中调用的一个子程序。
就是这个程序使之真正产生了输入信号X(n)。
*/
#include "math.h"
#include "gauss.c"
void sinwn(a,f,ph,m,fs,snr,seed,x,n)
int m,n;
long seed;
double fs,snr,a[],f[],ph[],x[];
{
  int i,k;
  double z,pi,nsr;
  pi=4.0*atan(1.0);
  z=snr/10.0;
  z=pow(10.0,z);
  z=1.0/(2*z);
  nsr=sqrt(z);
  for(i=0;i<m;i++)
  {
    f[i]=2*pi*f[i]/fs;
    ph[i]=ph[i]*pi/180.0;
  }
  for(k=0;k<n;k++)
  {
    x[k]=0.0;
    for(i=0;i<m;i++)
    {
      x[k]=x[k]+a[i]*sin(k*f[i]+ph[i]);
    }
    x[k]=x[k]+nsr*gauss(0.0,1.0,&seed);
  }
}