试验优化设计-数学建模非线性规划9-1_精品文档.ppt

试验优化设计-数学建模非线性规划9-1_精品文档.ppt

《试验优化设计-数学建模非线性规划9-1_精品文档.ppt》由会员分享，可在线阅读，更多相关《试验优化设计-数学建模非线性规划9-1_精品文档.ppt（33页珍藏版）》请在冰豆网上搜索。

数学建模数学建模非线性规划非线性规划1*非线性规划的基本解法非线性规划的基本解法非线性规划的基本概念非线性规划的基本概念非线性规划非线性规划2定义定义如果目标函数或约束条件中至少有一个是非线性函数时的最优化问题就叫做非线性规划问题非线性规划问题非现性规划的基本概念非现性规划的基本概念一般形式一般形式:

（1）其中，是定义在En上的实值函数，简记:

其它情况其它情况:

求目标函数的最大值或约束条件为小于等于零的情况，都可通过取其相反数化为上述一般形式3非线性规划的基本解法非线性规划的基本解法SUTM外点法外点法SUTM内点法（障碍罚函数法）内点法（障碍罚函数法）1、罚函数法、罚函数法2、近似规划法近似规划法4罚函数法罚函数法罚函数法罚函数法基本思想是通过构造罚函数把约束问题转化为一系列无约束最优化问题，进而用无约束最优化方法去求解这类方法称为序列无约束最小化方法序列无约束最小化方法简称为SUMTSUMT法法其一为SUMTSUMT外点法外点法，其二为SUMTSUMT内点法内点法5其中T（X,M）称为罚函数罚函数，M称为罚因子罚因子，带M的项称为罚项罚项，这里的罚函数只对不满足约束条件的点实行惩罚：

当时，满足各，故罚项=0，不受惩罚当时，必有的约束条件，故罚项0，要受惩罚SUTMSUTM外点法外点法6罚函数法的缺点缺点是：

每个近似最优解Xk往往不是容许解，而只能近似满足约束，在实际问题中这种结果可能不能使用；在解一系列无约束问题中，计算量太大，特别是随着Mk的增大，可能导致错误1、任意给定初始点X0，取M11，给定允许误差，令k=1；2、求无约束极值问题的最优解，设为Xk=X（Mk），即；3、若存在，使，则取MkM（）令k=k+1返回

（2），否则，停止迭代得最优解.计算时也可将收敛性判别准则改为.SUTMSUTM外点法外点法（罚函数法）的迭代步骤迭代步骤7SUTMSUTM内点法（内点法（障碍函数法）8内点法的迭代步骤内点法的迭代步骤9近似规划法的基本思想近似规划法的基本思想：

将问题（3）中的目标函数和约束条件近似为线性函数，并对变量的取值范围加以限制，从而得到一个近似线性规划问题，再用单纯形法求解之，把其符合原始条件的最优解作为（3）的解的近似近似规划法近似规划法每得到一个近似解后，都从这点出发，重复以上步骤这样，通过求解一系列线性规划问题，产生一个由线性规划最优解组成的序列，经验表明，这样的序列往往收敛于非线性规划问题的解。

10近似规划法的算法步骤如下算法步骤如下1112用MATLAB软件求解,其输入格式输入格式如下:

1.x=quadprog（H,C,A,b）;2.x=quadprog（H,C,A,b,Aeq,beq）;3.x=quadprog（H,C,A,b,Aeq,beq,VLB,VUB）;4.x=quadprog（H,C,A,b,Aeq,beq,VLB,VUB,X0）;5.x=quadprog（H,C,A,b,Aeq,beq,VLB,VUB,X0,options）;6.x,fval=quaprog（.）;7.x,fval,exitflag=quaprog（.）;8.x,fval,exitflag,output=quaprog（.）;1、二次规划、二次规划13例例11minf（x1,x2）=-2x1-6x2+x12-2x1x2+2x22s.t.x1+x22-x1+2x22x10,x20MATLAB（youh1）1、写成标准形式写成标准形式：

2、输入命令输入命令：

H=1-1;-12;c=-2;-6;A=11;-12;b=2;2;Aeq=;beq=;VLB=0;0;VUB=;x,z=quadprog（H,c,A,b,Aeq,beq,VLB,VUB）3、运算结果运算结果为：

x=0.66671.3333z=-8.2222s.t.141.首先建立M文件fun.m,定义目标函数F（X）:

functionf=fun（X）;f=F（X）;2、一般非线性规划、一般非线性规划其中X为n维变元向量，G（X）与Ceq（X）均为非线性函数组成的向量，其它变量的含义与线性规划、二次规划中相同.用Matlab求解上述问题，基本步骤分三步：

153.建立主程序.非线性规划求解的函数是fmincon,命令的基本格式如下：

（1）x=fmincon（fun,X0,A,b）

（2）x=fmincon（fun,X0,A,b,Aeq,beq）（3）x=fmincon（fun,X0,A,b,Aeq,beq,VLB,VUB）（4）x=fmincon（fun,X0,A,b,Aeq,beq,VLB,VUB,nonlcon）（5）x=fmincon（fun,X0,A,b,Aeq,beq,VLB,VUB,nonlcon,options）（6）x,fval=fmincon（.）（7）x,fval,exitflag=fmincon（.）（8）x,fval,exitflag,output=fmincon（.）输出极值点M文件迭代的初值参数说明变量上下限16注意：

注意：

1fmincon函数提供了大型优化算法和中型优化算法。

默认时，若在fun函数中提供了梯度（options参数的GradObj设置为on），并且只有上下界存在或只有等式约束，fmincon函数将选择大型算法。

当既有等式约束又有梯度约束时，使用中型算法。

2fmincon函数的中型算法使用的是序列二次规划法。

在每一步迭代中求解二次规划子问题，并用BFGS法更新拉格朗日Hessian矩阵。

3fmincon函数可能会给出局部最优解，这与初值X0的选取有关。

171、写成标准形式写成标准形式：

s.t.2x1+3x26s.tx1+4x25x1,x20例例2182、先建立先建立M-文件文件fun3.m:

functionf=fun3（x）;f=-x

（1）-2*x

（2）+（1/2）*x

（1）2+（1/2）*x

（2）2MATLAB（youh2）3、再建立主程序youh2.m：

x0=1;1;A=23;14;b=6;5;Aeq=;beq=;VLB=0;0;VUB=;x,fval=fmincon（fun3,x0,A,b,Aeq,beq,VLB,VUB）4、运算结果为：

运算结果为：

x=0.76471.0588fval=-2.0294191先建立先建立M文件文件fun4.m,定义目标函数定义目标函数:

functionf=fun4（x）;f=exp（x

（1）*（4*x

（1）2+2*x

（2）2+4*x

（1）*x

（2）+2*x

（2）+1）;x1+x2=0s.t.1.5+x1x2-x1-x20-x1x2100例例32再建立再建立M文件文件mycon.m定义非线性约束：

定义非线性约束：

functiong,ceq=mycon（x）g=x

（1）+x

（2）;1.5+x

（1）*x

（2）-x

（1）-x

（2）;-x

（1）*x

（2）-10;203主程序主程序youh3.m为为:

x0=-1;1;A=;b=;Aeq=11;beq=0;vlb=;vub=;x,fval=fmincon（fun4,x0,A,b,Aeq,beq,vlb,vub,mycon）MATLAB（youh3）3.运算结果为运算结果为：

x=-1.22501.2250fval=1.895121例41先建立先建立M-文件文件fun.m定义目标函数定义目标函数:

functionf=fun（x）;f=-2*x

（1）-x

（2）;2再建立再建立M文件文件mycon2.m定义非线性约束：

定义非线性约束：

functiong,ceq=mycon2（x）g=x

（1）2+x

（2）2-25;x

（1）2-x

（2）2-7;223.主程序主程序fxx.m为为:

x0=3;2.5;VLB=00;VUB=510;x,fval,exitflag,output=fmincon（fun,x0,VLB,VUB,mycon2）MATLAB（fxx（fun）234.运算结果为运算结果为:

x=4.00003.0000fval=-11.0000exitflag=1output=iterations:

4funcCount:

17stepsize:

1algorithm:

1x44charfirstorderopt:

cgiterations:

返回返回24应用实例：

应用实例：

供应与选址供应与选址某公司有6个建筑工地要开工，每个工地的位置（用平面坐标系a，b表示，距离单位：

千米）及水泥日用量d（吨）由下表给出。

目前有两个临时料场位于A（5,1），B（2,7），日储量各有20吨。

假设从料场到工地之间均有直线道路相连。

（1）试制定每天的供应计划，即从A，B两料场分别向各工地运送多少吨水泥，使总的吨千米数最小。

（2）为了进一步减少吨千米数，打算舍弃两个临时料场，改建两个新的，日储量各为20吨，问应建在何处，节省的吨千米数有多大？

25

（一）、建立模型

（一）、建立模型记工地的位置为记工地的位置为（ai，bi），水泥日用量为水泥日用量为di，i=1,6;料场位置料场位置为为（xj，yj），日储量为日储量为ej，j=1,2；从料场从料场j向工地向工地i的运送量为的运送量为Xij。

当用临时料场时决策变量为：

Xij，当不用临时料场时决策变量为：

Xij，xj，yj。

26

（二）使用临时料场的情形

（二）使用临时料场的情形使用两个临时料场A（5,1），B（2,7）.求从料场j向工地i的运送量为Xij，在各工地用量必须满足和各料场运送量不超过日储量的条件下，使总的吨千米数最小，这是线性规划问题.线性规划模型为：

设X11=X1,X21=X2,X31=X3,X41=X4,X51=X5,X61=X6X12=X7,X22=X8,X32=X9,X42=X10,X52=X11,X62=X12编写程序gying1.mMATLAB（gying1）27计算结果为：

计算结果为：

x=3.00005.00000.00007.00000.00001.00000.00000.00004.00000.00006.000010.0000fval=136.227528（三）改建两个新料场的情形（三）改建两个新料场的情形改建两个新料场，要同时确定料场的位置（xj,yj）和运送量Xij，在同样条件下使总吨千米数最小。

这是非线性规划问题。

非线性规划模型为：

29设X11=X1,X21=X2,X31=X3,X41=X4,X51=X5,X61=X6X12=X7,X22=X8,X32=X9,X42=X10,X52=X11,X62=X12x1=X13,y1=X14,x2=X15,y2=X16

（1）先编写M文件liaoch.m定义目标函数。

MATLAB（liaoch）

（2）取初值为线性规划的计算结果及临时料场的坐标:

x0=35070100406105127;编写主程序gying2.m.MATLAB（gying2）30（3）计算结果为：

x=3.00005.00000.07077.000000.9293003.929306.000010.07076.38754.39435.75117.1867fval=105.4626exitflag=131（4）若修改主程序gying2.m,取初值为上面的计算结果:

x0=3.00005.00000.07077.000000.9293003.929306.000010.07076.38754.39435.75117.1867得结果为:

x=3.00005.00000.30947.00000.01080.6798003.690605.989210.32025.53694.91945.82917.2852fval=103.4760exitflag=1总的吨千米数比上面