Discussion :

[stephan1609]

Bonjour tout le monde,

j'ai un petit problème dans matlab pour le code de lsim.
J'ai un programme qui simule la réponse pour une fonction de transfert du type:

(A*jw+B)/(C*jw+D) en utilsant tf puis lsim.

Je doit reproduire cette fonction en C++ mais j'aimerai d'abord comprendre ce que fait cette fonction exactement. Quelqu'un aurait il une idée?
merci

[tachmou]

Bonjour,

Tu t'attèles à un problème bien difficile là... à moins de compiler lsim depuis matlab en C++ directement, mais je ne pense pas que ça soit possible et je ne vois pas du tout comment faire.

lsim calcule la réponse temporelle à une excitation de commande en entrée. En clair :
- tu crées ta commande u(t) (pratiquement un tableau de valeurs, une valeur de ton entrée pour chaque pas de temps dans chacune des cases du tableau).
- tu envoies ta fonction de transfert et ton entrée u, ainsi que ta base temporelle, à lsim qui lui se charge de calculer la sortie s(t)

Pour effectuer un tel tour de force, théoriquement, il faut que tu transforme ta fonction de transfert qui est exprimée avec la variable de laplace "p" en fonction temporelle qui elle est exprimée avec la variable de temps "t" : c'est la transformation de Laplace inverse. La théorie est belle... la pratique : eh bien bon courage, du moins personnellement je ne vois pas trop comment faire, mais la solution m'intéresse.

Par contre si ta fonction de transfert ne change pas, sauf les valeurs de tes paramètres A, B, C, D et w, dans ce cas alors il te suffit de coder une bête fonction s(t_i) = f(A,B,C,D,w,u_i) sachant que u_i est la valeur de ton entrée u à l'instant t_i et tu fais une boucle sur i pour calculer s(t) complet. (sachant que f correspond à la transformation de laplace inverse de ta fonction de transfert : tu auras peut-être - au pire - des cas particulier à traiter lorsque tes paramètres sont nuls ou lorsque le dénominateur s'annule)

Mais si tu veux vraiment recoder lsim, il va falloir t'atteler à du calcul symbolique...

Cordialement