原创力文档- 1、本文档共20页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
最优化方法源程序1最优化方法源程序1
最优化方法源程序
一、线性规划问题
例 解线性规划函数linprog ()
线性规划规划矩阵形式：
格式为:
x= linprog (c,A,B,Aeq,Beq,LB,UB)
例 解线性规划
即
%线性规划求解
clear;
c=[-100,-200];
A=[1,1;1,0];
B=[500,200];
Aeq=[2,6];
Beq=1200;
LB=[0,0];UB=[];
[x,fval,exitflag,output]=linprog(c,A,B,Aeq,Beq,LB,UB)
执行结果：
Optimization terminated.
x =
200.0000
133.3333
fval =
-4.6667e+004
exitflag =
1
output =
iterations: 4
algorithm: large-scale: interior point
cgiterations: 0
message: Optimization terminated.
二、无约束一元非线性规划问题
例1 0.618法
function [s,phis,k,G,E]=golds(phi,a,b,delta,epsilon)
%功能: 0.618法精确线搜索
%输入: phi是目标函数, a, b 是搜索区间的两个端点
% delta, epsilon分别是自变量和函数值的容许误差
%输出: s, phis分别是近似极小点和极小值, G是nx4矩阵,
% 其第k行分别是a,p,q,b的第k次迭代值[ak,pk,qk,bk],
% E=[ds,dphi], 分别是s和phis的误差限.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
t=(sqrt(5)-1)/2; h=b-a; phia=feval(phi,a); phib=feval(phi,b);
p=a+(1-t)*h; q=a+t*h; phip=feval(phi,p); phiq=feval(phi,q);
k=1; G(k,:)=[a, p, q, b];
while(abs(phib-phia)epsilon)|(hdelta)
if(phipphiq)
b=q; phib=phiq; q=p; phiq=phip;
h=b-a; p=a+(1-t)*h; phip=feval(phi,p);
else
a=p; phia=phip; p=q; phip=phiq;
h=b-a; q=a+t*h; phiq=feval(phi,q);
end
k=k+1; G(k,:)=[a, p, q, b];
end
ds=abs(b-a); dphi=abs(phib-phia);
if(phip=phiq)
s=p; phis=phip;
else
s=q; phis=phiq;
end
E=[ds,dphi];
程序调用:
[s,phis,k,G,E]=golds(inline(2*x^2-x-1),-1,1,0.16,3)
s =
0.236067977499790
phis =
-1.124611797498107
k =
7
G =
-1.000000000000000 -0.236067977499790 0.236067977499790 1.000000000000000
-0.236067977499790 0.236067977499790 0.527864045000421 1.000000000000000
-0.236067977499790 0.055728090000841 0.236067977499790 0.527864045000421
0.055728090000841 0.236067977499790 0.347524157501472 0.527864045000421
0.055728090000841 0.167184270002524 0.236067977499790 0.347524157501472
0.167184270002524 0.236067977499790 0.278640450004206 0.347524157501
文档评论（0）