02_线性方程组的直接法.docVIP

习题课二 求解的直接法 摘要： 知道数值线性代数的四点内容： ①----规范线性方程组问题 ②----最小二乘问题 ③----代数特征值问题 ④----代数奇异值问题 理解求解的线性变换法的核心思想：把A三角化： 对，只要求P为非奇异 …列向量 对，要求P为正交阵 3.掌握每种主要线性变换：Guass变换与Householder变换，？？Givens变换. 4.Guass消元法——主元消元法 直接三角分解（LU分解法，选主元LU分解法） 对称正定阵有Cholesky分解. 三对角…… 这里运用了初等矩阵的性质. 5.向量与矩阵的范数 最大列和 最大行和 如，则为非奇异. 在Ax b问题下的矩阵条件数： 6．误差分析 一、基本方法：给定方程组，用Gauss消元法，主元消元法或直接三角分解法求解（数字题） 分别用顺序Gauss消元法，列主元Gauss消元法解下列方程组（取4位有效数字）比较计算结果与精度 解: 用matlab求得其精确解为 用顺序Gauss消元法 回代解: （下省略，请同学们自做） 2．用列主元消元法求矩阵A的行列式， 解：用Gauss列主元消元法，把A化为一个上三角阵 因为，的行列式均为1，故detA detU U对角线元素之积。 3. 用Doolittle法解方程组 。 解：用“西北角”分解法： ， ， 得： 用Crout分解法解方程组 。 解：设（要求解14+6个元素） 按“西北角”分解法，先后：求L的第一列 求U的第一行 求L的第二列 …… 由矩阵乘法，对比两边的元素，得Crout分解的计算公式为 i 1,2,3,4 j 2,3,4 k 2, 3, 4; i k, … , 4 k 2, 3, 4; j k+1, … , 4 具体的： 即其解为， 即， 然后分别解和， （下略，留给同学们） 5. 用平方根法（cholesky分解）求解下列方程组，A为对称正定阵 解：先做Cholesky分解 对角元： 时： Cholesky分解的优点是不必选主元，且数据稳定 解Ly b. 求出 ， ，， . 得 ，，， 不必选主元的原因是不用消元因子 数值稳定的原因： 因此： 故 ，即的数量级不会增长 6．用改进平方根法解 解：如用cholesky法求解时，要做一次开方运算，改进平方根法避免开方。 作矩阵乘法，并注意到， j k D为对角元，只有j k时元素非零，记 L为下三角阵，行数 列数时其元素为零 列k最多为j 1 ,j 1,2,…,n 对角元 i j i 1,2,…,n 显然. 于是得计算式 按列分解 j 1,2,…,n i j+1,…,n 对本题有j 1 j 2 j 3 7. 用追赶法解三对角方程组 解: 注意到A不是对称，也非对角占优。假如A的前n-1个顺序主子式都非零,那么A存在唯一的三角分解，且为保带的。 方法一：作初等行变换，使A上三角化 消元非常方便 于是有 ，；；。 方法二：作三角分解，实施“追赶法” 作Crout分解 作矩阵乘法： 故有 对i 1,2,3,… ，， U中 “追”步 Ly b 得 “赶”步 Ux y 得 方法三： 作Doolittle分解 略 8． 用列主元消元法解方程组系 Ax b Ay c Az d 其中，，， 解：作增广阵 第一个方程组 第二个 第三个 9．用平方根法求解Ax b，其中 解：设，则平方根法计算式为 j 1,2,3,…,n i j+1,…,n 因为A对称，且本题中， i j . 先用数学归纳法证明，所有 因为 ； ， i 2,3,…,n 作归纳假设 ， i j ，那么我们有 i j+1 归纳法完成。 因此有 由，即 得 由，即 得 10．举例说明之：非奇异矩阵不一定存在LU分解，但非奇异阵可作适当的行对调，可保证存在LU分解. 解: 设.A非奇异,若A存在LU分解,则 比较两边 ； 故不能LU分解 但A的行列对换，，则存在LU分解 11．设，求A的PLU分解。 解：由于A的一二阶主子式都为零，取使 ； 容易知，的一、二、三阶顺序主子式均非零，故有唯一的LU分解. 故 其实P不是唯一的，比如可取使得 12．设且各列线性无关，证明A可分解为，其中R为单位上三角方阵，为对角矩阵 证： A为列满秩 为对称正定矩阵 这样，存在分解（Cholesky分解），其中L为单位下三角矩阵，D为主对角元素为非零的对角阵. 记，即，注意到 故 令，则R为单位上三角阵， 于是 。 二、研究用Gauss消元法解Ax b时的性质 1、证明：单位下三角阵可逆，及两个单位下三角阵之积，仍为单位下三角阵。 证1：设L为单位下三角阵，显然它非奇异，有逆，记。 ，其中为L的伴随矩阵。 ，其中为L的元素的代数余