基于MATLAB的数值分析(3)课件展示.ppt

第三章 线性代数 3.2 矩阵的运算 二、矩阵的左除和右除 左除“ \ ”： 求矩阵方程AX=B的解；（ A 、B的行要保持一致） 解为 X=A\B； 当A为方阵且可逆时有X=A\B=inv(A)*B； 右除“ / ”： 求矩阵方程XA=B的解 （A 、B的列要保持一致） 解为 X=B/A ， 当A为方阵且可逆时有X=B/A=B*inv(A) 3.4 稀疏矩阵 3.5 用预条件共轭斜量法求解线性方程组 3.6 求解大型稀疏非对称正定的线性方程组 3.7 不可解问题 3.8 病态问题 1.不用预优矩阵的双共轭梯度法 x=bicg(a,b,tol,kmax) 2.用预优矩阵的双共轭梯度法 [l,u]=lu(m) x= bicg(a,b,tol,kmax,l,u,x0) 3.未给定预优矩阵的双共轭梯度法 [l,u]=luinc(sa,’0’) x= bicg(a,b,tol,kmax,l,u,x0) 一.双共轭梯度法 二.广义极小残差法 gmres 例：考察下面三个线性方程组的解 (A1) -x+y=1 (A2) -x+y=1 (A3) x+2y=-2 -2x+2y=2 -x+y=0 -x+y=1 2x-y=0 A1=[-1 1;-2 2], b1=[1;2],x1=A1\b1 A10=[-1 1], b10=[1],x10=A10\b10 A2=[-1 1;-1 1], b2=[1;0],x2=A2\b2 A3=[1 2;-1 1;2 -1], b3=[-2;1;0],x3=A3\b3 * * 3.1 常用矩阵函数 norm ： 矩阵或向量范数 For vectors... NORM(V,P) = sum(abs(V).^P)^(1/P). NORM(V) = norm(V,2). NORM(V,inf) = max(abs(V)). NORM(V,-inf) = min(abs(V)). 例 x=[1 2 3 4 5]; x=[1 20 300 40 5]; x=[10000 2000 300 40 5]; [norm(x,1), norm(x,2), norm(x,3), norm(x,inf)] ? 例: 不同范数意义下的单位圆 运行以下Matlab程序，文件名为：normpolt.m 描绘norm(x,1)=1;norm(x,2)=1;norm(x,inf)=1的图形 For matrices... NORM(X) is the largest singular value of X, max(svd(X)). NORM(X,2) is the same as NORM(X). NORM(X,1) is the 1-norm of X, the largest column sum, = max(sum(abs(X))). NORM(X,inf) is the infinity norm of X, the largest row sum, = max(sum(abs(X))). NORM(X,fro) is the Frobenius norm, sqrt(sum(diag(X*X))). NORM(X,P) is available for matrix X only if P is 1, 2, inf or fro. 例 x=[1 20 300 ;400 50 6;70 800 9]; [norm(x,1), norm(x,2), norm(x,inf) , norm(x, fro)] dot(x,y) 向量的内积 det(A) 方阵的行列式； rank(A) 矩阵的秩； trace(A) 矩阵的迹； rref(A) 初等变换化矩阵A为阶梯矩阵 Inv(A) 矩阵的逆；即 A-1 pinv(A) 矩阵的广义逆A+ null(A) 零空间的基阵 roth(A) 值空间的基阵 orth(A)