MATLAB矩阵分析与处理2012.pptVIP

Matlab矩阵分析与处理 3.1特殊矩阵 常见的特殊矩阵有零矩阵、幺矩阵、单位矩阵等，这类特 殊矩阵在应用中具有通用性。 1、通用的特殊矩阵 常用的产生通用特殊矩阵的函数有： zeros：产生全0矩阵（零矩阵）。 ones： 产生全1矩阵（幺矩阵）。 eye： 产生单位矩阵。 rand：产生0~1间均匀分布的随机矩阵。 randn：产生均值为0，方差为1的标准正态分布随机矩阵。 产生(0,1)区间均匀分布随机矩阵使用rand函数 产生均值为0，方差为1的标准正态分布随机矩阵使用randn函数 例3.1 建立随机矩阵： ① 在区间[20,50]内均匀分布的5阶随机矩阵。 ② 均值为0.6，方差为0.1的5阶正态分布随机矩阵。 命令如下： x = 20+(50-20)*rand(5) x = 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 y = 0.6 + sqrt(0.1)*randn(5) 3.2 矩阵结构变换 1、对角阵与三角阵 只有对角线上有非零元素的矩阵称为对角矩阵，在研究矩阵 时，有时候需要将矩阵的对角线上的元素提取出来形成一个列 向量，有时也需要用一个向量构造一个对角阵。 （1）提取矩阵的对角线元素函数：diag 例如： A = [1,2,3;4,5,6]; D = diag(A) D = 1 5 diag函数还有一种形式：diag(A,k)提取第k条对角线的元素。 例如： D1 = diag(A,1) D =2 6 3.2 矩阵结构变换 (2)构造对角矩阵 如果V是一个m个元素的向量，diag(V)将产生一个m×m对角矩阵， 其主对角线元素即为向量V的元素。 例如： 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) 例3.2 建立一个5×5矩阵A，然后将A的第一行元素乘以1，第二行乘以2，… 第五行乘以5。 解： 用一个对角矩阵左乘一个矩阵时，相当于用对角阵的第一个元素乘以 该矩阵的第一行，依次类推。 命令如下： A = ones(5); D = diag(1:5); D * A ans = 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 5 5 5 5 5 3.3 矩阵求逆与线性方程组求解 3.3.1 矩阵的逆 对于一个方阵A，如果存在一个与其同阶的方阵B，使得： (I为单位矩阵） 则称B为A的逆矩阵，当然，A也是B的逆矩阵。求方阵A的逆矩阵 可调用函数inv(A)。 例3.4 求方阵A的逆矩阵，且验证。 A = [1,-1,1;5,-4,3;2,1,1]; B = inv(A); A*B ans = 1.0000 0 0 -0.0000 1.0000 0 -0.0000 0 1.0000 3.3.2 用矩阵求逆方法求解线性方程组 将包含n个未知数，由n个方程构成的线性方程组表示为： 其矩阵表示形式为： 其中： 因此： 例3.3 用求逆矩阵A的方法解线性方程组 命令如下： A = [1,2,3;1,4,9;1,8,27]; b = [5,-2,