第5章线性代数相关运算.pptVIP

第5章 线性代数相关运算 * 5.1 矩阵 5.1.1 矩阵的修改 A(i,j) 表示矩阵A的第i行第j列元素. A(:,j) 表示矩阵A的第j列元素. A(i,:) 表示矩阵A的第i行元素. A(:) 表示矩阵A的所有元素以列优先排列成一个列矩阵. A(i) 表示矩阵A的第i个元素(以列优先排列). A(i:j) 表示矩阵A的第i个元素与第j个元素之间的所有元素. [] 表示空矩阵. 例5.1 已知矩阵 要求: 1)提取矩阵中第4个元素及第2行第3列的元素. A(4), A(2,3) 2)将原矩阵中第3行元素替换为(-1 -3 -5). A(3,:)=[-1 -3 -5],A 3)在以上基础上,再添加一行元素(-2 -4 -6). A(4,:)=[-2 -4 -6],A 4)在以上基础上,再删除第一列. A(:,1)=[] 例5.2 已知矩阵 解:MATLAB命令为: A=[1 -2;0 4]; B=[-1 0;0 2]; C=[A,B] D=[A;B] AB=[A,zeros(2);zeros(2),B] 5.1.2 矩阵的基本运算 MATLAB运算要满足如下矩阵运算法则： 1）矩阵加减法：C=A+B，C=A-B（要求A，B为同型矩阵. 矩阵的乘法: C=AB(要求A的列数要与B的行数相同). 若A为方阵,则A的乘幂为C=An. 3) 矩阵求逆:要求A为方阵,且A的行列式不等于零. 表5-1 矩阵的基本运算 A^n 方阵A的n次方 矩阵的乘幂 A*B 要求A的列数要与B的行数相同 矩阵与矩阵相乘 K*A 数与矩阵中的每个数相乘 数乘 A±B 同型矩阵相加减 矩阵的加减法 命令形式 功能 运算 rank(A) 求矩阵的秩 矩阵的秩 det(A) 求矩阵的行列式值 矩阵的行列式 B/A(B*inv(A)) XA=B， X=B/A，右除 矩阵的右除 rref(A) 求矩阵的行最简形式 矩阵行变换化简 A’ 求矩阵的转置AT 矩阵的转置 B\A(inv(A)*B) AX=B, X=B\A，左除 矩阵的左除 Inv(A)或A^(-1) 求方阵的逆 矩阵求逆 命令形式 功能 运算 例5.3 已知矩阵 求:A+B,A-B,5A,AB 解:MATLAB命令为: A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;2 2 0;3 3 3]; A+B A-B 5*A A*B 例5.4 已知矩阵 求：|B|，B-1. 解:MATLAB命令为: B=[3 0 -1 2;2 1 0 1;6 3 2 5;0 -1 -1 2]; det(B) inv(B)，B^(-1) 例5.5 已知矩阵 求：CT，CB，B3. 解:MATLAB命令为: B=[1 0 02;2 2 0 ;3 3 3 ]; C=[1 3 5;2 4 6]; C’ C*B B^3 例5.6 求矩阵 的秩和行最简型. 解:MATLAB命令为: A=[4 1 2 4;1 2 0 0;8 5 2 1;0 1 1 7 ] rank(A) rref(A) 5.1.3 矩阵的特征值和特征向量 p=poly(A):求A的特征多项式,给出的结果是多项式所对应的系数(降幂排列) 特征多项式形式：poly2str(p,’x’). 2)d=eig(A):返回A的特征值组成的列向量. 3)[V,D]=eig(A),D的主对角线为特征值,V为特征向量. 例5-8 已知矩阵 求A的特征多项式和全部特征值. 解:MATLAB命令为: A=[1 -1; 2 4 ]; p=poly(A) d=eig(A) 所以A的多项式为x2-5x+6;全部特征值为2和3. 例5-9 已知矩阵 求B的特征值及对应的特征向量. 解:MATLAB命令为: B=[1 2; 0 3 ]; [V,D]=eig(B) V= 1.0000 0.7071 0 0.7071 D= 1 0 0 3 所以B特征值为1和3及对应V的列特征向量. 例5-10略 5.3 线性方程