数值分析思考题4知识分享.docx

数 值 分 析 思 考 题 4 k 1 k 1 2 2 3 n n n n n n n 精品资料 数值分析思考题 4 1、 Gauss 消去法和 LU 三角分解法解线性方程组的工作量相同 吗？工作量为多少？平方根方法的工作量为多少？ 答： Gauss消去法所需的加，减，乘，除运算的次数为： n 1 (n k ) 2(n k )2 n(n 1) n(n 1)(2n 1) 2 n3 O( n2 ) 由于 Gauss消去法与 LU 三角分解法是等价的，因此， LU 三角分解 2 3 法的工作量也是 3 ，两者的工作量相同。 平方根法是 A 为实对称正定矩阵时，三角分解法的变形。其工作量 约为 LU 三角分解法的一半，该算法的计算量为 n3 2、求解一个线性方程的 LU 分解法什么条件下可以保障成功？选主 元的目的是什么？分别用列主元和全主元 Gauss 消去法求解下列方 程组： 12x1 3x2 3x3 15 18x1 3x2 x3 15 x1 x2 x3 6 答：（ 1）条件（三角分解的一个充分条件）：如果 A R 的顺序主子式 Ak Rk k (k 1,2, n 1)均非奇异，则存在唯一的单位 下三角阵 L R 和上三角阵 U R ，满足 A LU 。 并且，在三角分解过程中的除法运算要求分母不能太小， 否则将可能产生不稳定的情况。 仅供学习与交流，如有侵权请联系网站删除 谢谢 2 11 a11 11 a11 l 22 a22 l 1 l T 精品资料 （2） .选主元的目的就是为了完成消元并且避免不稳定情况 的发生。 3、用平方根方法（ 格式存储。 16 4 8 4 5 4 8 4 22 x1 x2 x3 答：系数矩阵 16 A= 4 5 8 4 22 4 1 2 2 3 Cholesky 分解法）求解下列方程组，并用紧凑 4 3 10 4 1 2 = 2 3 22 3 4 L= 1 2 2 3 3 求解方程 4 1 2 L = 2 - 3 b= 3 Ly=b y=( )T T 求解方程组 LTx=y 求得 x=( T 仅供学习与交流，如有侵权请联系网站删除 谢谢 3 1b2 10 1 b 2 10 4 T 精品资料 4、已知线性方程组 2.0002 1.9998 x1 4 1.9998 2.0002 x2 4 （ 1 ）求系数矩阵的逆 A 和条件数 （2 ）若方程组右端有微小扰动 cond ( A)； 2 10 4 , T ，不用求解方 程组，试利用解与系数扰动之间的关系式来估计解的相对变化率。 解：（ 1） A A 得到的结果：系数矩阵的逆 1 10001 8 9999 8 x1 4 x2 4 （2）线性方程组 2.0002 1.9998 1.9998 2.0002 其解为： x 1,1 若方程组右端有扰动，则原方程变为 : 2.0002 1.9998 x x x 1.9998 2.0002 1 2 x1 x 4.0002 3.9998 b b ，其解为： 1 20000 9999 8 10001 8 ,条件数为： 10000。 x