预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共28页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

第3章Matlab在最优化问题中的应用 优化理论是一门实践性很强的学科,广泛应用于生产管理、军事指挥和科学试验等各种领域,Matlab优化工具箱提供了对各种优化问题的一个完整的解决方案。 在数学上,所谓优化问题,就是求解如下形式的最优解: Minfun(x) Sub.to[C.E.] [B.C.] 其中fun(x)称为目标函数,“Sub.to”为“subjectto”的缩写,由其引导的部分称为约束条件。[C.E.]表示ConditionEquations,即条件方程,可为等式方程,也可为不等式方程。[B.C.]表示BoundaryConditions,即边界条件,用来约束自变量的求解域,以lb≤x≤ub的形式给出。当[C.E.]为空时,此优化问题称为自由优化或无约束优化问题;当[C.E.]不空时,称为有约束优化或强约束优化问题。 在优化问题中,根据变量、目标函数和约束函数的不同,可以将问题大致分为: ·线性优化目标函数和约束函数均为线性函数。 ·二次优化目标函数为二次函数,而约束条件为线性方程。线性优化和二次优化统称为简单优化。 ·非线性优化目标函数为非二次的非线性函数,或约束条件为非线性方程。 ·多目标优化目标函数并非一个时,称为多目标优化问题。 本章将对以上几类优化问题在Matlab中的实现作比较详细的讲解。另外还将介绍两个利用优化方法解非线性方程的函数。 通过本章的介绍,用户可以不必掌握艰涩的各种优化算法而轻易地解决一些常用的最优化问题了。 3.1线性规划问题 线性规划问题即目标函数和约束条件均为线性函数的问题。 其标准形式为: minC’x sub.ToAx=b x≥0 其中C,b,0∈Rn,A∈Rmn,均为数值矩阵,x∈Rn。 若目标函数为:maxC’x,则转换成:min–C’x。 标准形式的线性规划问题简称为LP(LinearProgramming)问题。其它形式的线性规划问题经过适当的变换均可以化为此种标准形。线性规划问题虽然简单,但在工农业及其他生产部门中应用十分广泛。 在Matlab中,线性规划问题由linprog函数求解。 函数:linprog%求解如下形式的线性规划问题: suchthat 其中f,x,b,beq,lb,ub为向量,A,Aeq为矩阵。 格式:x=linprog(f,A,b) x=linprog(f,A,b,Aeq,beq) x=linprog(f,A,b,Aeq,beq,lb,ub) x=linprog(f,A,b,Aeq,beq,lb,ub,x0) x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options) [x,fval]=linprog(...) [x,fval,exitflag]=linprog(...) [x,fval,exitflag,output]=linprog(...) [x,fval,exitflag,output,lambda]=linprog(...) 说明: x=linprog(f,A,b)求解问题minf’*x,约束条件为A*x<=b。 x=linprog(f,A,b,Aeq,beq)求解上面的问题,但增加等式约束,即Aeq*x=beq。若没有不等式存在,则令A=[]、b=[]。 x=linprog(f,A,b,Aeq,beq,lb,ub)定义设计变量x的下界lb和上界ub,使得x始终在该范围内。若没有等式约束,令Aeq=[]、beq=[]。 x=linprog(f,A,b,Aeq,beq,lb,ub,x0)设置初值为x0。该选项只适用于中型问题,默认时大型算法将忽略初值。 x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options)用options指定的优化参数进行最小化。 [x,fval]=linprog(...)返回解x处的目标函数值fval。 [x,fval,exitflag]=linprog(...)返回exitflag值,描述函数计算的退出条件。 [x,fval,exitflag,output]=linprog(...)返回包含优化信息的输出变量output。 [x,fval,exitflag,output,lambda]=linprog(...)将解x处的Lagrange乘子返回到lambda参数中。 exitflag参数 描述退出条件: ·>0表示目标函数收敛于解x处; ·=0表示已经达到函数评价或迭代的最大次数; ·<0表示目标函数不收敛。 output参数 该参数包含下列优化信息: ·output.iterations迭代次数; ·output.cgiterationsPCG迭代次数(只适用于大型规划问题); ·output.algorithm所