线性代数中的MATLAB的命令.pptVIP

第 4 章 线性代数中的MATLAB命令 ; ;矩阵运算;矩阵的基本运算 rank(A) 矩阵A的秩 det(A) 矩阵A的行列式 inv(A) 矩阵A的逆矩阵 % 或 A^(-1) norm(A) 矩阵A的范数 trace(A) 矩阵A的迹(对角线元素的和);特殊矩阵的生成 zeros(m, n) m行n列的零矩阵; ones(m, n) m行n列元素全为1的矩阵; eye(n) n阶单位矩阵; rand(m, n) m行n列[0,1]上均匀分布随机数矩阵 diag(A) 矩阵A的对角线元素构成的列向量 diag(x) 向量x的元素构成的对角矩阵 tril(A) 矩阵A的下三角部分，其余置0 triu(A) 矩阵A的上三角部分，其余置0;特殊矩阵的生成 flipud(A) 矩阵A的上下翻转 fliplr(A) 矩阵A的左右翻转 reshape(A, m, n) 矩阵A的元素重排成m行n列矩阵;矩阵的特征值、特征向量 eig(A) 求方阵A的特征值组成的列向量 [v, d]=eig(A) 其中，v是矩阵A的特征向量(列向量)构成的矩阵，d是矩阵A的特征值构成的对角阵。（每个特征向量与其特征值列号一致）;矩阵的初等行变换 rref(A) 对矩阵A进行初等行变换化简;参见教材第100-106页（自学）;求解线性方程组;求解线性方程组; 在MATLAB中，矩阵左除法可用于求线性方程组Ax=b的一个解 设未知数的个数为n ， (1)若rank(A)=rank([A,b])=n，则Ax=b存在唯一解， A\b 求得这个唯一解（用列向量表示，下同）； (2)若rank(A)=rank([A,b])n，则Ax=b有无穷多解， A\b 求得包含最多零元素的一个特解； (3)若rank(A)≠rank([A,b])，则Ax=b无解， A\b 求得一个最小二乘解;注意 求解线性方程组时，首先应该判断解的情况； 对于第(2)种情况，求得一个特解后，还应该求出齐次线性方程组Ax=0的基础解系，以便得到原方程组的通解。 在MATLAB中，求齐次线性方程组Ax=0的基础解系命令为 null(A) (用列向量表示);例 求解下列线性方程组 ;解: (1) A=[1 2 7; 3 -2 5; 6 9 -1], b=[1; -4; 3] rA=rank(A), rAb=rank([A, b]) x0=A\b (2) A=[5 2; -3 7; 1 -1], b=[2; 1; -3] rA=rank(A), rAb=rank([A, b]) x0=A\b (3) A=[1 -1 1 -1; -1 1 1 -1; 2 -2 -1 1], b=[1; 1; -1] rA=rank(A), rAb=rank([A, b]) x0=A\b x=null(A);例 对于任意一个线性方程组Ax=b，试编程判断解的情况，并求得相应的解。 演示（xianxingfangchengzu.m）