Ce bloc réalise un canal additif à bruit blanc gaussien. Il additionne des échantillons bruités au vecteur des entrées régulières et la sortie régulière peut s'exprimer par la formule suivante :
Ce bloc calcule toujours des échantillons de bruit complexes, mais comme montré par le modèle équivalent, il peut-être utilisé pour des vecteurs de signaux réels en utilisant seulement les premiers ports réguliers d'entrée/sortie.
/*
* 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 "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]);
}
}