[whxgj]

while(F>0)
b=m;
m=(a+b)/2;
omegam0=(1-m)*omegas/pairs; %omegam0重新赋值
vib02; %重新计算振动机参数
Tcapture=m0*r0^2*omegam0^2*abs(Wc);
Te01=Pe1*S0/(1+Tst1*S0^2);
Te02=Pe2*S0/(1+Tst2*S0^2);
Tdifference=abs(Te01-Te02-(f1-f2)*omegam0-m0*r0^2*omegam0^2*(Ws1-Ws2));
a1=[Te01-Te02-(f1-f2)*omegas*(1-S0)/pairs-Pm*(1-S0)^2*(Ws1-eta*Ws2)/(2*pairs^2)]/(Pm*(1-S0)^2*Wc/(2*pairs^2));
alpha1=asin(a1)-thetac; %2*alpha
F=Te01+Te02-(f1+f2)*omegas*(1-S0)/pairs-Pm*(1-S0)^2*(Ws1+eta*Ws2+2*Ws*sin(alpha1+thetas))/(2*pairs^2);%F(alpha,omegam0)
if abs(F)<delta
disp(' “Succeed: the values of 2alpha and omegm0 ”')
return
end
end