/* dsm1_c subroutine * First order Delta-Sigma Modulator * * Copyright (C) 2007-2009 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; }