- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
运
筹
学
实
习
报
告
姓 名: xxxxxxxxxx
学 号: xxxxxxxxxxx
专业班级: xxxxxxxxxxxx
2 0 1 3年 7 月 0 4 日
PAGE \* MERGEFORMAT
PAGE \* MERGEFORMAT 1
题目:用最速下降法求解无约束非线性规划问题
摘要:
无约束最优化问题的求解方法分为解析法和直接法两大类。解析法需要计算函数的梯度,其中最速下降法就属于解析法中的一种。对于一个无约束非线性规划利用最速下降法求解,首先需要确定其优化方向,此优化方向应该选择为f在当前点处的负梯度方向,利用一维搜索法找出沿此方向上的最小值及其对应点,此后将该点作为新的出发点重复上述过程,直到达到允许的误差为止。本文通过理论的计算方法,进一步分析,最后用c++编程实现求出允许误差内的最优解。此编程可用于计算符合下列形式的函数求最优解过程:
f(x)=a[0]x1*x1+a[1]x2*x2+a[2]x1*x2+a[3]x1+a[4]x2+a[5]其中:a[i] (i=0,1,2,3,4,5) 为函数的系数。
本文以“ 李占利 主编,中国矿业大学出版社出版”的《最优化理论与方法》 第五章 “无约束最优化方法,5.1 最速下降法 ”例5—1为实例,首先利用上述迭代的方法,计算出各迭代点的函数值,梯度及其模。然后应用c++语言编程,得到在精度范围内的精确最优解。
C++编程计算的最优解为 : 。
即转化为分数结果为: 。
满足精度要求的模为:。
关键词:无约束非线性规划 解析法 最速下降法 梯度 模 最优解
算法思想
无约束最优化方法中的最速下降法首先需要确定其优化方向,此优化方向应该选择为f在当前点处的负梯度方向,利用一维搜索法找出沿此方向上的最小值及其对应点,此后将该点作为新的出发点重复上述过程,直到达到允许的误差为止。主要依据解无约束非线性规划问题的最速下降法计算步骤进行设计算法。
具体步骤如下:
第1步 选取初始点,给定终止误差 0,令k=0;
第2步 计算,若,停止迭代,输出,否则进行第3步;
第3步 取;
第4步 进行一维搜索,求,使得,令,k=k+1。转第2步。
由以上计算步骤可知,最速下降法迭代终止时,求得的是目标函数驻点的一个近似点。 依据以上步骤就可以用C++编程实现最速下降法求解最优解的算法。
算法流程图
停
停
取,
k:=0
计算
是
否
求
令
k:=k+1
程序代码
#includeiostream.h
#includemath.h
double lamda(double x[2],double p[2],double a[2])
{
double lam1,lam2;
lam1=(pow(a[0],3)*x[0]*x[0]+pow(a[1],3)*x[1]*x[1]);
lam2=-(pow(a[0]*x[0],2)+pow(a[1]*x[1],2));
double s;
s=-lam2/(2*lam1);
return s;
}
void main()
{
cout@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@endl;
cout@@ 最速下降法求解最优解程序运行结果 @@endl;
cout@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@endl;
coutendl;
double lamd,x[3],a[6];
double p[2],g[2],e,y,m,n;
int i=0;
cout请输入精度eendl;
cine;
cout请输入初始点x[0],x[1]的值:\nendl;
cinm;
cinn;
x[0]=m;
x[1]=n;
cout函数通式为f(x)=a[0]x1*x1+a[1]x2*x2+a[2]x1*x2+a[3]x1+a[4]x2+a[5]endl;
cout请依次输入函数的系数:a[0]、a[1]、a[2]、a[3]、a[4]、a[5]:endl;
for(i=0;i6;i++)
cina[i];
p[0]=(2*a[0]*x[0]+a[2]*x[1]+a[3]);
p[1]=(2*a[1]*x[1]+a[2]*x[0]+a[4]);
g[0]=-p[0];
g[1]=-p[1];
i=0;
coutendl;
while(sqrt(g[0]*g[0]+g
文档评论(0)