1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
|
function f = valeur(X);
a1=X(1);
d2=X(2);
% Le polynome P1.
P1 = [a1 -3*a1 a1 0];
% Les differentes contraintes de P2, mise sous forme matricielle.
M2=[27 9 3 ; 216 36 6; 27 6 1];
B2=[ 3-d2; -d2; polyval(polyder(P1),3) ];
X2=M2/B2';
P2 = [X2(1) X2(2) X2(3) d2];
%Les differentes contraintes de P3, mise sous forme matricielle.
M3 = [6*6*6 6*6 6 0; 9*9*9 9*9 9 0; 3*9*9 2*9 1 0; 3*6*6 2*6 1 0];
B3 = [0; 0; polyval(polyder(P2),6) ; 1];
X3=M3/B3';
P3 = [X3(1) X3(2) X3(3) X3(4)];
% la fonction a optimiser :
g1 = @(t) sqrt(polyval(polyder(P1),t)^2+1);
g2 = @(t) sqrt(polyval(polyder(P2),t)^2+1);
g3 = @(t) sqrt(polyval(polyder(P3),t)^2+1);
f = quadl(g1,0,3) + quadl(g2,3,6)+ quadl(g3,6,9); |
Partager