针对60Hz/240Hz切换信号的Butterworth低通滤波方案咨询

背景 #

此前使用4阶低通Butterworth滤波器处理240Hz采样频率的信号，相关实现代码及滤波处理函数如下：

原有滤波器设置代码 #

library(signal) 
fo <- 4                   # 滤波器阶数
fs <- 240                 # 采样频率
fc10 <- 10                # 截止频率
w <- 2 * (fc10 / fs)      # 归一化临界频率（基于Nyquist频率）
ft <- "low"               # 滤波器类型：低通
 
bf <- butter(fo, w, type = ft)

原有滤波处理函数（首尾反射填充） #

pad_filt_unpad <- function(x, bf) {
  # 数据填充：首尾各复制fs个采样点
  padded_data <- c(rep(first(x), fs), x, rep(last(x), fs))
  # 双向滤波消除相位偏移
  filtered_data <- filtfilt(bf, padded_data)
  # 去除填充部分，还原原始长度
  filtered_data <- filtered_data[(fs + 1):(length(filtered_data) - fs)]
  return(filtered_data)
}

现有新信号的采样频率会在60Hz和240Hz之间切换（无中间值），若直接将60Hz数据复制4倍后沿用原有240Hz滤波器，会引入高频伪影等异常效果，以下是两种可行解决方案：

可行解决办法 #

方案1：按采样频率分段生成对应滤波器 #

核心逻辑：保持10Hz截止频率、4阶滤波器的参数不变，针对不同采样频率重新计算归一化临界频率，生成专属滤波器后分段处理信号。

实现代码：

# 通用滤波器生成函数
generate_butter_filter <- function(target_fs, cutoff = 10, order = 4) {
  w <- 2 * (cutoff / target_fs)
  butter(order, w, type = "low")
}

# 分段处理切换采样频率的信号
process_switching_signal <- function(signal_data, sampling_freqs) {
  filtered_result <- numeric(length(signal_data))
  current_fs <- sampling_freqs[1]
  current_bf <- generate_butter_filter(current_fs)
  segment_start <- 1
  
  # 遍历信号，识别采样频率切换点并分段处理
  for (i in 2:length(sampling_freqs)) {
    if (sampling_freqs[i] != current_fs) {
      # 处理当前频率段的信号
      segment_data <- signal_data[segment_start:(i-1)]
      filtered_result[segment_start:(i-1)] <- pad_filt_unpad(segment_data, current_bf)
      # 更新当前频率及对应滤波器
      current_fs <- sampling_freqs[i]
      current_bf <- generate_butter_filter(current_fs)
      segment_start <- i
    }
  }
  # 处理最后一段信号
  filtered_result[segment_start:length(signal_data)] <- pad_filt_unpad(signal_data[segment_start:length(signal_data)], current_bf)
  
  return(filtered_result)
}

说明：

• 每个采样频率段使用匹配的滤波器，保证10Hz截止频率在不同采样率下的物理意义一致
• 保留原有首尾填充+双向滤波的逻辑，避免边缘失真和相位偏移

方案2：将60Hz信号重采样至240Hz后统一滤波 #

若信号是分段切换（而非逐点切换），可先将60Hz段重采样到240Hz，再用原有240Hz滤波器处理，按需转回60Hz（若需要）。

实现代码：

library(signal)

# 将60Hz信号重采样至240Hz
resample_60to240 <- function(signal_60) {
  resample(signal_60, 240, 60)
}

# 示例：处理两段切换的信号
signal_segment1_60 <- ... # 60Hz采样的信号段
signal_segment2_240 <- ... # 240Hz采样的信号段

# 重采样60Hz段至240Hz
signal_segment1_240 <- resample_60to240(signal_segment1_60)
# 统一用原有240Hz滤波器处理
filtered_segment1 <- pad_filt_unpad(signal_segment1_240, bf)
filtered_segment2 <- pad_filt_unpad(signal_segment2_240, bf)

# 若需还原60Hz采样率，可反向重采样
filtered_segment1_60 <- resample(filtered_segment1, 60, 240)

说明：

• 使用signal包的resample函数重采样，比单纯复制采样点更合理，不会引入额外高频伪影
• 若最终输出不需要60Hz采样率，可直接保留240Hz的滤波结果

为什么不建议复制60Hz采样点？ #

直接复制60Hz数据4倍的方式存在明显缺陷：

• 重复采样点会在频谱中引入60Hz整数倍的谐波分量，导致滤波后信号失真
• 复制后的信号频谱与真实240Hz采样信号完全不同，无法保证10Hz截止频率的有效性

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