求解线性方程组.docVIP

《线性方程组求解》实验报告 实验名称：? 成绩：___________ 专业班级：? 姓名：学号： 实验日期 ：?2014年月日 实验报告日期：??年月日 一、实验目的 二、实验内容 ，给出系数矩阵A和b分别为： ， 例二:用LU分解法求解下面线性方程组，要求写出分解出的矩阵L和U. 例三：用QR分解法求解线性方程组，要求写出分解出的矩阵L和U. 例四：用Cholesky分解法求解线性方程组，给出A和b分别为： ， 2.2借用实例来区分四种方法的不同 例：用直接左除法、LU分解法、QR分解法求解下面的线性方程组 三、实验 该实验应用matlab2014来进行实验的验证和设计. 四、，给出系数矩阵A和b分别为： ， 给出matlab程序： A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; b=[13;-9;6;0]; x=A\b 得到线性方程组的解： x = -66.5556 25.6667 -18.7778 26.5556 所以，通过左除法得到的线性方程组的解为： （2）LU分解法 例二:用LU分解法求解下面线性方程组，要求写出分解出的矩阵L和U. 给出matlab程序： A=[10,-7,0,1;-3,2,6,2;5,-1,5,-1;2,1,0,2]; b=[8;6;5;1]; [L,U]=lu(A) x=U\(L\b) 得到方程组的解为： L = 1.0000 0 0 0 -0.3000 -0.0400 1.0000 0 0.5000 1.0000 0 0 0.2000 0.9600 -0.7742 1.0000 U = 10.0000 -7.0000 0 1.0000 0 2.5000 5.0000 -1.5000 0 0 6.2000 2.2400 0 0 0 4.9742 x = 0 -1 1 1 所以，通过左除法得到的线性方程组的解为： （3）QR分解法 例三：用QR分解法求解线性方程组，要求写出分解出的矩阵L和U. 求解：由方程组可以得到系数矩阵和初始向量分别为： ， 给出matlab程序： A=[1,-1,1;5,-4,3;2,7,10]; b=[4,7,5]; [Q,R]=qr(A) x=R\(Q\b) 得到结果： Q = -0.1826 -0.0956 -0.9785 -0.9129 -0.3532 0.2048 -0.3651 0.9307 -0.0228 R = -5.4772 1.2780 -6.5727 0 8.0229 8.1517 0 0 -0.5917 x = -4.6154 -4.2308 4.3846 所以，通过左除法得到的线性方程组的解为： （4）Cholesky分解法 例四：用Cholesky分解法求解线性方程组，给出A和b分别为： ， 给出matlab程序： A=[2,1,1;1,2,-1;1,-1,3]; b=[6;3;4]; R=chol(A) x=R\(R\b) 得到结果: R = 1.4142 0.7071 0.7071 0 1.2247 -1.2247 0 0 1.0000 x = 2.0000 1.0000 1.0000 所以，通过左除法得到的线性方程组的解为： 4.2借用实例来区分四种方法的不同 例1：用直接左除法、LU分解法、QR分解法求解下面的线性方程组 已知，该线性方程组的系数矩阵和初始向量分别是： ， 借用matlab求解该线性方程组，使用直接左除法、LU分解法、QR分解法来同时得到相对应的三种方法的结果。给出matlab的程序： format long A=[2,1,5,7;-1,3,-5,6;1,-1,3,9;-4,6,-1,5]; b=[16;7;19;41]; x1=A\b %左除法得到的解% [L,U]=lu(A); x2=U\(L\b) %LU分解法得到的结果% [Q,R]=qr(A); x3=R\(Q