This block realizes the sub-system H of Chua's circuit :
/* systemh Scicos Chua sub-systeme H block
* Type 2 simulation function ver 1.0 - scilab-2.6&2.7&3.0
* 19 novembre 2003 - IRCOM GROUP - Author : A.Layec
*/
/* REVISION HISTORY :
* $Log$
*/
#include "machine.h"
/*Cette fonction de simulation réalise le sous-systeme H
* du circuit de Chua :
*
* y=x1
*
* dx1/dt=B*(u-x1+x2)
* dx2/dt=C*x1
*
* où B et C sont des vecteurs de paramètres réels,
* E=1, u est le vecteurs des entrées, y le vecteur des sorties
* x1 et x2 les variables d'états
*
* Entrée régulière : u[0..nu-1] : vecteur des entrées
* Sortie régulière : y[0..nu-1] : vecteur des sorties
* Entrée évènementielle : néant (héritage)
* sortie évènementielle : néant
* paramètres : rpar[0..nu-1] vecteur B
* rpar[nu..2*nu-1] vecteur C
* etat continu : x[0..nu-1] vecteur de la variable d'état x1
* x[nu..2*nu-1] vecteur de la variable d'état x2
*/
/*prototype*/
void systemh(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 nu,i;
/*récupération des adresses des ports réguliers*/
y=(double *)outptr[0];
u=(double *)inptr[0];
/*récupération de la taille du port d'entrée u*/
nu = insz[0];
/*Le flag 1 place l'état continu x dans le reg y*/
if(*flag==1||*flag==6)
{
for(i=0;i<nu;i++) y[i]=x[i];
}
/*Le flag 0 calcule la dérivée de l'état continu*/
else if(*flag==0)
{
for(i=0;i<nu;i++)
{
xd[i]=rpar[i]*(u[i]-x[i]+x[nu+i]);
xd[nu+i]=rpar[nu+i]*x[i];
}
}
}