河北工业大学数值分析实验五实验六实验报告 (6000字).docVIP

河北工业大学数值分析实验五实验六实验报告 (6000字) 数值分析实验报告 姓 班 学 实验五 矩阵的LU分解法，雅可比迭代 1.实验名称 矩阵的LU分解法，雅可比迭代 2.实验目的 1）熟悉求解线性方程组的有关理论和方法； 2）会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序； 3）通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。 3.实验内容 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序，进一步了解各种方法的优缺点。 4.算法描述 1）列主元高斯消去法 算法：将方程用增广矩阵[A∣b]=(aij)n?(n?1)表示 1) 消元过程 对k=1,2,…,n-1 ①选主元，找ik??k,k?1,?,n?使得 aik,k = maxk?i?n aik ②如果a ik,k ?0 ，则矩阵A奇异，程序结束；否则执行③。 ③如果ik akj?aikj ?k ，则交换第k行与第ik行对应元素位置， j=k,┅,n+1 ④消元，对i=k+1, ┅,n计算 lik?aik/akk 对j=l+1, ┅,n+1计算 aij?aij?likakj 2) 回代过程 ①若ann②x n ?0 ，则矩阵A奇异，程序结束；否则执行②。 ；对i=n-1, ┅,2,1,计算 ?an,n?1/ann n ? xi??ai,n?1? ???/a?aijxj??ii j?i?1? 2）矩阵直接三角分解法 算法：将方程组Ax=b 中的A分解为A=LU，其中L为单位下三 角矩阵，U为上三角矩阵，则方程组Ax=b化为解2个方程组Ly=b，Ux=y，具体算法如下： ①对j=1,2,3,…,n计算 u1j?a1j 对i=2,3,…,n计算 li1?ai1/a11 ②对k=1,2,3,…,n: a. 对j=k,k+1,…,n计算 k?1 ukj?akj? ?l q?1 kq uqj b. 对i=k+1,k+2,…,n计算 k?1 lik?(aik? ?l q?1 iq uqk)/ukk ③y1?b1，对 k?1 k=2,3,…,n计算 yk?bk? ?l q?1 kq yq ④xn ?yn/unn ,对k=n-1,n-2,…,2,1计算 n xk?(yk? ?u q?k?1 kq xq)/ukk 注：由于计算u的公式于计算y的公式形式上一样，故可直接对增广矩阵 ?a11?a?21???an1 a1,n?1? ? a22?a2na2,n?1 ??????? ? an2?annan,n?1? ?a12 a1n [A∣b]= 施行算法②，③，此时U的第n+1列元素即为y。 3）雅可比迭代法 算法：设方程组Ax=b系数矩阵的对角线元素aii?0(i?1,2,?,n),M为 迭代次数容许的最大值,ε为容许误差。 T①取初始向量x=(x1(0),x(20),?,x(0))，令k=0。 n ②对i=1,2,…,n 计算 x (k?1) i ? 1aii n n (bi? ?a j?1j?i ij xj) (k) ③如果? i?1 xi (k?1) ?xi (k) ?? ，则输出x(k?1)，结束；否则执行④。 ?k?1,转②。 ④如果k≥M,则不收敛，终止程序；否则k 4）高斯-塞尔德迭代法 算法:设方程组Ax=b的系数矩阵的对角线元素,aii?0(i?1,2,?,n),M 为迭代次数容许的最大值,ε为容许误差 ①取初始向量x?(x (0)1 ,x2,?,xn) (0)(0)T