Fonction de calcul
fr


awgn

Contenu du fichier


/*
 * 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]);
  }
}