Matlab优化应用.docx
- 文档编号:30417682
- 上传时间:2023-08-14
- 格式:DOCX
- 页数:35
- 大小:78.25KB
Matlab优化应用.docx
《Matlab优化应用.docx》由会员分享,可在线阅读,更多相关《Matlab优化应用.docx(35页珍藏版)》请在冰豆网上搜索。
Matlab优化应用
MATLAB优化应用
§1 线性规划模型
一、线性规划课题:
实例1:
生产计划问题
假设某厂计划生产甲、乙两种产品,现库存主要材料有A类3600公斤,B类2000公斤,C类3000公斤。
每件甲产品需用材料A类9公斤,B类4公斤,C类3公斤。
每件乙产品,需用材料A类4公斤,B类5公斤,C类10公斤。
甲单位产品的利润70元,乙单位产品的利润120元。
问如何安排生产,才能使该厂所获的利润最大。
建立数学模型:
设x1、x2分别为生产甲、乙产品的件数。
f为该厂所获总润。
maxf=70x1+120x2
s.t 9x1+4x2≤3600
4x1+5x2≤2000
3x1+10x2≤3000
x1,x2≥0
实例2:
投资问题
某公司有一批资金用于4个工程项目的投资,其投资各项目时所得的净收益(投入资金锪百分比)如下表:
工程项目收益表
工程项目
A
B
C
D
收益(%)
15
10
8
12
由于某种原因,决定用于项目A的投资不大于其他各项投资之和而用于项目B和C的投资要大于项目D的投资。
试确定全文该公司收益最大的投资分配方案。
建立数学模型:
设x1、 x2 、x3 、x4分别代表用于项目A、B、C、D的投资百分数。
maxf=0.15x1+0.1x2+0.08x3+0.12x4
s.t x1-x2-x3-x4≤0
x2+x3-x4≥0
x1+x2+x3+x4=1
xj≥0 j=1,2,3,4
实例3:
运输问题
有A、B、C三个食品加工厂,负责供给甲、乙、丙、丁四个市场。
三个厂每天生产食品箱数上限如下表:
工厂
A
B
C
生产数
60
40
50
四个市场每天的需求量如下表:
市场
甲
乙
丙
丁
需求量
20
35
33
34
从各厂运到各市场的运输费(元/每箱)由下表给出:
市 场
甲
乙
丙
丁
工
厂
A
2
1
3
2
B
1
3
2
1
C
3
4
1
1
求在基本满足供需平衡的约束条件下使总运输费用最小。
建立数学模型:
设aij为由工厂i运到市场j的费用,xij 是由工厂i运到市场j的箱数。
bi是工厂i的产量,dj是市场j的需求量。
b=(604050) d=(20353334)
s.t
x ij≥0
当我们用MATLAB软件作优化问题时,所有求maxf 的问题化为求min(-f)来作。
约束g i (x)≥0,化为 –g i≤0来作。
上述实例去掉实际背景,归结出规划问题:
目标函数和约束条件都是变量x的线性函数。
形如:
(1) minf T X
s.t A X≤b
AeqX=beq
lb≤X≤ub
其中X为n维未知向量,f T=[f1,f2,…fn]为目标函数系数向量,小于等于约束系数矩阵A为m×n矩阵,b为其右端m维列向量,Aeq为等式约束系数矩阵,beq为等式约束右端常数列向量。
lb,ub为自变量取值上界与下界约束的n维常数向量。
二.线性规划问题求最优解函数:
调用格式:
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)返回值x为最优解向量。
x=linprog(f,A,b,Aeq,beq) 作有等式约束的问题。
若没有不等式约束,则令A=[]、b=[] 。
x=linprog(f,A,b,Aeq,beq,lb,ub,x0,options) 中lb,ub为变量x的下界和上界,x0为初值点,options为指定优化参数进行最小化。
Options的参数描述:
Display 显示水平。
选择’off’ 不显示输出;选择’iter’显示每一 步迭代过程的输出;选择’final’ 显示最终结果。
MaxFunEvals 函数评价的最大允许次数
Maxiter 最大允许迭代次数
TolX x处的终止容限
[x,fval]=linprog(…) 左端 fval 返回解x处的目标函数值。
[x,fval,exitflag,output,lambda]=linprog(f,A,b,Aeq,beq,lb,ub,x0) 的输出部分:
exitflag 描述函数计算的退出条件:
若为正值,表示目标函数收敛于解x处;若为负值,表示目标函数不收敛;若为零值,表示已经达到函数评价或迭代的最大次数。
output 返回优化信息:
output.iterations表示迭代次数;output.algorithm表示所采用的算法;outprt.funcCount表示函数评价次数。
lambda 返回x处的拉格朗日乘子。
它有以下属性:
lambda.lower-lambda的下界;
lambda.upper-lambda的上界;
lambda.ineqlin-lambda的线性不等式;
lambda.eqlin-lambda的线性等式。
三. 举例
例1:
求解线性规划问题:
maxf=2x1+5x2
s.t
先将目标函数转化成最小值问题:
min(-f)=-2x1-5x2
程序:
f=[-2-5];
A=[10;01;12];
b=[4;3;8];
[x,fval]=linprog(f,A,b)
f=fval*(-1)
结果:
x= 2
3
fval=-19.0000
maxf= 19
例2:
minf=5x1-x2+2x3+3x4-8x5
s.t –2x1+x2-x3+x4-3x5≤6
2x1+x2-x3+4x4+x5≤7
0≤xj≤15 j=1,2,3,4,5
程序:
f=[5-123-8];
A=[-21-11-3;21-141];
b=[6;7];
lb=[00000];
ub=[1515151515];
[x,fval]=linprog(f,A,b,[],[],lb,ub)
结果:
x=
0.0000
0.0000
8.0000
0.0000
15.0000
minf=
-104
例3:
求解线性规划问题:
minf=5x1+x2+2x3+3x4+x5
s.t –2x1+x2-x3+x4-3x5≤1
2x1+3x2-x3+2x4+x5≤-2
0≤xj≤1 j=1,2,3,4,5
程序:
f=[51231];
A=[-21-11-3;23-121];
b=[1;-2];
lb=[00000];
ub=[11111];
[x,fval,exitflag,output,lambda]=linprog(f,A,b,[],[],lb,ub) 运行结果:
Exiting:
Oneormoreoftheresiduals,dualitygap,ortotalrelativeerror
hasgrown100000timesgreaterthanitsminimumvaluesofar:
theprimalappearstobeinfeasible(andthedualunbounded).
(Thedualresidual x= 0.0000 0.0000 1.1987 0.0000 0.0000 fval= 2.3975 exitflag= -1 output= iterations: 7 cgiterations: 0 algorithm: 'lipsol' lambda= ineqlin: [2x1double] eqlin: [0x1double] upper: [5x1double] lower: [5x1double] 显示的信息表明该问题无可行解。 所给出的是对约束破坏最小的解。 例4: 求解实例1的生产计划问题 建立数学模型: 设x1、x2分别为生产甲、乙产品的件数。 f为该厂所获总润。 maxf=70x1+120x2 s.t 9x1+4x2≤3600 4x1+5x2≤2000 3x1+10x2≤3000 x1,x2≥0 将其转换为标准形式: minf=-70x1-120x2 s.t 9x1+4x2≤3600 4x1+5x2≤2000 3x1+10x2≤3000 x1,x2≥0 程序: f=[-70-120]; A=[94;45;310]; b=[3600;2000;3000]; lb=[00]; ub=[]; [x,fval,exitflag]=linprog(f,A,b,[],[],lb,ub) maxf=-fval 结果: x= 200.0000 240.0000 fval= -4.2800e+004 exitflag= 1 maxf= 4.2800e+004 例5: 求解实例2 建立数学模型: maxf=0.15x1+0.1x2+0.08x3+0.12x4 s.t x1-x2-x3-x4≤0 x2+x3-x4≥0 x1+x2+x3+x4=1 xj≥0 j=1,2,3,4 将其转换为标准形式: minz=-0.15x1-0.1x2-0.08x3-0.12x4 s.t x1-x2-x3-x4≤0 -x2-x3+x4≤0 x1+x2+x3+x4=1 xj≥0 j=1,2,3,4 程序: f=[-0.15;-0.1;-0.08;-0.12]; A= [1-1-1-1 0-1-11]; b=[0;0]; Aeq=[1111]; beq=[1]; lb=zeros(4,1); [x,fval,exitflag]=linprog(f,A,b,Aeq,beq,lb) f=-fval 结果: x= 0.5000 0.2500 0.0000 0.2500 fval= -0.1300 exitflag= 1 f= 0.1300 即4个项目的投资百分数分别为50%,25%,0, 25%时可使该公司获得最大的收益,其最大收益可到达13%。 过程正常收敛。 例6: 求解实例3 建立数学模型: 设aij为由工厂i运到市场j的费用,xij 是由工厂i运到市场j的箱数。 bi是工厂i的产量,dj是市场j的需求量。 b=(604050)T d=(20353334)T s.t x ij≥0 程序: A=[2132;1321;3411]; f=A(: ); B=[100100100100 010010010010 001001001001]; D=[111000000000 000111000000 000000111000 000000000111]; b=[60;40;50]; d=[20;35;33;34]; lb=zeros(12,1); [x,fval,exitflag]=linprog(f,B,b,D,d,lb) 结果: x= 0.0000 20.0000 0.0000 35.0000 0.0000 0.0000 0.0000 0.0000 33.0000 0.0000 18.4682 15.5318 fval= 122.0000 exitflag= 1 即运输方案为: 甲市场的货由B厂送20箱;乙市场的货由A厂送35箱;丙商场的货由C厂送33箱;丁市场的货由B厂送18箱,再由C厂送16箱。 最低总运费为: 122元。 §2 非线性规划模型 一.非线性规划课题 实例1 表面积为36平方米的最大长方体体积。 建立数学模型: 设x、y、z分别为长方体的三个棱长,f为长方体体积。 maxf=xy(36-2xy)/2(x+y) 实例2 投资决策问题 某公司准备用5000万元用于A、B两个项目的投资,设x1、x2分别表示配给项目A、B的投资。 预计项目A、B的年收益分别为20%和16%。 同时,投资后总的风险损失将随着总投资和单位投资的增加而增加,已知总的风险损失为2x12+x22+(x1+x2)2.问应如何分配资金,才能使期望的收益最大,同时使风险损失为最小。 建立数学模型: maxf=20x1+16x2-λ[2x12+x22+(x1+x2)2] s.t x1+x2≤5000 x 1≥0,x2≥0 目标函数中的λ≥0是权重系数。 由以上实例去掉实际背景,其目标函数与约束条件至少有一处是非线性的,称其为非线性问题。 非线性规划问题可分为无约束问题和有约束问题。 实例1为无约束问题,实例2为有约束问题。 二.无约束非线性规划问题: 求解无约束最优化问题的方法主要有两类: 直接搜索法(Searchmethod)和梯度法(Gradientmethod). 1.fminunc函数 调用格式: x=fminunc(fun,x0) x=fminunc(fun,x0,options) x=fminunc(fun,x0,options,P1,P2) [x,fval]=fminunc(…) [x,fval,exitflag]=fminunc(…) [x,fval,exitflag,output]=fminunc(…) [x,fval,exitflag,output,grad]=fminunc(…) [x,fval,exitflag,output,grad,hessian]=fminunc(…) 说明: fun为需最小化的目标函数,x0为给定的搜索的初始点。 options指定优化参数。 返回的x为最优解向量;fval为x处的目标函数值;exitflag描述函数的输出条件;output返回优化信息;grad返回目标函数在x处的梯度。 Hessian返回在x处目标函数的Hessian矩阵信息。 例1 : 求 程序: 编辑ff1.m文件 functionf=ff1(x) f=8*x (1)-4*x (2)+x (1)^2+3*x (2)^2; 通过绘图确定一个初始点: [x,y]=meshgrid(-10: .5: 10); z=8*x-4*y+x.^2+3*y.^2; surf(x,y,z) 选初始点: x0=(0,0) x0=[0,0]; [x,fval,exitflag]=fminunc(@ff1,x0) 结果: x= -4.0000 0.6667 fval= -17.3333 exitflag= 1 例2: 程序: 编辑ff2.m文件: functionf=ff2(x) f=4*x (1)^2+5*x (1)*x (2)+2*x (2)^2; 取初始点: x0=(1,1) x0=[1,1]; [x,fval,exitflag]=fminunc(@ff2,x0) 结果: x= 1.0e-007* -0.1721 0.1896 fval= 2.7239e-016 exitflag= 1 例3: 将上例用提供的梯度g最小化函数进行优化计算。 修改M文件为: function[f,g]=ff3(x) f=4*x (1)^2+5*x (1)*x (2)+2*x (2)^2; ifnargut>1 g (1)=8*x (1)+5*x (2); g (2)=5*x (1)+4*x (2); end 通过下面将优化选项结构options.GradObj设置为’on’来得到梯度值。 options=optimset(‘Gradobj’,’on’); x0=[1,1]; [x,fval,exitflag]=fminunc(@ff3,x0,options) 结果: x= 1.0e-015* -0.2220 -0.2220 fval= 5.4234e-031 exitflag=
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 优化 应用