线性方程组的最速下降法与共轭梯度法.doc_第1页
线性方程组的最速下降法与共轭梯度法.doc_第2页
线性方程组的最速下降法与共轭梯度法.doc_第3页
线性方程组的最速下降法与共轭梯度法.doc_第4页
免费预览已结束，剩余1页可下载查看

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

共轭梯度法一 共轭梯度法原理对于线性方程组，即： （1）其中，为对称正定矩阵，如何熟练地运用最速下降法与共轭梯度法的求解线性方程组。在求解线性方程组之前，首先用内积将问题转化为函数问题。1 最速下降法最速下降法是一种运用梯度与极值的性质，综合数值计算方法寻找局部极值。基本思想：任一点的负梯度方向是函数值在该点下降最快的方向。将维问题转化为一系列沿负梯度方向用一维搜索方法寻优的问题，利用负梯度作为搜索方向，故称最速下降法。具体步骤：1、搜索方向：，即最速下降方向。2、搜索步长：取最优步长，即满足：1 给定初始点，允许误差，令。2 计算搜索方向。3 若，则为所求的极值点，否则，求解最优步长，使得。4 令，最速下降方向是反映了目标函数的局部性质，它只是局部目标函数值下降最快的方向。2 共轭梯度法对于其中，是对称正定矩阵。基本思想：将共轭性与最速下降法相结合利用已知迭代点的梯度方向构造一组共轭方向，并沿此方向搜索，求出函数的极小值。具体步骤：1 取初始点，取第一次搜索方向为。2 设已求得，若，令，则下一个搜索方向 （1）由于与关于共轭，所以给（1）两边同时乘以，即：解得： （2）3 搜索步长的确定，已知迭代点，和搜索方向，确定步长，即：记，令既有：令 既有：解得：共轭梯度法是对最速下降法的一种改进，减少了迭代次数从而提高了程序运行效率。二 程序框图程序：%最速下降法%function x,k=fast(A,b)esp=input(esp=);N=input(N=);x0=input(x0=);k=0;tol=1;while tol=esp r=b-A*x0; q=dot(r,r)/dot(A*r,r); x=x0+q*r; k=k+1; tol=norm(x-x0); x0=x if k=N disp(); return; endendxk%共轭梯度法%function k,x=C_G(A,b)esp=input(请输入最大误差=);x0=input(请输入初值x0=);k = 0 ;r0 = b-A*x0; %dangqianwhile norm(r0)esp r0 = b -A*x0; k = k + 1 ; if k=1 p0 = r0 ; else lamda=(r0*r0)/(p0*A*p0); r1 = r0 - lamda*A*p0 ; p0=r0+(r0*r0)/(r1*r1)*p
人人文库> 全部分类> 生活休闲 > 科普知识

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论