线性共轭梯度法解决二次正定优化问题.docVIP

最优化理论与算法实验报告 四 实验名称 线性共轭梯度法解决二次正定优化问题 实验时间 姓名 专业班级 学号 成绩 实验目的和内容 实验目的：通过实验, 让学生掌握解决二次优化问题的具体实现, 同时对于具体的问题设计, 让学生根据在实验中出现的数值计算结果, 分析并发现存在的不足. 理解的基本思想. 实验内容：请采用求解 终止的准则, 相关背景知识介绍 其中为正定矩阵，则. 对二次函数总有: 1）设为初始点. 首先取，令 （为精确步长因子） 则有：（精确一维搜索性质）. 2）令，适当选择，使， 得 （从而得到） 且有正交性：，, ， 3）再令，适当选择，，使得 （），由此得： ， 4） 一般地，在第次迭代中，令 适当选取，使 （）， 可得到 （） 利用归纳法可以证明: （），再由与的关系得： （） 所以当时，，而 . 共轭梯度法的迭代公式为： （为共轭方向，为最佳步长因子） 对二次函数: ；算法如下描述: Step 1: 取初始点, , 则停机; 否则, Step 2: 计算, 如果, 则停机; 否则 go to S；Step 3: , , go to Step 2 代码 function f getd x0,t x x0; syms xi yi a f xi^2+4*yi^2-4*xi-8*yi; fx diff f,xi ; fy diff f,yi ; fx subs fx, xi,yi ,x0 ; fy subs fy, xi,yi ,x0 ; fi [fx,fy]; c 0; while double sqrt fx^2+fy^2 t s -fi; if c 0 s -fi; else s s1; end x x+a*s; f subs f, xi,yi ,x ; f1 diff f ; f1 solve f1 ; if f1~ 0 ai double f1 ; else break x,f subs f, xi,yi ,x ,c end x subs x,a,ai ; f xi-xi^2+2*xi*yi+yi^2; fxi diff f,xi ; fyi diff f,yi ; fxi subs fxi, xi,yi ,x ; fyi subs fyi, xi,yi ,x ; fii [fxi,fyi]; d fxi^2+fyi^2 / fx^2+fy^2 ; s1 -fii+d*s; c c+1; fx fxi; fy fyi; end x,f subs f, xi,yi ,x ,c 数值结果f getd [0,0],0.000001 运算次数 最优解 X0.2500, x2 -0.2500 函数最小值 0.1250 计算结果的分析 计算中出现的问题，解决方法及体会 教 师 评 语 指导教师： 年 月 日 实验原始数据记录