MATLAB实验第四讲.pptVIP

* 第四讲 行列式、矩阵、线性方程组、特征值和特征多项式的求解 一、 行列式 行列式的计算命令 D=[ , ; , ; , ; , ;]; Det(D) 当行列式是符号表示时，一定要注明 Simple() 化简结果 zeros(m,n) 生成0矩阵 ones(m,n) 生成单位矩阵 rand（m，n） 0到1之间的均匀分布的随机矩阵 randn(m,n) 标准正态分布的随机元素矩阵 magic(n) 每行每列和对角线上元素之和相等矩阵 diag（M） 从矩阵M对角线元素生成一列数组或将一 列数组M生成一个对角矩阵 triu(M) 取矩阵M的对应元素生成上三角矩阵 tril(M) 取矩阵M的对应元素生成下三角矩阵 length(M) 返回向量M的长度 size(M) 返回矩阵M的行数和列数 eye(n) 生成n阶单位阵 hilb(n) 生成n阶希尔伯特病态矩阵 pascal(n) 帕斯卡矩阵 二、 矩阵的生成 三、 矩阵的取块和变换 A(i,:) 提取A的第i行 A(:,j) 提取A的第j列 A(:) 将矩阵A的各列元素依次排成一列向量 A(i:j) 将矩阵A(:) 中的第i个到第j个的元素一次排成一行 向量 A(i:j,:) 提取矩阵A的第i行到第j行的所有元素所成的矩阵 A(:,i:j) 提取矩阵A的第i列到第j列的所有元素所成的矩阵 A(i:j,k:l) 提取矩阵A的第i行到第j行以及第k列到第l列的 所有元素所成的矩阵 B=reshape(A,m,n) 将矩阵A中的元素按次序组成一个m 行n列的矩阵 B= rot90(A) 将矩阵逆时针旋转90度得到新矩阵 B= fliplr(A) 将矩阵A中的元素左右对称得到新矩阵 B= flipud(A) 将矩阵A中的元素上下对称得到新矩阵 A(2,:)=[ ] 将矩阵A的第2行删除 A(A5)=1 将矩阵中所有大于6的元素赋值为1 三、 矩阵的基本运算 A+B A-B A*B A/B 加， 减 ，乘， 除 A+B A-B A.*B A./B 对应元素的运算 A.’ 矩阵A的转置 Inv(A) 或(A)^(-1) 矩阵A的逆矩阵 A^k 矩阵的k次幂 Squrtm(A)或A^(1/2) 矩阵A的开方 四、 线性方程组的求解 常用的命令; rank(A) 矩阵的秩 rref(A) 矩阵A的行最简形式 null(A) 系数矩阵为A的齐次方程组 基础解系 null(A,’r’) 系数矩阵为A的齐次方程组 有理数形式的基础解系 五、 特征值和特征多项式 常用的命令; trace(A) 得到矩阵A的迹 poly(A) 矩阵A的特征多项式形式 [a,b]=eig(A) 得到矩阵A的特征列向量矩阵a和对应 特征值组成的对角阵b B=orth(A) 正交化空间，即矩阵B的列向量正交 且生成的线性空间与矩阵A的列向量 生成的线性空间等价 A=[-1,-2,1,-2,1;-1,1,1,-1,1;2,-2,-1,1,-1] A=[0.1,0.1^2;0.2,0.2^2;0.3,0.3^2;0.4,0.4^2] B==[0.1,0.1^2,0.045; 0.2,0.2^2,0.12; 0.3,0.3^2,0.20; 0.4,0.4^2,0.33] A=[1,2,3;2,-1,-1;1,-2,-2;1,-1,-1] b=[2;1;-1;0] B=[A,b]; *