如何在R中精准识别信号的最大峰值及持续时长？

针对峰值两侧最低点识别不准确、无法计算信号持续时长的问题，可通过阈值匹配+线性插值的方式实现精准定位，以下是具体实现步骤：

步骤1：加载并预处理信号数据 #

先读取信号数据（假设数据包含Time和Peak两列）：

# 读取数据（替换为你的文件读取方式）
signal <- read.csv("你的信号文件路径.csv")
# 确保时间序列按顺序排列
signal <- signal[order(signal$Time), ]

步骤2：定位最大峰值点 #

先锁定信号中的最大峰值位置与对应值：

# 找到最大峰值的索引、时间和数值
peak_idx <- which.max(signal$Peak)
peak_time <- signal$Time[peak_idx]
peak_value <- signal$Peak[peak_idx]

步骤3：定位峰值左侧的最低点 #

遍历峰值左侧数据，找到最小值对应的位置：

# 提取峰值左侧的所有数据
left_data <- signal[1:peak_idx, ]
# 找到左侧最低点的索引、时间和数值
left_min_idx <- which.min(left_data$Peak)
left_min_time <- left_data$Time[left_min_idx]
left_min_value <- left_data$Peak[left_min_idx]

步骤4：计算右侧与左侧最低点值的交点 #

若右侧无恰好等于左侧最小值的点，用线性插值计算交点时间：

# 提取峰值右侧的所有数据
right_data <- signal[peak_idx:nrow(signal), ]
# 找到右侧第一个小于等于左侧最小值的位置及前一个位置
cross_idx <- which(right_data$Peak <= left_min_value)[1]

if (!is.na(cross_idx)) {
  # 取交点前后的两个数据点
  point1 <- right_data[cross_idx - 1, ]
  point2 <- right_data[cross_idx, ]
  # 线性插值计算交点对应的时间
  right_cross_time <- point1$Time + (left_min_value - point1$Peak) * (point2$Time - point1$Time) / (point2$Peak - point1$Peak)
} else {
  # 若右侧无低于左侧最小值的点，取最后一个点作为端点
  right_cross_time <- right_data$Time[nrow(right_data)]
}

步骤5：标记起止点并计算持续时长 #

# 标记起止点与峰值点
signal$startEndPeak <- NA
signal$startEndPeak[left_min_idx] <- "start"
signal$startEndPeak[peak_idx] <- "peak"

# 找到右侧交点在原数据中的近似索引
right_cross_idx <- which.min(abs(signal$Time - right_cross_time))
signal$startEndPeak[right_cross_idx] <- "end"

# 计算信号持续时长
duration <- right_cross_time - left_min_time
cat("信号持续时长：", duration, "秒\n")

# 查看最终标记结果
result <- signal[!is.na(signal$startEndPeak), c("Time", "Peak", "startEndPeak")]
print(result)

补充：用专业包简化操作 #

也可使用pracma包的findpeaks函数，结合自定义阈值快速识别峰值范围：

library(pracma)
# 识别所有峰值，设置最小高度（峰值的10%）和最小宽度
peaks <- findpeaks(signal$Peak, minpeakheight = peak_value * 0.1, minpeakwidth = 2)
# 提取最大峰值的起止索引（findpeaks返回的第4、5列为起止位置）
max_peak_row <- which(peaks[, 1] == max(peaks[, 1]))
start_idx <- peaks[max_peak_row, 4]
end_idx <- peaks[max_peak_row, 5]

# 计算持续时长
start_time <- signal$Time[start_idx]
end_time <- signal$Time[end_idx]
duration <- end_time - start_time

这种方法通过左侧最低点的阈值匹配+线性插值，能精准定位峰值的起止端点，解决原有代码识别不准确的问题。

内容的提问来源于stack exchange，提问作者M Sandoval