Matlab矩阵分析与处理1.pptVIP

第三章 Matlab矩阵分析与处理 * 矩阵是Matlab的基本特征，也是Matlab重要特性，它的运算功能丰富而方便，前一章介绍的矩阵的建立及基本运算，本章介绍矩阵分析与处理。 3.1 特殊矩阵 Matlab中产生特殊矩阵的函数有： zeros: 产生全0矩阵（零矩阵） ones：产生全1矩阵 (幺矩阵） eye： 产生单位矩阵 rand：%产生0-1间均匀分布的随机矩阵； randn %产生均值为0，方差为1的标准正态分布随机矩阵 调用格式： zeros(m): % 产生m×m零矩阵； zreos(m,n): %产生m×n零矩阵； zreos(size(A)): %产生于矩阵A同样大小的零矩阵 其它的几个函数的调用格式相似 例：分别建立3×3，3×2，2×3零矩阵 zeros(3) zeros(3,2) zeros(2,3) 例：建立随机矩阵 （1）在区间[20，50]内均匀分布的5阶矩阵 （2）均值为0.6、方差为0.1的5阶正态分布随机矩阵 解（1） z=20+(50-20)*rand(5) 输出结果： z = 48.5039 42.8629 38.4630 32.1712 21.7367 26.9342 33.6940 43.7581 48.0641 30.5860 38.2053 20.5551 47.6544 47.5071 44.3950 34.5795 44.6422 42.1462 32.3081 20.2958 46.7390 33.3411 25.2880 46.8095 24.1667 解（2） y=0.6+sqrt(0.1)*randn(5) 输出结果（略） 3.2 矩阵结构变换 一、对角阵 对角线上有非0元素的矩阵为对角阵，对角线上元素相等的对角阵称为数量矩阵，对角线上的元素全为1称为单位矩阵。 （1）提取矩阵对角线元素 A=[1 2 3;4 5 6]; D=diag(A) diag(A)函数还有一种形式diag(A,k),可提取第k条对角线元素，主对角向上为1，向下为-1，类推。 （2）构造对角矩阵 diag([1 2 -1 4]) ans = 1 0 0 0 0 2 0 0 0 0 -1 0 0 0 0 4 输出结果： diag(1:3,1) diag(1:3,-1) 例：建立一个5×5矩阵，然后将第1行乘1第2行乘2，…… A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15 16 17 18 19 20;21 22 23 24 25] d=diag(1:5) B=d*A 二、三角阵 三角阵又分上三角阵、下三角阵 A=[7 13 -28;2 -9 8;0 34 5] B=triu(A) triu(A,k)表示第k条对角线以上保留，其它置0 A=[7 13 -28;2 -9 8;0 34 5] tril(A) %取下三角 tril(A,k)与triu(A,k)一样使用 3.2.2 矩阵的转置和旋转 A=[7 13 -28;2 -9 8;0 34 5] B=A 旋转 A=[7 13 -28;2 -9 8;0 34 5] B=rot90(A) 转置 格式：fliplr(A) %左右翻转 flipud(A) %上下翻转 3.3 矩阵求逆于线性方程组求解 称B为A矩阵的逆矩阵 A=[7 13 -28;2 -9 8;0 34 5] B=inv(A) C=A*B 故A与A-1是互逆的，满秩矩阵才可能互逆。 二、用矩阵求逆的方法求解线性方程组 线性方程组： 其矩阵表达式： 线性方程Ax=b的解 例：用求逆矩阵的方法解线性方程组 Matlab程序： A=[1 2 3;1 4 9;1 8 27]; b=[5;-2;6]; x=inv(A)*b 输出结果： x = 23.0000 -14.5000 3.6667 x=A\b 也可用矩阵左除 det(A) 求A(方阵）行列式的值 A=rand(5) B=det(A) 三、矩阵的秩和迹 求秩函数格式：rank(A)：求矩阵A的秩 1、矩阵的秩—矩阵线形无关的行或列数称为矩阵的秩。 何谓矩阵线形无关的行或列？ 求迹函数格式： trace(A)：求矩阵A的迹 2、矩阵的迹—矩阵的