r=0.145*t*(sqrt(1+d^2*f/10)+sqrt(f/5)*d^2/(8*t))/(d^2*P);
m=0.000294*t*(1+1/sqrt(9+d^2*f/20)+0.85*d/t);
a=4*r/((1+r)^2+(2*pi*f*m-cot(2*pi*f*m))^2);
     确定决策变量和约束条件
其中孔径，板厚、空腔深度的单位为mm，穿孔率以%为单位，确定四个待优化参数范围分别为d[O.1,1],t[O.1,1] ,
P[o.5,5]，D[10,200]。