数学建模第二章非线性规划分析.pptVIP

%半无穷约束 k1=sin(w1*x(1)).*cos(w1*x(2))-1/1000*(w1-50).^2-sin(w1*x(3))-x(3)-1; k2=sin(w2*x(2)).*cos(w2*x(1))-1/1000*(w2-50).^2-sin(w2*x(3))-x(3)-1; %画出半无穷约束的图形,此处可以不用 plot(w1,k1,-,w2,k2,+); （3）调用函数fseminf 在Matlab 的命令窗口输入 [x,y]=fseminf(@fun6,rand(3,1),2,@fun7) 即可。 3 fminimax 函数 上述问题的Matlab 命令为 X=FMINIMAX(FUN,X0,A,B,Aeq,Beq,LB,UB,NONLCON) 例3.13 解 （1）编写M 文件fun8.m 定义向量函数如下： function f=fun8(x); f=[2*x(1)^2+x(2)^2-48*x(1)-40*x(2)+304 -x(1)^2-3*x(2)^2 x(1)+3*x(2)-18 -x(1)-x(2) x(1)+x(2)-8]; （2）调用函数fminimax [x,y]=fminimax(@fun8,rand(2,1)) 4 利用梯度求解约束优化问题 分析：当使用梯度求解上述问题时，效率更高并且结果更准确。 题目中目标函数的梯度为(分别对x1,x2