)/4-pi*x(1)*(i+1)/8; a2=x(1)^2*(i-1)^2/8; a=a1+sqrt(a1^2-a2); g(3)=120-180*(1-x(1)*(i-1)/a/pi); %小带轮包角>=alpmin g(4)=0.7*x(1)*(i+1)-a; %中心距范围a>=amin ceq=[];6.1 函数fgoalattain三、例题 [x,fval,exitflag,output, grad,hessian]= fgoalattain(@fun,x0,goal,w,A,b,Aeq,beq,Lb,Ub,’Nlc’,options,P1,P2…)(3)编制优化设计的M文件%V带传动多目标优化设计的调用命令P=4;i=3;n1=1440;KA=1.1;%已知条件:功率,传动比,转速,工况系数x0=[100;1250];%初始点(小带轮直径,V带基准长度)goal=[75,280,2];%分目标w=[10^(-2),40^(-2),1.5^(-2)];%分目标加权系数lb=[80,630];%最小带轮直径和A型V带的基准长度ub=[100;4000];%最大带轮直径和A型V带基准长度[xopt,fopt]=fgoalattain(@VDCD_3mb_MB,x0,goal,w,[],[],[],[],lb,ub,@VDCD_3mb_YS)6.1 函数fgoalattain三、例题 [x,fval,exitflag,output, grad,hessian]= fgoalattain(@fun,x0,goal,w,A,b,Aeq,beq,Lb,Ub,’Nlc’,options,P1,P2…)[ ],[ ],[ ],[ ],