Calcul d'un coefficient de puissance

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é

Pièce jointe 326891

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.

Pièce jointe 326896

Je ne comprend pas pourquoi j'obtiens cette courbe.

Pièce jointe 326897

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

Merci

Timothée

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:

---

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

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 ?

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)

C'est ce que j'ai fais mais j'obtiens le même résultat:
Pièce jointe 327457
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...

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

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

Pièce jointe 327662

donc je comprend vraiment pas !

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:

---

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

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é ;)