问题标签 [toeplitz]

我正在尝试制作一个如下所示的 numpy 数组：

所以这涉及更新主对角线和它上面的两条对角线。

这样做的有效方法是什么？

任何人都知道在matlab中构建以下方阵的快速有效方法

除了对角线的第一个和最后一个元素以及两个相邻对角线上的 -1 之外，它在对角线上有 2。

这是一个 6x6 的示例，我想在 Matlab nxn 中尽可能快速高效地生成一个示例。我尝试使用内置函数 kron 但无法摆脱它。

非常感谢

在 Octave 或 Matlab 中，有一种简洁、紧凑的方法可以创建大型 Toeplitz 矩阵，例如：

这样可以节省大量时间，否则这些时间会通过使用额外的代码行来用数十个或数百个零填充矩阵。

然而，我似乎无法对 scipy 或 numpy 做同样的事情，尽管这两个库都有 toeplitz() 和 zeros() 函数。有没有类似的方法可以做到这一点，还是我必须自己制定一个例程才能做到这一点（不是一个大问题，但仍然很麻烦）？

谢谢，

F。

离散卷积可以通过 Toeplitz 矩阵执行，如下所示（Wiki 文章）：

请注意，这与一般 Toeplitz 矩阵的形式并不完全相同，但它经历了各种移位和零填充。

有没有办法在 numpy 中实现这一点，纯粹基于等roll，hstack即不使用任何for循环？我已经尝试了各种各样的转变，但我无法真正将它变成上面显示的形式。

考虑一个离散的动态系统，其中x[0]=rand()表示系统的初始条件。

我通过以下步骤生成了一个 m x n 矩阵——生成m具有m不同初始条件的向量，每个向量都有维度N（N 表示样本或元素的数量）。该矩阵称为 R。如何使用 R 创建 Toeplitz 矩阵T？数学上，

托普利兹矩阵T=

我尝试使用，toeplitz(R)但尺寸发生了变化。从数学上看，尺寸应该没有变化。

我想使用 toeplitz 命令在 Matlab 中创建一个包含复杂元素的对称矩阵。但是，如果我为 toeplitz 命令提供复杂的条目，它会返回一个 Hermitian 矩阵，也就是说，虚部的符号在对角线的上方和下方颠倒。

如何创建一个对称矩阵，其中虚部的符号在对角线上方和下方相同，如下所示：

我正在尝试在 java 中初始化一个 Toeplitz 矩阵。我希望它有这种形式
6 -4 1 0 0 ... 0 -4 6 -4 1 0 ... 0 1 -4 6 -4 1 ...0 ................ 0 ... 1 -4 6 -4 1 0 ... ...1 -4 6-4 0 .. ... 0 1 -4 6

我意识到问题出if(j>i)在data[i-j-1]. 我试图改变它，但我得到了IndexOutOfBounds错误。这是我到目前为止写的代码

输出是

我有下一个问题，如何在知道矩阵大小为 NxN 的情况下构建具有特定值的矩阵。

这是我的问题

我正在尝试使用下一个代码：

但是正如您所看到的，(1,1) 和 (N,N) 中没有所需的值，那么我该如何构建这个特定的矩阵呢？MATLAB 中的代码如何？

我有 anxn symmetrix toeplitz matrix T，一个v长度为 n 的向量，我想快速计算矩阵向量乘积T%*%v。R 中是否有一个包可以使用快速傅立叶变换计算方法T%*%v（或其他方法，如果存在的话）？例如，Matlab 有 Toeplitzmult 包。

模块 scikits.talkbox 包含一些用于 Levinson-Durbin 递归的 C 代码。不幸的是，这段代码在最新版本的 Python 中不起作用，我想用纯 Python 实现替换它。（速度不是问题，只要它有效。）

损坏的 C 函数的文档字符串如下：

我看到有一个scipy.linalg.solve_toeplitz看起来像我想要的功能。但是，它无法指定顺序，而是将数组元组作为输入。

我承认我在这里有点超出我的深度，并且不完全理解这段代码应该做什么。有没有一种简单的方法可以用 Numpy's 替换对损坏的 C 函数的调用solve_toeplitz？