Function gear_opt1
clc；
clear all；
close all；
n1=960；
p1=10；
n2=850；
p2=8.6；
x0=[2,3,27,27,4.13,pi/18,3.15,0.8]；
[x,fval]=fmincon(@gearoptfun2,x0,[],[],[],[],[],[],
@gearoptcons2,[],p1,p2,n1,n2)
function f=gearoptfun2(x,p1,p2,n1,n2)
f=((x(1)*x(3))*(1+x(5))+x(2)*x(4)*(1+13/x(5)))/(2*cos(x(6)));
function[C,Ceq]=gearoptcons2(x,p1,p2,n1,n2)
c(1)=2-x(1)；
c(2)=x(1)-5；
c(3)=3.5-x(2)；
c(4)=x(2)-6；
c(5)=18-x(3)；
c(6)=x(3)-32；
c(7)=20-x(4)；
c(8)=x(4)-32；
c(9)=1.3*x(5)-x(7)；
c(10)=x(7)-1.5*x(7)；
x(5)=13/x(7)；k1=1.9；
k2=1.92；
x(6)=3.14/18；
z1=189.8^(1/2)*2.5*sqrt((1+x(5))/x(5))；
z2=189.8^(1/2)*2.5*sqrt((1+x(7))/x(7))；
T1=9550*1000*p1/n1；
T2=9550*1000*p2/n2；
c(11)=z1*sqrt(2*k1*T2/(x(1)*x(3))^3*x(8)*1.7)-570；
c(12)=z2*sqrt(2*k2*T2/(x(1)*x(3))^3*x(8)*1.64)-400；
[YFa1,Ysa1]=yfasa(x(3)/(cos(6))^3)；
[YFa2,Ysa2]=yfasa(x(5)*x(3)/(cos(6))^3)；
[YFa3,Ysa3]=yfasa(x(4)/(cos(6))^3)；
[YFa4,Ysa4]=yfasa(x(7)*x(4)/(cos(6))^3)；
c(13)=2*k1*T1/(x(1)^3*x(3)^2*x(8)*1.7)*YFa1*Ysa1*0.88-440；
c(14)=2*k1*T1/(x(1)^3*x(3)^2*x(8)*1.7)*YFa2*Ysa2*0.88-330；
c(15)=2*k2*T2/(x(2)^3*x(4)^2*x(8)*1.64)*YFa1*Ysa1*0.88-440；
c(16)=2*k2*T2/(x(2)^3*x(4)^2*x(8)*1.64)*YFa2*Ysa2*0.88-330；
ceq=[]；
function[yfa,ysa]=yfasa(z)
x=[17 18 19 20 21 22 23 24 25 26 27 28 29 30 35 40 45 50 60 70 80 90 100 150]；
y1=[2.97 2.91 2.85 2.80 2.76 2.72 2.69 2.65 2.62 2.60 2.57 2.55 2.53 2.52 2.45 2.40 2.35 2.32 2.28 2.24 2.22 2.20 2.18 2.14]；
y2=[1.52 1.53 1.54 1.55 1.56 1.57 1.575 1.58 1.59 1.595 1.60 1.61 1.62 1.625 1.65 1.67 1.68 1.70 1.73 1.75 1.77 1.78 1.79 1.83]；
yfa=interp1(x,y1,z)；
ysa=interp1(x,y2,z)；