写在前面: 这个帖子主要是整理我常用到的matlab矩阵运算,给自己作为笔记,也提供给需要的同学们看。里面的内容在很多别的博客也都有介绍,我开始用的时候也是跟他们学习的,在此表示感谢。能附上链接的我都会尽量附上。
1. 矩阵合并
a. 矩阵沿行或列合并
按行合并:
C=[A,B];
按列合并:
C=[A;B];
b. 沿对角线合并,维度由原来两矩阵的2*2变为4*4.
C=blkdiag(A,B);
2. 矩阵的Kronecher 积
注:定义来自维基百科 https://en.wikipedia.org/wiki/Kronecker_product
C=kron(A,B);
3. 三对角矩阵
三对角矩阵就是对角线、邻近对角线的上下次对角线上有元素,其他位置均为0的矩阵。形式如下:matlab实现三对角矩阵的输入_Flipped_新浪博客
主对角线元素生成:
A0=5*diag(ones(5,1));
注:ones(m, n)生成m*n维的全1矩阵;zeros(m, n) 生成 m*n维的全0矩阵;
上次对角线元素生成
A1=diag(ones(4,1),1);
下次对角线元素生成
A2=diag(ones(4,1),-1);
三对角矩阵
A=5*diag(ones(5,1))+diag(ones(4,1),1)+diag(ones(4,1),-1);
通过将上面的基本矩阵生成语句组合就可以得到我们所需要的一些特殊形式的矩阵。
4. 应用实例
我在实际应用中遇到了如下形式的矩阵
这在广义上来说是一个三对角矩阵,三对角线上的元素又分别由2*2维的矩阵构成
但是特殊的是,对角线上的矩阵unit在奇数行与偶数行又不相同,因此这种只能用编程语句来实现。
首先构造矩阵元素即单位2*2矩阵,这里共有四种元素:
i_unit=[i w; w i];
u_unit=[0 u; u 0];
l_unit=[l 0; 0 l];
c_unit=[c 0; 0 c];
对于主对角元素可以把i_unit 和i_unit 沿对角线扩展为一个4*4的矩阵,再将其作为单位矩阵;
ui_unit=blkdiag(i_unit,u_unit);
N=2;
A0=kron(diag(ones(N,1)),ui_unit);
对于次对角元素,奇数行与偶数行要分开处理。不妨构造四个4*4矩阵,如下所示,再用Kronecker积将相应元素填充到相应位置。
N=2;
A1=zeros(2*N,2*N);
A2=zeros(2*N,2*N);
A3=zeros(2*N,2*N);
A4=zeros(2*N,2*N);
for i=1:N
A1(2*i-1,2*i)=-1; //上三角奇数行
A3(2*i,2*i-1)=1; //下三角偶数行
end
for i=1:N-1
A2(2*i,2*i+1)=-1; //上三角偶数行,注意最后一行是没有上三角元素的,所以到2*(N-1)截止
A4(2*i+1,2*i)=1; // 下三角奇数行
end
A=A0+kron(A1,l_unit)+kron(A2,c_unit)+kron(A3,c_unit)+kron(A4,l_unit);
对于更高维度的矩阵也是同样的方法构造,只要把N值增大即可,这个矩阵可以看做是4N*4N维。
希望对你有用!