南理工高等数值分析课程设计张军.docxVIP

MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 1 \h \* MERGEFORMAT 南京理工大学 课程考核论文 课程名称： 高等数值分析 论文题目： 求解线性方程组Ax=b的极 小化方法比较 姓 名： 学 号： 成 绩： 任课教师评语： 任课教师评语： 签名： 年 月 日 摘要 本文对求解线性方程组Ax=b的极小化方法进行了理论上的阐述，并且选取三种具有代表性的方法：最速下降法、共轭梯度法（CG）、预处理共轭梯度法（PCG），使用Matlab编程并分别求解相同的线性方程组，在准确性和收敛速度方面进行比较。结果表明，如果预处理矩阵选择得当，使用预处理共轭梯度法（PCG）效果最好。 极小化方法 极小化方法的基本原理是变分原理。 设A对称正定，求解的方程组为 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 其中，，。考虑2次函数，定义为 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) 有如下性质 = 1 \* GB2 ⑴ 对一切， MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) = 2 \* GB2 ⑵ 对一切， MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4) = 3 \* GB2 ⑶ 设为(1.1)的解，则 而且对一切，有 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 5) 则有定理：设A对称正定，则为（1.1）解的充分必要条件是满足 求，使取最小，这就是求解等价于方程组（1.1）的变分问题。求解的方法一般是构造一个向量序列，使。 最速下降法 可通过求泛函的极小点来获得式（1.1）的解。为此，可以从任一出发，沿着泛函在处下降最快的方向搜索下一个近似点，使得在该方向上达到极小值。 由多元微积分知识，在处下降最快的方向是在该点的负梯度方向，经过简单计算，得： , 此处，也称为近似点对应的残量。取， 确定的值使取得极小值。由式（1.4），令 , 得 从上式求出并记为： MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 1. SEQ MTEqn \c \* Arabic \* MERGEFORMAT 6) 综合上述推导过程，可得最速下降法的算法描述： = 1 \* GB2 ⑴ 给定初始点，容许误差，置。 = 2 \* GB2 ⑵ 计算，若，停算，输出作为近似解。 = 3 \* GB2 ⑶ 按式（1.6）计算步长因子，置，，转步骤 = 2 \* GB2