[工学]第三章MATLAB线性方程组及矩阵特征值.pptVIP

然后按公式 依顺序 进行回代，求得方程组的解，该过程称为 “赶”。 具体算法如下： 三对角方程组的求解程序tri_diag.m如下： % tri_diag(a,b,c,d,n) solves a tridiagonal equation. function f = tri_diag(a,b,c,d,n) for i=2:n r=a(i)/b(i-1); b(i)=b(i)-r*c(i-1); d(i)=d(i)-r*d(i-1); end d(n)=d(n)/b(n); for i=n-1:-1:1 d(i)=(d(i)-c(i)*d(i+1))/b(i); end f=d; 一.简单迭代法 1.迭代法建立. 考虑 (矩阵B不唯一) 对应写出 §3.2 解线性方程组的迭代法 产生向量序列 若收敛,记 则于(3.4)两端取极限有： 上式说明： 是解向量 ,从而当k充分大时 注意: 迭代阵B不唯一,影响收敛性。 解向量 (3.4)叫简单迭代法,B叫迭代矩阵。 2.收敛性. 定义3.3 称 为矩阵B的谱半径。 定理3.4 定理3.3 简单迭代法 设有方程组( 其中　　　 ) Ax = b,即 (3.5)