Discussion :

[sara.ber]

Bonjour,

s_v_p, j'ai utilisé la fonction diff pour calculer la dérivée (vitesse), elle est bien exécutée en MATLAB, mais lorsque je l’appelle dans une S-Function, le résultat est toujours zéro, et vraiment je ne comprends pas pourquoi, est-ce que vous pouvez m'aider et merci d'avance.

[duf42]

Bonjour,

Pourrais-tu nous montrer le code que tu utilises dans ta S-function?

Duf

[sara.ber]

Merci bien
voici le code de la S-Function sachant que la sortie doit être la simulation de vitesse.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
function [sys,x0,str,ts]=essaiii(T,xi,X,flag) 
switch flag
 
    case 0  % Initialisation
 
        str = [];
        ts = [0 0]; 
        s = simsizes;
          s.NumContStates = 1;
          s.NumDiscStates = 0;
          s.NumOutputs = 1;
          s.NumInputs = 1;
          s.DirFeedthrough = 0;
          s.NumSampleTimes = 1;
          sys = simsizes(s);
          x0 = [0] % valeur initiale des états
 
     case 1 % calcul de dx/dt
         disp 1
 
 
VX = vitessex(T,X)
 
 sys = accelerationx(VX,T)
 
     case 3 % calcul des sorties
         disp 3
 
       sys=xi
 
     case{2 4 9}
         disp 4
 
         sys=[]
 
     otherwise
 
          error(['unhandled flag = ',num2str(flag)]);
  end

La fonction de vitesse est la suivante:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
function VX=vitessex(T,X)
%calcul de la vitesse suivant X
 
% les entrées sont la composante suivant X et le temps t
% la sortie est la composante de vitesse suivant X, VX
 
 
dD=diff(X)
 
dt=diff(T)
 
VX=(dD)./(dt)
 
 
 
%VX=[0 ;VX(1:length(VX)-1); 0];% on ajoute le zéro comme vitesse(état) initial
%VX=[VX; 0]% de plus le vitesse à la fin du mouvement toujours égale à zéro 
 
VX=[0,VX];
 
return ;

fonction accélération:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
% la fonction accelerationx calcule la composante en accélération suivant l'axe des X 
% les paramètres d'entrée sont le signal de vitesse suivant l'axe des X et la composante temps 
function AX=accelerationx(VX,T)
 
%-----------différentiel du temps
dt=diff(T);
%-----------différentiel du signal vitesse suivant X
dDVX=diff(VX);
AX=(dDVX)./(dt);
AX=[0,AX];
 
return;