Page 175 - 数学建模算法与应用
P. 175

第六章  非线性规划的研究


                      s=[0.2,0;0.2 0];
                   end

                   % 取样值
                   w1=1:s(1,1):100;
                   w2=1:s(2,1):100;
                   % 半无穷约束
                   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,’+’);

                   调用函数 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)
                   (三) 利用梯度求解约束优化问题

                   例  已知函数                                    且满足非线性约束:




                   求
                   分析:当使用梯度求解上述问题时,效率更高并且结果更准确。



                                                                                      165
   170   171   172   173   174   175   176   177   178   179   180