2.2 矩阵变换

写在前面：作者本人是纯纯的菜鸟，学习的内容来自于 中国大学MOOC 中南大学 《科学计算与MATLAB语言》，欢迎各位大佬或新手在这里和平讨论，如果我有错误请各位不吝赐教，提前感谢各位捧场！

---

何为矩阵变化？

指对矩阵进行某种运算与操作，其结果还是一个矩阵。

今天要讲到的矩阵变换：

• 对角阵
• 三角阵
• 矩阵的转置
• 矩阵的翻转
• 矩阵的旋转
• 矩阵求逆

1.对角阵

对角线上有非零元素的矩阵称为对角阵。

数量矩阵：对角线上数字相同的对角矩阵。

单位矩阵：对角线上数字全为1的对角矩阵。

MATLAB可能用到的操作：提取一个矩阵对角线上的元素或以一个向量构造一个对角矩阵。分别讨论：

（1）提取对角线上的元素：

diag(A,k)
%若提取主对角线上的元素，k可以省略
%若提取与主对角线平行且在主对角线之上的对角线的元素，k为正数；反之，k为负数。
%主对角线默认k为0

（2）利用已有向量V建立对角矩阵

diag(V,k)
%与提取时用法类似，不做说明

例题：构造一个5阶方阵，若要使该方阵每行元素乘以对应的行数应如何写命令？若要使该方阵每列元素乘以对应列数该如何写命令？

思路：利用对角矩阵，通过对原矩阵左乘或右乘达到要求。

A=rand(5);
%生成一个5阶随机方阵
D=diag(1:5);
%构造对角矩阵，对角线元素分别为1~5
B=D*A
%第1题：使A矩阵每行与对应行数相乘
C=A*D
%第2题：使A矩阵每列与对应列数相乘

2.三角阵

三角针的分类：

• 上三角阵：对角线以下全为0的矩阵
• 下三角阵：对角线以上全为0的矩阵

求矩阵的上三角阵和下三角阵

triu(A,k)
%u代表up，求上三角矩阵
%该命令代表求矩阵A第k条对角线及以上的元素
%k=0可以省略
tril(A,k)
%l代表low，求下三角矩阵
%其余参数使用与triu完全相同

3.矩阵的转置

矩阵的转置：将原矩阵的同一行列对应元素互换，即第一行转为第一列，或者（个人理解）将原矩阵关于主对角线对称就可得到该矩阵的转置矩阵。

转置的运算符：

（1）小数点+单引号  .'

（2）共轭转置：单引号 '，他在转置的基础上还要将原矩阵每个元素求共轭

4.矩阵的旋转

所需函数：rot90()

rot90(A,k)
%个人理解：rot -> rotate 旋转
%将A矩阵逆时针旋转90°的k倍，k=1时可省略

5.矩阵的翻转

翻转，一定有两种方向：（1）与翻书类似，左右翻转；（2）与翻作业本类似，上下翻转

fliplr(A)
%flip -> 翻转，l->left r->right ，该命令表示左右翻转
flipud(A)
%u->up d->down ，该命令表示上下翻转

例题：验证魔方矩阵主对角线和副对角线元素之和相等（以5阶魔方矩阵为例）

思路：主对角线自不必说，用diag命令可提出，对于副对角线，可以令原矩阵上下翻转形成一个新的矩阵，提出该矩阵的主对角线即为原矩阵的副对角线

A=magic(5);
D1=diag(A)
B=flipud(A);
D2=diag(B);
(sum(D1))==(sum(D2));

6.矩阵求逆

逆矩阵的定义：对于方阵A，若存在一个与其同阶的方阵B，使得A*B=B*A=I（I为单位矩阵），则称B为A的逆矩阵，A也为B的逆矩阵。

计算逆矩阵很复杂，但在MATLAB中仅需一条命令inv(A)即可得到。

例题：用求逆的方法解线性方程组

思路：线性方程组简化为Ax=b，其中，b=[5 -2 6]

方程两边左乘即可得出x的结果.也可以利用左除运算符求出x

A=[1 2 3;1 4 9;1 8 27];
b=[5 -2 6];
x=inv(A)*b；
%思路1：求逆矩阵
x=A\b;
%思路2：利用左除运算符

---

以上就是今天的全部内容了，再次感谢各位耐心观看！