最优化实验六.docVIP

实验六 应用MATLAB求解约束非线性规划 实验目的： 掌握应用MATLAB求解约束非线性规划的指令； 灵活运用非线性规划进行建模并用MATLAB求解，进行数值结果分析，提交实验报告。 实验内容： 用MATLAB优化工具箱解约束非线性规划： 模型：min f(x) s.t. Axb Aeq x=beq ceq(x)=0 VLBxVUB 命令为： [x,fval]=fmincon(‘fun’, x0,A,b,Aeq,beq,VLB,VUB, ‘mycon’) 注：(1) fun表示目标函数，x0 为初始点，mycon 表示非线性约束函数，[x,fval]中的 x表示最优解，fval表示最优值。若约束非线性规划的模型中无等式约束条件，则Aeq=[];beq=[];且ceq(x)=[]。 (2) 约束非线性规划中的约束条件： 等式约束：Aeq x=beq表示的是关于决策变量的线性（一次）等式约束； ceq(x)=0表示的是关于决策变量的线性或非线性（一次或高次）等式约束。 不等式约束：Axb表示的是关于决策变量的线性（一次）不等式约束； 表示的是关于决策变量的线性或非线性（一次或高次）不等式约束。 i) 若题目的约束条件都是线性约束（包括等式约束与不等式约束），则将约束条件化成Axb，Aeq x=beq或，ceq(x)=0的形式均可，从而确定A，b或在定义mycon函数时确定和ceq。 ii) 若题目的约束条件包含非线性约束（包括等式约束与不等式约束），则将约束条件化成，ceq(x)=0的形式，从而在定义mycon函数时确定和ceq。 ceq指construction equation等式约束。 求解下列约束非线性规划 （1） min , s.t , 与模型相比，需将其变形为 , . （2）课本P3的例1.2资金使用问题。 数学模型为： 变形为： 变形为： 变形为： 变形为： 变形为： （3）, 实验方法步骤： （1）首先通过建立M文件来对目标函数进行定义 function f=fun(x); f=x(1)^2+x(2)^2+8; 由于该题目的约束条件包含非线性约束，所以需要建立M文件来对非线性约束函数进行定义。 function [g,ceq]=mycon(x); g=[-x(1)^2+x(2)]; ceq=[-x(1)-x(2)^2+2]; 于是在MATLAB命令窗口中输入： x0=[1;1];A=[];b=[]; Aeq=[];beq=[]; VLB=[0;0];VUB=[]; [x,fval]=fmincon(fun,x0,A,b,Aeq,beq,VLB,VUB,mycon) 运行结果为： x = 1 1 fval = 10 （2）首先通过建立M文件来对目标函数进行定义 function f=fun1(x); f=-sqrt(x(1))-sqrt(x(2))-sqrt(x(3))-sqrt(x(4)); 由于该题目的约束条件均为线性约束，所以可用以下两种方法来求解。 法一：建立M文件来对非线性约束函数进行定义。 function [g,ceq]=mycon1(x); g=[x(1)-400;1.1*x(1)+x(2)-440;1.21*x(1)+1.1*x(2)+x(3)-484;1.331*x(1)+1.21*x(2)+1.1*x(3)+x(4)-532.4]; ceq=[]; 于是在MATLAB命令窗口中输入： x0=[400;400;400;400]; A=[];b=[]; Aeq=[];beq=[]; VLB=[0;0;0;0];VUB=[]; [x,fval]=fmincon(fun1,x0,A,b,Aeq,beq,VLB,VUB,mycon1) 运行结果为: x = 86.2295 102.7654 127.1028 153.4692 fval = -43.0856 或在MATLAB命令窗口中输入： [x,fval]=fmincon(fun1,400*ones(4,1),[],[],[],[],zeros(4,1),[],mycon1) 运行结果同上。 即当第一年、第二年、第三年和第四年使用的资金数分别为86.2295万元、 102.7654万元、 127.10