如何在R语言中创建(t-2)×t维度的三对角矩阵？

要生成这种每行包含1, -2, 1三个非零元素的三对角矩阵（滑动窗口式的三对角线结构），可以用两种实用方法：基础R手动构造，或者用Matrix包的高效函数。

方法1：基础R手动构造 #

假设t是目标列数，行数固定为t-2，我们可以先创建零矩阵，再分别填充三条对角线：

t <- 7
n_rows <- t - 2
mat <- matrix(0, nrow = n_rows, ncol = t)

# 填充第一条对角线（行i，列i）
diag(mat) <- 1
# 填充第二条对角线（行i，列i+1）
diag(mat[, -1]) <- -2
# 填充第三条对角线（行i，列i+2）
diag(mat[, -c(1,2)]) <- 1

mat

运行结果：

[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,]    1   -2    1    0    0    0    0
[2,]    0    1   -2    1    0    0    0
[3,]    0    0    1   -2    1    0    0
[4,]    0    0    0    1   -2    1    0
[5,]    0    0    0    0    1   -2    1

这就是你需要的5×7（t=7）结构矩阵，三条对角线分别由重复的1、-2、1构成。

方法2：使用Matrix包的bandSparse函数（适合大矩阵） #

如果t取值很大，手动构造效率较低，可借助Matrix包的稀疏矩阵函数，节省内存且运行更快：

library(Matrix)
t <- 7
n_rows <- t - 2

mat_sparse <- bandSparse(n_rows, t, 
                         k = c(0, 1, 2),  # 对角线偏移：0=主对角线，1=右移1位，2=右移2位
                         diagonals = list(rep(1, n_rows), 
                                          rep(-2, n_rows), 
                                          rep(1, n_rows)))

# 转换为普通矩阵查看结果
as.matrix(mat_sparse)

输出结果与方法1完全一致。

原代码无效的原因 #

你之前使用的diag(c(1,-2,1), nrow=5, ncol=7)只会循环将向量填充到主对角线位置，无法识别并填充其他对角线，因此只能得到主对角线的循环值，而非三对角结构。

内容的提问来源于stack exchange，提问作者Jorge Paredes