Discussion :

[timot91]

Bonjour,
J'essaye de calculer un coefficient de puissance à partir d'une équation prise dans une thèse mais je n'arrive pas à la courbe espérer.
Voila l'équation ainsi que la courbe demandé



Et voila ce que j'ai fais sur simulink avec dans f(u) l'équation du sujet, et ramp qui simule des valeurs entre 0 et 9.



Je ne comprend pas pourquoi j'obtiens cette courbe.



Si vous pouvez m'aider je vous en serais très reconnaissant !

Merci

Timothée

[le fab]

salut

tu as quoi dans f(u) ? (parce que l'équation du sujet, c'est flou!)
ca devrait être lambda1*u(1)+lambda2*u(2)+lambda3*u(3)+lambda4*u(4) +lambda5.

tu sais que tu aurais pu éviter le truc illisible avec les * et le mux à gauche de f(u) et mettre qu'une entrée et directement l''équation lambda1*u^4+lambda2*u^3+lambda3*u^2+lambda4*u +lambda5 ?

sinon verifie bien que tu as rentré les bon coefficients, parce que c'est sensé marcher :

Code :

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

ezplot(@(x) 1.92e-4*x^4 - 5.2e-3*x^3 + 2.44e-2*x^2 +5.79e-2*x +1.38e-4,[0 9]);

Fabien

[timot91]

Salut Fabien, merci pour ta réponse!
Dans f(u) j'avais mis ça: a1*u[1]+a2*u[2]+a3*u[3]+a4*u[4]+a5 parce que je croyais avoir compris que les u justement venait de ramp. Et j'ai vérifié j'ai bien mis les bons coeffs.
Je vais essayer d'enlever les trucs illisibles avec ton équation mais c'est sur l'entrée que je bloque en fait, il faut bien que ça soit plusieurs valeurs aléatoire non ?

[le fab]

tu mets une seule entrée (qui viens de la rampe) et tu élèves les puissances dans l'équation
ca me parait quand même plus lisible

ou alors tu développes tout en blocs simulink (bloc puissance, addition, gain ...)

mais le mélange des deux est peu esthétique
(et tu dois avoir une erreur d'implémentation quelque part, qui ne saute pas aux yeux)

[timot91]

C'est ce que j'ai fais mais j'obtiens le même résultat:

j'ai mis 1.92e-4*u^4 - 5.2e-3*u^3 + 2.44e-2*u^2 +5.79e-2*u +1.38e-4 dans f(u) mais ca me donne exactement la même courbe...

[le fab]

bizarre bizarre, je viens de copier coller ton équation et je l'ai mise dans un bloc interpreted function comme toi : et ca marche !!

tu as mis quoi dans ta ramp ? tu aurais pas mis 10 dans la pente des fois ??
met un scope ne sortie de ta rampe, je parie qu'elle ne va pas de 0 à 9

[timot91]

Si si j'ai bien mis de 0 à 9, le scope après la ramp donne ca



donc je comprend vraiment pas !

[le fab]

bein voila , tu a mis une pente de 10, du coup la sortie de la rampe varie de 0 à 90 comme le montre ton graphique (faut regarder les ordonnées qui reprensnete la sortie du bloc rampe en fonction du temps en abcisse)
et du coup c'est normal !!

essaye ça :

Code :

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

1
2
3
4
5
figure
subplot 211
ezplot(@(x) 1.92e-4*x^4 - 5.2e-3*x^3 + 2.44e-2*x^2 +5.79e-2*x +1.38e-4,[0 9]);
subplot 212
ezplot(@(x) 1.92e-4*x^4 - 5.2e-3*x^3 + 2.44e-2*x^2 +5.79e-2*x +1.38e-4,[0 90]);

met une pente de 1 dans le bloc

[timot91]

Ah mais oui d'accord j'ai compris ! Erreur bête, ca fonctionne parfaitement!
Merci beaucoup Fabien!
Je continue sur mon projet, y'a moyen que je revienne ici quand je serais bloqué