Matlab笔记数值计算—线代篇016.docxVIP

16. 数值计算—线代篇一、行列式det(A)——矩阵A的行列式； inv(A)——矩阵A的逆；rank(A)——矩阵A的秩；B(: , i)=b——将向量b赋给矩阵B的第i行；[A, eye(5)]——在矩阵A右端，拼接5阶单位矩阵；[U,s]=rref(A)——对矩阵A作行变换，U返回A的最简行阶梯形矩阵，s为行向量存储U的各行首个非0元所在列号，length(s)即为A的秩；例1 用初等行变换法求矩阵的逆。代码：format short g % 省略小数位多余的0A=[1 2 3; 2 2 1; 3 4 3];B=rref([A,eye(3)])% 对矩阵[A,I]进行初等行变换,得到最简行阶梯矩阵Bif(rank(B(:,1:3))==3)% 判断B的前3列是否为单位阵，若是取出后3列，即A逆 A1=B(:,4:6)elsedisp(A不可逆);end运行结果：B =1 0 0 1 3 -20 1 0 -1.5 -3 2.50 0 1 1 1 -1A1 = 1 3 -2 -1.5 -3 2.51 1 -1例2解方程代码：syms x;A=[3 2 1 1;3 2 2-x^2 1;5 1 3 2;7-x^2 1 3 2];D=det(A) f=factor(D) % 对行列式D进行因式分解X=solve(D) % 求