计算物理 第2章 习题及答案.docVIP

第2章：方程数值解法 2.1线性代数方程组的数值解法 ===【2.1.1】=== 用高斯消去法解方程组 【解】 对增广矩阵进行初等变换，化为上三角方程组然后回代求解 所以有： 得到： ===【2.1.2】 用列主元消元法的方程组（采用列选主元：每次消元时主元的选取是各列中系数最大的） 【解】 第1列主元为3，交换第1、2方程位置后消元得: 余下方程列主元为，交换第2、3方程位置后消元得: 回代解得 ===【2.1.3】 将矩阵A进行LU三角分解（采用Doolittle分解） 其中 【解】 （1）则矩阵的Doolittle分解为： （2）解方程组 （3）解方程组 所以 ===【2.1.4】 用LU分解方法给出下列代数方程组的解及L,U矩阵 【解】 ，，，， ===【2.1.5】 用消元法将三对角系数矩阵变成上三角阵，推出变换公式， 并用matlab编程计算矩阵方程 【解】 （1）三对角矩阵方程的形式为 消去系数矩阵主对角元下面的矩阵元，就变成上三角矩阵， 第二行变换公式为：， 同样：第i行变换为：， 则其解： ，， 由此可容易实现编程计算。 （2） % exa_2_1_5.m clc;clear;format long; a =