Discussion :

[vtyrode]

Bonjour,
Je suis en train de réaliser une simulation d'un oscillateur harmonique excité par un Dirac et je me retrouve face à un problème.
En effet j'ai fait une modélisation sous SimuLink et une simulation sous Matlab.
La simulation sous Simulink se présente sous la forme suivante:
La fonction de transfert est excité par un dirac.

Sous matlab j'ai utilisé la fonction impulse comme suit:


Ce qui est bizarre est que je ne trouve pas les mêmes résultats. En effet la période est bien la même dans les deux cas mais l'amplitude est complètement différente (facteur d'au moins 1000).

Qq a t'il déjà rencontré ce genre de problème?
Merci

[phryte]

Bonjour,
Je trouve les mêmes résultats.
Qu'as-tu comme coefficients ?
Qu'as-tu mis comme signal Dirac ?

[vtyrode]

Bonjour,
J'ai pris M1=1, K1=1.10^5 et C1=32.
Pour le dirac j'utilise la fonction signal builder et j'essaye de créer le signal le plus court possible.
J'obtiens alors la réponse suivante
Et avec Matlab j'obtiens cette réponse
Le problème vient peut être du dirac qui ne ressemble pas à la fonction impulse de Matlab?

merci

[phryte]

Bonjour,
Tu ne respecte pas la définition de Dirac :

l'amplitude de l'impulsion de Dirac est infinie (délicat), son support ou largeur est nulle et son aire est égale à 1.

Donc tu prends un signal 1x1 ou 0.5x2.....
Pour Matlab j'ai cette réponse :

[vtyrode]

Ok merci pour la précision sur le dirac.
Quelle valeurs de m,c et k utilise tu pour ta réponse?

[phryte]

Bonjour,

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
clear
M1=1;C1=32.;K1=1.10^5;
num1=[1];
den1=[M1 C1 K1];
H1=tf(num1,den1);
figure
impulse(H1)
grid
[y,t,x,ysd]=impulse(H1);
figure
plot(t,y)
grid

[vtyrode]

Bonjour,
Attention car en l'écrivant ainsi la valeur de K1 n'est pas la bonne.
Ainsi cela fausse les résultats...

[vtyrode]

Mais cela n'explique pas le fait que tu ais trouvé la même réponse aux deux modèles...
Est ce que tu l'as fait avec Simulink?

[phryte]

J'avais pris K1=1.10^5 comme tu l'avais écrit. Avec K1=10^5 je trouve comme toi.
Les coefficients me paraissent mal adaptés.
K1 pour un passe-bas d'ordre deux vaut : Omega0^2, de même que le numérateur !