最优控制程序.docxVIP

上机实验报告 课程名称：静态优化算法及软件实现 课程时间：2014年夏 授课教师：李传江 学生姓名： 学生班级：1236104 学生学号： 所在组别：Group 62 联系电话： 选作附加题：是 哈尔滨工业大学航天学院 2014年6月 问题 = 1 \* CHINESENUM3 一 MATLAB解答过程 程序代码： c=[-2,1,-5,3,-4]; A=[3,-2,7,-5,4;1,-1,2,-4,2]; b=[6;0]; [x,fval,exitflag,output]=bintprog(c,A,b) 运行结果： Optimization terminated. x = 0 0 1 1 1 fval = -6 exitflag = 1 output = iterations: 14 nodes: 5 time: 0.0936 algorithm: LP-based branch-and-bound branchStrategy: maximum integer infeasibility nodeSrchStrategy: best node search message: Optimization terminated. MATHEMATICA解答过程 Maximize程序代码： Maximize2x-y+5z-3u+4v,3x-2y+7z-5u+4v?6x-y+2z-4u+2v?00?x?10?y?10?z?10?u?10?v?1x,y,z,u,v∈Integers,x,y,z,u,v 运行结果： 6,x→0,y→0,z→1,u→1,v→1 Minimize程序代码： Minimize-2x+y-5z+3u-4v,3x-2y+7z-5u+4v?6x-y+2z-4u+2v?00?x?10?y?10?z?10?u?10?v?1x,y,z,u,v∈Integers,x,y,z,u,v 运行结果： -6,x→0,y→0,z→1,u→1,v→1 NMaximize程序代码： NMaximize2x-y+5z-3u+4v,3x-2y+7z-5u+4v?6x-y+2z-4u+2v?00?x?10?y?10?z?10?u?10?v?1x,y,z,u,v∈Integers,x,y,z,u,v 运行结果： 6.,x→0,y→0,z→1,u→1,v→1 NMinimize程序代码： NMinimize-2x+y-5z+3u-4v,3x-2y+7z-5u+4v?6x-y+2z-4u+2v?00?x?10?y?10?z?10?u?10?v?1x,y,z,u,v∈Integers,x,y,z,u,v 运行结果： -6.,x→0,y→0,z→1,u→1,v→1 总结：由上面的计算过程我们可以发现MATHEMATICA相对于MATLAB的使用更加简洁，但是MATLAB可以根据需要返回详细的计算迭代过程，为我们进一步理解优化过程提供了便利。MATLAB只能解决0-1整数规划问题，无法解决纯整数规划问题或混合整数规划问题。而MATHEMATICA则可以解决以上问题，并且MATHEMATICA中提供了多种函数供我们选择：Minimize，Maximize，NMaximize，NMinimize，FindMinimum等。当输入给定为精确条件时，Minimize，Maximize返回精确解，NMaximize，NMinimize则会保留一定位的有效数字，返回近似解，FindMinimum同样会对结果保留一定有效数字，并且还可以指定初始点。因为题目中函数以及约束条件均为线性，所以返回结果为全局最优量。 问题 = 2 \* CHINESENUM3 二 MATLAB解答过程：给定初始条件x0=[1;1] A=[-1,-1];b=[-2]; x0=[1;1]; options=optimset(Display,iter); fun=@(x)(1-x(1))*(1-x(1))-10*(x(2)-x(1)*x(1))*(x(2)-x(1)*x(1))+x(1)*x(1)-2*x(1)*x(2); [x, fmin, exitflag, output, lambda] = fmincon (fun,x0,A,b,[],[],[],[],@confun,options) M文件定义函数： function [c ceq] = confun(x) c=[x(1)*x(1)+x(2)*x(2)-16;(x(1)-x(2))*(