请问怎么用R语言做lanczos滤波呢？

小跳蛙吃奥利奥

想对某一年的数据进行10~20d的滤波，但是R语言里面好像没有这个滤波器的函数，所以自己想着模仿写了一个，但是感觉公式好像有点问题？然后参数啥的设置不对，有大佬帮忙看看吗？Lanczos_filter <- function(n, fs, f1, f2) {
  a <- (n - 1) / 2
  h <- sin(2 * pi * f2 * (-(a):(a)) * fs) / (-(a):(a)) / pi - sin(2 * pi * f1 * (-(a):(a)) * fs) / (-(a):(a)) / pi
  h[a+1] <- (2 * (f2 - f1)) / fs
  return(h)
}

# 定义滤波器参数
fs <- 1/24  # 采样频率
f1 <- 10     # 低通截止频率
f2 <- 20     # 高通截止频率
n <- 10    # 滤波器长度

# 创建Lanczos带通滤波器
h <- Lanczos_filter(n, fs, f1, f2)

# 应用滤波器，x是需要滤波的数据，日资料来的
y <- signal:: filter(h, 1, x)

一大碗年糕

https://rdrr.io/github/jrevenaugh/TSAUMN/man/lanczos.html

小跳蛙吃奥利奥

一大碗年糕 发表于 2023-2-17 23:59
https://rdrr.io/github/jrevenaugh/TSAUMN/man/lanczos.html

按照这页的代码创建了滤波器之后怎么应用呢？是用filter函数吗？

一大碗年糕

小跳蛙吃奥利奥 发表于 2023-2-20 23:09
按照这页的代码创建了滤波器之后怎么应用呢？是用filter函数吗？

不好意思我也没怎么了解，只是搜了一下有

小跳蛙吃奥利奥

一大碗年糕 发表于 2023-2-20 23:50
不好意思我也没怎么了解，只是搜了一下有

好的也谢谢你！