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
| % Optimisation
F=100;
omega=[1];
m=10;
k=100;
p1min=0.001;
p3min=0.001;
p1max=100;
p3max=100;
p1=omega/(sqrt(k/m));
p3=[0.01];
resultats=zeros(7,11);
for alpha=0:0.1:1
i = alpha*10 + 1;
g=f([p1 p3],alpha);
option = optimset('Algorithm','interior-point','display','iter');
[s,fval]=fmincon(@g,[p1 p3],[ ],[ ],[ ],[ ],[p1min p3min],[p1max p3max],[],option);
p1_cal=s(1);
p3_cal=s(2);
H(i)=1/sqrt(1-2*p1_cal^2+p1_cal^4+4*p3_cal^2*p1_cal^2);
T(i)=sqrt(4*p3_cal^2*p1_cal^2+1)/sqrt(1-2*p1_cal^2+p1_cal^4+4*p3_cal^2*p1_cal^2);
W(i)=(p3_cal*p1_cal^2)/(1-2*p1_cal^2+p1_cal^4+4*p3_cal^2*p1_cal^2);
resultats(:,i)=[alpha;s(1);s(2);fval;H(i);T(i);W(i)];
end; |
Partager