/* * awgn Scicos Additive White Gaussian Noise channel block * Type 4 simulation function ver 1.0 - scilab-3.0 * 23 décembre 2004 - IRCOM GROUP - Author : A.Layec */ /* REVISION HISTORY : * $Log$ */ #include "modnum_lib.h" #include "scicos_block.h" /* Cette fonction de simulation calcule un vecteur aléatoire gaussien * par la méthode "Muler Box" et l'ajoute au vecteur complexe d'entrée * * y1[i]=u1[i]+mean+sigma*sqrt(-log(rand1))*cos(2*%pi*rand2) * y2[i]=u2[i]+mean+sigma*sqrt(-log(rand1))*sin(2*%pi*rand2) * * où mean et sigma sont des paramètres et rand1 et rand2 sont * des variables aléatoires uniformément réparties * * entrée régulière : u1[0..nu-1] vecteur des parties rélles des échantillons d'entrée * u2[0..nu-1] vecteur des parties imaginaires des échantillons d'entrée * sortie régulière : y1[0..ny-1] vecteur des parties rélles des échantillons de sortie * y2[0..ny-1] vecteur de parties imaginaires des échantillons de sortie * entrée évènementielle : Dates de déclenchement * sortie évènementeille : néant * * paramètre entier : néant * paramètre réel : rpar[0] : valeur du sigma * rpar[1] : valeur de la moyenne */ /*prototype*/ void awgn(scicos_block *block,int flag) { /*déclaration des variables*/ double *u1,*u2; double *y1,*y2; int k; int ny; /*récupération de l'adresses des ports réguliers*/ u1=(double *)block->inptr[0]; u2=(double *)block->inptr[1]; y1=(double *)block->outptr[0]; y2=(double *)block->outptr[1]; /*récupère taille de sortie*/ ny=block->outsz[0]; /*Le flag 1 calcule le registre de sortie y*/ if(flag==1) { /*Appel noiseiq_c*/ noiseiq_c(&ny,&block->rpar[0],&block->rpar[1],&y1[0],&y2[0]); /*Appel complxa_c*/ cmplxa_c(&ny,(k=1,&k),&u1[0],&u2[0],&y1[0],&y2[0],&y1[0],&y2[0]); } }