如何在R中计算满足v2[i]=v1[i-1]-2v1[i]+v1[i+1]的向量？

核心思路：利用二阶差分简化计算 #

你需要的公式 v2[i] = v1[i-1] - 2*v1[i] + v1[i+1] 本质就是二阶差分的变形：

• 一阶差分diff(v1)得到的是 v1[2]-v1[1], v1[3]-v1[2], ..., v1[20]-v1[19]
• 二阶差分diff(diff(v1))计算的是 (v1[3]-v1[2]) - (v1[2]-v1[1]) = v1[3] - 2*v1[2] + v1[1]，和你的公式完全对应（仅索引偏移一位）

用二阶差分就能快速实现需求，无需循环或rollingapply。

具体实现代码 #

v1 <- 1:20
# 直接计算二阶差分，得到v2的有效值（长度为length(v1)-2）
v2 <- diff(diff(v1))

# 若需要v2和原向量长度一致（首尾补NA）
v2_full <- c(NA, v2, NA)

逻辑验证 #

拿前几个元素手动核对：

• v1[1]=1, v1[2]=2, v1[3]=3 → v2[2] = 1 - 2*2 +3 = 0，diff(diff(v1))[1]结果为(3-2)-(2-1)=0，完全匹配
• v1[2]=2, v1[3]=3, v1[4]=4 → v2[3] =2 -2*3 +4=0，二阶差分结果同样为0，计算正确

时间序列库替代方案 #

如果想用时间序列相关包（比如zoo），底层依然基于差分逻辑，示例如下：

library(zoo)
v1_zoo <- zoo(v1)
v2_zoo <- diff(diff(v1_zoo))
# 转换回普通向量
v2 <- as.vector(v2_zoo)

内容的提问来源于stack exchange，提问作者Anthony Tan