我会使用线性插值。在base R中有一个approx()函数的实现。语法如下所示
y2 <- approx(x1, y1, x2)$y
其中x1和y1是一个频率的时间和值向量，x2是另一个频率的时间向量。我想这取决于你想插值到哪个频率。

library(tidyverse)
library(lubridate)

# simulate data
origin <- as_date("1970-01-01")

high_freq <- tibble(
  freq = "high",
  tme = Sys.time() + seconds(seq(0, 1000, by = 10)),
  x = as.double(tme, origin = origin),
  y = sin(x * 1e-2)
)

low_freq <- tibble(
  freq = "low",
  tme = Sys.time() + seconds(seq(0, 1000, by = 100)),
  x = as.double(tme, origin = origin),
  y = as.numeric(scale(x))
)

bind_rows(high_freq, low_freq) %>%
  ggplot() + aes(tme, y) + geom_point() + facet_wrap(~freq, ncol = 1, labeller = label_both) + 
  ggtitle("both low and high freq")

# low freq values interpolated to high freq axis
high_freq$z <- approx(low_freq$x, low_freq$y, high_freq$x)$y

high_freq %>%
  rename(high = y, low = z) %>%
  pivot_longer(c(high, low), names_to = "frequ") %>%
  ggplot() + aes(tme, value) + geom_point(na.rm = TRUE) + facet_wrap(~frequ, ncol = 1, labeller = label_both) + 
  ggtitle("low freq values interpolated to high freq axis")

创建于2023-04-21使用reprex v2.0.2