Add here a paragraph of the function description.
/* noiseblk_c subroutine * Gaussian Noise generator * with Box Muller Law method * IRCOM GROUP - Author : A.Layec */ /* REVISION HISTORY : * $Log$ */ #include "mod_num_lib.h" /* noiseblk_c routine de calcul d'échantillons bruités par la méthode "Box Muller Law" * * Entrées : * n : taille des vecteurs * typ : type de sortie (0:cos/1:sin) * sig : variance (scalaire) * mean : moyenne (scalaire) * Sorties : * y : vecteur des sorties * * dépendances * math.h * */ void noiseblk_c(int *n,int *typ,double *sig,double *mean,double *y) { /*déclaration des variables*/ int i; double rand1, rand2; double rand_m; double ampl, phase; /*récupération de la valeur de RAND_MAX*/ rand_m=RAND_MAX; for(i=0;i<(*n);i++) { /*calcul rand1*/ rand1=rand()/rand_m; /*test rand1*/ while((rand1<=0)||(rand1>=1)) rand1=rand()/rand_m; /*calcul rand2*/ rand2=rand()/rand_m; /*test rand2*/ while((rand2<=0)||(rand2>=1)) rand2=rand()/rand_m; /*Calcul amplitude et phase*/ ampl=(*sig)*sqrt(-log(rand1)); phase=2*M_PI*rand2; /*Calcul y*/ if((*typ)==0) y[i]=(*mean)+ampl*cos(phase); else if((*typ)==1) y[i]=(*mean)+ampl*sin(phase); } return; }