最优化课程设计Word格式文档下载.docx

最优化课程设计Word格式文档下载.docx

《最优化课程设计Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《最优化课程设计Word格式文档下载.docx（10页珍藏版）》请在冰豆网上搜索。

报告上交时间：

课程设计结束时上交报告；

7.严禁抄袭行为，一旦发现，课程设计成绩为不及格。

一、运用共轭梯度法求解无约束最优化问题

要求：

1）了解求解无约束最优化问题的共轭梯度法；

2）绘出程序流程图；

3）编写求解无约束最优化问题的共轭梯度法MATLAB?

序；

4）利用编写文件求解某无约束最优化问题；

5）给出程序注释。

指导教师年月曰

负责教师年月曰

学生签字年月曰

沈阳航空航天大学

课程设计成绩评定单

课程名称最优化理论与算法课程设计院（系）理学院专业信息与计算科学课程设计题目共轭梯度法算法分析与实现学号2011041401011姓名黄中武指导教师评语：

课程设计成绩

指导教师签字

年月曰

最优化方法课程设计沈阳航空航天大学课程设计用纸目录

一、正

文1

二、总结

8参考文

献9附

录10

第I页

最优化方法课程设计沈阳航空航天大学课程设计用纸正文

一、正文

一无约束最优化问题的共轭梯度法

共轭梯度法最初是由Hesteness和Stiefel于1952年为求解线形方程组而提出的。

后来，人们把这种方法用于求解无约束最优化问题，使之成为一种重要的最优化方法。

下面，重点介绍Fletcher-Reeves共轭梯度法，简称FR法。

共轭梯度法的基本思想是把共轭性与最速下降法相结合，利用已知点处的梯度构造一组共轭方向，并沿这组方向进行搜索，求出目标函数的极小点。

根据共轭梯度方向的基本性质，这种方法具有二次终止性。

首先讨论对于二次函数的共轭梯度法，然后再把这种方法推广到极小化一般函数的情形。

考虑问题

TTmifx（）,，，0.5xAxbxcn

其中A是对称正定矩阵，c是常数。

具体求解方法如下：

⑴首先，任意给定一个初始点，计算出目变函数在这点的梯度，fx（）x

g若=0,则停止计算，否则，令1

（1）

（1）dfxg,,,,,（）1

（1）

（2）

（2）,0沿方向搜索，得到点，计算在处的梯度，若g，贝闲用dxx2

（1）

（2）

（2）和构造第二个搜索方向，再沿搜索。

gddd2

（）k（）k（）k（）k一般地，若已知点和搜索方向，贝以出发，沿进行xdxd搜索,

得到

（1）（）（）kkk，xxd,，,k

其中步长满足,k

（）（）kk（）（）kkfxd（），,,fxd（），,mink

此时可求出的显式表达。

令,k

（）（）kk,,,（）（）,，fxd

,（）求得极小点，令

'

（1）（）kTk,,,（）（）O,,,fxd

第1页

最优化方法课程设计沈阳航空航天大学课程设计用纸正文根据二次函数的梯度表达式，即

（1）（）kTk，（）OAxbd，,

T（）（）kkkAxdbd（+）O,，,,,k

T（）（）kk（10.3.16）gAdd，,,0,,（）kk

由（10.3.16）式得到

Tk（）gd（）（）kTkk（10.3.17）,,,dAdk

（1）k，计算在处的梯度，若，则停止计算；

否则，用fx（）g,0x,gk，1k，1

（）k

（1）k，

（1）k，（）k和构造下一个搜索方向，并使和关于A共轭。

按此设dddd

想，令

（1）（）kk，dgd,,，,（10.3.18）kk，1

（）kT上式两端左乘，并令dA

（）

（1）（）（）（）kTkkTkTk，dAddAgdAd,,，,,0kk，1

由此得到

（）kT（）（）kTk,,dAg（10.3.19）dAd，kk1

（1）k，

（1）k，再从出发，沿方向搜索xd

综上分析，在第一个搜索方向取负梯度的前提下，重复使用公式（10.3.14），（10.3.17），（10.3.18）和（10.3.19），就能伴随计算点的增加，构造出一组搜索方向。

第2页

最优化方法课程设计沈阳航空航天大学课程设计用纸正文二程序流程图

开始

输入初始点X0

精度e

继续迭代

满足精度要N

求

Y

输出迭代结果

结束

图一共轭梯度法程序流程图

三共轭梯度法的MATLAB?

序

function[x,val,k]=frcg（fun,gfun,xO）

maxk=5000;

rho=0.6;

sigma=0.4;

k=0;

epsilon=1e-7;

n=length（x0）;

while（kvmaxk）

g=feval（gfun,x0）;

第3页

itern=k-（n+1）*floor（k/（n+1））;

itern=iterin+1;

if（itern==1）

d=-g;

else

beta=（g'

*g）/（gO'

*gO）;

d=-g+beta*dO;

gd=g'

*d;

if（gd>

=0）

end

if（norm（g）<

epsilon）

break;

m=0;

mk=0;

while（m<

20）

if

（feval（fun,x0+rh°

Am*d）vfeval（fun,x0）+sigma*rh°

Am*g'

*d）

mk=m;

m=m+1;

x0=x0+rhoAmk*d;

val=feval（fun,x0）;

g0=g;

dO=d;

k=k+1;

第4页

最优化方法课程设计沈阳航空航天大学课程设计用纸正文x=xO;

val=feval（fun,x）;

functionf=fun（x）

f=100*（x

（1）A2-x

（2））A2+（x

（1）-1）A2;

functiong=gfun（x）

g=[400*x

（1）*（x

（1）A2-x

（2））+2*（x

（1）-1）,-200*（x（1F2-x

（2））]'

四利用所编程序求解实际问题

222mifxxxx（）100*（）

（1）,,，,例n121

**Txfx,,（1,1）,（）0该问题有精确解

利用所编程序求解此问题，在命令窗口输入：

x0=[-1.2叮‘；

[工，iral,k]=frce■fun'

'

rl'

unTxO）

运行结果为：

1.DCOO

I.0000

val=

乩7224^017

it-

第5页

最优化方法课程设计沈阳航空航天大学课程设计用纸正文五程序注释

function[x,val,k]=frcg（fun,gfun,x0）

%用FR共轭梯度法求解无约束问题:

minf（x）%输入:

x0时初始点，fun，gfun分别是目标函数和梯度％输出:

x，val分别是近似最优点和最优值，k是迭代次数。

%最大迭代次数

%计算梯度

itern=itern+1;

%计算搜索方向

epsilon）%检验终止条件

第6页

（feval（fun,x0+rhoAm*d）<

feval（fun,x0）+sigma*rhoAm*g'

val=feval（fun,xO）;

go=g;

x=xO;

f=100*（x

（1）A2-x

（2））A2+（x

（1）-1F2;

%需要求解的函数

g=[400*x

（1）*（x（1F2-x

（2））+2*（x

（1）-1）,-200*（x

（1）A2-x

（2））]'

%梯度

第7页

二、总结

不知不觉，一个学期就这样过去了，随着期末尾声的到来，《最优化》这

门课也已经结课了，接着而来的是为期两个星期的课程设计。

作为信息与计算科学系的一员，我们在课堂上接触最多的就是理论部分，刚好

课程设计可以给我们验证理论，联系实际的机会，所以，这次课程设计，我一直都在认真完成。

在拿到自己的课程设计题目的时候，我对它进行了一定时间的审题，题目中要求我用共轭梯度法求解无约束优化问题，所以我拿到课本，仔细地看了

共轭梯度法的原理与例题求解，以便帮助我更好的理解共轭梯度法的应用。

我记得，再一次研究生面试中，有的老师问过考生：

你在你的大学生涯中学到

了什么，如果让我回答，我会说：

我学到了一种发现问题与解决问题的方法，这种方法是在我的学习中慢慢形成的一种潜移默化的思想，关于学到的内容，我们不会一直保留太清晰的记忆，也不会在以后的生活中用到太多，就像买菜不会用到微积分，不会求导，但是，数学方法和思想给我们思考问题所带来的方法论却是大有裨益。

这次课程设计，我发现我的很多不足，比如在MATLA上的应用上，我还不熟

练，这与我在平时的练习中没有加以重视有关，而且在《最优化》的很多知识点上，自己还没有很好地掌握。

最后，我要感谢我们的任课教师王吉波老师，感谢他在给我们上课中的带来的幽默感和知识的传授。

在此，祝愿老师身体健康，工作顺利。

完稿日期：

2014年6月21日

第8页

最优化方法课程设计沈阳航空航天大学课程设计用纸参考文献

参考文献

[1]陈宝林•《最优化理论与算法》[M].北京,清华大学出版社，2013.[2]刘卫国.《MATLABg序设计教程》[M].北京，中国水利水电出版社,2010.[3]马昌凤.

《最优化方法及其MATLAB?

序设计》[M].北京,科学出版社,2010.

第9页

最优化方法课程设计沈阳航空航天大学课程设计用纸附录

附录

源程序：

n=length（xO）;

g=feval（gfun,xO）;

=O）

m=O;

mk=O;

2O）

（feval（fun,xO+rh°

Am*d）vfeval（fun,xO）+sigma*rh°

第10页

g=[400*x

（1）*（x（1F2-x

（2））+2*（x

（1）-1）,-200*（x（1F2-x

（2））]'

第11页