天津大学《数学模型》课件数学模型第三章约束规划.pptVIP

例：求函数的极小值。解：编写M文件fun5.mfunctiong=fun5(x)f=(x-3)^2-1;主程序[x,fval]=fminbnd(‘fun5’,0,5);1.fminbnd函数：求单变量非线性函数在区间上的极小值Matlab命令为[x,fval]=fminbnd(fun,x1,x2,options).x=3fval=-12.fseminf函数：求其中，w是附加的向量变量，其每个分量都限定在某个区间内Matlab命令为x=fseminf(fun,x0,ntheta,seminfcon,A,b,Aeq,beq).其中，ntheta是半无穷约束phi(x,w)的个数;seminfcon用于定义非线性不等式约束c(x),非线性等式约束ceq(x)和半无穷约束phi(x,w)的每一个分量函数，函数seminfcon有两个输入参量x和s,s是推荐的取样步长，也可以不使用.3.fminimax函数：求其中F(x)={F1(x),F2(x),…,Fm(x)}.Matlab命令为x=fminimax(fun,x0,A,b,Aeq,beq,LB,UB,nonlcon).例：求下列函数构成的函数族取极大极小值时的x值。解：编写fun8.mfunctionf=fun8(x)F=[2*x(1)^2+x(2)^2-48*x(1)-40*x(2)+304-x(1)^2-3*x(2)^x(1)+3*x(2)-18x(1)-x(2)x(1)+x(2)-8];调用[x,fval]=fminimax(@fun8,rand(2,1));x=4.00004.0000fval=0.0000-64.0000-2.0000-0.0000-0.00004.fmincon函数：求约束极小化问题Matlab命令为[x,fval]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options);其中nonlcon表示非线性等式与不等式约束.例：求约束规划解：编写fun9.mfunction[f,g]=fun9(x)f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1);g=[exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+8*x(1)+6*x(2)+1exp(x(1))*(4*x(2)+4*x(1)+2];编写fun10.mfucntion[c,ceq,dc,dceq]=fun10(x)c=[x(1)*x(2)-x(1)-x(2)+1.5;-x(1)*x(2)-10];dc=[x(2)-1,x(1)-1;-x(2),-x(1)];ceq=[];deq=[];调用opts=optimset(‘largescale’,’off’);Opts=optimset(opts,GradObj,on,GradConstr,on);[x,fval]=fmincon(@fun9,rand(2,1),[],[],[],[],[],[],@fun10,opts);x=-9.13991.0570fval=0.0323四.飞行管理问题在约10000m高空的某边长160km的正方形区域内，经常有若干架飞机做水平飞行。区域内每架飞机的位置和速度向量均由计算机记录其数据，以便进行飞行管理。当一架欲进入该区域的飞机到达区域边缘时，记录其数据后，要立即计算并判断是否会与区域内的飞机发生碰撞。如果会碰撞，则应计算如何调整各架(包括新进入的)飞机飞行的方向角，以避免碰撞。先假定条件如下：1)不碰撞的标准为任意两架飞机的距离大于8km;2)飞机飞行方向角调整到幅度不应超过30度；3)所有飞机飞行速度均为每小时800km;4)进入该区域的飞机在到达区域边缘时，与区域内飞机的距离应在60km以上；5)最多需考虑6架飞机；6)不必考虑飞机离开此区域后的状况。请你对这个避免碰撞的飞机管理问题建立数学模型，列出计算步骤，对以下数据进行计算(方向角误差不超过0.01度)，要求飞机飞行方向角调整的幅度尽量小。设该区域