This is my program, it is connected with a function:
m_s=10*ones(3,1); %[Kg]
Mtot=sum(m_s);
k_s=1000*ones(3,1); %[N/m]
C_s=[1.246,0.176,0.097;0.176,1.343,0.272;0.097,0.272,1.519];
omega_s=4.4504; %first mode
m_t_1_1=0.05*Mtot;
beta=1;
k_t_1_1=(beta*omega_s)^2*m_t_1_1;
ksi_t=0;
c_t_1_1=2*ksi_t*sqrt(m_t_1_1*k_t_1_1);
v0=[k_t_1_1,...
c_t_1_1];
LB=[0.8*(omega_s^2*m_t_1_1) 0.01*(2*sqrt(k_t_1_1*m_t_1_1))];
UB=[1.2*(omega_s^2*m_t_1_1) 0.35*(2*sqrt(k_t_1_1*m_t_1_1))];
[v,h]=fminsearchbnd(@(v) hinf11tmd(v,m_s,k_s,C_s,m_t_1_1),v0,LB,UB);
beta_1_1=sqrt(v(1)/m_t_1_1)/omega_s
ksi_1_1=v(2)/(2*sqrt(m_t_1_1*v(1)))
I need to put in the right side of LB and UB:
0.01*(2*sqrt(k_t_1_1*m_t_1_1))
0.35*(2*sqrt(k_t_1_1*m_t_1_1))
instead of k_t_1_1 the second component of v, v(2).
In other words v(2) is ok only if it respects this constraint.
How I could do it?
Thank you
