求C或C++的白噪声产生程序!

要满足白噪声的各种指标哟。 问题点数：50、回复次数：6Top

1 楼wqs6（竹山）回复于 2004-07-01 12:35:33 得分 50

#include    
  #include    
  #include    
  #include   "msp.h"    
  float   randnu(long   *iseed)    
  {    
  float   z;    
  *iseed=2045*(*iseed)+1;    
  *iseed=*iseed-(*iseed/1048576)*1048576;    
  z=(float)((*iseed+1)/1048577.0);    
  return(z);    
  }    
  /*--------------------------------------------------------------------*/    
  void   meavar(float   u[],int   *n,float   *pum,float   *puv)    
  {    
  int   i,k;    
  *pum=0.0;    
  for(k=0;k<*n;k++)    
  *pum=*pum+u[k];    
  *pum=*pum/(*n);    
  *puv=0.0;    
  for(i=0;i<*n;i++)    
  *puv=*puv+pow((u[i]-*pum),2);    
  *puv=*puv+pow((u[i]-*pum),2);    
  *puv=*puv/(*n-1.);    
  return;    
  }    
  /*---------------------------------------------------------------------    
  Routine   mrandom   :   To   generate   the   random   number(pseudo-white   noise).    
  input   Parameters:    
  n   :   the   random   data   number   requested;   integer   .    
  iseed:   the   seed   for   pseudo-random   data   generation.it   must   be    
  initialized   by   main   program(suggested   value   is   ISEED=12357),    
  and   the   random   number   is   cycled,the   cycle   length=1,048,576    
  itype:   random   data   distribution   type,   see   below:    
  itype=1:   Uniform   distributed,from   0.0   to   1.0    
  itype=2:   Uniform   distributed,Mean=0.0,   Variance(Power)   p=1.0    
  itype=3:   Uniform   distributed,Mean=0.0,   Variance(Power)   p=p.    
  itype=4:   Gaussian   distributed,Mean=0.0,   Variance(Power)   p=1.0    
  itype=5:   Gaussian   distributed,Mean=0.0,   Variance(Power)   p=p.    
  p   :variance(Power)   of   random,   only   used   when   itype=3   or   itype=5.    
  out   parameters:    
  u   :n   dimensioned   real   array,   data   is   stored   in   u(0)   to   u(n-1).    
  in   Chapter   1    
  ---------------------------------------------------------------------*/    
  void   mrandom(float   u[],int   *n,long   *piseed,int   itype,float   p)    
  {    
  int   k,ns,ksection,ks,j;    
  float   a,v,umean,uvari;    
  float   *pum,*puv;    
  pum=&umean;    
  puv=&uvari;    
  if(itype   >6   |r   |r   itype   <1)    
  return;    
  for(k=0;k<*n;k++)    
  u[k]=randnu(piseed);    
  if(itype==2   |r   |r   itype==3)    
  {    
  meavar(u,n,pum,puv);    
  /*   to   obtain   a   zero   mean   and   P-power   random   sequence   u[k]).*/    
  a=12.;    
  if(itype==2)    
  p=sqrt(a);    
  if(itype==3)    
  p=sqrt(p*a);    
  for(k=0;k<*n;k++)    
  u[k]=(u[k]-umean)*p;    
   
  }    
  if(itype==4   |r   |ritype==5)    
  {    
  /*   to   generate   the   Gaussian   randow   sequence   u[k],k=0,1,2,...,ns-1*/    
  ksection=12;    
  ns=*n/ksection;    
  ks=0;    
  if   (itype==4)   p=1;    
  p=sqrt(p);    
  for(k=0;k   {v=0.0;    
  for(j=0;j   v+=p*(u[j+ks]-.5);    
  u[k]=v;    
  ks=ks+ksection;    
  }    
  *n=ns;    
  }    
  meavar(u,n,pum,puv);    
  printf("   The   mean   of   u[n]=%f\n",umean);    
  printf("   The   variance   of   u[n]=%f\n",uvari);    
  return;    
  }    
   
   
  其中msp.h头文件：    
  #define   abs_error   1.e-10    
   
  #ifndef   _MSP_H_    
  #define   _MSP_H_    
   
  typedef   struct   {float   real,imag;}   complex;    
  /*-------------------------------------------------------------------*/    
  float   mabs(complex   a)    
  {    
  float   m;    
  m=a.real*a.real+a.imag*a.imag;    
  m=sqrt(m);    
  return(m);    
  }    
  /*-------------------------------------------------------------------*/    
  float   msign(float   a,float   b)    
  {    
  float   z;    
  if(b>=0)   z=sqrt(pow(a,2));    
  else   z=-sqrt(pow(a,2));    
  return(z);    
  }    
  /*-------------------------------------------------------------------*/    
  complex   cexp(complex   a)    
  {    
  complex   z;    
  z.real=exp(a.real)*cos(a.imag);    
  z.imag=exp(a.real)*sin(a.imag);    
  return(z);    
  }    
  /*-------------------------------------------------------------------*/    
  void   maftodf(float   d[],float   c[],int   ln,int   iband,float   fln,float   fhn,    
  float   b[],float   a[],int   *ierror);    
  void   mampres(complex   h[],float   amp[],int   n,float   fs,int   iamp,char   filename[]);    
  void   mar1psd(complex   a[],int   ip,int   mfre,float   *ep,float   ts);    
  void   marburg(complex   x[],complex   a[],complex   ef[],complex   eb[],    
  int   n,int   ip,float   *ep,int   *ierror);    
  void   marmach(complex   x[],complex   ef[],int   n,complex   a[],    
  complex   b[],int   ip,int   iq,int   m,float   *ep,float   ts);    
  void   maryuwa(complex   x[],complex   a[],complex   r[],int   n,int   ip,    
  float   *ep,   int   *ierror);    
  void   mbiline(float   d[],float   c[],int   ln,float   b[],float   a[],int   *ierror);    
  void   mbutwcf(int   l,int   k,int   ln,float   d[],float   c[],int   *ierror);    
  void   mchebsh(int   l,int   k,int   ln,float   d[],float   c[],float   phi2,    
  int   *ierror);    
  void   mcholsk(complex   a[],complex   b[],int   n,float   eps,int   *iflag);    
  void   mcmpdft(complex   x[],complex   y[],int   n,int   isign);    
  void   mcmpfft(complex   x[],int   n,int   isign);    
  void   mconvo1(float   x[],float   h[],float   y[],int   n,int   m,int   L);    
  void   mconvo2(complex   x[],complex   h[],complex   y[],int   n1,int   n2,int   n);    
  void   mcorpsd(complex   x[],complex   r[],int   n,int   lag,int   iwindow,float   t);    
  void   mcorre1(complex   x[],complex   y[],complex   r[],int   n,int   lag);    
  void   mcorre2(complex   x[],complex   y[],int   m,int   n,int   icorre);    
  void   mcztfft(complex   x[],int   n,int   m,int   maxnm,float   dltomg,    
  float   omg0,float   fs,int   *ierror);    
  void   mdecint(float   x[],float   h[],float   y[],int   nh,int   ny,int   m,    
  int   l,int   *k);    
  void   mdefir1(int   l,int   iband,float   fl,float   fh,float   fs,int   iwindow,    
  float   b[],float   w[],int   *ierror);    
  void   mdefir2(int   l,int   iband,float   fl,float   fh,complex   b[],    
  float   trans,float   fs,int   *ierror);    
  void   mdefir3(int   nfilt,int   nbands,float   edge[],float   fx[],    
  float   wtx[],float   h[]);    
  void   mdesiir(float   *f1,float   *f2,float   *f3,float   *f4,float   fs,    
  float   alpha1,float   alpha2,int   iband,int   itype);    
  void   mfirres(float   b[],int   lb,int   n,complex   h[]);    
  void   mfitout(float   b[],float   a[],int   lb,int   la,float   x[],    
  int   n,float   y[]);    
  void   miirres(float   a[],float   b[],int   lb,int   la,complex   h[],int   n);    
  void   mlattic(float   b[],float   a[],int   l,float   k[],    
  float   c[],int   itype   ,int   *ierror);    
  void   mmayuwa(complex   x[],int   n,complex   a[],int   ip,complex   b[],int   iq,    
  complex   r[],float   *ep,   float   ts,int   *ierror);    
  void   mmvseps(complex   x[],complex   ef[],complex   eb[],int   n,complex   a[],    
  int   ip,int   *ierror,float   ts);    
  void   morderb(float   *f1,float   *f2,float   *f3,float   *f4,float   fs,float   alpha1,    
  float   alpha2,int   *l,int   iband,int   itype,int   *ierror);    
  void   mperpsd(complex   x[],int   n,int   nshift,int   nsamp,int   iwidow,float   ts);    
  void   mphares(complex   h[],float   phase[],int   n,float   fs,char   filename[]);    
  void   mprgfft(complex   x[],int   n,int   l,int   lf,int   k1,int   isign);    
  void   mpsplot(float   psdr[],float   psdi[],int   mfre,float   ts);    
  float   randnu(long   *iseed);    
  void   meavar(float   u[],int   *n,float   *pum,float   *puv);    
  void   mrandom(float   u[],int   *n,long   *piseed,int   ITYPE,float   p);    
  void   mrelfft(float   xr[],float   xi[],int   n,int   isign);    
  float   d(int   k,int   n,int   m);    
  float   gee(int   k,int   n);    
  void   mremez1();    
  void   msplfft(complex   x[],int   n,int   isign);    
  void   munwrap(float   phase[],int   n);    
  void   mwindow(float   w[],int   n,int   iwindow,int   *ierror);    
  int   mspbfct(int   i1,int   i2);    
   
  /*-------------------------------------------------------------------*/    
  #endif  
  Top

2 楼olivine1（）回复于 2004-07-01 16:34:00 得分 0

兄弟,你给的代码有部分没清楚哟,帮忙帮到底嘛,分还是少不了给的。  
  你看下面的代码，循环有问题，我又不好乱改。  
  for(k=0;k   {v=0.0;    
  for(j=0;j   v+=p*(u[j+ks]-.5);Top

3 楼olivine1（）回复于 2004-07-01 17:11:22 得分 0

兄弟,你给的代码有部分没清楚哟,帮忙帮到底嘛,分还是少不了给的。  
  你看下面的代码，循环有问题，我又不好乱改。  
  for(k=0;k   {v=0.0;    
  for(j=0;j   v+=p*(u[j+ks]-.5);  
  Top

4 楼wqs6（竹山）回复于 2004-07-01 22:06:16 得分 0

for(k=0;k   {v=0.0;    
  for(j=0;j   v+=p*(u[j+ks]-.5)  
   
  ???  
   
  看不懂，是不是弄错了？Top

5 楼olivine1（）回复于 2004-07-02 11:28:20 得分 0

是你上面给的代码有问题哟,给个全的代码嘛.Top

6 楼wqs6（竹山）回复于 2004-07-02 12:46:56 得分 0

for(k=0;k   {v=0.0;    
  for(j=0;j   v+=p*(u[j+ks]-.5);    
  u[k]=v;    
  ks=ks+ksection;    
  }    
  替换成：  
  for(k=0;k   <ns;k++)  
  {v=0.0;    
  for(j=0;j<k;j++)    
  v+=p*(u[j+ks]-.5);    
  u[k]=v;    
  ks=ks+ksection;    
  }    
     
  Top

相关问题

• iseed
• float
• include

• 帖主：olivine1
• wqs6

• Visual C++类图书
• Visual C++类源码下载