MATLAB矩阵处理（二）

矩阵变换

1．对角阵

 对角阵：只有对角线上有非零元素的矩阵。

 数量矩阵：对角线上的元素相等的对角矩阵。

 单位矩阵：对角线上的元素都为1的对角矩阵。

(1) 提取矩阵的对角线元素

 diag(A)：提取矩阵A主对角线元素，产生一个列向量。

 diag(A,k)：提取矩阵A第k条对角线的元素，产生一个列向量。

矩阵的对角线：与主对角线平行，往上为第1条、第2条、一直到第n条对角线，往下为第-1条、-2条、一直到-n条对角线。主对角线为第0条对角线。

(2) 构造对角阵

 diag(V)：以向量 V为主对角线元素，产生对角矩阵。

 diag(V,k)：以向量 V为第k条对角线元素，产生对角矩阵。

2．三角阵

 上三角阵：矩阵的对角线以下的元素全为零的矩阵。

 下三角阵：对角线以上的元素全为零的矩阵。

（1）上三角矩阵

 triu(A)：提取矩阵A的主对角线及以上的元素。

 triu(A,k)：提取矩阵A的第k条对角线及以上的元素。

(2) 下三角矩阵

在MATLAB中，提取矩阵A的下三角矩阵的函数是tril，其用法与triu函

数完全相同。

 

3．矩阵的转置

 转置运算符是小数点后面接单引号（.'）。

 共轭转置，其运算符是单引号（'），它在转置的基础上还要取每个数的复共轭。

若矩阵元素为实数，则转置与共轭转置结果相同

4．矩阵的旋转

rot90(A,k)：将矩阵A逆时针方向旋转90º的k倍，当k为1时可省略。

5．矩阵的翻转

对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第 二列调换，…，依此类推。

 fliplr(A)：对矩阵A实施左右翻转

 flipud(A)：对矩阵A实施上下翻转。

6．矩阵的求逆

 对于一个方阵A，如果存在一个与其同阶的方阵B，使得AB=BA=I (I为单位矩阵)，则称B为A的逆矩阵，当然，A也是B的逆矩阵。

 inv(A)：求方阵A的逆矩阵。

矩阵求值

1．方阵的行列式

 把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为所对应的行列式的值。

 det(A)：求方阵A所对应的行列式的值。

2．矩阵的秩

 矩阵线性无关的行数或列数称为矩阵的秩。

 rank(A)：求矩阵A的秩。

 

 奇数阶魔方阵秩为n，即奇数阶魔方阵是满秩矩阵。

 一重偶数阶魔方阵秩为n/2+2（n是2的倍 数，但非4的倍数）。

 双重偶数阶魔方阵秩均为3 （阶数是4的 倍数） 。

3．矩阵的迹

 矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。

 trace(A)：求矩阵A的迹。

直接求：trace(A)

求对角线元素之和：sum(diag(A))

求特征值的和：sum(eig(A))

4．向量和矩阵的范数

矩阵或向量的范数用来度量矩阵或向量在某种意义下的长度。

在MATLAB中，求向量范数的函数为：

 norm(V)或norm(V,2)：计算向量V的2—范数。

 norm(V,1)：计算向量V的1—范数。

 norm(V,inf)：计算向量V的∞—范数。

5．矩阵的条件数

 矩阵A的条件数等于A的范数与A的逆矩阵的范数的乘积。

 条件数越接近于1，矩阵的性能越好，反之，矩阵的性能越差。

在MATLAB中，求向量条件数的函数为：

 cond(A)或cond(A,2)：计算A的2—范数下的条件数

 norm(V,1)：计算A的1—范数下的条件数。

 norm(V,inf)：计算A的∞—范数下的条件数。