Discussion :

[lecteur1001]

Bonjour,

je possède les résultats d'une simulation et je souhaite comparer un modèle du 1er ordre retardé avec ces résultats de simulation pour savoir si le modèle est fidèle à ma courbe.
Mes données de simu sont sous Excel et je les ai donc importé avec xlsread, pas de problème. Je peux donc "ploter" ces résultats de simu.

J'ai déclaré mon modèle de la manière suivante :

Code :

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

1
2
s = tf('s');
H = exp(-0.4*s)*0.06/(0.2*s+1);

puis je fais un

Code :

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

step(H)

pour connaître sa réponse indicielle.

Ma question est : comment ploter sur la même figure le step(H) et les résultats de simulation pour faire la comparaison ?

Merci d'avance pour les réponses

[duf42]

Bonjour,

Le plus simple c'est de faire un hold on après avoir tracé ton step et tracer ta courbe (avec plot j'imagine...)

Duf

[lecteur1001]

J'avais déjà essayé, mais ça ne fonctionne pas.
Par contre, j'ai tenté de faire un

Code :

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

step(H,t)

avec t étant la variable de temps dans mon fichier de simulation.
Le problème est que t varie entre 160 et 175 et à ce moment là, ben la réponse indicielle de mon modèle est déjà à son maximum depuis longtemps... Comment faire pour que la réponse indicielle commence à t=160 plutôt qu'à 0 ?

Edit : d'ailleurs il est dit dans l'aide de step :

STEP(SYS,T) uses the user-supplied time vector T for simulation.
For discrete-time models, T should be of the form Ti:Ts:Tf
where Ts is the sample time. For continuous-time models,
T should be of the form Ti:dt:Tf where dt will become the sample
time for the discrete approximation to the continuous system. The
step input is always assumed to start at t=0 (regardless of Ti)

[duf42]

J'avais déjà essayé, mais ça ne fonctionne pas.

C'est à dire? Qu'est-ce qui ne fonctionne pas?

[lecteur1001]

C'est à dire? Qu'est-ce qui ne fonctionne pas?

Ben le plot après le hold on n'était pas pris en compte.

De tout façon, le mieux est de faire un

Code :

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

[t,Y] = step(H,t)

Et puis après je pourrais faire un

Code :

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

plot(t,donnee_simu,t,Y)

Le problème est, comme je l'ai dit, que la fonction step débute à t=0 au lieu de --comme je le souhaiterai-- de débuter à 160.

[duf42]

Ben le plot après le hold on n'était pas pris en compte.

Etrange, chez moi ca marche

De tout façon, le mieux est de faire un

Code :

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

[t,Y] = step(H,t)

Et puis après je pourrais faire un

Code :

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

plot(t,donnee_simu,t,Y)

Je suis tout à fait d'accord avec toi

Le problème est, comme je l'ai dit, que la fonction step débute à t=0 au lieu de --comme je le souhaiterai-- de débuter à 160.

Si c'est juste pour la représentation, tu peux peut-être décaler l'un de tes graphes :

Code :

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

1
2
3
4
 
plot(t+160,donnee_simu);
hold on;
plot(t,Y);

Si je ne me suis pas trompé de sens

[lecteur1001]

Si c'est juste pour la représentation, tu peux peut-être décaler l'un de tes graphes :

Code :

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

1
2
3
4
 
plot(t+160,donnee_simu);
hold on;
plot(t,Y);

Si je ne me suis pas trompé de sens

Tu vas me trouver chiant, mais au lieu de décaler mes simus, je veux décaler le step

[duf42]

Ah ben après tu fais bien comme tu l'entends, c'est pas génant

[lecteur1001]

Ok, je m'en sors bêtement (grrr) avec un :

Code :

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

1
2
3
4
s = tf('s');  
H = exp(-0.4*s)*0.06/(0.2*s+1);
[Y,T] = step(H);
plot(T+t(1),Y,t,donnee_simu)

Merci pour l'aide