Routine bas-niveau
fr -
eng
dsm_c - routine de calcul de modulateur Sigma-Delta
/* dsm1_c subroutine
* First order Delta-Sigma Modulator
*
* Copyright (C) 2007-2011 Alan Layec
*
* This file is part of modnumlib.
*
* modnumlib is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* modnumlib is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with modnumlib; if not, write to the Free Software
* Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
*
*/
/* REVISION HISTORY :
* $Log$
*/
#include "modnum_lib.h"
/* dsm1_c routine de calcul d'un modulateur sigma-delta du 1er ordre
*
* entrées :
* n : longueur du vecteur d'entrée (scalaire)
* m : gain du modulateur (scalaire)
* u : adresse de départ du vecteur à convertir
* entrées/sorties
* e : intégral du signal d'erreur du dernier élément (scalaire)
* z : etat de sortie du dernier élément +m;-m (scalaire)
* y : adresse de départ du vecteur de sortie (+1;-1)
* sorties :
* q : adresse de départ du vecteur de sortie de l'erreur de quantification
*/
void dsm1_c(int *n,double *m,double *u,double *e,double *z,double *y,double *q)
{
/*déclaration*/
int i;
for(i=0;i<(*n);i++)
{
/*Calcul de l'intégrale du signal d'erreur*/
*e+=u[i]-(*z);
/*Réalisation de la comparaison de signe*/
if((*e)>0) y[i]=1;
else y[i]=-1;
/*mise en mémoire de la sortie*/
*z = *m*y[i];
/*Calcul du bruit de quantification*/
q[i]=(*z)-(*e);
}
return;
}
/* dsm2_c routine de calcul d'un modulateur sigma-delta du 2eme ordre
*
* entrées :
* n : longueur du vecteur d'entrée (scalaire)
* m : gain du modulateur (scalaire)
* u : adresse de départ du vecteur à convertir
* entrées/sorties
* e : intégral du signal d'erreur du dernier élément (vecteur)
* e[0] : modulateur 1
* e[1] : modulateur 2
* z : etat de sortie du dernier élément +m;-m (scalaire)
* z[0] : modulateur 1
* z[1] : modulateur 2
* z[2] : etat mémoire du dsm2
* y : adresse de départ du vecteur de sortie (+1;-1)
* q : adresse de départ du vecteur de sortie de l'erreur de quantification
* w : adresse de départ du vecteur de travail
*
* dépendances :
* dsm1_c
*/
void dsm2_c(int *n,double *m,double *u,double *e,double *z,double *y,double *q,double *w)
{
/*déclaration*/
int i;
/*Appel dsm1*/
dsm1_c(n,m,&u[0],&e[0],&z[0],&y[0],&q[0]);
/*Appel dsm1*/
dsm1_c(n,m,&q[0],&e[1],&z[1],&w[0],&q[0]);
/*réalise opération sur bruit*/
for(i=0;i<(*n);i++)
{
y[i]=(y[i]+z[2]-w[i]);
/*mise en mémoire de l'etat mémoire du dsm2*/
z[2]=w[i];
}
return;
}
/* dsm3_c routine de calcul d'un modulateur sigma-delta du 3eme ordre
*
* entrées :
* n : longueur du vecteur d'entrée (scalaire)
* m : gain du modulateur (scalaire)
* u : adresse de départ du vecteur à convertir
* entrées/sorties
* e : intégral du signal d'erreur du dernier élément (vecteur)
* e[0] : modulateur 1
* e[1] : modulateur 2
* e[2] : modulateur 3
* z : etat de sortie du dernier élément +m;-m (scalaire)
* z[0] : modulateur 1
* z[1] : modulateur 2
* z[2] : etat mémoire du dsm2
* z[3] : modulateur 3
* z[4],z[5]: etats mémoire du dsm3
* y : adresse de départ du vecteur de sortie (+1;-1)
* q : adresse de départ du vecteur de sortie de l'erreur de quantification
* w : adresse de départ du vecteur de travail
*
* dépendances :
* dsm1_c
* dsm2_c
*/
void dsm3_c(int *n,double *m,double *u,double *e,double *z,double *y,double *q,double *w)
{
/*déclaration*/
int i;
/*Appel dsm2*/
dsm2_c(n,m,&u[0],&e[0],&z[0],&y[0],&q[0],&w[0]);
/*Appel dsm1*/
dsm1_c(n,m,&q[0],&e[2],&z[3],&w[0],&q[0]);
/*réalise opération sur bruit*/
for(i=0;i<(*n);i++)
{
y[i]=(y[i]+(w[i]-2*z[4]+z[5]));
/*mise en mémoire des états mémoires du dsm3*/
z[5]=z[4];
z[4]=w[i];
}
return;
}
A. Layec