最速下降法求解线性代数方程组.docxVIP

《MATLAB语言及应用》大作业 姓名： 学号： 学院： 班级： 题目编号： 2013 年 10 月 13 最速下降法求解线性代数方程组。 最速下降法的数学理论 从某个初始点X?出发，沿/(X)在点x(°)处的负梯度方向 严=—巧(灯))-卅 求得/(X)的极小值点X⑴，即 min/(X(0) + 2r(0)) 20 7 然后从X⑴出发，重复上而的过程得到X⑵。如此下去，得到序列｛X? /(x(0))/(x(,)).../(x(A)) 可以证明，从任一初始点X?出发，用最速下降法所得到的序列｛X)｝均收敛 于问题使X最小化/(X)的解，也就是方程组AX=h的解。其收敛速度取决于 牛二芈，其屮人，入分别为A的最小，最大特征值。最速下降法迭代格式：给定 血+人 初值x(0), X⑹按如下方法决定： rM=-\7f(X{k)) = b-AXa) r(k)T,Ar(k) r(k)T,Ar(k) 最速下降法的算法和流程图 X(k+i)= x)+W) 算法: ⑴任取 ^a(0)gRah ,计算 rA(0) = h-ArA(O) (2)对 0, 1,2,…,计算 a(k)=(rA(k),rA(k))/(A*rA(k),rA(k)) xA(k+1 )=xA(k)+a(k)*rA(k) rA(k+1 )=b-A*xA(k+ l)=rA(k)-a(k)*A*rA(k) 若(||(k+l)|