sinwn.txt

来自「一些信号分析的C 语言程序」· 文本 代码 · 共 24 行

#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);
}
}