模型04：规划模型.pptxVIP

数 学 建 模规划模型通知放假：清明补课：4月9日（星期六）上午1234节 G101？一、数学规划模型1、规划模型的一般形式例1：选址问题某公司有6个建筑工地，位置坐标为(ai,bi)(单位：公里),水泥日用量ri (单位：吨）i１２３４５６a1.258.750.55.7537.25b1.250.754.7556.57.75r3547611假设：料场和工地之间有直线道路（1）现有2料场，位于A(5,1),B(2,7),记(xj,yj),j=1,2, 日储量qj各有20吨。目标：制定每天的供应计划，即从A, B两料场分别向各工地运送多少吨水泥，使总的吨公里数最小。料场到施工点距离运输吨公里数容量需求解:设(xj,yj) 表示 n 个料场的位置坐标, wij 表示第 j 料场向第 i 施工点的材料运量目标函数 约束条件决策变量?线性规划模型?（2）改建两个新料场，需要确定新料场位置(xj,yj)和运量cij ，在其它条件不变下使总吨公里数最小。目标约束条件wi j，(xj,yj)~16维决策变量?线性规划模型?非线性规划模型Mathematic Modeling规划模型的一般形式决策变量 x =(x1, x2, …, xn )目标函数 Min Z = f (x)约束条件 s.t x ? A (? Rn )约束条件x ? A 一般用等式或不等式方程表示 hi (x1, x2, …, xn ) ? 0 , i=1,2,…,m gj (x1, x2, …, xn ) = 0 , j=1,2,…,l无约束条件数学规划非线性规划NLP线性规划LP整数规划IP连续规划纯整数规划PIP混合规划MIP0-1规划2、数学规划类型划分：表达式变量取值二、线性规划目标函数：决策变量的线性函数 ——代数和数乘 约束条件：决策变量的线性等式或不等式线性规划的一般形式求解图解法: 二元单纯形法,灵敏度分析: 20世纪大型优化算法：Lipsol法数学软件Matlab ……Lindo Lingo : 解规划问题的数学软件Cx2min-2x1+ x2 =2f = x1-x2Dx1+x2 =5x1 - 2x2=2B0x1Ax1-x2=1x1-x2=21、图解法——二元概念:可行解 可行基在图中可以看出C(1,4)∴ min f =-3 x*= C = (1,4)另: min f= x1+x2 例min f = x1-x2 s.t -2x1+ x2 ≤2 x1 - 2x2≤2 x1+ x2 ≤5 x1,x2 ≥02、线性规划： Matlab求解fxAb线性规划基本格式x=linprog(f,A,b)例 min f = x1-x2 s.t: -2x1+ x2 ≤2 x1 - 2x2≤2 x1+ x2 ≤5f=[1 -1];A=[-2 1 1 -2 1 1];b=[2 2 5];x=linprog(f,A,b)l01.m优化信息基本拉格朗日乘子目标退出条件基本变量上下界指定优化参数等式约束初值 相关命令线性规划fAeqxbeqAlbbub格式[x,fval,exitflag,output,lambda]=linprog(f,A,b,Aeq,beq,lb,ub,x0,options)[x,fval] =linprog(f,A,b,Aeq,beq,lb,ub )例min f =-5 x1-4x2 -6x3 s.t x1 - x2 + x3≤20 3x1+2x2 +4x3 ≤42 3x1+2x2 ≤30 x1,x2 ,x3 ≥0 f=[-5 -4 -6];A=[1 -1 1 3 2 4 3 2 0];b=[20 42 30];lb=zeros(3,1);[x,fval]=linprog(f,A,b,[ ],[ ],lb)l02.ml03.ml04.m解例1：选址问题 其中i１２３４５６a1.258.750.55.7537.25b1.250.754.7556.57.75r3547611jABx52y17q2020 注：wij三、非线性规划模型非线性规划有约束、无约束规划连续、离散优化最优解局部、全局Mathematic Modelingyf(x)x求解无约束优化寻找 x* 使得 f(x*)最小局部最优解必要条件充分条件Hessian阵求解方法搜索法：黄金分割法、Fibonacci、随机搜索法、最速下降法梯度法、牛顿法、拟牛顿法、方向加速法、单纯形法、模拟退火法遗传算法、神经网络……Matlab求解1、一元无约束优化模型基本格式全局最优讨论：函数讨论：区间例：[x, fval]= fminbnd (f,x1,x2)字符串 内联函数 m文件 句柄函数f01.m2、多元无约束优化模型基本格式局部最优 例：单变量符号[x, fva