Bloc passage à zéro discret
- Step :
Type 'vec' de taille 1. La description du paramètre 1.
- Value of transition :
Type 'vec' de taille 1. La description du paramètre 2.
- Initital Value of transition :
Type 'vec' de taille 1. La description du paramètre 3.
- Date of initial output event :
Type 'vec' de taille 1. La description du paramètre 4.
- Dynamical Value(0:No/1:Yes) :
Type 'vec' de taille 1. La description du paramètre 5.
- toujours actif: non
- direct-feedthrough: non
- détection de passage à zéro: non
- mode: non
- nombre/taille des entrées régulières: 1 / 1
- nombre/taille des sorties sorties régulières: 0 /
- nombre/taille des entrées évènementielles: 1 / 1
- nombre/taille des sorties évènementielles: 1 / 1
- possède un état continu: non
- possède un état discret: oui
- nom de la fonction de calcul: vcoevt
VCOEVT_f.sci
/* vcoevt Scicos Time transition calculator block
* Type 4 simulation function ver 1.1 - scilab-3
* 7 Avril 2004 - IRCOM GROUP - Author : A.Layec
* 21 Novembre 2004 : rajout option ipar[0]
* 10 janvier 2005 : passage au type 4
*/
/* REVISION HISTORY :
* $Log$
*/
#include "machine.h"
#include "scicos_block.h"
/* Cette fonction réalise le calcul approximé de la date d'évenement
* correspondant au passage du signal d'entrée aux valeurs rpar[1]
* (rpar[1]*u2[0] si ipar[0]=1).
* Rmq : L'entrée est considérée strictement croissante et monotone.
*
* entrées régulières : u1[0] : wot+phi(t)
* u2[0] : entrée dynamique(rapport de division)
* sortie régulière : néant
* entrée d'évenement : Pas de calcul : Tsampl
* sortie d'évenement : date du passage à rpar[1] (ou rpar[1]*u2[0])
* paramètre entier : ipar[0] : 0 : passage aux valeurs rpar[1]
* 1 : passage aux valeurs rpar[1]*u2[0]
* paramètres réels : rpar[0] : Tsampl
* rpar[1] : valeur de passage en radian
* état discret : z[0] : u(k-1)
* z[1] : valeur de u à la date de passage
*/
/*prototype*/
void vcoevt(scicos_block *block,int flag)
{
/*déclaration des variables*/
double *u1;
double v; /*valeur de passage*/
/*Récupération des adresses des ports d'entrées*/
u1=(double *)block->inptr[0];
/*Test si entrée 2*/
if (block->ipar[0]==1)
{
/*déclaration*/
double *u2;
/*récupération adresse*/
u2=(double *)block->inptr[1];
v=block->rpar[1]*u2[0];
}
else v=block->rpar[1];
/*le flag 3 calcul la date de transition*/
if(flag==3)
{
/*Test si la valeur est >=rpar[1]*/
if((u1[0]-block->z[1])>=v)
{
/*Calcul la date par extrapolation*/
/*tvec[0] = *t + (((z[1]+v)-z[0])*rpar[0])/(u1[0]-z[0]);*/
block->evout[0]=(((block->z[1]+v)-block->z[0])*block->rpar[0])/(u1[0]-block->z[0]);
/*mémorise valeur à la date*/
block->z[1] = block->z[1] + v;
}
/*mémorise u*/
block->z[0]=u1[0];
}
}
IRCOM Group
Alan Layec