Matlab 矩阵运算
说明：这一段时间用Matlab做了LDPC码的性能仿真，过程中涉及了大量的矩阵运算，本文记录了Matlab中矩阵的相关知识，特别的说明了稀疏矩阵和有限域中的矩阵。Matlab的运算是在矩阵意义下进行的，这里所提到的是狭义上的矩阵，即通常意义上的矩阵。
目录
目录
第一部分：矩阵基本知识
一、矩阵的创建
一、矩阵的创建
1. 直接输入法
2. 利用Matlab函数创建矩阵
3. 利用文件创建矩阵
二、矩阵的拆分
1. 矩阵元素
2. 矩阵拆分
3. 特殊矩阵
三、矩阵的运算
1. 算术运算
2. 关系运算
3. 逻辑运算
四、矩阵分析
1. 对角阵
linspace函数调用的格式为2. 三角阵
3. 矩阵的转置与旋转
4. 矩阵的翻转
5. 矩阵的逆与伪逆
6. 方阵的行列式
7. 矩阵的秩与迹
8. 向量和矩阵的范数
9. 矩阵的特征值与特征向量
五、字符串
六、其他
六、其他
第二部分 矩阵的应用
一、稀疏矩阵
一、稀疏矩阵
1. 稀疏矩阵的创建
2. 稀疏矩阵的运算
3. 其他
二、有限域中的矩阵
内容
第一部分：矩阵基本知识（只作基本介绍，详细说明请参考Matlab帮助文档）
矩阵是进行数据处理和运算的基本元素。在MATLAB中
a、通常意义上的数量（标量）可看成是”1*1″的矩阵；
b、n维矢量可看成是”n*1″的矩阵；
c、多项式可由它的系数矩阵完全确定。
a、通常意义上的数量（标量）可看成是”1*1″的矩阵；
b、n维矢量可看成是”n*1″的矩阵；
c、多项式可由它的系数矩阵完全确定。
一、矩阵的创建
在MATLAB中创建矩阵有以下规则：
a、矩阵元素必须在”[ ]“内；
b、矩阵的同行元素之间用空格（或”,”）隔开；
a、矩阵元素必须在”[ ]“内；
b、矩阵的同行元素之间用空格（或”,”）隔开；
c、矩阵的行与行之间用”;”（或回车符）隔开；
d、矩阵的元素可以是数值、变量、表达式或函数；
e、矩阵的尺寸不必预先定义。
d、矩阵的元素可以是数值、变量、表达式或函数；
e、矩阵的尺寸不必预先定义。
下面介绍四种矩阵的创建方法：
1、直接输入法
最简单的建立矩阵的方法是从键盘直接输入矩阵的元素，输入的方法按照上面的规则。建立向量的时候可以利用冒号表达式，冒号表达式可以产生一个行向量，一般格式是： e1:e2:e3，其中e1为初始值，e2为步长，e3为终止值。 还可以用linspace函数产生行向量，其调用格式为：linspace(a,b,n) ，其中a和b是生成向量的第一个和最后一个元素，n是元素总数。 可以看出来linspace(a,b,n)与a:(b-a)/(n-1):b等价。
最简单的建立矩阵的方法是从键盘直接输入矩阵的元素，输入的方法按照上面的规则。建立向量的时候可以利用冒号表达式，冒号表达式可以产生一个行向量，一般格式是： e1:e2:e3，其中e1为初始值，e2为步长，e3为终止值。 还可以用linspace函数产生行向量，其调用格式为：linspace(a,b,n) ，其中a和b是生成向量的第一个和最后一个元素，n是元素总数。 可以看出来linspace(a,b,n)与a:(b-a)/(n-1):b等价。
2、利用MATLAB函数创建矩阵
基本矩阵函数如下：
(1) ones()函数：产生全为1的矩阵，ones(n)：产生n*n维的全1矩阵，ones(m,n)：产生m*n维的全1矩阵；
基本矩阵函数如下：
(1) ones()函数：产生全为1的矩阵，ones(n)：产生n*n维的全1矩阵，ones(m,n)：产生m*n维的全1矩阵；
(2) zeros()函数：产生全为0的矩阵；
(3) rand()函数：产生在（0，1）区间均匀分布的随机阵；
(4) eye()函数：产生单位阵；
(5) randn()函数：产生均值为0，方差为1的标准正态分布随机矩阵。
(3) rand()函数：产生在（0，1）区间均匀分布的随机阵；
(4) eye()函数：产生单位阵；
(5) randn()函数：产生均值为0，方差为1的标准正态分布随机矩阵。
3、利用文件建立矩阵
当矩阵尺寸较大或为经常使用的数据矩阵，则可以将此矩阵保存为文件，在需要时直接将文件利用load命令调入工作环境中使用即可。同时可以利用命令reshape对调入的矩阵进行重排。reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m*n的二维矩阵。
当矩阵尺寸较大或为经常使用的数据矩阵，则可以将此矩阵保存为文件，在需要时直接将文件利用load命令调入工作环境中使用即可。同时可以利用命令reshape对调入的矩阵进行重排。reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m*n的二维矩阵。
二、矩阵的拆分
1．矩阵元素
可以通过下标（行列索引）引用矩阵的元素，如 Matrix(m,n)。也可以采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在MATLAB中，矩阵元素按列存储，先第一列，再第二列，依次类推。序号(Index)与下标(Subscript )是一一
可以通过下标（行列索引）引用矩阵的元素，如 Matrix(m,n)。也可以采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在MATLAB中，矩阵元素按列存储，先第一列，再第二列，依次类推。序号(Index)与下标(Subscript )是一一
对应的，以m*n矩阵A为例，矩阵元素A(i,j)的序号为(j-1)*m+i。其相互转换关系也可利用sub2ind和ind2sub函数求得。
2．矩阵拆分
利用冒号表达式获得子矩阵：
(1) A(:,j)表示取A矩阵的第j列全部元素；A(i,:)表示A矩阵第i行的全部元素；A(i,j)表示取A矩阵第i行、第j列的元素。
(2) A(i:i+m,:)表示取A矩阵第i~i+m行的全部元素；A(:,k:k+m)表示取A矩阵第k~k+m列的全部元素，A(i:i+m,k:k+m)表示取A矩阵第i~i+m行内，并在第k~k+m列中的所有元素。此外，还可利用一般向量和end运算符来表示矩阵下标，从而获得子矩阵。end表示某一维的末尾元素下标。
利用空矩阵删除矩阵的元素：
在MATLAB中，定义[]为空矩阵。给变量X赋空矩阵的语句为X=[]。注意，X=[]与clear X不同，clear是将X从工作空间中删除，而空矩阵则存在于工作空间中，只是维数为0。
利用冒号表达式获得子矩阵：
(1) A(:,j)表示取A矩阵的第j列全部元素；A(i,:)表示A矩阵第i行的全部元素；A(i,j)表示取A矩阵第i行、第j列的元素。
(2) A(i:i+m,:)表示取A矩阵第i~i+m行的全部元素；A(:,k:k+m)表示取A矩阵第k~k+m列的全部元素，A(i:i+m,k:k+m)表示取A矩阵第i~i+m行内，并在第k~k+m列中的所有元素。此外，还可利用一般向量和end运算符来表示矩阵下标，从而获得子矩阵。end表示某一维的末尾元素下标。
利用空矩阵删除矩阵的元素：
在MATLAB中，定义[]为空矩阵。给变量X赋空矩阵的语句为X=[]。注意，X=[]与clear X不同，clear是将X从工作空间中删除，而空矩阵则存在于工作空间中，只是维数为0。
3、特殊矩阵
(1) 魔方矩阵 魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵，其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n)，其功能是生成一个n阶魔方阵。
(2) 范得蒙矩阵 范得蒙(Vandermonde)矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中，函数vander(V)生成以向量V为基础向量的范得蒙矩阵。
(3) 希尔伯特矩阵 在MATLAB中，生成希尔伯特矩阵的函数是hilb(n)。 使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中，有一个专门求希尔伯特矩阵的逆的函数invhilb(n)，其功能是求n阶的希尔伯特矩阵的逆矩阵。
(4) 托普利兹矩阵 托普利兹(Toeplitz)矩阵除第一行第一列外，其他每个元素都与左上角的元素相同。生成托普利兹矩阵的函数是toeplitz(x,y)，它生成一个以x为第一列，y为第一行的托普利兹矩阵。这里x, y均为向量，两者不必等长。toeplitz(x)用向量x生成一个对称的托普利兹矩阵。
(5) 伴随矩阵 MATLAB生成伴随矩阵的函数是compan(p)，其中p是一个多项式的系数向量，高次幂系数排在前，低次幂排在后。
(2) 范得蒙矩阵 范得蒙(Vandermonde)矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中，函数vander(V)生成以向量V为基础向量的范得蒙矩阵。
(3) 希尔伯特矩阵 在MATLAB中，生成希尔伯特矩阵的函数是hilb(n)。 使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中，有一个专门求希尔伯特矩阵的逆的函数invhilb(n)，其功能是求n阶的希尔伯特矩阵的逆矩阵。
(4) 托普利兹矩阵 托普利兹(Toeplitz)矩阵除第一行第一列外，其他每个元素都与左上角的元素相同。生成托普利兹矩阵的函数是toeplitz(x,y)，它生成一个以x为第一列，y为第一行的托普利兹矩阵。这里x, y均为向量，两者不必等长。toeplitz(x)用向量x生成一个对称的托普利兹矩阵。
(5) 伴随矩阵 MATLAB生成伴随矩阵的函数是compan(p)，其中p是一个多项式的系数向量，高次幂系数排在前，低次幂排在后。
(6) 帕斯卡矩阵 我们知道，二次项(x+y)n展开后的系数随n的增大组成一个三角形表，称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵。函数pascal(n)生成一个n阶帕斯卡矩阵。
三、矩阵的运算
1、算术运算
MATLAB的基本算术运算有：＋(加)、－(减)、*(乘)、/(右除)、\(左除)、^(乘方)、’(转置)。 运算是在矩阵意义下进行的，单个数据的算术运算只是一种特例。
(1) 矩阵加减运算 假定有两个矩阵A和B，则可以由A+B和A-B实现矩阵的加减运算。运算规则是：若A和B矩阵的维数相同，则可以执行矩阵的加减运算，A和B矩阵的相应元素相加减。如果A与B的维数不相同，则MATLAB将给出错误信息，提示用户两个矩阵的维数不匹配。
(2) 矩阵乘法 假定有两个矩阵A和B，若A为m*n矩阵，B为n*p矩阵，则C=A*B为m*p矩阵。
(3) 矩阵除法 在MATLAB中，有两种矩阵除法运算：\和/，分别表示左除和右除。如果A矩阵是非奇异方阵，则A\B和B/A运算可以实现。A\B等效于A的逆左乘B矩阵，也就是inv(A)*
MATLAB的基本算术运算有：＋(加)、－(减)、*(乘)、/(右除)、\(左除)、^(乘方)、’(转置)。 运算是在矩阵意义下进行的，单个数据的算术运算只是一种特例。
(1) 矩阵加减运算 假定有两个矩阵A和B，则可以由A+B和A-B实现矩阵的加减运算。运算规则是：若A和B矩阵的维数相同，则可以执行矩阵的加减运算，A和B矩阵的相应元素相加减。如果A与B的维数不相同，则MATLAB将给出错误信息，提示用户两个矩阵的维数不匹配。
(2) 矩阵乘法 假定有两个矩阵A和B，若A为m*n矩阵，B为n*p矩阵，则C=A*B为m*p矩阵。
(3) 矩阵除法 在MATLAB中，有两种矩阵除法运算：\和/，分别表示左除和右除。如果A矩阵是非奇异方阵，则A\B和B/A运算可以实现。A\B等效于A的逆左乘B矩阵，也就是inv(A)*
B，而B/A等效于A矩阵的逆右乘B矩阵，也就是B*inv(A)。 对于含有标量的运算，两种除法运算的结果相同。对于矩阵来说，左除和右除表示两种不同的除数矩阵和被除数矩阵的关系，一般A\B≠B/A。
(4) 矩阵的乘方 一个矩阵的乘方运算可以表示成A^x，要求A为方阵，x为标量。
(5) 矩阵的转置 对实数矩阵进行行列互换，对复数矩阵，共轭转置，特殊的，操作符.’共轭不转置（见点运算）；
(6) 点运算 在MATLAB中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。点运算符有.*、./、.\和.^。两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。
(4) 矩阵的乘方 一个矩阵的乘方运算可以表示成A^x，要求A为方阵，x为标量。
(5) 矩阵的转置 对实数矩阵进行行列互换，对复数矩阵，共轭转置，特殊的，操作符.’共轭不转置（见点运算）；
(6) 点运算 在MATLAB中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。点运算符有.*、./、.\和.^。两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。
2、关系运算
MATLAB提供了6种关系运算符：<(小于)、<=(小于或等于)、>(大于)、>=(大于或等于)、==(等于)、~=(不等于)。关系运算符的运算法则为：
(1) 当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为1，否则为0；
(2) 当参与比较的量是两个维数相同的矩阵时，比较是对两矩阵相同位置的元素按标量关系
MATLAB提供了6种关系运算符：<(小于)、<=(小于或等于)、>(大于)、>=(大于或等于)、==(等于)、~=(不等于)。关系运算符的运算法则为：
(1) 当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为1，否则为0；
(2) 当参与比较的量是两个维数相同的矩阵时，比较是对两矩阵相同位置的元素按标量关系
运算规则逐个进行，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成；
(3) 当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。
(3) 当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。
3、逻辑运算
MATLAB提供了3种逻辑运算符：&(与)、|(或)和~(非)。 逻辑运算的运算法则为：
(1) 在逻辑运算中，确认非零元素为真，用1表示，零元素为假，用0表示；
(2) 设参与逻辑运算的是两个标量a和b，那么，a&b a,b全为非零时，运算结果为1，否则为0。 a|b a,b中只要有一个非零，运算结果为1。~a 当a是零时，运算结果为1；当a非零时，运算结果为0。
(3) 若参与逻辑运算的是两个同维矩阵，那么运算将对矩阵相同位置上的元素按标量规则逐个进行。最终运算结果是一个与原矩阵同维的矩阵，其元素由1或0组成；
(4) 若参与逻辑运算的一个是标量，一个是矩阵，那么运算将在标量与矩阵中的每个元素之间按标量规则逐个进行。最终运算结果是一个与矩阵同维的矩阵，其元素由1或0组成；
MATLAB提供了3种逻辑运算符：&(与)、|(或)和~(非)。 逻辑运算的运算法则为：
(1) 在逻辑运算中，确认非零元素为真，用1表示，零元素为假，用0表示；
(2) 设参与逻辑运算的是两个标量a和b，那么，a&b a,b全为非零时，运算结果为1，否则为0。 a|b a,b中只要有一个非零，运算结果为1。~a 当a是零时，运算结果为1；当a非零时，运算结果为0。
(3) 若参与逻辑运算的是两个同维矩阵，那么运算将对矩阵相同位置上的元素按标量规则逐个进行。最终运算结果是一个与原矩阵同维的矩阵，其元素由1或0组成；
(4) 若参与逻辑运算的一个是标量，一个是矩阵，那么运算将在标量与矩阵中的每个元素之间按标量规则逐个进行。最终运算结果是一个与矩阵同维的矩阵，其元素由1或0组成；
(5) 逻辑非是单目运算符，也服从矩阵运算规则；
(6) 在算术、关系、逻辑运算中，算术运算优先级最高，逻辑运算优先级最低。
(6) 在算术、关系、逻辑运算中，算术运算优先级最高，逻辑运算优先级最低。
四、矩阵分析
1、对角阵
(1) 对角阵 只有对角线上有非0元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为1的对角矩阵称为单位矩阵。
(1) 提取矩阵的对角线元素 设A为m*n矩阵，diag(A)函数用于提取矩阵A主对角线元素，产生一个具有min(m,n)个元素的列向量。diag(A)函数还有一种形式diag(A,k)，其功能是提取第k条对角线的元素。
(2) 构造对角矩阵 设V为具有m个元素的向量，diag(V)将产生一个m*m对角矩阵，其主对角线元素即为向量V的元素。diag(V)函数也有另一种形式diag(V,k)，其功能是产生一个n*n(n=m+k)对角阵，其第m条对角线的元素即为向量V的元素。
(1) 对角阵 只有对角线上有非0元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为1的对角矩阵称为单位矩阵。
(1) 提取矩阵的对角线元素 设A为m*n矩阵，diag(A)函数用于提取矩阵A主对角线元素，产生一个具有min(m,n)个元素的列向量。diag(A)函数还有一种形式diag(A,k)，其功能是提取第k条对角线的元素。
(2) 构造对角矩阵 设V为具有m个元素的向量，diag(V)将产生一个m*m对角矩阵，其主对角线元素即为向量V的元素。diag(V)函数也有另一种形式diag(V,k)，其功能是产生一个n*n(n=m+k)对角阵，其第m条对角线的元素即为向量V的元素。
2、三角阵
三角阵又进一步分为上三角阵和下三角阵，所谓上三角阵，即矩阵的对角线以下的元素全
三角阵又进一步分为上三角阵和下三角阵，所谓上三角阵，即矩阵的对角线以下的元素全
为0的一种矩阵，而下三角阵则是对角线以上的元素全为0的一种矩阵。
(1) 上三角矩阵 求矩阵A的上三角阵的MATLAB函数是triu(A)。 triu(A)函数也有另一种形式triu(A,k)，其功能是求矩阵A的第k条对角线以上的元素。
(2) 下三角矩阵 在MATLAB中，提取矩阵A的下三角矩阵的函数是tril(A)和tril(A,k)，其用法与提取上三角矩阵的函数triu(A)和triu(A,k)完全相同。
(1) 上三角矩阵 求矩阵A的上三角阵的MATLAB函数是triu(A)。 triu(A)函数也有另一种形式triu(A,k)，其功能是求矩阵A的第k条对角线以上的元素。
(2) 下三角矩阵 在MATLAB中，提取矩阵A的下三角矩阵的函数是tril(A)和tril(A,k)，其用法与提取上三角矩阵的函数triu(A)和triu(A,k)完全相同。
3、矩阵的转置与旋转
(1) 矩阵的转置 转置运算符是单撇号(’)。
(2) 矩阵的旋转 利用函数rot90(A,k)将矩阵A旋转90o的k倍，当k为1时可省略。
(1) 矩阵的转置 转置运算符是单撇号(’)。
(2) 矩阵的旋转 利用函数rot90(A,k)将矩阵A旋转90o的k倍，当k为1时可省略。
4、矩阵的翻转
对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第二列调换，…，依次类推。矩阵A实施左右翻转的函数是fliplr(A)，对矩阵A实施上下翻转的函数是flipud(A)。
对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第二列调换，…，依次类推。矩阵A实施左右翻转的函数是fliplr(A)，对矩阵A实施上下翻转的函数是flipud(A)。
5、矩阵的逆与伪逆
(1) 矩阵的逆 对于一个方阵A，如果存在一个与其同阶的方阵B，使得：AB=BA=I (I为单位
(1) 矩阵的逆 对于一个方阵A，如果存在一个与其同阶的方阵B，使得：AB=BA=I (I为单位
矩阵) 则称B为A的逆矩阵，当然，A也是B的逆矩阵。求方阵A的逆矩阵可调用函数inv(A)。
(2) 矩阵的伪逆 如果矩阵A不是一个方阵，或者A是一个非满秩的方阵时，矩阵A没有逆矩阵，但可以到一个与A的转置矩阵A’同型的矩阵B，使得：ABA=A，BAB=B 此时称矩阵B为矩阵A的伪逆，也称为广义逆矩阵。在MATLAB中，求一个矩阵伪逆的函数是pinv(A)。
(2) 矩阵的伪逆 如果矩阵A不是一个方阵，或者A是一个非满秩的方阵时，矩阵A没有逆矩阵，但可以到一个与A的转置矩阵A’同型的矩阵B，使得：ABA=A，BAB=B 此时称矩阵B为矩阵A的伪逆，也称为广义逆矩阵。在MATLAB中，求一个矩阵伪逆的函数是pinv(A)。
6、方阵的行列式
把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为矩阵所对应的行列式的值。在MATLAB中，求方阵A所对应的行列式的值的函数是det(A)。
把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为矩阵所对应的行列式的值。在MATLAB中，求方阵A所对应的行列式的值的函数是det(A)。
7、矩阵的秩与迹
(1) 矩阵的秩 矩阵线性无关的行数与列数称为矩阵的秩。在MATLAB中，求矩阵秩的函数是rank(A)。
(2) 矩阵的迹 矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。在MATLAB中，求矩阵的迹的函数是trace(A)。
(1) 矩阵的秩 矩阵线性无关的行数与列数称为矩阵的秩。在MATLAB中，求矩阵秩的函数是rank(A)。
(2) 矩阵的迹 矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。在MATLAB中，求矩阵的迹的函数是trace(A)。
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论