wilson法设计_wilson_Wilson法_叶片wilson_

function [c,ceq]=mycon(x) D=2.3;%直径 R=0.5*D; %半径 r=1.1;%截面到中心距离 B=5;%叶片数 Y0=3.5;%尖速比 Y=r/R*Y0;%周速比 Q=atan((1-x(1))/((1+x(2))*Y)); k=(-B/2)*(R-r)/(r*sin(Q)); F=(2/pi)*acos(exp(k)); c=[]; ceq=x(2)*(1+x(2))*Y^2-x(1)*(1-x(1)*F);