《MATLAB求解线性规划问题.pptx》由会员分享，可在线阅读，更多相关《MATLAB求解线性规划问题.pptx（25页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、一、实验目的：第一节 线性方程组的应用1、了解线性规划问题及可行解、最优解的概念；2、掌握Matlab软件关于求解线性规划的语句和方法。二、实验原理和方法：在生活实践中，很多重要的实际问题都是线性的（至少能够用线性函数很好的近似表示），所以我们一般把这些问题化为线性的目标函数和约束条件进行分析，通常将目标函数和约束都是线性表达式的规划问题称为线性规划。第1页/共25页它的一般形式是：也可以用矩阵形式来表示：第2页/共25页线性规划的可行解是满足约束条件的解；线性规划的最优解是使目标函数达到最优的可行解。线性规划关于解的情况可以是：1、无可行解，即不存在满足约束条件的解；2、有唯一最优解，即在可
2、行解中有唯一的最有解；4、有可行解，但由于目标函数值无界而无最优解。3、有无穷最优解，即在可行解中有无穷个解都可使目 标函数达到最优；第3页/共25页 一般求解线性规划的常用方法是单纯形法和改进的单纯形法，这类方法的基本思路是先求得一个可行解，检验是否为最优解；若不是，可用迭代的方法找到另一个更优的可行解，经过有限次迭代后，可以找到可行解中的最优解或者判定无最优解。第4页/共25页三、内容与步骤：在Matlab优化工具箱中，linprog函数是使用单纯形法求解下述线性规划问题的函数。第5页/共25页它的命令格式为：其中：A为约束条件矩阵，b,c分别为目标函数的系数向量和约束条件中最右边的数值向
3、量；也可设置解向量的上界vlb和下界vub，即解向量必须满足vlb=x=vub；还可预先设置初始解向量x0。如没有不等式，而只有等式时，A=,b=;输出的结果：x表示最优解向量；fval表示最优值。第6页/共25页【例 1】求解线性规划问题：解：考虑到linprog函数只解决形如的线性规划。所以先要将线性规划变为如下形式：第7页/共25页然后建立M文件如下：c=-3;1;1;A=1-21;4-1-2;b=11;-3;aeq=20-1;beq=-1;vlb=0;0;0;x,fval=linprog(c,A,b,aeq,beq,vlb)Matlab程序：第8页/共25页以ch701作为文件名保存此
4、M文件后，在命令窗口输入ch701后即可得到结果：同时返回fval=-2对应到原来的线性规划中即知目标函数的最大值为2，此时x1=4,x2=1,x3=9。第9页/共25页第二节 无约束规划计算方法一、实验目的1、了解无约束规划问题的求解原理与方法；2、会用Matlab软件求解无约束规划问题。二、实验原理和方法无约束规划问题的解法一般按目标函数的形式分为两大类：一类是一元函数的一维搜索法，如黄金分割法、插值法等；另一类是求解多元函数的下降迭代法。第10页/共25页迭代的基本思想和步骤大致可分为以下四步：第11页/共25页三、实验内容与步骤在Matlab软件中，求解无约束规划的常用命令是：x=fm
5、inunc(fun,x0)其中，fun函数应预先定义到M文件中，并设置初始解向量为x0。第12页/共25页【例 2】求解取解：首先建立函数文件第13页/共25页 以fun702为文件名保存此函数文件。在命令窗口输入：x0=-2;4;x=fminunc(fun702,x0)结果显示：f=x=即极小值为-1，是x1=1,x2=1时取得。Matlab程序：第14页/共25页【例 3】解非线性方程组解：解此非线性方程组等价于求解无约束非线性规划问题：然后建立函数文件fun703.m 第15页/共25页在命令窗口输入：x0=0;0;x=fminunc(fun703,x0)结果显示：则非线性方程组的解为。
6、Matlab程序：第16页/共25页第三节 约束非线性规划计算方法 一、实验目的1、了解约束非线性规划问题的求解原理与方法；2、会用Matlab软件求解约束非线性规划问题。二、实验原理和方法对于约束非线性规划，随着目标函数和约束条件的不同，解法也不同，一般来说，有两类方法：（1）、将约束问题化为无约束问题的求解方法；（2）、用线性规划来逼近非线性规划；第17页/共25页三、实验内容与步骤 约束非线性规划的一般形式为：其中，f(x)为多元实值函数;g(x)为向量函数,并且f(x),g(x)中至少有一个函数是非线性函数的（否则成为线性规划问题）。第18页/共25页x=fmincon(fun,x0,
7、A,b)x=fmincon(fun,x0,A,b,Aeq,beq)x=fmincon(fun,x0,A,b,Aeq,beq,lb,ub)x=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)在Matlab优化工具箱中，fmincon函数是用SQP算法来解决一般的约束非线性规划的函数，它的命令格式为：第19页/共25页【例 4】求解约束非线性规划：(初值为1;1)解:首先建立一个m文件functiony=fun7041(x)y=-exp(x(1)*x(2)2*(3-exp(x(1)-x(2)2);存储为首先将问题转化为matlab要求的格式;即求出fun,A,b,
8、Aeq,Beq,X0,Lb,Ub第20页/共25页然后建立一个m文件functionc,cep=fun7042(x)c=;%c为非线性不等式,且为c=0cep=exp(x(1)+x(2)2-3;%cep为非线性等式然后存储为最后在命令窗口中输入：A=;b=;Aeq=;Beq=;Lb=;Ub=;x,f=fmincon(fun7041,1;1,fun7042)-f因题目中有非线性约束条件，所以建立非线性约束m-文件。Matlab程序：第21页/共25页结果为:ans最后的结果为:e-016第22页/共25页【例 5】求解约束非线性规划：解：首先建立一个m文件functiony=fun705(x)y=(x(1)-1)2+(x(2)-2)2+(x(3)-3)2+(x(4)-4)2;存储为文件.第23页/共25页x0=1;1;1;1;A=1111;3321;B=5;10;Aeq=;Beq=;Lb=0;0;0;0;x,g=fmincon(fun705,x0,A,B,Aeq,Beq,Lb)答案为:x=g=Matlab程序：第24页/共25页感谢您的观看！第25页/共25页