原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
最优化大作业 张俊飞
最优化大作业
班级：021051
学号
姓名：张俊飞
时间：2013 年1 月1 日
最优化大作业 张俊飞
最优化上机作业
一、分别用梯度下降法和牛顿法求解下列问题：
min f (x) (x 10x )2 5(x x )2 (x 2x )4 10(x x )4
1 2 3 4 2 3 1 4
其中一维搜索可用0.618 法
1.1 最速下降法
解无约束问题的的最速下降法计算步骤
0
第1 步 选取初始点 ，给定终止误差 ，令k ：=0 ；
x 0
k
f (xk ) f (xk ) x
第2 步 计算 ，若‖ ‖ ，停止迭代,输出 ,否则进行
第三步；
k k
p f (x )
第3 步 取 =− ；
第4 步 进行一维搜索，求 ，使得
k
f ( k k k k
x p ) m i nf (x p )
k
0
令xk 1 xk pk ，k:=k+1 ，转第2 步.
k
算法流程图如下：
第 1 页
最优化大作业 张俊飞
Matlab 代码如下
tic
clc
e=1.0e-6;
syms x1 x2 x3 x4 a1
f=(x1+10*x2)^2+5*(x3-x4)^2+(x2-2*x3)^4+10*(x1-x4)^4;
v=[x1,x2,x3,x4];
g=jacobian(f,v);%求梯度
y = input([x1,x2,x3,x4] = );
x1 = y(1);
x2 = y(2);
x3 = y(3);
x4 = y(4);
第 2 页
最优化大作业 张俊飞
x=[x1,x2,x3,x4];
i=0;%迭代次数
p=-g;
while norm(subs(p,v,x))=0.04
norm(subs(p,v,x))
i=i+1;
x1=x+a1*subs(p,v,x);
h=subs(f,v,x1);
a=gold(h,0,100);
x=x+a*subs(p,v,x);
end
disp(sprintf(迭代次数为: %d次,i));
disp(最优的x为);
x
disp(最优解为:);
subs(f,v,x)
toc
其中一维搜索用0.618 法
function min_1=gold(h,a,b)
f
文档评论（0）