Scicos Block
eng - fr

Single step integrator block by trapezoidal method

Contents

• Single step integrator block by trapezoidal method
  • Palette
  • Description
  • Dialog box
  • Default properties
  • Interfacing function
  • Computational function (type 2)
  • See also
  • Authors

Palette

• Tools.cosf - Tools Palette

Description

This block is a single step discrete integrator. The output vector is the integral of the input vector with the trapezoidal method define by the following equation :

$$y_{k}=\frac{h}{2}\left(u_{k}+u_{k-1}\right)+y{k-1}$$

where h is the time step, $u_{k}$ the input vector and $y_{k}$ the output vector.

Dialog box

• Step :
  Type 'vec' of size 1. Integral time step.
• Initial input state :
  Type 'vec' of size -1. Initial value of input state.
• Initial output sate :
  Type 'vec' of size -1. Initial value of output state.

Default properties

• always active: no
• direct-feedthrough: no
• zero-crossing: no
• mode: no
• number/sizes of inputs: 1 / 1
• number/sizes of outputs: 1 / 1
• number/sizes of activation inputs: 1 / 1
• number/sizes of activation outputs: 0 /
• continuous-time state: no
• discrete-time state: yes
• name of computational function: int_trap

Interfacing function

TRAPINTEGRAL_f.sci

Computational function (type 2)

/* int_trap Scicos discrete integral by trapezoidal method block
 * Type 2 simulation function ver 1.0 - scilab-2.6&2.7
 * 17 novembre 2003 - IRCOM GROUP - Author : A.Layec
 */
 
/* REVISION HISTORY :
 * $Log$
 */
 
#include "machine.h"

/* Cette fonction de simulation réalise une intégration discrète
 * du vecteur d'entrée par la méthode des trapèzes :
 * y[k]=h/2*(u[k]+u[k-1])+y[k-1]
 * où h est le pas d'intégration, y[k] le vecteur de sortie,
 * et u[k] le vecteur d'entrée
 *
 * entrée régulières : u[0..nu-1] vecteur d'entrée
 * sorties régulières : y[0..nu-1] registre de sortie
 * entrée évènementielle : Instants de déclenchement
 * sortie évènementielle : néant.
 * état discret : z[0..nu-1]    vecteur des états de sortie discrets précédents
 *                z[nu..2*nu-1] vecteur des états d'entrée discrets précédents
 * 
 * paramètre réel : rpar[0] pas d'intégration
 *
 */
 
/*prototype*/
void int_trap(flag,nevprt,t,xd,x,nx,z,nz,tvec,ntvec,rpar,nrpar,
              ipar,nipar,inptr,insz,nin,outptr,outsz,nout)
integer *flag,*nevprt,*nx,*nz,*ntvec,*nrpar,ipar[],*nipar,insz[],*nin,outsz[],*nout;
double x[],xd[],z[],tvec[],rpar[];
double *inptr[],*outptr[],*t;
{
  /*déclaration des variables*/
  double *y;
  double *u;
  int i,nu;

  /*récupération de l'adresses des ports réguliers*/
  y=(double *)outptr[0];
  u=(double *)inptr[0];

  /*Récupération de la taille du port d'entrée*/
  nu=insz[0];

  /*Le flag 1 calcule l'integrale et place le résultat dans le registre y[]*/
  if(*flag==1)
  {
   for(i=0;i<nu;i++) y[i]=(rpar[0]/2)*(u[i]+z[nu+i])+z[i];
  }

  /*le flag 2 place l'état de sortie y et l'état d'entrée dans z*/
  else if(*flag==2)
  {
   for(i=0;i<nu;i++)
   {
    z[i]=y[i];
    z[nu+i]=u[i];
   }
  }
}

See also

• EULERINTEGRAL_f - Single step integrator block by Euler method (Scicos Block)

Authors

IRCOM Group Alan Layec