How do I optimize my code with respect to 4 variables (x_1 x_3 k1 k2 k3)?

1 次查看（过去 30 天）

Alex M 2017-11-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/368260-how-do-i-optimize-my-code-with-respect-to-4-variables-x_1-x_3-k1-k2-k3

评论： Walter Roberson 2017-11-21

%%Parameters
format short
% Constants
rho=1.225;
U=170;
s=1;
c=0.5;
c_2=1;
b=0.5*c;
m=5;
I_CG=0.05;
m_p=0.5;
% x_p=(0:0.01:c);
zero=[0 0; 0 0];
% Variables
x_g=0.5*c;
x_ac=0.25*c;
x_1=0.1*c;
x_3=0.7*c;
q=0.5*rho*u^2;
c_2=1;
%x_1 x_3 k1 k2 k3 have to changed in order to result in the highest u;
%k1+ k2+ k3 <30000
for u = [0:1:75]
M=[ m   -m*(x_ac - x_g);  -m*(x_ac - x_g)       I_CG + m*(x_ac - x_g)^2  ];
E=[k1 + k2 + k3     k1*(x_1 - x_ac) + k3*(x_3 - x_ac); k1*(x_1 - x_ac) + k3*(x_3 - x_ac)  k1*(x_1 - x_ac)^2 + k3*(x_3 - x_ac)^2];
B=[1+c_2/(pi*s*c*rho) c/2; 0 c^2/16];
Bbar=pi*s*c*rho*[B];
K=[0 1; 0 0];
Kbar=pi*s*c*rho*[K];
A=[[M] u*[Bbar];[zero] [M]];
C=[[zero] [E]+u^2*[Kbar] ; -[M] [zero] ];
c_2;
u;
lambda=eig(C,-A);
   if any(lambda>0) 
   break 
   end 
u

4 个评论
显示 2更早的评论隐藏 2更早的评论

Alex M 2017-11-21

u velocities are being tested, and the best combination of the x_1 x_3 k1 k2 k3 will give me the greatest u velocity before the code reads a positive value in the eigenvector.