计算物理学chp5摘要.ppt

有唯一解（非奇异方程组） 无解 有无穷多解 计算物理 (Computational Physics) 几个要注意的问题 目录 线性代数方程组的解法 1.奇异方程组 对任意一组线性方程式，它的解一般有三种可能 奇异方程组 矩阵除法: X=A\B 矩阵方程 AX=B 的解 X=A/B 矩阵方程 XA=B 的解 A为 m?n阶矩阵, m=n, 求精确解 mn, 求最小二乘解 mn, 求基本解 A=randn(3); B=randn(3,1); A\B ans = -0.3583 0.6884 -0.1500 A=randn(3,2); A\B ans = -0.4697 0.7990 A=randn(2,3); B=randn(2,1); A\B ans = 0 -1.2119 -0.4689 计算物理 (Computational Physics) 奇异方程组 目录 线性代数方程组的解法 A=[1 2;1 -3]; B=[13;3]; A\B ans = 9 2 A=[1 2;2 4]; B=[4;3]; A\B Warning: Matrix is singular to working precision. (Type warning off MATLAB:singularMatrix to suppress this warning.) ans = Inf Inf B=[4;8]; A\B Warning: Matrix is singular to working precision. (Type warning off MATLAB:singularMatrix to suppress this warning.) ans = Inf Inf seidel(A,B) x = -6 5 k = 2 计算物理 (Computational Physics) 病态方程组 目录 线性代数方程组的解法 2.病态方程组 方程组系数的微小变化导致方程组解的巨大变化的方程组 A1=[1 10;10 101]; B=[11;111]; x1=gauss2(A1,B) x1 = 1 1 A2=[1 10;10.1 100]; x2=gauss2(A2,B) x2 = 10.0000 0.1000 可用矩阵A的条件数度量： cond() 矩阵条件数函数，当矩阵条件数接近1时，说明矩阵非病态 cond(A1) ans = 1.0402e+004 cond(A2) ans = 1.0203e+004 A=[1 2;1 -3]; cond(A) ans = 2.6180 计算物理 (Computational Physics) 迭代解法的收敛性 目录 线性代数方程组的解法 3.迭代解法的收敛性 迭代解法的前提条件是迭代解出的近似解序列必须具有收敛性。如果近似解序列是发散的，则迭代解法不能获得其解。 求线性方程组的迭代法公式： 对任意初值x0和任意d都收敛的充分必要条件： λi 表示M的特征值，ρ(M) 称为M的谱半径 函数eig() 求矩阵特征值 函数det() 求矩阵行列式 函数inv() 求矩阵逆 计算物理 (Computational Physics) 迭代解法的收敛性 目录 线性代数方程组的解法 A=[10 -1 -2;1 -10 2;1 1 -5]; v=diag(A); D=diag(v); L=tril(A,-1) L = 0 0 0 1 0 0 1 1 0 U=triu(A,1) U = 0 -1 -2 0 0 2 0 0 0 A=[2 3 1;1 1 1;1 -2 -1]; v=diag(A); D=diag(v); L=tril(A,-1); U=triu(A,1); M1=-inv(D)*(L+U); eig(M1) ans = 1.8100 -0.1683 -1.6418 M1=-inv(D)*(L+U); eig(M1) ans = -0.2372 -0.1000 0.3372 M2=-inv(D+L)*U; eig(M2) ans = 0 -0.0318 0.1258 计算物理 (