托普利兹矩阵(T矩阵)及其应用(Matlab demo测试)

托普利兹矩阵(T矩阵)及其应用(Matlab demo测试)

  • [1. 概念](#1. 概念)
  • [2. Matlab简单测试](#2. Matlab简单测试)
    • [2.1 生成测试](#2.1 生成测试)
    • [2.2 基本性质及原理](#2.2 基本性质及原理)
    • [2.3 性质验证](#2.3 性质验证)
  • [3. 其他应用总结](#3. 其他应用总结)
    • [3.1 其他性质](#3.1 其他性质)
    • [3.2 文献阅读看到的](#3.2 文献阅读看到的)
  • 参考资料

1. 概念

托普利兹矩阵,简称为T型矩阵,托普利兹矩阵的主对角线上的元素相等,平行于主对角线的线上的元素也相等;矩阵中的各元素关于次对角线对称,即T型矩阵为次对称矩阵。即 a i j = a j i a_{ij}=a_{ji} aij=aji

2. Matlab简单测试

2.1 生成测试

Matlab中可以用toeplitz(x,y)。它生成一个以 x 为第一列,y 为第一行的托普利兹矩阵。

函数中x=(x1,x2,...,xk) y=(y1,y2,...,yj)为向量形式,代表托普利兹矩阵的第一行、第一列。

matlab 复制代码
x=[1, 2, 3, 3, 4, 4];
y=[1, 3, 3, 2, 3, 4];
T=toeplitz(x,y)

生成结果如下:

ans =

     1     3     3     2     3     4
     2     1     3     3     2     3
     3     2     1     3     3     2
     3     3     2     1     3     3
     4     3     3     2     1     3
     4     4     3     3     2     1

2.2 基本性质及原理

其中,最基础的性质,是托普利兹矩阵可以表示为前向位移矩阵和后向位移矩阵之和。

  • 前向位移矩阵
    F = ( 0 1 . . . 0 0 . . . . . . . . . . . . . . . . . . 1 0 . . . 0 0 ) ∈ R n × n F=\left( \begin{matrix} 0& 1& ...& 0\\ 0& ...& ...& ...\\ ...& ...& ...& 1\\ 0& ...& 0& 0\\ \end{matrix} \right) \in \mathbb{R} ^{n\times n} F= 00...01..................00...10 ∈Rn×n
  • 后向位移矩阵
    B = ( 0 0 . . . 0 1 . . . . . . . . . . . . . . . . . . 0 0 . . . 1 0 ) ∈ R n × n B=\left( \begin{matrix} 0& 0& ...& 0\\ 1& ...& ...& ...\\ ...& ...& ...& 0\\ 0& ...& 1& 0\\ \end{matrix} \right) \in \mathbb{R} ^{n\times n} B= 01...00..................10...00 ∈Rn×n
  • 基于性质 前向、后向矩阵幂次和

T = ∑ k − 1 n − 1 t − k B k + ∑ k = 0 n − 1 t k F k    T=\sum_{k-1}^{n-1}{t_{-k}B^k+\sum_{k=0}^{n-1}{t_kF^k}}\,\, T=k−1∑n−1t−kBk+k=0∑n−1tkFk

2.3 性质验证

  • 简单前向后向矩阵 后向矩阵 的幂次性质
matlab 复制代码
n = 5; % Define the size of the matrix
F = diag(ones(1, n-1), 1); % Create the forward matrix
B = F'

这性质确实有点意思... 位置变化了

matlab 复制代码
>> B^2

ans =

     0     0     0     0     0
     0     0     0     0     0
     1     0     0     0     0
     0     1     0     0     0
     0     0     1     0     0

>> B^3

ans =

     0     0     0     0     0
     0     0     0     0     0
     0     0     0     0     0
     1     0     0     0     0
     0     1     0     0     0

>> B^4

ans =

     0     0     0     0     0
     0     0     0     0     0
     0     0     0     0     0
     0     0     0     0     0
     1     0     0     0     0

>> F^2

ans =

     0     0     1     0     0
     0     0     0     1     0
     0     0     0     0     1
     0     0     0     0     0
     0     0     0     0     0

>> F^3

ans =

     0     0     0     1     0
     0     0     0     0     1
     0     0     0     0     0
     0     0     0     0     0
     0     0     0     0     0

>> F^4

ans =

     0     0     0     0     1
     0     0     0     0     0
     0     0     0     0     0
     0     0     0     0     0
     0     0     0     0     0
  • 生成 托普利兹矩阵
matlab 复制代码
n = 5; % Define the size of the matrix
F = diag(ones(1, n-1), 1); % Create the forward matrix
B = F';

% Define the coefficients t_{-k} and t_k
t_neg = [1, 2, 3, 4, 5]; % Example coefficients for t_{-k}
t_pos = [1, 3, 3, 2, 1]; % Example coefficients for t_k

T = zeros(n); % Initialize the Toeplitz matrix

for k = 1:n
    T = T + t_neg(k) * (B^(k-1));
end

for k = 2:n
    T = T + t_pos(k) * (F^(k-1));
end

定义的信息如下:

t_neg = [1, 2, 3, 4, 5]; % Example coefficients for t_{-k}

t_pos = [1, 3, 3, 2, 1]; % Example coefficients for t_k

matlab 复制代码
T =

     1     3     3     2     1
     2     1     3     3     2
     3     2     1     3     3
     4     3     2     1     3
     5     4     3     2     1

3. 其他应用总结

3.1 其他性质

  • Python实现版本可以参考哈工大 赵老师的博客。

  • 其他的一些性质,

    • 包括可以高效率的计算卷积...
    • 对于Ax=b的系统(线性代数中),当A为托普利兹矩阵时,可以称其为托普利兹系统, 且此时的系统自由度为2-1而不是n^2, (究其原因,和托普利兹矩阵的形式有关), 因此,可以用Levinson求解方法快速计算
    • 托普利兹矩阵可以被分解,如LU分解中的Bareiss算法

PS: LU分解,顾名思义,L 是单位下三角矩阵, U 是单位上三角矩阵。 LU分解有两种实现,分别是. Gauss消去法. 待定系数法.

    • 关于对称块矩阵(Block Toepliz)和对称矩阵(Toepliz) 虽然托普利茨矩阵具有与对角线恒定性相关的特定特征,但对称块矩阵的特征在于其子矩阵的对称性。

这些具体的性质,等到需要用的时候,再推导吧...

3.2 文献阅读看到的

对于一些工程应用,最近在一篇论文中,就用到了这个性质,需要分析一个能量传播矩阵,这个能量传播矩阵可以表示为一个近似的对称块托普利兹矩阵,因此,可以利用其卷积性质,得到不变卷积核:

参考资料

【1】-csdn 托普利兹矩阵

相关推荐
顶呱呱程序7 分钟前
2-143 基于matlab-GUI的脉冲响应不变法实现音频滤波功能
算法·matlab·音视频·matlab-gui·音频滤波·脉冲响应不变法
爱吃生蚝的于勒28 分钟前
深入学习指针(5)!!!!!!!!!!!!!!!
c语言·开发语言·数据结构·学习·计算机网络·算法
羊小猪~~32 分钟前
数据结构C语言描述2(图文结合)--有头单链表,无头单链表(两种方法),链表反转、有序链表构建、排序等操作,考研可看
c语言·数据结构·c++·考研·算法·链表·visual studio
王哈哈^_^1 小时前
【数据集】【YOLO】【VOC】目标检测数据集,查找数据集,yolo目标检测算法详细实战训练步骤!
人工智能·深度学习·算法·yolo·目标检测·计算机视觉·pyqt
星沁城1 小时前
240. 搜索二维矩阵 II
java·线性代数·算法·leetcode·矩阵
脉牛杂德1 小时前
多项式加法——C语言
数据结构·c++·算法
legend_jz1 小时前
STL--哈希
c++·算法·哈希算法
kingmax542120081 小时前
初三数学,最优解问题
算法
一直学习永不止步2 小时前
LeetCode题练习与总结:赎金信--383
java·数据结构·算法·leetcode·字符串·哈希表·计数
小刘|2 小时前
《Java 实现希尔排序:原理剖析与代码详解》
java·算法·排序算法